The Data Exploration and Integration Lab (DEIL) and the Urban Data Lab (UDL) at the Center of Special Business Projects (CSBP) are grateful to Statistics Canada, particularly Christian Wolfe, Shujaat Ansari and Serge Godbout, for their knowledge of the Business Register, Dr. Mahamat Hamit-Haggar for coordinating the editorial process and internal reviews, Chris Li for the institutional review, Dr. Bjenk Ellefsen, for elaborating the vision and supporting the feasibility around the project, Dr. Ala’a Al-Habashna for initiating the data processing, and Zheng Yu and Wafa Ashraf for the final technical revision of the document. We are also grateful to Dr. Stephen Tapp and Patrick Gill from the Canadian Chamber of Commerce (CCC) for their review, comments, and suggestions.
Executive summary
Businesses don’t just choose a city to locate their operations; they choose a neighbourhood.Note In turn, the clustering of businesses in a neighbourhood shapes the economic opportunities and quality of life in the area. Despite the relevance of these local-level dynamics documented in the literature, research on business clusters in Canada has focused primarily on the regional or metropolitan scale. This focus has limited the possible applications of cluster analysis for urban planning, infrastructure development, and local development by actors operating or delivering programs at the local scale.
The steady improvement in the geolocations of business data is providing new analytical opportunities. This paper presents a method to define business clusters at a granular sub-metropolitan level. Using data from Statistics Canada’s Business Register (BR) for selected industries, employment locations at the establishment level are spatially distributed within their respective dissemination block (DB) (a block in urban and rural areas). A spatial kernel density estimation (KDE) approach is performed on these employment locations to define the boundaries of business clusters. A new approach to define the kernel bandwidth is detailed since the traditional Silverman’s rule bandwidth method fails, in the case of our applications, to directly recognize the configuration of the DB structure within the cities of interest. Results are developed for three industries (manufacturing, retail trade, and accommodations and food services), as well as some industry clusters as defined by Delgado et al. (2014), for four major metropolitan areas (Montreal, Toronto, Winnipeg, Vancouver).
The results are mapped for each type of cluster and metropolitan area showing different spatial configurations for different industry sectors. As expected, retail trade and accommodation and food services clusters are relatively more scattered across metropolitan areas, compared to manufacturing clusters. However, simple statistics on establishments and employment counts show that the geographic boundaries of neighborhood clusters generated by this analysis capture most of the employment and establishment counts located in the Census Metropolitan Area (CMA) of reference and associated with the industry cluster. For instance, the manufacturing sector cluster contains 89.7% (Montreal), 94.4% (Toronto), 91.9% (Winnipeg), and 90.8% (Vancouver) of total manufacturing employment in the respective CMAs. The results also point to greater co-location of specific types of business, such as retail trade and accommodation and food services. This preliminary analysis appears promising in revealing patterns of business co-location in defined neighborhood areas.
The methods used to define these neighbourhood-level clusters open new opportunities for timely analysis of business conditions at the local level, as well as broader types of analyses at the neighborhood-level (e.g., social disparities, and quality of life), accounting for the business composition of the area. The use of geographic boundary files of specific business clusters, as a geofencing tool, can be employed to monitor local business performance and trends in combination with other Statistics Canada data holdings or alternative data sources, such as mobility flows.
This project proposes an experimental methodology and a set of experimental industry clusters. People’s feedback is valued. Commentaries, and suggestions can be communicated to the lead author, Jérôme Blanchet (819-576-5502), Unit Head at the Data Exploration, and Integration Laboratory (DEIL) and Urban Data Laboratory (UDL), of the Centre for Special Business Projects (CSBP), Statistics Canada.
Introduction
The academic and policy debate on business clusters has spanned over three decades. In most of the ensuing literature and policy documents, business clusters are defined as a geographic concentration of interconnected enterprises, organizations, and institutions within a specific industry or sector (Wolfe and Gertler, 2004). Theory and evidence suggest that spatial proximity and agglomeration facilitate collaboration, linkages, resource sharing, and other synergies. Thus, policy support for clusters draws on the idea that the proximity of businesses and supporting organizations in a particular geographic area fosters synergies, innovation, and competitive advantages (Bekar and Lipsey, 2001).
The cluster literature in Canada has delved into the effect of clusters on business performance, employment, and wages (Lucas et al., 2009; Niosi and Bas, 2001; Steiner and Ali, 2011; Spencer et al., 2010). It has examined policies in support of clusters, assessed their effectiveness (Niosi and Bas, 2001), and it has explored methodologies to identify and map clusters across space (Spencer, 2014). Overall, the insights generated by this literature are that spatial clustering of businesses creates an environment conducive to innovation, resource efficiency, collaboration, and overall economic development. These clusters contribute to the growth and success of individual businesses while enhancing the competitiveness and resilience of regions and metropolitan areas.
Most of the Canadian literature has focused on a regional scale or on individual metropolitan areas, and for many applications and policy purposes, an analysis at the regional level (city, metropolitan area, or labor market area) will remain adequate. Nevertheless, there is a growing number of applications which require neighborhood-level analysis and provide unique insights to both local-level actors (municipalities and other local business organizations) as well as to provincial or federal actors. The demand for geographically granular data on business conditions and trends is constantly growing from federal agencies and local stakeholders, such as municipalities, business organizations, and the business community.
This paper brings the analysis of business clusters to a more granular geographic scale by developing a methodology for identifying business clusters at the neighborhood level. The proposed method identifies clusters of businesses at the DB level, which is one of the most granular spatial units of analysis defined by Statistics Canada.Note The method is developed with an application to four census metropolitan areas (CMAs) of different sizes and for different industry cluster specifications, including simple 2-digit North American Industry Classification System (NAICS) groups as well as industry clusters resulting from groupings of NAICS codes, as defined by Delgado et al. (2014).
The increasing accuracy of geolocation for establishments in the BR of Statistics Canada enables this analysis, while the possibility offered by business microdata linkages creates opportunities to explore a multitude of business performance dimensions.Note In this context, one of the main challenges in this type of analysis is to preserve the confidentiality of sensitive business information while providing valuable insights to the business community and policymakers. Hence, this paper remains a first exploration at the intersection of the highest level of granularity and confidentiality preservation of business information.
This paper is organized into five main sections. The next section presents a cursory review of the existing literature on cluster analysis in Canada, intended to highlight the information gap at the neighborhood level and the main motivational aspects for this research, as well as examples of neighborhood-level analysis from other countries. This is followed by a detailed presentation of the data and methodology applied in this analysis, which involves a new approach for the kernel bandwidth calculation. The next section presents selected results and validates the findings, followed by a discussion on further development and possible applications of the cluster delineations. Finally, the annex includes a large supplementary content of 23 high-resolution cluster maps.
Why a neighborhood dimension of clusters
Most studies on business clusters in Canada use regions or cities as geographic units of analysis. For many applications and policy analyses, that geographic scale provides an adequate level of granularity. Business performance indicators have become relatively abundant at the municipal or regional scale. Within cities or regions, it is assumed that business proximity is sufficient to allow for the interactions that underpin the very concept of clusters and the benefits that derive from it. Some of the well-known clusters, for instance, Silicon Valley in Northern California, and the Research Triangle Park in North Carolina, are spread over several municipalities and a regional scale appears appropriate to study the dynamics and development of these clusters.
Nevertheless, there is evidence of clusters, such as financial, cultural, retail, or manufacturing clusters, concentrating in specific neighborhoods within a metropolitan area. Similarly, there is evidence that spatial disparities and differences in business performance may be as pronounced at neighborhood scale as they are at the regional scale (Wheeler, 2006; OECD, 2018). Several stakeholders and policymakers, such as local business organizations and municipalities, are operating with a neighborhood lens and are developing or delivering policies that impact businesses in specific areas within a municipality. Hence, these stakeholders require a more geographically granular analysis of clusters and cluster performance.
To respond to these information needs, different streams of research have analyzed business clusters at the neighborhood level, looking at location choices within metropolitan areas, the impact of municipal regulations and policies on the formation and growth of clusters, the role of local associations, and the impact and economic spillover of clusters on the surrounding neighborhoods. This research comes from different disciplinary perspectives. Ranging from the more traditional analysis of business clusters to urban planning studies, and applied research and analysis generated by local associations, boards of trade, and municipal planning departments. The remainder of this section provides an overview of this literature, highlighting key insights and pointing to the current information gap in the Canadian context.
Spatial disparities within cities, at neighborhood level, are a well-known and researched phenomenon (OECD, 2018). These disparities are reflected and exacerbated by business location choices and clustering in different neighborhoods. For example, Wheeler (2006) shows that positive business growth in the St. Louis metropolitan area, in the US, is, in fact, the result of substantial business growth in one neighborhood combined with a decline in another neighborhood. These dynamics are driven by both neighborhood and business characteristics. These insights are relevant from both a business perspective (for location choices), as well as from a municipal perspective (for policy to support business clustering and reduce neighborhood socio-economic disparities). The importance of a neighborhood dimension of spatial clusters is also highlighted by Gabaix (2011), who uses satellite imagery data to define clusters of population density. These results suggest that data generated from cluster density-based approaches for city boundaries are easier to integrate within a spatial model compared to regions defined by administrative boundaries directly.
Insights from analysis on business clusters at the neighborhood level are related to the location choices of businesses. Wheeler (2006) notes that prospective investors do not choose just a region or city, they also choose a neighborhood, which in some contexts may refer to a specific industrial park or business district. Therefore, understanding the dynamics of industrial parksNote within a metropolitan area, and their surroundings, is of specific relevance. Similar insights come from Arauzo-Carod (2021), who examines the location choices of high-tech firms at the neighbourhood level in Barcelona. This analysis shows that both neighbourhood characteristics and amenities matter in business location choices and that spatial spillovers are relevant for some high-tech industries.
Business clusters do not just determine economic composition and job availability in a neighborhood. Several studies highlighted that they may also shape the quality of life of neighborhoods and their attractiveness for residential use (Shybalkina, 2022; Stern and Seifert, 2010). Research has studied the impact of specific clusters on neighborhoods, particularly arts and culture clusters within metropolitan areas. Grodach et al. (2014) identify arts clusters at the regional and neighborhood levels, using Zip Codes, for US metropolitan areas of different sizes. Moreover, their findings show that arts industries exhibit distinct metropolitan area and neighborhood-level location patterns. It is also revealed that while many of the characteristics of arts clusters are place specific, the arts are associated with broad measures of local innovation and development, suggesting that these business clusters can play a larger role in economic development for metropolitan areas.
Another stream of literature on intra-urban business clusters focuses on the delineation and role of business districts. This literature has often been categorized under the heading of central business districts (CBD) analysis (Meltzer, 2012; Yu et al., 2015); it has both methodological and policy relevance, with part of this literature focusing on the modeling aspects, while other parts delving into the role and dynamics of organizations related to the management, development, and promotion of these business districts. The methods used in delineating CBD range from use of remote sensing data (Taubenböck et al., 2013) to census data on job density (Yu et al., 2015).
The approach taken in the present analysis is inspired by this literature and, in particular, by the analysis of Sergerie et al. (2021), which aimed to develop a method to identify the geographic boundaries of Canada’s downtown neighbourhoods. Sergerie et al. (2021) use spatial KDE to calculate a density surface of job location data at the dissemination area (DA) level, allowing for comparisons of these areas across Canada.
Business cluster formation, at the neighborhood level, like at the regional level, is not just a spontaneous process or a cumulative result of historical accidents. That is, it is not a random phenomenon, and it can be explained. Municipalities play a key role in shaping, developing, and supporting the clustering of businesses in specific areas (Zhang, 2019). Typically, this is done through zoning, which is the regulatory method used by municipalities or local governments to set rules that define the activities and buildings that can occur in a certain location. In this way, municipalities provide the space, infrastructure, and services for the development of industrial parks, as areas dedicated to industrial uses.
Like municipalities, local business associations are key actors in supporting the development of local business clusters (Dhamo et al., 2023). In Toronto, for instance, the Toronto Board of Trade (2021) released a study mapping five types of district areas across the broader metropolitan region, which include a metro center, goods production and distribution areas, services and mixed-use areas, regional centers, and knowledge creation centers. In parallel, that municipality is home to 84 Business Improvement Areas (BIAs). These are local associations of businesses aiming to support competitive and attractive business areas for consumers and new businesses.
The proactive role that local actors can take in shaping economic development and prosperity of their area explains why neighborhood-level analysis is becoming increasingly relevant. Analyzing local clusters or concentrations of businesses within a neighborhood can provide insights into consumer behavior, local employment, and the overall economic health of the neighborhood. These local actors operate in geographically defined ecosystems for which data may be generated from local sources or analyzed at the local level. What is missing, particularly in the Canadian context, is the broader and comparative framework that would allow neighborhood cluster analysis nationwide, with standardized definitions across jurisdictions, which is the information gap that this paper is intended to fill. Improvements in georeferencing of business microdata and advancements in spatial data analytics with large databases are making it increasingly feasible to conduct fine-grained analyses at the neighborhood scale.
A proposed methodology
The general approach used in this analysis is drawn from the work of Sergerie et al. (2021) on the definition of downtown areas for metropolitan areas in Canada. That analysis applies spatial KDEs on the geolocation of total jobs derived from the place-of-work status variable of the Census of Population; the geographic unit of analysis in that application is the DA.
In the present study, business clusters are defined with an analogous method, using spatial KDEs applied to the geolocation of employment recorded at the establishment level. Data on establishments are extracted from the BR for selected industries, and the unit of analysis is the DB, a more granular unit than the DA. Given the substantially greater level of granularity (both for selected industries and geography), the methodology used to define business clusters presents several additional steps compared to the workflow outlined by Sergerie et al. (2021). The following sections describe the proposed methodology in detail.
Study Areas
Four study areas were selected for the development of the methodology, representing different levels of urban density in Canada. These study areas are the CMA of Montréal, Toronto, Winnipeg, and Vancouver. Each CMA comprises of a different number of municipalities (census subdivisions), with neighborhoods reporting substantially different population and employment densities and degrees of urbanization.
The geographic unit of analysis used for the geolocation of businesses is the DB. A DB is an area bounded on all sides by roads and/or boundaries; that is, in urban areas, it is what is commonly referred to as a block. DBs are part of Statistics Canada’s standard geographic areas for dissemination, and they are the smallest geographic area for which population and dwelling counts are disseminated with coverage for all the territory of Canada.Note
Business Register
The data used in the analysis come from Statistics Canada’s BR, which is the continuously maintained central repository of baseline information on businesses and institutions operating in Canada. For this analysis, data are for the reference period of December 2023.
In the BR, industry sectors are defined by NAICS codes. The use of BR data presents advantages when compared to other possible data sources, such as place-of-work data from the Census of Population. Establishment-level data from the BR is classified with more detailed NAICS (6-digit codes), which allows for detailed custom clusters. Employment data in the BR are also updated with greater frequency. Although they are not as accurate as dedicated employment statistics, they are a timelier alternative to census data.
The NAICS codes included in this analysis represent six different industry clusters. Three industry clusters are composed of 2-digit NAICS codes, and the remaining are defined according to the work by Delgado et al. (2014). These industry clusters and their corresponding NAICS codes are summarized in Table 1.
The industry clusters generated at the 2-digit NAICS level, include the Manufacturing Sector (NAICS codes 31, 32, and 33), the Retail Trade Sector (NAICS codes 44 and 45), and the Accommodations and Food Services Sector (NAICS code 72).
The three industry clusters defined according to Delgado et al. (2014), are as follows. First, the Distribution and Electronic Commerce (cluster 10); this cluster consists primarily of traditional wholesalers as well as mail-order housesNote and electronic merchants. The companies in this cluster mostly buy, hold, and distribute a wide range of products such as apparel, food, chemicals, gases, minerals, farm materials, machinery, and other merchandise. The cluster also contains firms that support distribution and electronic commerce operations, including packaging, labelling, and equipment rental and leasing. The second cluster is Financial Services (cluster 16); this cluster contains establishments involved in aiding the transaction and growth of financial assets for businesses and individuals. These firms include securities brokers, dealers, and exchanges; credit institutions; and financial investment support. Insurance firms are located in a separate Insurance Services cluster. Finally, the Hospitality and Tourism (Cluster 22) contains establishments related to hospitality and tourism services and venues. This includes sports venues, casinos, museums, and other attractions. It also includes hotels and other accommodations, transportation, and services related to recreational travel, such as reservation services and tour operators.
The hierarchical structure between the 3 industry clusters generated at the 2-digit NAICS level, and the 3-industry cluster defined according to Delgado et al. (2014) is not straightforward. That is, the several 4-digit NAICS level used according to Delgado et al. (2014) are not necessarily composed of the first 2-digit 31, 32, 33, 44, 45 and 72. Furthermore, for space convenience and simplicity, the recurrence of the 6 industry clusters is not necessarily consistent across this paper. That is, some results analysis and methodology steps focus on the 3 industry clusters generated at the 2-digit NAICS level alone, while some other focuses on the 6 industry clusters.
Table 1
Industry clusters with associated NAICS codes
Table summary The information is grouped by Industry cluster (appearing as row headers), , calculated using (appearing as column headers).
Industry cluster
NAICS codes included
Source: authors’ computations from the BR database.
The spatial concentration of industries was determined using a spatial KDE method. The spatial KDE method is a non-parametric technique that estimates the probability density function of a random variable over a spatial domain. For each metropolitan area, clusters were identified using KDE results by aggregating adjacent DBs with a minimum density of employment in each industry or combination of industries.
Employment counts from establishment-level data of the BR were geocoded to the DB spatial boundaries files. The total number of employees in each DB was then calculated. Job locations representing each employee were randomly and uniformly distributed within the DB boundary.Note This processing step must be acknowledged, not as a lack of accuracy and weakness of the data, but as a way to leverage the methodology forward. That is, this approach smooths the relatively sparse spatial distribution of jobs within the DB, and eventually facilitates the KDE process by generating a relatively more continuous distribution.Note The uniformity argument was made to not prioritize any sub-region of the DB during the randomization process. That is, the concentration of establishment density within some specific spots of the DB doesn’t affect the job randomization process. Some sections of this report below describe in more detail the randomization process.
The density estimation section of this paper is made of two sub-sections. First, we document partial information about the Polynomial Kernel Density Function, with a focus on the necessary arguments to understand the use and significance of the kernel bandwidth. We also explain why the Silverman kernel bandwidth approach is not a good fit for our application, and then, we describe a new method for the computation of the bandwidth. Second, we describe the full information about the Polynomial Kernel Density Function. That is, all remaining arguments of the function not described so far. We also document the random processes used.
Kernel Density Bandwidth Methodology and Other Parameters
Before proceeding to the full KDE description, we specify partial details of the Polynomial Kernel Density Function,
, which is the main instrument of the KDE approach. The function involves a minimum of 3 arguments: the grid output cell id,
, the grid output cell centroid,
, and finally, the bandwidth or radius, . The polynomial functional form was selected according to Sergerie et al. (2021), even though a normal and uniform form was also available. was defined geometrically, and not weighted according to the establishment spatial density concentration of the DB. Sets of specifications and testing for,
and
were considered for economic interpretability of the results and computational efficiency matter. Those are briefly documented below.
The spatial resolution, or dimension of a cell (or tile) for a grid
is a parameter of choice. For the purpose of this analysis, a squared grid with a cell dimension of 50x50 metres was adopted. Following the testing of different specifications (e.g., a range of 10x10 metres to 100x100 metres), the 50x50 metres grid was adopted to strike a balance between spatial detail and computational efficiency. Figure 1 below shows the distribution of the number of 50x50 metres cell per DB, in the four metropolitan areas. The Figure confirm that a minimum and reasonable number of cells are available in most of the DB of each CMA. Furthermore, the median number of cells per DB ranges from 7.5 to 12.5. Finally, the Figure also shows variations across CMAs for the range of total number of cells that can fit within a single DB. This variation is particularly important, since it captures the characteristics of both large and small CMAs within Canada. That is, smaller CMAs generally involves a larger range of DB dimension, while the largest CMAs of our study are made of a smaller and lower range of DB dimension. It is worth noting that Montreal is the only CMA to get a symmetric distribution of DBs, as the opposite of the 3 other CMAs where the distribution is relatively more skewed.
Description for Figure 1
Table figure 1
Table summary This table displays the results of Table figure 1 , calculated using (appearing as column headers).
"Montréal CMA"
"Toronto CMA"
"Winnipeg CMA"
"Vancouver CMA"
Source : authors’ computations from the BR database.
Min
0
0
0
0
1st Qu
4
5
6
6
Median
8
9
12
9
Mean
53.12
76.26
249.4
77.41
3rd Qu
14
22
48
21
Max
11,357
38,258
37,816
162,076
Similarly to the choice of grid cell dimension, the choice of the kernel bandwidth,
, also has implications for the results. To understand the relevance of the bandwidth, consider the analogy to a simple histogram. Excessively large bandwidths, i.e., a histogram with a few numbers of bars, will mask the underlying distribution. On the contrary, a bandwidth that is too small may result in a frequency of one unit for every outcome of the distribution, hindering a clear understanding of the data distribution with a flat histogram.
To define
, various specifications were tested, including the well-known B. W. Silverman’s Rule of Thumb’s Silverman, (1986) as applied in Sergerie et al. (2021). Given the grid cell size and DB configuration used in our analysis, the Silverman’s bandwidth did not succeed to include enough data content in its local surroundings. Consequently, the Silverman’s process was not effective in smoothing the spatial discrete distribution of job locations, making the KDE ineffective and resulting in fragmented clusters. The following paragraph proposes an interpretation about the reasons for the Silverman’s process to fails in providing a satisfactory bandwidth for the case of our applications in this research.
The Silverman equation captures the dispersion of the job locations around a point of reference within the CMA. From the previous section, we know that the job locations are originally geo-coded within their respective establishment fixed spatial locations and are then distributed randomly and uniformly within the boundaries of their respective DB, without prioritizing the location of the establishment or any sub-region of the DB. Therefore, we can state that the Silverman equation only encapsulates partial information about the distribution of DB dimension within the CMA. However, the equation doesn’t necessarily acknowledge the dimension of a typical DB, that is, an average or median DB. Consequently, in the specific examples of our applications, the Silverman bandwidth fails to aggregate the data across DBs and only propose kernel data transformation within the boundaries of the DB of reference, leaving the final DB level results unchanged compared to the original BR data.
We quickly provide an interpretation for this problem of why the bandwidth is too small. The Silverman bandwidth , where
is the sample size and , where is the standard dispersion of the spatial job data points distribution
to the unique mean center point of the CMA, and where μ is the median dispersion statistical moment (50th percentile) in the distribution of dispersion of
to the mean center. Note, another version of the Silverman equation substitutes the median
for the interquartile range
where
percentile minus the 25% percentile.
Formerly, the unweighted equation for the standard distance is,
where
and
are the longitude and latitude numerical coordinates of the elements of , respectively, and where
and
are the longitude and latitude numerical coordinates of the mean CMA center, respectively. Intuitively, both
and
cannot be contained within a typical or median DB of the CMA because the job locations are largely available spatially across the whole superficies of the CMA. Consequently, we assume the minimum (min) of the 2 quantities to be reasonably large enough and not able to explain a too small Silverman bandwidth. On the other hand, the sample size
is at the denominator of the Silverman formula and
can provide a very small bandwidth if
is too large. BR distribution of employees’ number per establishment can be extremely right skewed, including positive outliers reaching very large values and adding up to a very large due to the potential heaviness of the right tail. Consequently, if
is large enough to keep
large enough, then the bandwidth will be unrealistically small. This research doesn’t dedicate time to an outlier removal process for the BR data, but could be the objective of a future study. Specifically, the application of papers published by Statistics Canada, such as, Outliers in Sample Surveys by Lee et al. (1992), A Cautionary Note on Clark Winsorization by Mulry et al. (2016), A Method of Determining the Winsorization Threshold, with an Application to Domain Estimation by Martinoz et al. (2015), and On Searls’ Winsorized Mean for Skewed Populations by Rivest et al. (1995). It is also worth noting that the Silverman equation assumes normally distributed kernel density values, which is not the case for our highly right skewed density distribution presented above. Finally, we point out that Mathematica uses a modification of the Silverman’s equation, that is, , for large sample size , reaching beyond 100,000 observations. However, we did not use analytical software Mathematica for this research and decided to not leverage this version of the Silverman equation. Intuitively, taking off the large sample size of the BR at the denominator of the Silverman ratio would contribute to resolve the issue of a too large denominator and a too small bandwidth, if substituting a coefficient of 0.9 with 0.09 (around 90% reduction) is not reducing the bandwidth too much.
To address the problem related to the Silverman’s rule of Thumb in our applications, this research develops its own bandwidth methodology that focuses on the distribution of DB superficies provided by the data. In other words, our custom bandwidth accounts for the dimension of the median DB in the CMA of reference.Note This new approach ensures that the bandwidth would be sufficient to span across both a typical DB, and the neighbouring DBs. This condition is fundamental, because the final product of this project is at the DB level. Consequently, if density for all grid cells represents the aggregation of information located within the boundaries of their respective DB only, there would be no exclusive spatial information generated from this research.
Figure 2 illustrates the bandwidth approach developed for the present analysis. The grid
populated by DBs, not to be confused with the grid of output cells introduced above, assumes a local environment with equal squared DBs whose dimensions are those of the empirical median DB of the reference CMA. The median DB is derived from a projected area in the coordinate reference system (CRS) EPSG:3347.
is the side length of a median DB. Consequently,
is the volume or superficies of a median DB. The red points are the geometric centroids of grid cells within a median DB of , and should not be mistaken as the DB’s centroid.
Based on the left-sided example of Figure 2, bandwidth #2, represented as a circle, succeeds in reaching out to the totality of the requested information, which is highlighted in grey and yellow. That is, the full reference DB of the cell and all 8 direct neighborhood DBs of the reference DB, respectively. It is worth noting that it is the smallest possible bandwidth to meet these conditions. Specifically, it covers
times the requested superficies. The calculation of the bandwidth length uses a simple Euclidean distance logic; , therefore, . This logic makes sure that a grid output cell centroid overlapping at the DB’s geometric centroid will be able to capture all needed data. Furthermore, based on the right-sided example of Figure 2, the grid output cell is now located at the top-right corner of the same reference DB. Bandwidth #2 succeed to reach out to all request DBs, except for 1 DB, at the bottom left, whose superficies is reachable at 50%. This unreachable area is acknowledged as acceptable in our methodology considering new DBs now become reachable at the top-right, even though they are not a direct neighborhood of the reference DB.
Cluster heat maps produced using this new bandwidth approach yielded a smooth surface while preserving neighborhood details (see details in the annex). For clarity, this paper acknowledges the bandwidth notation as , since the custom radius is now dependent on a single argument, the Median DB (MBD) configuration, and the MDB is itself specific to the CMA. That is, the bandwidth changes across CMAs but it is fixed within CMAs, no matter the variance of concentration of job locations in the neighborhoods of the CMA. However, notation
will remain in use for space efficiency matter. Finally, it is worth noting that our custom bandwidth,
, is not dependent on the dispersion of the job locations. Consequently, if CMA A and B have the same MDB and CMA B has twice the dispersion of job locations than CMA A, then both CMAs will get the same , which is improper because the dispersion of job location matters. On the other hand, our custom bandwidth is robust to large outliers in the distribution of BR’s number of employees per establishment. An improvement of our bandwidth would include both a dependence on the MDB and job location dispersion.
Description for Figure 2
Source: authors’ design and methodology.
Legend:
,
,
, and
Uniform Multinomial Random Processes and Polynomial Kernel Density Function
Following the specification of arguments
and , this section documents the remaining details of the spatial polynomial kernel density function, PKDF, applied in the analysis.Note This kernel density model is partially inspired by Sergerie et al. (2021) regarding the way to apply a kernel density model to a granular population of businesses in Canada, however it takes the fundamental structure from statisticians Emanuel Parzen (1962) and Murray Rosenblatt (1956) who independently created the kernel density theoretical form. The notation used in this paper is our very own. This paper is about an applied contribution, with the generation of the cluster heat maps. This paper has no theoretical contribution, outside of the kernel bandwidth design (presented above). This paper uses several existing results from the literature for the explanation of the theoretical kernel density mathematical terms and concepts, and for the explanation of the random processes applied before the kernel density model (explained below), for the perspective of a spatial population of DB polygons and job locations. The remaining terms to be explained below are;
,
,
, ,
,
, and .
represents the finite non-negative number of job locations available within the circular and symmetric surrounding
(kernel area) of the grid output cell of interest
and included in the calculation of the cell total density. Therefore,
is dependent on the grid output cell id being assessed
, the cell geometric centroid location
, and the length of the kernel bandwidth
. Consequently, the number of DB it covers will vary depending on where
is located, and it can cover more than one DB or less than one DB.
is the job index for all unique jobs included in
around the grid output cell of interest
and can only range from 1 to
.
The next few pages document the random processes used in our methodology. This documentation is important to understand the structure of such processes.
is the job geographic location
for the purpose of the kernel density estimation, it is uniformly and randomly distributed
on its respective DB polygon superficies DBP, and it is not necessarily located at the location of the original establishment the job belong to. Also, job location is not allowed to be randomly distributed outside of the boundaries of its own DB, not even within the boundaries of its own DA. Furthermore, all job locations 's, included in to
and belonging to same DB, are i.i.d., that is, independent and identically distributed from the same spatial distribution. We call this random process RP.
is i.i.d. in the sense that a unique job location can be spatially and randomly allocated independently of the spatial random location of another unique job. However, like it will be explained later in this paper, the random process must also be seen as the number of unique jobs randomly allocated on a specific spatial spot of the DB. In this perspective, the larger the number of unique jobs randomly allocated on a specific spatial spot, the less likely another spatial spot will get many jobs, because the number of unique jobs in a DB is finite and not infinite. For this reason, the key insight is to visualize the phenomena as data points (a large set of unique jobs) randomly allocated to the multi-categories (several spatial spots on a DB) of a multinomial random variable (a DB).
Some perspectives about the random process RP are worth mentioning. Even though the data allocation is intended to be uniform, the random process cannot always generate a uniform coverage for a DB, and some neighborhood of the DB could be more covered than others. More specifically, the most likely event of our process RP is a spatially uniform allocation of data points within the DB, where each spatial spot gets the same number of jobs. This is our expected event . On the other hand, the most unlikely, but still possible event, is when all distinct data points overlap at the same unique spot in the DB. This is our rare or tail event .
Example for a small DB using a uniform binomial random process
To better understand the random process RP of our methodology, an analogy can be provided with a simple uniform binomial random process (BRP) with
unique very small DBs of two spatial spots and
unique jobs per DB, where is a large and even number, and is a large number. To build the intuition in a simple manner, we analyse mainly the expected and rarest event of the distribution, and some other events surrounding the expected and rarest event, and not any other details of the distribution. Analogically, an unbiased coin proposes a uniform spatial process made of two spots since both sides of the coin are of equal probability and there is no priority for any of the two sides. However, there is no guarantee that the allocation of head and tail will be even all the time. The most likely event is a perfect mixed bag of heads and tails. That is, the expected event
is a situation where each side of the coin gets an equal allocation of
unique jobs. This situation will happen to a very large proportion of the
DBs. On the other hand, the most unlikely event, but still possible event, is to get all
jobs on one side of the coin only (either head or tail). That is, all unique jobs get randomly allocated on only one of the two available spatial spots. This situation is our
event, and it will happen to a negligible proportion of the
DBs. Intuitively, all possible combinations of
jobs allocated on the spatial spot of head and tail gets the same probability of realization, but since the event of head or tail is one possible combination, respectively, and the event of a perfect mixed bag of head and tail is about a very large number of distinct combinations, then the perfect mixed bag event is more likely to become the expected event.
Formally, we have the following structure for the number of possibilities,
where , is the number of ways of picking
item among
item without order,Note the right inequality () is for the rarest events, the middle inequality (
) is for the expected event, the left inequality (
) is for the total number of possible combinations, and the
index of value 0 and value
is for tail and head, respectively. As a supplementary analysis, let’s note that
represents the number of combinations for events that are one step away from a perfectly balanced spatial allocation, respectively, and are equally less likely to happen than the expected event related to number of combination
but still very likely to happen in probability and quite close to the expected event probability. The same way around,
can be seen as the number of combinations for events that are one step forward a more uniform spatial allocation, respectively, and still very unlikely to happen in probability, but quite more likely to happen than the
event (Grimaldi, 2003).
We should also visualize the analysis in terms of the summation of equally weighted random variables (non-weighted average). That is, if we label each of the two spatial spots of the DB with numerical values +1 and -1, then for the expected event , the number of +1 is equal to the number of -1 and the average is,
If we move one step away from a perfectly balanced spatial allocation, we have an average of,
which are different than an average of 0, respectively, but still close to 0, respectively.
For
, we get either the numerical value -1, v time in a row or the numerical value +1, v time in a row, and the average is,
respectively. If we move one step forward a more uniform spatial allocation, we have an average of,
which are different than an average of -1 and +1, respectively, but still close to -1 and +1, respectively.
We now understand that the distribution of DBs is not only about a perfectly uniform spatial allocation for each DB of the CMA and a symmetric and smooth degradation of uniformity is present on both side (left and right) of the average event, if the number of jobs in each DB is large and if the number of DBs in the CMA is large. Repeating the exercise for the whole support of the distribution (and not just the average and rarest events) will generate the shape of an approximated bell curve distribution. It is now time to bridge the existing univariate central limit theorem for the binomial distribution and its finite variance to our explanation.Note That is, Asymptotically, if , then BRP converge in distribution to the univariate normal distribution
.Note In a finite context, if
is large enough, then we can approximate BRPNote such that,Note
where the term 0.5 represents the equal chance of allocating a job on one of the two spatial spots of the DB, the
term represents the expected number of jobs to be allocated on one of the two spots, which is half of the total number of unique jobs available from all establishments of the DB, and where
represents the variance of number of jobs allocated on one of the two spots of the DB (Severini, 2005). Last but not least, in a finite context, we can expect to observe a realised approximated normal distribution over the set of DBs of a large CMA if
is large enough in the CMA and
is large enough in each DB of the CMA, or in a theoretical context, a normal distribution if
in the CMA and
in each DB. The paragraph closes our binomial random process example.
Uniform multinomial random processes for large DBs
The mathematical calculations for our main random process RP are a similar idea, however, uses a uniform multinomial process (with finite variance) instead of a uniform binomial process (with finite variance) because the number of spatial spots per DB is finiteNote and large, and not only equal to 2. Furthermore, from the existing Multivariate Central Limit Theorem (MCLT) for the multinomial process, if
then the convergence in distribution of the uniform multinomial process is a multivariate normal distribution instead of a univariate normal distribution (Severini, 2005) (. In a finite context, if is very large, then the process RP can be approximated such as,
where
is our usual large number of job notation, and
is an
-dimensionNote vector made of
probabilities summing to 1 for the large number of spots
of a typical DB of a CMA. For convenience matter, and without losing too much generality,
is also divisible by
to allow an exact uniform spatial allocation.Note In the special case of our applications,Note
and has an equal probability for each of the element of the vector because the process RP is a spatially uniform multinomial random process and doesn’t prioritize any spatial spot among the full set of unique spots of the DB in the CMA. Consequently,
is the vector of
entries representing the expected number of jobs allocated to each of the
spatial spot of the DB, which is
for each spot of the DB, because
. In other words,
, and
is the variance-covariance matrix of dimension
of the multivariate normal distribution of our converging uniform multinomial random process. It is expressed in a finite context as,
and it include only two unique components among the
non-zero components due to the simplification related to the uniformity. That is, the component on the diagonal and the one off the diagonal. The matrix notation of
follows the notation of Ericson, 1969 (Appendix, equation A2) for the proper way of notating a matrix with equal terms on the diagonal, equal terms off diagonal and when the number of rows and columns is even and large.
, where
, and where
is the identity matrix of dimension
. That is,
is a diagonal matrix whose diagonal elements are the items of vector .
Also,
The variance-covariance matrix
informs that the variance of job allocation to any spatial spot is
and the covariance between any two distinct spatial spots of the DB is . The minus term of the covariance is due to the less likely chance of getting more job counts on one of the spatial spots as the job counts of the other spatial spot increases (Aitkin, 2022).
The multivariate normal distribution of our converging spatially uniform multinomial random process can be expressed in a finite context in the following manner,
All terms of the S-variateNote normal distribution
are already defined above, with the exception of
Note , which is a row vector of S-dimension, and made of
entries for the respective count of job allocation for the
unique spatial spots of our DB. Even though multivariate and made of vectors and matrices,
outputs a scalar value
and takes the shape of a univariate normal distribution (Aitkin, 2022).
reaches its maximum densityNote when , which is the expected job allocation vector presented above. That is, when . This statement was pointed out by Thompson and al, 2022 (page 78 (8) equation 2.9) for the general case of a multivariate normal distribution with log transformation. Also,
reaches its minimum density for , which is a situation when the total of unique jobs of the establishments of the DB are allocated exclusively to the first spatial spot only. The same minimum density will be reached for any other single spatial spot of the DB. Finally,
represent the determinant of the
variance-covariance matrix.
is a non-zero scalar value, because
is non-singular and the inverse matrix of ,
, does exist. This paragraph completes our explanation about the stability of the random processes involved in our methodology. In other words, we documented the idea that our random processes relate to the normal distribution and the use of random variables converges toward a well centered result and are not a synonym for inaccuracy.
It is worth noting that, since the process RP for job allocation is a random process, then, consequently, any other metric generated on top of this pre-processing is also a random variable. This basically includes all steps of the KDE program. The quantification of uncertainty due to the random pre-processes of our KDE program is not the focus of this paper, but could be the topic of future research. This topic is relevant because the local neighborhood spatial accuracy matters and could vary depending solely on some random realization of rare events or other less uniform events. However, if large number of jobs are available then the magnitude of the uncertainty is expected to be reasonable. Furthermore, as documented above theoretically, if the random processes relate to the normal distribution (that is, Student’s t-distribution with degree of freedom
), then most of the random events will be quite centered toward the perfectly balanced expected event and the two tailsNote of the distribution will be thin in density for the rare, and extremely unbalanced events, which is the opposite of what a Cauchy distribution would be (that is, Student’s t-distribution with 1 degree of freedom) (Fisher, 1925) and (Hurst, 2010).
Polynomial Kernel Density Function
We now proceed with the remaining parameters of the methodology.
is the job weight variable for job , where
is positive and equal or below , and
is equal to unity, that is, a value of 1 for all unique jobs. The alternative option would have been the extremum scenario. That is, to spatially allocate a number of job points equal to the number of unique establishments within the DB and select a weight equal to the number of unique jobs within the respective establishment for each job point. This scenario makes better use of the
weight variable because it is not equal to unity. However, this scenario is not in favor of pre-processing and cleaning the data for the purpose of the KDE, as it makes the spatial discrete distribution of jobs sparser and harder to smooth. The
setting is preferred because it is about a very large number of equal weighted entities spanning across space and pre-smoothing the original spatial distribution of jobs before applying the KDE.
For a more compact notation of the latest terms introduced in the previous paragraphs, we define the set
of dimension Note
as
and we include it as a dependent argument of the . Furthermore, is the Euclidean distance between the grid output cell geometric centroid and the random job location . It is dependent on , , and the bandwidth . The ratio
is a simple normalizer term outside of the summation and constant for all grid output cells within a CMA and that vary only across CMAs. Finally, if we substitute the two squared exponents with a repetition of their terms, the ratio,
is the key term of the kernel density function . The
ratio properly allocate a job location
closer to the grid output cell’s geometric centroid with a larger individual contribution to the total density of the cell of reference. For its part, a job location
farther from the cell geometric centroid gets a smaller individual contribution. Formally, the polynomial kernel density function generating the total density value
and can be noted as,
is dependent of , the full CMA spatial job distribution, and not just the subset of job elements of
indexed in for a specific grid cell of interest . For this paper, to simplify the methodology,
is acknowledged as a fixed distribution, and not subject to the uncertainty related to a super-population model generating a complete random realization of a CMA population, nor to a random or non-random sample aiming at representing the CMA population.
For simplicity matter, we note the ratio , included within the ratio , as the ratio . We also note
. A quick function study is essential to understand the shape of .
can be expressed as . The first derivative of
is
and
is equal to zero for
. The second derivative of
is
and
equal zero for . Consequently, is a decreasing polynomial function for , that is, .
is an increasing polynomial function for , that is, .
is concave up for , that is .
is concave down for , that is, .
Furthermore, the original domain of
is not limited to a specific subset of the real numbers . That is, . The original image of
is the set of the non-negative real numbers . That is, .
Finally, even though polynomial, the shape of
is similar to a normal distribution for domain . That is,
is similar to a Taylor Series Approximation (infinite sum of polynomial terms) for a normal distribution in a neighborhood centered at zero.
Formally,
, by definition of a univariate normal distribution of mean 0 and variance 1, and
by exact equality of experimental Taylor Series result,Note and
, by definition of an infinite summation, and
, by simple distributivity property, and
, after removing an infinity of terms of negligible respective magnitude compared to the first 3 terms, and
, after arithmetic simplification and a tough rounding of the
terms
The last approximation above, , is similar enoughNote to the expression of
and for 2 reasons: 1) the sign +/- structure is the same, 2) the exponential terms structure is the same. However, the coefficients magnitude structure is not the same. That is,
is quite stable while the other equation is sharply decreasing, which is typical for Taylor Series. Also, the third coefficient of magnitude 0.05 is enough to understand the relative respective negligence of the remaining infinity of removed terms compared to the first 3. The similarities between
and
are relevant to understanding the theoretical connection of
with the normal distribution, especially since the normal was an optional parameter in Sergerie et al, 2021. Figure 3 presents the 3 curves of interest to understand the full picture of , that is, the blue function , its red first derivative function and its green second derivative function. The x-axis is for the ratio (Bartle and Sherbert, 2011).
Description for Figure 3
Figure 3
Polynomial kernel density function (zoom out for f(x), f'(x) and f''(x))
Table summary The information is grouped by Distance Between Cell Centroid and
Job Location to Bandwidth Size Ratio (x) (appearing as row headers), , calculated using (appearing as column headers).
Distance Between Cell Centroid and
Job Location to Bandwidth Size Ratio (x)
Original Function (blue)
f(x)
First Derivative (red)
f'(x)
Second Derivative (green)
f''(x)
Source : authors’ computations from the theoretical kernel density estimation function of Parzen (1962) and Rosenblatt (1956). Legend: Blue = function of interest
, red = first derivative of
, and green = second derivative of
-2.50
26.32024871
-50.13380707
67.80000576
-2.49
25.8222911
-49.45866799
67.22819388
-2.48
25.33105633
-48.78923556
66.65867383
-2.47
24.84648745
-48.12548687
66.09144561
-2.46
24.36852773
-47.46739901
65.52650923
-2.45
23.89712067
-46.81494905
64.96386467
-2.44
23.43221003
-46.16811407
64.40351195
-2.43
22.97373975
-45.52687117
63.84545106
-2.42
22.52165404
-44.89119741
63.28968199
-2.41
22.07589732
-44.26106989
62.73620476
-2.40
21.63641423
-43.63646568
62.18501936
-2.39
21.20314967
-43.01736186
61.6361258
-2.38
20.77604874
-42.40373552
61.08952406
-2.37
20.35505678
-41.79556374
60.54521416
-2.36
19.94011936
-41.1928236
60.00319608
-2.35
19.53118227
-40.59549218
59.46346984
-2.34
19.12819156
-40.00354656
58.92603543
-2.33
18.73109347
-39.41696383
58.39089285
-2.32
18.33983448
-38.83572107
57.8580421
-2.31
17.95436132
-38.25979535
57.32748318
-2.30
17.57462093
-37.68916376
56.79921609
-2.29
17.20056048
-37.12380339
56.27324083
-2.28
16.83212737
-36.56369131
55.74955741
-2.27
16.46926923
-36.0088046
55.22816582
-2.26
16.11193393
-35.45912035
54.70906605
-2.25
15.76006956
-34.91461564
54.19225812
-2.24
15.41362443
-34.37526755
53.67774202
-2.23
15.0725471
-33.84105316
53.16551775
-2.22
14.73678633
-33.31194956
52.65558532
-2.21
14.40629115
-32.78793382
52.14794471
-2.20
14.08101077
-32.26898302
51.64259593
-2.19
13.76089468
-31.75507426
51.13953899
-2.18
13.44589256
-31.2461846
50.63877388
-2.17
13.13595433
-30.74229114
50.1403006
-2.16
12.83103016
-30.24337095
49.64411915
-2.15
12.53107041
-29.74940112
49.15022953
-2.14
12.23602571
-29.26035872
48.65863174
-2.13
11.94584689
-28.77622084
48.16932578
-2.12
11.66048502
-28.29696457
47.68231165
-2.11
11.3798914
-27.82256697
47.19758936
-2.10
11.10401756
-27.35300514
46.7151589
-2.09
10.83281526
-26.88825615
46.23502026
-2.08
10.56623647
-26.42829709
45.75717346
-2.07
10.30423343
-25.97310504
45.28161849
-2.06
10.04675856
-25.52265709
44.80835535
-2.05
9.793764546
-25.0769303
44.33738405
-2.04
9.545204292
-24.63590177
43.86870457
-2.03
9.301030927
-24.19954857
43.40231692
-2.02
9.061197812
-23.76784779
42.93822111
-2.01
8.825658539
-23.34077651
42.47641713
-2.00
8.594366927
-22.91831181
42.01690498
-1.99
8.367277024
-22.50043077
41.55968466
-1.98
8.144343109
-22.08711047
41.10475617
-1.97
7.925519689
-21.678328
40.65211951
-1.96
7.710761499
-21.27406044
40.20177468
-1.95
7.500023507
-20.87428487
39.75372169
-1.94
7.293260905
-20.47897837
39.30796052
-1.93
7.090429119
-20.08811802
38.86449119
-1.92
6.891483801
-19.70168091
38.42331369
-1.91
6.696380833
-19.31964411
37.98442801
-1.90
6.505076327
-18.94198471
37.54783417
-1.89
6.317526624
-18.56867979
37.11353217
-1.88
6.133688293
-18.19970642
36.68152199
-1.87
5.953518133
-17.83504171
36.25180364
-1.86
5.776973173
-17.47466271
35.82437713
-1.85
5.60401067
-17.11854652
35.39924244
-1.84
5.434588109
-16.76667022
34.97639959
-1.83
5.268663208
-16.41901089
34.55584857
-1.82
5.106193911
-16.07554561
34.13758938
-1.81
4.947138392
-15.73625147
33.72162202
-1.80
4.791455055
-15.40110553
33.30794649
-1.79
4.639102531
-15.0700849
32.89656279
-1.78
4.490039682
-14.74316664
32.48747093
-1.77
4.3442256
-14.42032784
32.08067089
-1.76
4.201619604
-14.10154558
31.67616269
-1.75
4.062181243
-13.78679695
31.27394632
-1.74
3.925870296
-13.47605901
30.87402178
-1.73
3.79264677
-13.16930887
30.47638907
-1.72
3.662470902
-12.86652359
30.08104819
-1.71
3.535303158
-12.56768027
29.68799914
-1.70
3.411104233
-12.27275597
29.29724192
-1.69
3.289835052
-11.98172779
28.90877654
-1.68
3.171456767
-11.6945728
28.52260299
-1.67
3.055930761
-11.41126809
28.13872126
-1.66
2.943218647
-11.13179074
27.75713137
-1.65
2.833282265
-10.85611782
27.37783331
-1.64
2.726083686
-10.58422643
27.00082708
-1.63
2.621585208
-10.31609364
26.62611268
-1.62
2.519749361
-10.05169654
26.25369012
-1.61
2.420538901
-9.7910122
25.88355938
-1.60
2.323916817
-9.534017711
25.51572048
-1.59
2.229846324
-9.280690151
25.1501734
-1.58
2.138290867
-9.031006603
24.78691816
-1.57
2.049214122
-8.784944149
24.42595475
-1.56
1.96257999
-8.542479869
24.06728317
-1.55
1.878352607
-8.303590846
23.71090342
-1.54
1.796496332
-8.068254161
23.3568155
-1.53
1.716975759
-7.836446896
23.00501942
-1.52
1.639755706
-7.608146133
22.65551516
-1.51
1.564801224
-7.383328954
22.30830274
-1.50
1.492077591
-7.161972439
21.96338215
-1.49
1.421550316
-6.944053671
21.62075339
-1.48
1.353185135
-6.729549732
21.28041645
-1.47
1.286948015
-6.518437703
20.94237136
-1.46
1.222805151
-6.310694665
20.60661809
-1.45
1.160722968
-6.106297702
20.27315665
-1.44
1.10066812
-5.905223893
19.94198705
-1.43
1.04260749
-5.707450321
19.61310927
-1.42
0.986508189
-5.512954068
19.28652333
-1.41
0.93233756
-5.321712215
18.96222922
-1.40
0.880063173
-5.133701844
18.64022693
-1.39
0.829652828
-4.948900037
18.32051649
-1.38
0.781074554
-4.767283875
18.00309787
-1.37
0.734296608
-4.58883044
17.68797108
-1.36
0.689287479
-4.413516814
17.37513612
-1.35
0.646015882
-4.241320078
17.064593
-1.34
0.604450764
-4.072217315
16.7563417
-1.33
0.564561299
-3.906185605
16.45038224
-1.32
0.526316892
-3.743202031
16.14671461
-1.31
0.489687175
-3.583243673
15.84533881
-1.30
0.45464201
-3.426287615
15.54625484
-1.29
0.421151491
-3.272310937
15.2494627
-1.28
0.389185937
-3.121290721
14.9549624
-1.27
0.358715898
-2.97320405
14.66275392
-1.26
0.329712154
-2.828028004
14.37283728
-1.25
0.302145712
-2.685739665
14.08521246
-1.24
0.275987811
-2.546316115
13.79987948
-1.23
0.251209917
-2.409734436
13.51683833
-1.22
0.227783726
-2.275971709
13.23608901
-1.21
0.205681163
-2.145005016
12.95763152
-1.20
0.184874382
-2.016811439
12.68146587
-1.19
0.165335767
-1.891368059
12.40759204
-1.18
0.14703793
-1.768651959
12.13601004
-1.17
0.129953713
-1.648640219
11.86671988
-1.16
0.114056187
-1.531309922
11.59972155
-1.15
0.099318653
-1.416638148
11.33501505
-1.14
0.085714639
-1.304601981
11.07260038
-1.13
0.073217904
-1.195178501
10.81247754
-1.12
0.061802436
-1.088344791
10.55464653
-1.11
0.051442452
-0.984077931
10.29910735
-1.10
0.042112398
-0.882355005
10.04586001
-1.09
0.033786949
-0.783153092
9.794904494
-1.08
0.026441009
-0.686449275
9.546240811
-1.07
0.020049713
-0.592220636
9.299868959
-1.06
0.014588422
-0.500444257
9.055788938
-1.05
0.01003273
-0.411097218
8.814000748
-1.04
0.006358456
-0.324156602
8.57450439
-1.03
0.003541653
-0.239599491
8.337299863
-1.02
0.001558598
-0.157402965
8.102387167
-1.01
0.000385801
-0.077544108
7.869766302
-1.00
0
-7.4607E-14
7.639437268
-0.99
0.000378162
0.075252277
7.411400066
-0.98
0.001497482
0.148235641
7.185654695
-0.97
0.003335388
0.21897301
6.962201155
-0.96
0.005869532
0.287487303
6.741039446
-0.95
0.0090778
0.353801438
6.522169568
-0.94
0.012938304
0.417938334
6.305591521
-0.93
0.017429386
0.479920908
6.091305306
-0.92
0.022529617
0.53977208
5.879310922
-0.91
0.028217799
0.597514766
5.669608369
-0.90
0.034472961
0.653171886
5.462197647
-0.89
0.041274361
0.706766359
5.257078756
-0.88
0.048601489
0.758321101
5.054251697
-0.87
0.056434061
0.807859032
4.853716468
-0.86
0.064752023
0.85540307
4.655473071
-0.85
0.073535552
0.900976133
4.459521505
-0.84
0.082765052
0.944601139
4.265861771
-0.83
0.092421158
0.986301008
4.074493867
-0.82
0.102484732
1.026098656
3.885417795
-0.81
0.112936866
1.064017003
3.698633554
-0.80
0.123758884
1.100078967
3.514141143
-0.79
0.134932334
1.134307465
3.331940565
-0.78
0.146438998
1.166725417
3.152031817
-0.77
0.158260884
1.197355741
2.9744149
-0.76
0.17038023
1.226221355
2.799089815
-0.75
0.182779505
1.253345177
2.626056561
-0.74
0.195441404
1.278750125
2.455315138
-0.73
0.208348854
1.302459119
2.286865546
-0.72
0.22148501
1.324495076
2.120707786
-0.71
0.234833255
1.344880914
1.956841856
-0.70
0.248377204
1.363639552
1.795267758
-0.69
0.262100699
1.380793909
1.635985491
-0.68
0.275987811
1.396366902
1.478995055
-0.67
0.290022842
1.410381449
1.32429645
-0.66
0.304190322
1.42286047
1.171889677
-0.65
0.318475009
1.433826882
1.021774735
-0.64
0.332861894
1.443303604
0.873951624
-0.63
0.347336192
1.451313554
0.728420344
-0.62
0.361883352
1.457879651
0.585180895
-0.61
0.376489049
1.463024812
0.444233277
-0.60
0.391139188
1.466771956
0.305577491
-0.59
0.405819904
1.469144001
0.169213535
-0.58
0.420517561
1.470163866
0.035141411
-0.57
0.435218751
1.469854468
-0.096638881
-0.56
0.449910296
1.468238727
-0.226127343
-0.55
0.464579247
1.465339561
-0.353323974
-0.54
0.479212885
1.461179887
-0.478228773
-0.53
0.49379872
1.455782625
-0.600841741
-0.52
0.508324489
1.449170692
-0.721162878
-0.51
0.522778161
1.441367007
-0.839192184
-0.50
0.537147933
1.432394488
-0.954929659
-0.49
0.551422231
1.422276053
-1.068375302
-0.48
0.565589711
1.411034621
-1.179529114
-0.47
0.579639257
1.39869311
-1.288391095
-0.46
0.593559982
1.385274439
-1.394961245
-0.45
0.607341231
1.370801525
-1.499239564
-0.44
0.620972575
1.355297287
-1.601226051
-0.43
0.634443816
1.338784643
-1.700920708
-0.42
0.647744983
1.321286512
-1.798323533
-0.41
0.660866337
1.302825812
-1.893434527
-0.40
0.673798367
1.283425461
-1.98625369
-0.39
0.686531791
1.263108378
-2.076781021
-0.38
0.699057555
1.24189748
-2.165016522
-0.37
0.711366837
1.219815687
-2.250960191
-0.36
0.723451042
1.196885916
-2.334612029
-0.35
0.735301805
1.173131086
-2.415972036
-0.34
0.74691099
1.148574114
-2.495040212
-0.33
0.75827069
1.123237921
-2.571816556
-0.32
0.769373228
1.097145423
-2.64630107
-0.31
0.780211154
1.070319539
-2.718493752
-0.30
0.79077725
1.042783187
-2.788394603
-0.29
0.801064526
1.014559286
-2.856003623
-0.28
0.811066221
0.985670754
-2.921320811
-0.27
0.820775802
0.956140509
-2.984346169
-0.26
0.830186968
0.92599147
-3.045079695
-0.25
0.839293645
0.895246555
-3.10352139
-0.24
0.848089989
0.863928682
-3.159671254
-0.23
0.856570385
0.832060769
-3.213529287
-0.22
0.864729448
0.799665736
-3.265095489
-0.21
0.872562019
0.766766499
-3.314369859
-0.20
0.880063173
0.733385978
-3.361352398
-0.19
0.887228211
0.69954709
-3.406043106
-0.18
0.894052664
0.665272755
-3.448441983
-0.17
0.900532291
0.63058589
-3.488549029
-0.16
0.906663083
0.595509414
-3.526364243
-0.15
0.912441257
0.560066245
-3.561887626
-0.14
0.917863262
0.524279301
-3.595119179
-0.13
0.922925774
0.488171501
-3.626058899
-0.12
0.927625699
0.451765762
-3.654706789
-0.11
0.931960172
0.415085004
-3.681062848
-0.10
0.935926558
0.378152145
-3.705127075
-0.09
0.939522451
0.340990102
-3.726899471
-0.08
0.942745673
0.303621795
-3.746380036
-0.07
0.945594276
0.266070141
-3.76356877
-0.06
0.948066541
0.228358059
-3.778465673
-0.05
0.950160979
0.190508467
-3.791070744
-0.04
0.951876328
0.152544283
-3.801383985
-0.03
0.953211559
0.114488427
-3.809405394
-0.02
0.954165868
0.076363815
-3.815134972
-0.01
0.954738682
0.038193367
-3.818572719
0.00
0.954929659
3.44428E-14
-3.819718634
0.01
0.954738682
-0.038193367
-3.818572719
0.02
0.954165868
-0.076363815
-3.815134972
0.03
0.953211559
-0.114488427
-3.809405394
0.04
0.951876328
-0.152544283
-3.801383985
0.05
0.950160979
-0.190508467
-3.791070744
0.06
0.948066541
-0.228358059
-3.778465673
0.07
0.945594276
-0.266070141
-3.76356877
0.08
0.942745673
-0.303621795
-3.746380036
0.09
0.939522451
-0.340990102
-3.726899471
0.10
0.935926558
-0.378152145
-3.705127075
0.11
0.931960172
-0.415085004
-3.681062848
0.12
0.927625699
-0.451765762
-3.654706789
0.13
0.922925774
-0.488171501
-3.626058899
0.14
0.917863262
-0.524279301
-3.595119179
0.15
0.912441257
-0.560066245
-3.561887626
0.16
0.906663083
-0.595509414
-3.526364243
0.17
0.900532291
-0.63058589
-3.488549029
0.18
0.894052664
-0.665272755
-3.448441983
0.19
0.887228211
-0.69954709
-3.406043106
0.20
0.880063173
-0.733385978
-3.361352398
0.21
0.872562019
-0.766766499
-3.314369859
0.22
0.864729448
-0.799665736
-3.265095489
0.23
0.856570385
-0.832060769
-3.213529287
0.24
0.848089989
-0.863928682
-3.159671254
0.25
0.839293645
-0.895246555
-3.10352139
0.26
0.830186968
-0.92599147
-3.045079695
0.27
0.820775802
-0.956140509
-2.984346169
0.28
0.811066221
-0.985670754
-2.921320811
0.29
0.801064526
-1.014559286
-2.856003623
0.30
0.79077725
-1.042783187
-2.788394603
0.31
0.780211154
-1.070319539
-2.718493752
0.32
0.769373228
-1.097145423
-2.64630107
0.33
0.75827069
-1.123237921
-2.571816556
0.34
0.74691099
-1.148574114
-2.495040212
0.35
0.735301805
-1.173131086
-2.415972036
0.36
0.723451042
-1.196885916
-2.334612029
0.37
0.711366837
-1.219815687
-2.250960191
0.38
0.699057555
-1.24189748
-2.165016522
0.39
0.686531791
-1.263108378
-2.076781021
0.40
0.673798367
-1.283425461
-1.98625369
0.41
0.660866337
-1.302825812
-1.893434527
0.42
0.647744983
-1.321286512
-1.798323533
0.43
0.634443816
-1.338784643
-1.700920708
0.44
0.620972575
-1.355297287
-1.601226051
0.45
0.607341231
-1.370801525
-1.499239564
0.46
0.593559982
-1.385274439
-1.394961245
0.47
0.579639257
-1.39869311
-1.288391095
0.48
0.565589711
-1.411034621
-1.179529114
0.49
0.551422231
-1.422276053
-1.068375302
0.50
0.537147933
-1.432394488
-0.954929659
0.51
0.522778161
-1.441367007
-0.839192184
0.52
0.508324489
-1.449170692
-0.721162878
0.53
0.49379872
-1.455782625
-0.600841741
0.54
0.479212885
-1.461179887
-0.478228773
0.55
0.464579247
-1.465339561
-0.353323974
0.56
0.449910296
-1.468238727
-0.226127343
0.57
0.435218751
-1.469854468
-0.096638881
0.58
0.420517561
-1.470163866
0.035141411
0.59
0.405819904
-1.469144001
0.169213535
0.60
0.391139188
-1.466771956
0.305577491
0.61
0.376489049
-1.463024812
0.444233277
0.62
0.361883352
-1.457879651
0.585180895
0.63
0.347336192
-1.451313554
0.728420344
0.64
0.332861894
-1.443303604
0.873951624
0.65
0.318475009
-1.433826882
1.021774735
0.66
0.304190322
-1.42286047
1.171889677
0.67
0.290022842
-1.410381449
1.32429645
0.68
0.275987811
-1.396366902
1.478995055
0.69
0.262100699
-1.380793909
1.635985491
0.70
0.248377204
-1.363639552
1.795267758
0.71
0.234833255
-1.344880914
1.956841856
0.72
0.22148501
-1.324495076
2.120707786
0.73
0.208348854
-1.302459119
2.286865546
0.74
0.195441404
-1.278750125
2.455315138
0.75
0.182779505
-1.253345177
2.626056561
0.76
0.17038023
-1.226221355
2.799089815
0.77
0.158260884
-1.197355741
2.9744149
0.78
0.146438998
-1.166725417
3.152031817
0.79
0.134932334
-1.134307465
3.331940565
0.80
0.123758884
-1.100078967
3.514141143
0.81
0.112936866
-1.064017003
3.698633554
0.82
0.102484732
-1.026098656
3.885417795
0.83
0.092421158
-0.986301008
4.074493867
0.84
0.082765052
-0.944601139
4.265861771
0.85
0.073535552
-0.900976133
4.459521505
0.86
0.064752023
-0.85540307
4.655473071
0.87
0.056434061
-0.807859032
4.853716468
0.88
0.048601489
-0.758321101
5.054251697
0.89
0.041274361
-0.706766359
5.257078756
0.90
0.034472961
-0.653171886
5.462197647
0.91
0.028217799
-0.597514766
5.669608369
0.92
0.022529617
-0.53977208
5.879310922
0.93
0.017429386
-0.479920908
6.091305306
0.94
0.012938304
-0.417938334
6.305591521
0.95
0.0090778
-0.353801438
6.522169568
0.96
0.005869532
-0.287487303
6.741039446
0.97
0.003335388
-0.21897301
6.962201155
0.98
0.001497482
-0.148235641
7.185654695
0.99
0.000378162
-0.075252277
7.411400066
1.00
0
-6.35048E-14
7.639437268
1.01
0.000385801
0.077544108
7.869766302
1.02
0.001558598
0.157402965
8.102387167
1.03
0.003541653
0.239599491
8.337299863
1.04
0.006358456
0.324156602
8.57450439
1.05
0.01003273
0.411097218
8.814000748
1.06
0.014588422
0.500444257
9.055788938
1.07
0.020049713
0.592220636
9.299868959
1.08
0.026441009
0.686449275
9.546240811
1.09
0.033786949
0.783153092
9.794904494
1.10
0.042112398
0.882355005
10.04586001
1.11
0.051442452
0.984077931
10.29910735
1.12
0.061802436
1.088344791
10.55464653
1.13
0.073217904
1.195178501
10.81247754
1.14
0.085714639
1.304601981
11.07260038
1.15
0.099318653
1.416638148
11.33501505
1.16
0.114056187
1.531309922
11.59972155
1.17
0.129953713
1.648640219
11.86671988
1.18
0.14703793
1.768651959
12.13601004
1.19
0.165335767
1.891368059
12.40759204
1.20
0.184874382
2.016811439
12.68146587
1.21
0.205681163
2.145005016
12.95763152
1.22
0.227783726
2.275971709
13.23608901
1.23
0.251209917
2.409734436
13.51683833
1.24
0.275987811
2.546316115
13.79987948
1.25
0.302145712
2.685739665
14.08521246
1.26
0.329712154
2.828028004
14.37283728
1.27
0.358715898
2.97320405
14.66275392
1.28
0.389185937
3.121290721
14.9549624
1.29
0.421151491
3.272310937
15.2494627
1.30
0.45464201
3.426287615
15.54625484
1.31
0.489687175
3.583243673
15.84533881
1.32
0.526316892
3.743202031
16.14671461
1.33
0.564561299
3.906185605
16.45038224
1.34
0.604450764
4.072217315
16.7563417
1.35
0.646015882
4.241320078
17.064593
1.36
0.689287479
4.413516814
17.37513612
1.37
0.734296608
4.58883044
17.68797108
1.38
0.781074554
4.767283875
18.00309787
1.39
0.829652828
4.948900037
18.32051649
1.40
0.880063173
5.133701844
18.64022693
1.41
0.93233756
5.321712215
18.96222922
1.42
0.986508189
5.512954068
19.28652333
1.43
1.04260749
5.707450321
19.61310927
1.44
1.10066812
5.905223893
19.94198705
1.45
1.160722968
6.106297702
20.27315665
1.46
1.222805151
6.310694665
20.60661809
1.47
1.286948015
6.518437703
20.94237136
1.48
1.353185135
6.729549732
21.28041645
1.49
1.421550316
6.944053671
21.62075339
1.50
1.492077591
7.161972439
21.96338215
1.51
1.564801224
7.383328954
22.30830274
1.52
1.639755706
7.608146133
22.65551516
1.53
1.716975759
7.836446896
23.00501942
1.54
1.796496332
8.068254161
23.3568155
1.55
1.878352607
8.303590846
23.71090342
1.56
1.96257999
8.542479869
24.06728317
1.57
2.049214122
8.784944149
24.42595475
1.58
2.138290867
9.031006603
24.78691816
1.59
2.229846324
9.280690151
25.1501734
1.60
2.323916817
9.534017711
25.51572048
1.61
2.420538901
9.7910122
25.88355938
1.62
2.519749361
10.05169654
26.25369012
1.63
2.621585208
10.31609364
26.62611268
1.64
2.726083686
10.58422643
27.00082708
1.65
2.833282265
10.85611782
27.37783331
1.66
2.943218647
11.13179074
27.75713137
1.67
3.055930761
11.41126809
28.13872126
1.68
3.171456767
11.6945728
28.52260299
1.69
3.289835052
11.98172779
28.90877654
1.70
3.411104233
12.27275597
29.29724192
1.71
3.535303158
12.56768027
29.68799914
1.72
3.662470902
12.86652359
30.08104819
1.73
3.79264677
13.16930887
30.47638907
1.74
3.925870296
13.47605901
30.87402178
1.75
4.062181243
13.78679695
31.27394632
1.76
4.201619604
14.10154558
31.67616269
1.77
4.3442256
14.42032784
32.08067089
1.78
4.490039682
14.74316664
32.48747093
1.79
4.639102531
15.0700849
32.89656279
1.80
4.791455055
15.40110553
33.30794649
1.81
4.947138392
15.73625147
33.72162202
1.82
5.106193911
16.07554561
34.13758938
1.83
5.268663208
16.41901089
34.55584857
1.84
5.434588109
16.76667022
34.97639959
1.85
5.60401067
17.11854652
35.39924244
1.86
5.776973173
17.47466271
35.82437713
1.87
5.953518133
17.83504171
36.25180364
1.88
6.133688293
18.19970642
36.68152199
1.89
6.317526624
18.56867979
37.11353217
1.90
6.505076327
18.94198471
37.54783417
1.91
6.696380833
19.31964411
37.98442801
1.92
6.891483801
19.70168091
38.42331369
1.93
7.090429119
20.08811802
38.86449119
1.94
7.293260905
20.47897837
39.30796052
1.95
7.500023507
20.87428487
39.75372169
1.96
7.710761499
21.27406044
40.20177468
1.97
7.925519689
21.678328
40.65211951
1.98
8.144343109
22.08711047
41.10475617
1.99
8.367277024
22.50043077
41.55968466
2.00
8.594366927
22.91831181
42.01690498
2.01
8.825658539
23.34077651
42.47641713
2.02
9.061197812
23.76784779
42.93822111
2.03
9.301030927
24.19954857
43.40231692
2.04
9.545204292
24.63590177
43.86870457
2.05
9.793764546
25.0769303
44.33738405
2.06
10.04675856
25.52265709
44.80835535
2.07
10.30423343
25.97310504
45.28161849
2.08
10.56623647
26.42829709
45.75717346
2.09
10.83281526
26.88825615
46.23502026
2.10
11.10401756
27.35300514
46.7151589
2.11
11.3798914
27.82256697
47.19758936
2.12
11.66048502
28.29696457
47.68231165
2.13
11.94584689
28.77622084
48.16932578
2.14
12.23602571
29.26035872
48.65863174
2.15
12.53107041
29.74940112
49.15022953
2.16
12.83103016
30.24337095
49.64411915
2.17
13.13595433
30.74229114
50.1403006
2.18
13.44589256
31.2461846
50.63877388
2.19
13.76089468
31.75507426
51.13953899
2.20
14.08101077
32.26898302
51.64259593
2.21
14.40629115
32.78793382
52.14794471
2.22
14.73678633
33.31194956
52.65558532
2.23
15.0725471
33.84105316
53.16551775
2.24
15.41362443
34.37526755
53.67774202
2.25
15.76006956
34.91461564
54.19225812
2.26
16.11193393
35.45912035
54.70906605
2.27
16.46926923
36.0088046
55.22816582
2.28
16.83212737
36.56369131
55.74955741
2.29
17.20056048
37.12380339
56.27324083
2.30
17.57462093
37.68916376
56.79921609
2.31
17.95436132
38.25979535
57.32748318
2.32
18.33983448
38.83572107
57.8580421
2.33
18.73109347
39.41696383
58.39089285
2.34
19.12819156
40.00354656
58.92603543
2.35
19.53118227
40.59549218
59.46346984
2.36
19.94011936
41.1928236
60.00319608
2.37
20.35505678
41.79556374
60.54521416
2.38
20.77604874
42.40373552
61.08952406
2.39
21.20314967
43.01736186
61.6361258
2.40
21.63641423
43.63646568
62.18501936
2.41
22.07589732
44.26106989
62.73620476
2.42
22.52165404
44.89119741
63.28968199
2.43
22.97373975
45.52687117
63.84545106
2.44
23.43221003
46.16811407
64.40351195
2.45
23.89712067
46.81494905
64.96386467
2.46
24.36852773
47.46739901
65.52650923
2.47
24.84648745
48.12548687
66.09144561
2.48
25.33105633
48.78923556
66.65867383
2.49
25.8222911
49.45866799
67.22819388
2.50
26.32024871
50.13380707
67.80000576
Now that the polynomial kernel density function is well described mathematically, we focus on the relevant sub-domain of the function , which is , and not . Indeed, the ratio
has no purpose below zero and above 1 in the context of a density estimation exercise, since
cannot be below zero (that is, a negative distance from the spatial point of interest is irrelevant) and cannot be beyond
(that is, reaching out above the circular area of interest is also irrelevant). Consequently, the relevant image of
is , because
and
and is a decreasing function between 0 and 1. More importantly, the distinct contribution of a job location
within the radius of interest
is approximately normalized between 0 and 1 (that is, between values 0 and
and
is a finite sum of approximately normalized terms, since the number of job locations within the circular and symmetric surrounding
is finite. Figure 4 is the same as the previous Figure 3. However, it focuses on the relevant domain and image of
and therefore, enable us to better visualize the value of
when , which is not 1.0, but 0.95. The contribution of job location
decreases as
increases between . The decrease is initially less drastic around zero, but eventually becomes sharper and then smoother again just before reaching a ratio
of 1. This structure is due to the second derivative (green curve) reaching a value of zero when
reaches +0.577. A square highlighted yellow bound the region of interest (Bartle and Sherbert, 2011).
Description for Figure 4
Figure 4 Polynomial kernel density function (zoom in at domain [0,1] for f(x), f'(x) and f''(x))
Table summary The information is grouped by Distance Between Cell Centroid and
Job Location to Bandwidth Size Ratio (x) (appearing as row headers), , calculated using (appearing as column headers).
Distance Between Cell Centroid and
Job Location to Bandwidth Size Ratio (x)
Original Function (blue)
f(x)
First Derivative (red)
f'(x)
Second Derivative (green)
f''(x)
Source : authors’ computations from the theoretical kernel density estimation function of Parzen (1962) and Rosenblatt (1956).
Legend: Blue = function of interest
, red = first derivative of
, and green = second derivative of
-2.50
26.32024871
-50.13380707
67.80000576
-2.49
25.8222911
-49.45866799
67.22819388
-2.48
25.33105633
-48.78923556
66.65867383
-2.47
24.84648745
-48.12548687
66.09144561
-2.46
24.36852773
-47.46739901
65.52650923
-2.45
23.89712067
-46.81494905
64.96386467
-2.44
23.43221003
-46.16811407
64.40351195
-2.43
22.97373975
-45.52687117
63.84545106
-2.42
22.52165404
-44.89119741
63.28968199
-2.41
22.07589732
-44.26106989
62.73620476
-2.40
21.63641423
-43.63646568
62.18501936
-2.39
21.20314967
-43.01736186
61.6361258
-2.38
20.77604874
-42.40373552
61.08952406
-2.37
20.35505678
-41.79556374
60.54521416
-2.36
19.94011936
-41.1928236
60.00319608
-2.35
19.53118227
-40.59549218
59.46346984
-2.34
19.12819156
-40.00354656
58.92603543
-2.33
18.73109347
-39.41696383
58.39089285
-2.32
18.33983448
-38.83572107
57.8580421
-2.31
17.95436132
-38.25979535
57.32748318
-2.30
17.57462093
-37.68916376
56.79921609
-2.29
17.20056048
-37.12380339
56.27324083
-2.28
16.83212737
-36.56369131
55.74955741
-2.27
16.46926923
-36.0088046
55.22816582
-2.26
16.11193393
-35.45912035
54.70906605
-2.25
15.76006956
-34.91461564
54.19225812
-2.24
15.41362443
-34.37526755
53.67774202
-2.23
15.0725471
-33.84105316
53.16551775
-2.22
14.73678633
-33.31194956
52.65558532
-2.21
14.40629115
-32.78793382
52.14794471
-2.20
14.08101077
-32.26898302
51.64259593
-2.19
13.76089468
-31.75507426
51.13953899
-2.18
13.44589256
-31.2461846
50.63877388
-2.17
13.13595433
-30.74229114
50.1403006
-2.16
12.83103016
-30.24337095
49.64411915
-2.15
12.53107041
-29.74940112
49.15022953
-2.14
12.23602571
-29.26035872
48.65863174
-2.13
11.94584689
-28.77622084
48.16932578
-2.12
11.66048502
-28.29696457
47.68231165
-2.11
11.3798914
-27.82256697
47.19758936
-2.10
11.10401756
-27.35300514
46.7151589
-2.09
10.83281526
-26.88825615
46.23502026
-2.08
10.56623647
-26.42829709
45.75717346
-2.07
10.30423343
-25.97310504
45.28161849
-2.06
10.04675856
-25.52265709
44.80835535
-2.05
9.793764546
-25.0769303
44.33738405
-2.04
9.545204292
-24.63590177
43.86870457
-2.03
9.301030927
-24.19954857
43.40231692
-2.02
9.061197812
-23.76784779
42.93822111
-2.01
8.825658539
-23.34077651
42.47641713
-2.00
8.594366927
-22.91831181
42.01690498
-1.99
8.367277024
-22.50043077
41.55968466
-1.98
8.144343109
-22.08711047
41.10475617
-1.97
7.925519689
-21.678328
40.65211951
-1.96
7.710761499
-21.27406044
40.20177468
-1.95
7.500023507
-20.87428487
39.75372169
-1.94
7.293260905
-20.47897837
39.30796052
-1.93
7.090429119
-20.08811802
38.86449119
-1.92
6.891483801
-19.70168091
38.42331369
-1.91
6.696380833
-19.31964411
37.98442801
-1.90
6.505076327
-18.94198471
37.54783417
-1.89
6.317526624
-18.56867979
37.11353217
-1.88
6.133688293
-18.19970642
36.68152199
-1.87
5.953518133
-17.83504171
36.25180364
-1.86
5.776973173
-17.47466271
35.82437713
-1.85
5.60401067
-17.11854652
35.39924244
-1.84
5.434588109
-16.76667022
34.97639959
-1.83
5.268663208
-16.41901089
34.55584857
-1.82
5.106193911
-16.07554561
34.13758938
-1.81
4.947138392
-15.73625147
33.72162202
-1.80
4.791455055
-15.40110553
33.30794649
-1.79
4.639102531
-15.0700849
32.89656279
-1.78
4.490039682
-14.74316664
32.48747093
-1.77
4.3442256
-14.42032784
32.08067089
-1.76
4.201619604
-14.10154558
31.67616269
-1.75
4.062181243
-13.78679695
31.27394632
-1.74
3.925870296
-13.47605901
30.87402178
-1.73
3.79264677
-13.16930887
30.47638907
-1.72
3.662470902
-12.86652359
30.08104819
-1.71
3.535303158
-12.56768027
29.68799914
-1.70
3.411104233
-12.27275597
29.29724192
-1.69
3.289835052
-11.98172779
28.90877654
-1.68
3.171456767
-11.6945728
28.52260299
-1.67
3.055930761
-11.41126809
28.13872126
-1.66
2.943218647
-11.13179074
27.75713137
-1.65
2.833282265
-10.85611782
27.37783331
-1.64
2.726083686
-10.58422643
27.00082708
-1.63
2.621585208
-10.31609364
26.62611268
-1.62
2.519749361
-10.05169654
26.25369012
-1.61
2.420538901
-9.7910122
25.88355938
-1.60
2.323916817
-9.534017711
25.51572048
-1.59
2.229846324
-9.280690151
25.1501734
-1.58
2.138290867
-9.031006603
24.78691816
-1.57
2.049214122
-8.784944149
24.42595475
-1.56
1.96257999
-8.542479869
24.06728317
-1.55
1.878352607
-8.303590846
23.71090342
-1.54
1.796496332
-8.068254161
23.3568155
-1.53
1.716975759
-7.836446896
23.00501942
-1.52
1.639755706
-7.608146133
22.65551516
-1.51
1.564801224
-7.383328954
22.30830274
-1.50
1.492077591
-7.161972439
21.96338215
-1.49
1.421550316
-6.944053671
21.62075339
-1.48
1.353185135
-6.729549732
21.28041645
-1.47
1.286948015
-6.518437703
20.94237136
-1.46
1.222805151
-6.310694665
20.60661809
-1.45
1.160722968
-6.106297702
20.27315665
-1.44
1.10066812
-5.905223893
19.94198705
-1.43
1.04260749
-5.707450321
19.61310927
-1.42
0.986508189
-5.512954068
19.28652333
-1.41
0.93233756
-5.321712215
18.96222922
-1.40
0.880063173
-5.133701844
18.64022693
-1.39
0.829652828
-4.948900037
18.32051649
-1.38
0.781074554
-4.767283875
18.00309787
-1.37
0.734296608
-4.58883044
17.68797108
-1.36
0.689287479
-4.413516814
17.37513612
-1.35
0.646015882
-4.241320078
17.064593
-1.34
0.604450764
-4.072217315
16.7563417
-1.33
0.564561299
-3.906185605
16.45038224
-1.32
0.526316892
-3.743202031
16.14671461
-1.31
0.489687175
-3.583243673
15.84533881
-1.30
0.45464201
-3.426287615
15.54625484
-1.29
0.421151491
-3.272310937
15.2494627
-1.28
0.389185937
-3.121290721
14.9549624
-1.27
0.358715898
-2.97320405
14.66275392
-1.26
0.329712154
-2.828028004
14.37283728
-1.25
0.302145712
-2.685739665
14.08521246
-1.24
0.275987811
-2.546316115
13.79987948
-1.23
0.251209917
-2.409734436
13.51683833
-1.22
0.227783726
-2.275971709
13.23608901
-1.21
0.205681163
-2.145005016
12.95763152
-1.20
0.184874382
-2.016811439
12.68146587
-1.19
0.165335767
-1.891368059
12.40759204
-1.18
0.14703793
-1.768651959
12.13601004
-1.17
0.129953713
-1.648640219
11.86671988
-1.16
0.114056187
-1.531309922
11.59972155
-1.15
0.099318653
-1.416638148
11.33501505
-1.14
0.085714639
-1.304601981
11.07260038
-1.13
0.073217904
-1.195178501
10.81247754
-1.12
0.061802436
-1.088344791
10.55464653
-1.11
0.051442452
-0.984077931
10.29910735
-1.10
0.042112398
-0.882355005
10.04586001
-1.09
0.033786949
-0.783153092
9.794904494
-1.08
0.026441009
-0.686449275
9.546240811
-1.07
0.020049713
-0.592220636
9.299868959
-1.06
0.014588422
-0.500444257
9.055788938
-1.05
0.01003273
-0.411097218
8.814000748
-1.04
0.006358456
-0.324156602
8.57450439
-1.03
0.003541653
-0.239599491
8.337299863
-1.02
0.001558598
-0.157402965
8.102387167
-1.01
0.000385801
-0.077544108
7.869766302
-1.00
0
-7.46E-14
7.639437268
-0.99
0.000378162
0.075252277
7.411400066
-0.98
0.001497482
0.148235641
7.185654695
-0.97
0.003335388
0.21897301
6.962201155
-0.96
0.005869532
0.287487303
6.741039446
-0.95
0.0090778
0.353801438
6.522169568
-0.94
0.012938304
0.417938334
6.305591521
-0.93
0.017429386
0.479920908
6.091305306
-0.92
0.022529617
0.53977208
5.879310922
-0.91
0.028217799
0.597514766
5.669608369
-0.90
0.034472961
0.653171886
5.462197647
-0.89
0.041274361
0.706766359
5.257078756
-0.88
0.048601489
0.758321101
5.054251697
-0.87
0.056434061
0.807859032
4.853716468
-0.86
0.064752023
0.85540307
4.655473071
-0.85
0.073535552
0.900976133
4.459521505
-0.84
0.082765052
0.944601139
4.265861771
-0.83
0.092421158
0.986301008
4.074493867
-0.82
0.102484732
1.026098656
3.885417795
-0.81
0.112936866
1.064017003
3.698633554
-0.80
0.123758884
1.100078967
3.514141143
-0.79
0.134932334
1.134307465
3.331940565
-0.78
0.146438998
1.166725417
3.152031817
-0.77
0.158260884
1.197355741
2.9744149
-0.76
0.17038023
1.226221355
2.799089815
-0.75
0.182779505
1.253345177
2.626056561
-0.74
0.195441404
1.278750125
2.455315138
-0.73
0.208348854
1.302459119
2.286865546
-0.72
0.22148501
1.324495076
2.120707786
-0.71
0.234833255
1.344880914
1.956841856
-0.70
0.248377204
1.363639552
1.795267758
-0.69
0.262100699
1.380793909
1.635985491
-0.68
0.275987811
1.396366902
1.478995055
-0.67
0.290022842
1.410381449
1.32429645
-0.66
0.304190322
1.42286047
1.171889677
-0.65
0.318475009
1.433826882
1.021774735
-0.64
0.332861894
1.443303604
0.873951624
-0.63
0.347336192
1.451313554
0.728420344
-0.62
0.361883352
1.457879651
0.585180895
-0.61
0.376489049
1.463024812
0.444233277
-0.60
0.391139188
1.466771956
0.305577491
-0.59
0.405819904
1.469144001
0.169213535
-0.58
0.420517561
1.470163866
0.035141411
-0.57
0.435218751
1.469854468
-0.096638881
-0.56
0.449910296
1.468238727
-0.226127343
-0.55
0.464579247
1.465339561
-0.353323974
-0.54
0.479212885
1.461179887
-0.478228773
-0.53
0.49379872
1.455782625
-0.600841741
-0.52
0.508324489
1.449170692
-0.721162878
-0.51
0.522778161
1.441367007
-0.839192184
-0.50
0.537147933
1.432394488
-0.954929659
-0.49
0.551422231
1.422276053
-1.068375302
-0.48
0.565589711
1.411034621
-1.179529114
-0.47
0.579639257
1.39869311
-1.288391095
-0.46
0.593559982
1.385274439
-1.394961245
-0.45
0.607341231
1.370801525
-1.499239564
-0.44
0.620972575
1.355297287
-1.601226051
-0.43
0.634443816
1.338784643
-1.700920708
-0.42
0.647744983
1.321286512
-1.798323533
-0.41
0.660866337
1.302825812
-1.893434527
-0.40
0.673798367
1.283425461
-1.98625369
-0.39
0.686531791
1.263108378
-2.076781021
-0.38
0.699057555
1.24189748
-2.165016522
-0.37
0.711366837
1.219815687
-2.250960191
-0.36
0.723451042
1.196885916
-2.334612029
-0.35
0.735301805
1.173131086
-2.415972036
-0.34
0.74691099
1.148574114
-2.495040212
-0.33
0.75827069
1.123237921
-2.571816556
-0.32
0.769373228
1.097145423
-2.64630107
-0.31
0.780211154
1.070319539
-2.718493752
-0.30
0.79077725
1.042783187
-2.788394603
-0.29
0.801064526
1.014559286
-2.856003623
-0.28
0.811066221
0.985670754
-2.921320811
-0.27
0.820775802
0.956140509
-2.984346169
-0.26
0.830186968
0.92599147
-3.045079695
-0.25
0.839293645
0.895246555
-3.10352139
-0.24
0.848089989
0.863928682
-3.159671254
-0.23
0.856570385
0.832060769
-3.213529287
-0.22
0.864729448
0.799665736
-3.265095489
-0.21
0.872562019
0.766766499
-3.314369859
-0.20
0.880063173
0.733385978
-3.361352398
-0.19
0.887228211
0.69954709
-3.406043106
-0.18
0.894052664
0.665272755
-3.448441983
-0.17
0.900532291
0.63058589
-3.488549029
-0.16
0.906663083
0.595509414
-3.526364243
-0.15
0.912441257
0.560066245
-3.561887626
-0.14
0.917863262
0.524279301
-3.595119179
-0.13
0.922925774
0.488171501
-3.626058899
-0.12
0.927625699
0.451765762
-3.654706789
-0.11
0.931960172
0.415085004
-3.681062848
-0.10
0.935926558
0.378152145
-3.705127075
-0.09
0.939522451
0.340990102
-3.726899471
-0.08
0.942745673
0.303621795
-3.746380036
-0.07
0.945594276
0.266070141
-3.76356877
-0.06
0.948066541
0.228358059
-3.778465673
-0.05
0.950160979
0.190508467
-3.791070744
-0.04
0.951876328
0.152544283
-3.801383985
-0.03
0.953211559
0.114488427
-3.809405394
-0.02
0.954165868
0.076363815
-3.815134972
-0.01
0.954738682
0.038193367
-3.818572719
0.00
0.954929659
3.44E-14
-3.819718634
0.01
0.954738682
-0.038193367
-3.818572719
0.02
0.954165868
-0.076363815
-3.815134972
0.03
0.953211559
-0.114488427
-3.809405394
0.04
0.951876328
-0.152544283
-3.801383985
0.05
0.950160979
-0.190508467
-3.791070744
0.06
0.948066541
-0.228358059
-3.778465673
0.07
0.945594276
-0.266070141
-3.76356877
0.08
0.942745673
-0.303621795
-3.746380036
0.09
0.939522451
-0.340990102
-3.726899471
0.10
0.935926558
-0.378152145
-3.705127075
0.11
0.931960172
-0.415085004
-3.681062848
0.12
0.927625699
-0.451765762
-3.654706789
0.13
0.922925774
-0.488171501
-3.626058899
0.14
0.917863262
-0.524279301
-3.595119179
0.15
0.912441257
-0.560066245
-3.561887626
0.16
0.906663083
-0.595509414
-3.526364243
0.17
0.900532291
-0.63058589
-3.488549029
0.18
0.894052664
-0.665272755
-3.448441983
0.19
0.887228211
-0.69954709
-3.406043106
0.20
0.880063173
-0.733385978
-3.361352398
0.21
0.872562019
-0.766766499
-3.314369859
0.22
0.864729448
-0.799665736
-3.265095489
0.23
0.856570385
-0.832060769
-3.213529287
0.24
0.848089989
-0.863928682
-3.159671254
0.25
0.839293645
-0.895246555
-3.10352139
0.26
0.830186968
-0.92599147
-3.045079695
0.27
0.820775802
-0.956140509
-2.984346169
0.28
0.811066221
-0.985670754
-2.921320811
0.29
0.801064526
-1.014559286
-2.856003623
0.30
0.79077725
-1.042783187
-2.788394603
0.31
0.780211154
-1.070319539
-2.718493752
0.32
0.769373228
-1.097145423
-2.64630107
0.33
0.75827069
-1.123237921
-2.571816556
0.34
0.74691099
-1.148574114
-2.495040212
0.35
0.735301805
-1.173131086
-2.415972036
0.36
0.723451042
-1.196885916
-2.334612029
0.37
0.711366837
-1.219815687
-2.250960191
0.38
0.699057555
-1.24189748
-2.165016522
0.39
0.686531791
-1.263108378
-2.076781021
0.40
0.673798367
-1.283425461
-1.98625369
0.41
0.660866337
-1.302825812
-1.893434527
0.42
0.647744983
-1.321286512
-1.798323533
0.43
0.634443816
-1.338784643
-1.700920708
0.44
0.620972575
-1.355297287
-1.601226051
0.45
0.607341231
-1.370801525
-1.499239564
0.46
0.593559982
-1.385274439
-1.394961245
0.47
0.579639257
-1.39869311
-1.288391095
0.48
0.565589711
-1.411034621
-1.179529114
0.49
0.551422231
-1.422276053
-1.068375302
0.50
0.537147933
-1.432394488
-0.954929659
0.51
0.522778161
-1.441367007
-0.839192184
0.52
0.508324489
-1.449170692
-0.721162878
0.53
0.49379872
-1.455782625
-0.600841741
0.54
0.479212885
-1.461179887
-0.478228773
0.55
0.464579247
-1.465339561
-0.353323974
0.56
0.449910296
-1.468238727
-0.226127343
0.57
0.435218751
-1.469854468
-0.096638881
0.58
0.420517561
-1.470163866
0.035141411
0.59
0.405819904
-1.469144001
0.169213535
0.60
0.391139188
-1.466771956
0.305577491
0.61
0.376489049
-1.463024812
0.444233277
0.62
0.361883352
-1.457879651
0.585180895
0.63
0.347336192
-1.451313554
0.728420344
0.64
0.332861894
-1.443303604
0.873951624
0.65
0.318475009
-1.433826882
1.021774735
0.66
0.304190322
-1.42286047
1.171889677
0.67
0.290022842
-1.410381449
1.32429645
0.68
0.275987811
-1.396366902
1.478995055
0.69
0.262100699
-1.380793909
1.635985491
0.70
0.248377204
-1.363639552
1.795267758
0.71
0.234833255
-1.344880914
1.956841856
0.72
0.22148501
-1.324495076
2.120707786
0.73
0.208348854
-1.302459119
2.286865546
0.74
0.195441404
-1.278750125
2.455315138
0.75
0.182779505
-1.253345177
2.626056561
0.76
0.17038023
-1.226221355
2.799089815
0.77
0.158260884
-1.197355741
2.9744149
0.78
0.146438998
-1.166725417
3.152031817
0.79
0.134932334
-1.134307465
3.331940565
0.80
0.123758884
-1.100078967
3.514141143
0.81
0.112936866
-1.064017003
3.698633554
0.82
0.102484732
-1.026098656
3.885417795
0.83
0.092421158
-0.986301008
4.074493867
0.84
0.082765052
-0.944601139
4.265861771
0.85
0.073535552
-0.900976133
4.459521505
0.86
0.064752023
-0.85540307
4.655473071
0.87
0.056434061
-0.807859032
4.853716468
0.88
0.048601489
-0.758321101
5.054251697
0.89
0.041274361
-0.706766359
5.257078756
0.90
0.034472961
-0.653171886
5.462197647
0.91
0.028217799
-0.597514766
5.669608369
0.92
0.022529617
-0.53977208
5.879310922
0.93
0.017429386
-0.479920908
6.091305306
0.94
0.012938304
-0.417938334
6.305591521
0.95
0.0090778
-0.353801438
6.522169568
0.96
0.005869532
-0.287487303
6.741039446
0.97
0.003335388
-0.21897301
6.962201155
0.98
0.001497482
-0.148235641
7.185654695
0.99
0.000378162
-0.075252277
7.411400066
1.00
0
-6.35E-14
7.639437268
1.01
0.000385801
0.077544108
7.869766302
1.02
0.001558598
0.157402965
8.102387167
1.03
0.003541653
0.239599491
8.337299863
1.04
0.006358456
0.324156602
8.57450439
1.05
0.01003273
0.411097218
8.814000748
1.06
0.014588422
0.500444257
9.055788938
1.07
0.020049713
0.592220636
9.299868959
1.08
0.026441009
0.686449275
9.546240811
1.09
0.033786949
0.783153092
9.794904494
1.10
0.042112398
0.882355005
10.04586001
1.11
0.051442452
0.984077931
10.29910735
1.12
0.061802436
1.088344791
10.55464653
1.13
0.073217904
1.195178501
10.81247754
1.14
0.085714639
1.304601981
11.07260038
1.15
0.099318653
1.416638148
11.33501505
1.16
0.114056187
1.531309922
11.59972155
1.17
0.129953713
1.648640219
11.86671988
1.18
0.14703793
1.768651959
12.13601004
1.19
0.165335767
1.891368059
12.40759204
1.20
0.184874382
2.016811439
12.68146587
1.21
0.205681163
2.145005016
12.95763152
1.22
0.227783726
2.275971709
13.23608901
1.23
0.251209917
2.409734436
13.51683833
1.24
0.275987811
2.546316115
13.79987948
1.25
0.302145712
2.685739665
14.08521246
1.26
0.329712154
2.828028004
14.37283728
1.27
0.358715898
2.97320405
14.66275392
1.28
0.389185937
3.121290721
14.9549624
1.29
0.421151491
3.272310937
15.2494627
1.30
0.45464201
3.426287615
15.54625484
1.31
0.489687175
3.583243673
15.84533881
1.32
0.526316892
3.743202031
16.14671461
1.33
0.564561299
3.906185605
16.45038224
1.34
0.604450764
4.072217315
16.7563417
1.35
0.646015882
4.241320078
17.064593
1.36
0.689287479
4.413516814
17.37513612
1.37
0.734296608
4.58883044
17.68797108
1.38
0.781074554
4.767283875
18.00309787
1.39
0.829652828
4.948900037
18.32051649
1.40
0.880063173
5.133701844
18.64022693
1.41
0.93233756
5.321712215
18.96222922
1.42
0.986508189
5.512954068
19.28652333
1.43
1.04260749
5.707450321
19.61310927
1.44
1.10066812
5.905223893
19.94198705
1.45
1.160722968
6.106297702
20.27315665
1.46
1.222805151
6.310694665
20.60661809
1.47
1.286948015
6.518437703
20.94237136
1.48
1.353185135
6.729549732
21.28041645
1.49
1.421550316
6.944053671
21.62075339
1.50
1.492077591
7.161972439
21.96338215
1.51
1.564801224
7.383328954
22.30830274
1.52
1.639755706
7.608146133
22.65551516
1.53
1.716975759
7.836446896
23.00501942
1.54
1.796496332
8.068254161
23.3568155
1.55
1.878352607
8.303590846
23.71090342
1.56
1.96257999
8.542479869
24.06728317
1.57
2.049214122
8.784944149
24.42595475
1.58
2.138290867
9.031006603
24.78691816
1.59
2.229846324
9.280690151
25.1501734
1.60
2.323916817
9.534017711
25.51572048
1.61
2.420538901
9.7910122
25.88355938
1.62
2.519749361
10.05169654
26.25369012
1.63
2.621585208
10.31609364
26.62611268
1.64
2.726083686
10.58422643
27.00082708
1.65
2.833282265
10.85611782
27.37783331
1.66
2.943218647
11.13179074
27.75713137
1.67
3.055930761
11.41126809
28.13872126
1.68
3.171456767
11.6945728
28.52260299
1.69
3.289835052
11.98172779
28.90877654
1.70
3.411104233
12.27275597
29.29724192
1.71
3.535303158
12.56768027
29.68799914
1.72
3.662470902
12.86652359
30.08104819
1.73
3.79264677
13.16930887
30.47638907
1.74
3.925870296
13.47605901
30.87402178
1.75
4.062181243
13.78679695
31.27394632
1.76
4.201619604
14.10154558
31.67616269
1.77
4.3442256
14.42032784
32.08067089
1.78
4.490039682
14.74316664
32.48747093
1.79
4.639102531
15.0700849
32.89656279
1.80
4.791455055
15.40110553
33.30794649
1.81
4.947138392
15.73625147
33.72162202
1.82
5.106193911
16.07554561
34.13758938
1.83
5.268663208
16.41901089
34.55584857
1.84
5.434588109
16.76667022
34.97639959
1.85
5.60401067
17.11854652
35.39924244
1.86
5.776973173
17.47466271
35.82437713
1.87
5.953518133
17.83504171
36.25180364
1.88
6.133688293
18.19970642
36.68152199
1.89
6.317526624
18.56867979
37.11353217
1.90
6.505076327
18.94198471
37.54783417
1.91
6.696380833
19.31964411
37.98442801
1.92
6.891483801
19.70168091
38.42331369
1.93
7.090429119
20.08811802
38.86449119
1.94
7.293260905
20.47897837
39.30796052
1.95
7.500023507
20.87428487
39.75372169
1.96
7.710761499
21.27406044
40.20177468
1.97
7.925519689
21.678328
40.65211951
1.98
8.144343109
22.08711047
41.10475617
1.99
8.367277024
22.50043077
41.55968466
2.00
8.594366927
22.91831181
42.01690498
2.01
8.825658539
23.34077651
42.47641713
2.02
9.061197812
23.76784779
42.93822111
2.03
9.301030927
24.19954857
43.40231692
2.04
9.545204292
24.63590177
43.86870457
2.05
9.793764546
25.0769303
44.33738405
2.06
10.04675856
25.52265709
44.80835535
2.07
10.30423343
25.97310504
45.28161849
2.08
10.56623647
26.42829709
45.75717346
2.09
10.83281526
26.88825615
46.23502026
2.10
11.10401756
27.35300514
46.7151589
2.11
11.3798914
27.82256697
47.19758936
2.12
11.66048502
28.29696457
47.68231165
2.13
11.94584689
28.77622084
48.16932578
2.14
12.23602571
29.26035872
48.65863174
2.15
12.53107041
29.74940112
49.15022953
2.16
12.83103016
30.24337095
49.64411915
2.17
13.13595433
30.74229114
50.1403006
2.18
13.44589256
31.2461846
50.63877388
2.19
13.76089468
31.75507426
51.13953899
2.20
14.08101077
32.26898302
51.64259593
2.21
14.40629115
32.78793382
52.14794471
2.22
14.73678633
33.31194956
52.65558532
2.23
15.0725471
33.84105316
53.16551775
2.24
15.41362443
34.37526755
53.67774202
2.25
15.76006956
34.91461564
54.19225812
2.26
16.11193393
35.45912035
54.70906605
2.27
16.46926923
36.0088046
55.22816582
2.28
16.83212737
36.56369131
55.74955741
2.29
17.20056048
37.12380339
56.27324083
2.30
17.57462093
37.68916376
56.79921609
2.31
17.95436132
38.25979535
57.32748318
2.32
18.33983448
38.83572107
57.8580421
2.33
18.73109347
39.41696383
58.39089285
2.34
19.12819156
40.00354656
58.92603543
2.35
19.53118227
40.59549218
59.46346984
2.36
19.94011936
41.1928236
60.00319608
2.37
20.35505678
41.79556374
60.54521416
2.38
20.77604874
42.40373552
61.08952406
2.39
21.20314967
43.01736186
61.6361258
2.40
21.63641423
43.63646568
62.18501936
2.41
22.07589732
44.26106989
62.73620476
2.42
22.52165404
44.89119741
63.28968199
2.43
22.97373975
45.52687117
63.84545106
2.44
23.43221003
46.16811407
64.40351195
2.45
23.89712067
46.81494905
64.96386467
2.46
24.36852773
47.46739901
65.52650923
2.47
24.84648745
48.12548687
66.09144561
2.48
25.33105633
48.78923556
66.65867383
2.49
25.8222911
49.45866799
67.22819388
2.50
26.32024871
50.13380707
67.80000576
Spatial connection between our study and Sergerie et al. (2021)
Before moving to the last part of the methodology where processing steps are done on top of the estimated kernel density distribution for the CMAs and NAICS of interest, we quickly provide an intuitive explanation about the dynamic between spatial density of job locations and spatial distance of<p>Source: authors’ computations from job locations. This explanation is essential to properly understand the very nature of our data and its limits. It also provides some connections between our study and Sergerie et al. (2021) in terms of spatial mapping. A paper from the US Census Bureau (Biemer and Stokes, 1984), highlights the following ‘’Justification: The average distance between randomly distributed points in a plane is increased by
when the density of those points is decreased by a factor of . This explanation is intuitively straightforward, however, to better understand their idea, we generate the function
and perform a quick study of its first and second derivative properties, domain, and image. In Figure 5, the horizontal axis is the -axis, and a positive number represents the decreasing factor in spatial density. The blue curve is the function and it represents the increase in average distance between spatial points. The domain and image of the function
is the set of the non-negative real numbers . That is,
and . The red curve is the first derivative
and it is equal to
and there is no value of
such that it is equal to zero even though it converge to zero, as tend to infinity . Both its domain and image are
since no zero can be at the denominator. The green curve is the second derivative
and it is equal to
and again, there is no value of
such that it is equal to zero even though it converge to zero, as
tend to infinity
and, at a faster rate of convergence than the first derivative function due to the larger absolute magnitude of parameters at the denominator (that is , and ). Its domain and image are
and , respectively (Bartle and Sherbert, 2011). Therefore, we explain the justification of (Biemer and Stokes, 1984) in the following way for our own applications about a large population of job locations spatially distributed within a population of DB in a CMA: For a large and fixed number of uniform job locations within the DB of a CMA, as the DB polygon boundaries tend to infinity , and as the spatial density of job’s location tend to infinity within the
,Note and as the average spatial distance between job’s location tend to infinity , and as the spatial random uniformity of job locations keep holding within the DB, the variation in average spatial distance between the job’s location converge in probability toward 0 between increment. In a finite context, for two distinct large expansions,
and , where , and each expansion are implemented independently and from the initial DB polygon boundaries, then . The information provided in this paragraph connects the theoretical dots between our study and Sergerie et al. (2021). In other words, the material enables the understanding of how the data of a DB would behave if the surface of randomization expands from a DB to a DA (Dissemination Area) or even an ADA (Aggregate Dissemination Area).
Description for Figure 5
Figure 5
Dynamics between spatial density and spatial distance between points ( g(k), g'(k) & g''(k))
Table summary The information is grouped by Dynamics Between Spatial Density and Spatial Distance Between Points (k) (appearing as row headers), , calculated using (appearing as column headers).
Dynamics Between Spatial Density and Spatial Distance Between Points (k)
Original Function (blue)
g(k)
First Derivative (red)
g'(k)
Second Derivative (green)
g''(k)
Note ...
not applicable
Source: authors’ computations based on a textual justification provided in Biemer and Stokes, 1984.
Legend: Blue = function of interest
, red = first derivative of
, and green = second derivative of
0.00
0.00
... not applicable
... not applicable
0.01
0.10
5
-250
0.02
0.14
3.535533906
-88.38834765
0.03
0.17
2.886751346
-48.11252243
0.04
0.20
2.5
-31.25
0.05
0.22
2.236067977
-22.36067977
0.06
0.24
2.041241452
-17.01034544
0.07
0.26
1.889822365
-13.49873118
0.08
0.28
1.767766953
-11.04854346
0.09
0.30
1.666666667
-9.259259259
0.10
0.32
1.58113883
-7.90569415
0.11
0.33
1.507556723
-6.852530559
0.12
0.35
1.443375673
-6.014065304
0.13
0.36
1.386750491
-5.333655733
0.14
0.37
1.33630621
-4.772522177
0.15
0.39
1.290994449
-4.303314829
0.16
0.40
1.25
-3.90625
0.17
0.41
1.212678125
-3.566700368
0.18
0.42
1.178511302
-3.273642505
0.19
0.44
1.147078669
-3.018628077
0.20
0.45
1.118033989
-2.795084972
0.21
0.46
1.091089451
-2.597832027
0.22
0.47
1.066003582
-2.422735413
0.23
0.48
1.04257207
-2.266461022
0.24
0.49
1.020620726
-2.126293179
0.25
0.50
1
-2
0.26
0.51
0.980580676
-1.885732069
0.27
0.52
0.962250449
-1.781945275
0.28
0.53
0.944911183
-1.687341397
0.29
0.54
0.928476691
-1.600821881
0.30
0.55
0.912870929
-1.521451549
0.31
0.56
0.89802651
-1.448429855
0.32
0.57
0.883883476
-1.381067932
0.33
0.57
0.87038828
-1.318770121
0.34
0.58
0.857492926
-1.261019008
0.35
0.59
0.845154255
-1.207363221
0.36
0.60
0.833333333
-1.157407407
0.37
0.61
0.821994937
-1.110803968
0.38
0.62
0.811107106
-1.067246192
0.39
0.62
0.800640769
-1.026462524
0.40
0.63
0.790569415
-0.988211769
0.41
0.64
0.780868809
-0.952279036
0.42
0.65
0.77151675
-0.918472321
0.43
0.66
0.762492852
-0.886619595
0.44
0.66
0.753778361
-0.85656632
0.45
0.67
0.745355992
-0.828173325
0.46
0.68
0.737209781
-0.801314979
0.47
0.69
0.729324957
-0.775877614
0.48
0.69
0.721687836
-0.751758163
0.49
0.70
0.714285714
-0.728862974
0.50
0.71
0.707106781
-0.707106781
0.51
0.71
0.700140042
-0.686411806
0.52
0.72
0.693375245
-0.666706967
0.53
0.73
0.68680282
-0.647927188
0.54
0.73
0.680413817
-0.630012794
0.55
0.74
0.674199862
-0.612908966
0.56
0.75
0.668153105
-0.596565272
0.57
0.75
0.662266179
-0.580935244
0.58
0.76
0.656532164
-0.565976004
0.59
0.77
0.650944555
-0.551647928
0.60
0.77
0.645497224
-0.537914354
0.61
0.78
0.6401844
-0.524741311
0.62
0.79
0.635000635
-0.512097286
0.63
0.79
0.629940788
-0.499953007
0.64
0.80
0.625
-0.48828125
0.65
0.81
0.620173673
-0.477056671
0.66
0.81
0.615457455
-0.466255648
0.67
0.82
0.610847222
-0.455856136
0.68
0.82
0.606339063
-0.445837546
0.69
0.83
0.601929265
-0.436180627
0.70
0.84
0.597614305
-0.42686736
0.71
0.84
0.593390829
-0.417880866
0.72
0.85
0.589255651
-0.409205313
0.73
0.85
0.585205736
-0.400825847
0.74
0.86
0.581238194
-0.392728509
0.75
0.87
0.577350269
-0.384900179
0.76
0.87
0.573539335
-0.37732851
0.77
0.88
0.569802882
-0.370001872
0.78
0.88
0.566138517
-0.362909306
0.79
0.89
0.56254395
-0.356040475
0.80
0.89
0.559016994
-0.349385621
0.81
0.90
0.555555556
-0.342935528
0.82
0.91
0.55215763
-0.336681482
0.83
0.91
0.5488213
-0.330615241
0.84
0.92
0.545544726
-0.324729003
0.85
0.92
0.542326145
-0.319015379
0.86
0.93
0.539163866
-0.313467364
0.87
0.93
0.536056267
-0.308078315
0.88
0.94
0.533001791
-0.302841927
0.89
0.94
0.52999894
-0.297752213
0.90
0.95
0.527046277
-0.292803487
0.91
0.95
0.524142418
-0.28799034
0.92
0.96
0.521286035
-0.283307628
0.93
0.96
0.518475847
-0.278750456
0.94
0.97
0.515710623
-0.274314161
0.95
0.97
0.512989176
-0.269994303
0.96
0.98
0.510310363
-0.265786647
0.97
0.98
0.507673083
-0.261687156
0.98
0.99
0.505076272
-0.257691976
0.99
0.99
0.502518908
-0.253797428
1.00
1.00
0.5
-0.25
1.01
1.00
0.497518595
-0.246296334
1.02
1.01
0.495073771
-0.242683221
1.03
1.01
0.492664639
-0.239157592
1.04
1.02
0.490290338
-0.235716509
1.05
1.02
0.487950036
-0.23235716
1.06
1.03
0.485642931
-0.229076854
1.07
1.03
0.483368245
-0.225873011
1.08
1.04
0.481125224
-0.222743159
1.09
1.04
0.478913143
-0.219684928
1.10
1.05
0.476731295
-0.216696043
1.11
1.05
0.474578998
-0.213774323
1.12
1.06
0.472455591
-0.210917675
1.13
1.06
0.470360434
-0.208124086
1.14
1.07
0.468292906
-0.205391625
1.15
1.07
0.466252404
-0.202718437
1.16
1.08
0.464238345
-0.200102735
1.17
1.08
0.462250164
-0.197542805
1.18
1.09
0.460287309
-0.195036995
1.19
1.09
0.458349249
-0.192583718
1.20
1.10
0.456435465
-0.190181444
1.21
1.10
0.454545455
-0.1878287
1.22
1.10
0.45267873
-0.18552407
1.23
1.11
0.450834817
-0.183266186
1.24
1.11
0.449013255
-0.181053732
1.25
1.12
0.447213595
-0.178885438
1.26
1.12
0.445435403
-0.176760081
1.27
1.13
0.443678255
-0.174676478
1.28
1.13
0.441941738
-0.172633492
1.29
1.14
0.440225453
-0.170630021
1.30
1.14
0.43852901
-0.168665004
1.31
1.14
0.436852028
-0.166737415
1.32
1.15
0.43519414
-0.164846265
1.33
1.15
0.433554985
-0.162990596
1.34
1.16
0.431934213
-0.161169482
1.35
1.16
0.430331483
-0.159382031
1.36
1.17
0.428746463
-0.157627376
1.37
1.17
0.427178829
-0.155904682
1.38
1.17
0.425628265
-0.15421314
1.39
1.18
0.424094465
-0.152551966
1.40
1.18
0.422577127
-0.150920403
1.41
1.19
0.421075961
-0.149317717
1.42
1.19
0.419590679
-0.147743197
1.43
1.20
0.418121005
-0.146196156
1.44
1.20
0.416666667
-0.144675926
1.45
1.20
0.415227399
-0.143181862
1.46
1.21
0.413802944
-0.141713337
1.47
1.21
0.412393049
-0.140269745
1.48
1.22
0.410997468
-0.138850496
1.49
1.22
0.40961596
-0.13745502
1.50
1.22
0.40824829
-0.136082763
1.51
1.23
0.406894229
-0.134733189
1.52
1.23
0.405553553
-0.133405774
1.53
1.24
0.404226042
-0.132100014
1.54
1.24
0.402911482
-0.130815416
1.55
1.24
0.401609664
-0.129551505
1.56
1.25
0.400320385
-0.128307816
1.57
1.25
0.399043442
-0.127083899
1.58
1.26
0.397778642
-0.125879317
1.59
1.26
0.396525793
-0.124693646
1.60
1.26
0.395284708
-0.123526471
1.61
1.27
0.394055203
-0.122377392
1.62
1.27
0.392837101
-0.121246019
1.63
1.28
0.391630225
-0.120131971
1.64
1.28
0.390434405
-0.119034879
1.65
1.28
0.389249472
-0.117954385
1.66
1.29
0.388075263
-0.116890139
1.67
1.29
0.386911616
-0.115841801
1.68
1.30
0.385758375
-0.11480904
1.69
1.30
0.384615385
-0.113791534
1.70
1.30
0.383482494
-0.112788969
1.71
1.31
0.382359556
-0.11180104
1.72
1.31
0.381246426
-0.110827449
1.73
1.32
0.380142961
-0.109867908
1.74
1.32
0.379049022
-0.108922133
1.75
1.32
0.377964473
-0.107989849
1.76
1.33
0.376889181
-0.10707079
1.77
1.33
0.375823014
-0.106164693
1.78
1.33
0.374765844
-0.105271305
1.79
1.34
0.373717546
-0.104390376
1.80
1.34
0.372677996
-0.103521666
1.81
1.35
0.371647073
-0.102664937
1.82
1.35
0.370624658
-0.101819961
1.83
1.35
0.369610635
-0.100986512
1.84
1.36
0.36860489
-0.100164372
1.85
1.36
0.367607311
-0.099353327
1.86
1.36
0.366617788
-0.098553169
1.87
1.37
0.365636212
-0.097763693
1.88
1.37
0.364662479
-0.096984702
1.89
1.37
0.363696484
-0.096216001
1.90
1.38
0.362738125
-0.095457401
1.91
1.38
0.361787303
-0.094708718
1.92
1.39
0.360843918
-0.09396977
1.93
1.39
0.359907875
-0.093240382
1.94
1.39
0.358979079
-0.092520381
1.95
1.40
0.358057437
-0.091809599
1.96
1.40
0.357142857
-0.091107872
1.97
1.40
0.35623525
-0.090415038
1.98
1.41
0.355334527
-0.089730941
1.99
1.41
0.354440603
-0.089055428
2.00
1.41
0.353553391
-0.088388348
2.01
1.42
0.352672808
-0.087729554
2.02
1.42
0.351798772
-0.087078904
2.03
1.42
0.350931203
-0.086436257
2.04
1.43
0.350070021
-0.085801476
2.05
1.43
0.349215148
-0.085174426
2.06
1.44
0.348366507
-0.084554977
2.07
1.44
0.347524023
-0.083943001
2.08
1.44
0.346687623
-0.083338371
2.09
1.45
0.345857232
-0.082740965
2.10
1.45
0.34503278
-0.082150662
2.11
1.45
0.344214195
-0.081567345
2.12
1.46
0.34340141
-0.080990899
2.13
1.46
0.342594355
-0.08042121
2.14
1.46
0.341792964
-0.079858169
2.15
1.47
0.34099717
-0.079301667
2.16
1.47
0.340206909
-0.078751599
2.17
1.47
0.339422117
-0.078207861
2.18
1.48
0.338642731
-0.077670351
2.19
1.48
0.337868689
-0.07713897
2.20
1.48
0.337099931
-0.076613621
2.21
1.49
0.336336397
-0.076094207
2.22
1.49
0.335578028
-0.075580637
2.23
1.49
0.334824765
-0.075072817
2.24
1.50
0.334076552
-0.074570659
2.25
1.50
0.333333333
-0.074074074
2.26
1.50
0.332595053
-0.073582976
2.27
1.51
0.331861656
-0.073097281
2.28
1.51
0.331133089
-0.072616906
2.29
1.51
0.3304093
-0.072141769
2.30
1.52
0.329690237
-0.071671791
2.31
1.52
0.328975847
-0.071206893
2.32
1.52
0.328266082
-0.070747
2.33
1.53
0.327560891
-0.070292037
2.34
1.53
0.326860225
-0.069841928
2.35
1.53
0.326164037
-0.069396604
2.36
1.54
0.325472277
-0.068955991
2.37
1.54
0.324784901
-0.068520021
2.38
1.54
0.324101862
-0.068088626
2.39
1.55
0.323423114
-0.067661739
2.40
1.55
0.322748612
-0.067239294
2.41
1.55
0.322078313
-0.066821227
2.42
1.56
0.321412173
-0.066407474
2.43
1.56
0.32075015
-0.065997973
2.44
1.56
0.3200922
-0.065592664
2.45
1.57
0.319438282
-0.065191486
2.46
1.57
0.318788357
-0.064794381
2.47
1.57
0.318142381
-0.064401292
2.48
1.57
0.317500318
-0.064012161
2.49
1.58
0.316862125
-0.063626933
2.50
1.58
0.316227766
-0.063245553
2.51
1.58
0.315597202
-0.062867968
2.52
1.59
0.314970394
-0.062494126
2.53
1.59
0.314347307
-0.062123974
2.54
1.59
0.313727903
-0.061757461
2.55
1.60
0.313112146
-0.061394538
2.56
1.60
0.3125
-0.061035156
2.57
1.60
0.311891431
-0.060679267
2.58
1.61
0.311286403
-0.060326822
2.59
1.61
0.310684883
-0.059977777
2.60
1.61
0.310086836
-0.059632084
2.61
1.62
0.30949223
-0.059289699
2.62
1.62
0.308901032
-0.058950579
2.63
1.62
0.308313208
-0.058614678
2.64
1.62
0.307728727
-0.058281956
2.65
1.63
0.307147558
-0.05795237
2.66
1.63
0.30656967
-0.057625878
2.67
1.63
0.305995031
-0.05730244
2.68
1.64
0.305423611
-0.056982017
2.69
1.64
0.30485538
-0.056664569
2.70
1.64
0.30429031
-0.056350057
2.71
1.65
0.30372837
-0.056038445
2.72
1.65
0.303169531
-0.055729693
2.73
1.65
0.302613766
-0.055423767
2.74
1.66
0.302061047
-0.055120629
2.75
1.66
0.301511345
-0.054820244
2.76
1.66
0.300964633
-0.054522578
2.77
1.66
0.300420884
-0.054227596
2.78
1.67
0.299880072
-0.053935265
2.79
1.67
0.29934217
-0.05364555
2.80
1.67
0.298807152
-0.05335842
2.81
1.68
0.298274993
-0.053073842
2.82
1.68
0.297745667
-0.052791785
2.83
1.68
0.297219149
-0.052512217
2.84
1.69
0.296695415
-0.052235108
2.85
1.69
0.296174439
-0.051960428
2.86
1.69
0.295656198
-0.051688147
2.87
1.69
0.295140668
-0.051418235
2.88
1.70
0.294627825
-0.051150664
2.89
1.70
0.294117647
-0.050885406
2.90
1.70
0.29361011
-0.050622433
2.91
1.71
0.293105191
-0.050361717
2.92
1.71
0.292602868
-0.050103231
2.93
1.71
0.292103119
-0.049846949
2.94
1.71
0.291605922
-0.049592844
2.95
1.72
0.291111255
-0.049340891
2.96
1.72
0.290619097
-0.049091064
2.97
1.72
0.290129427
-0.048843338
2.98
1.73
0.289642223
-0.048597688
2.99
1.73
0.289157466
-0.048354091
3.00
1.73
0.288675135
-0.048112522
3.01
1.73
0.288195209
-0.047872958
3.02
1.74
0.287717669
-0.047635376
3.03
1.74
0.287242495
-0.047399752
3.04
1.74
0.286769667
-0.047166064
3.05
1.75
0.286299167
-0.04693429
3.06
1.75
0.285830975
-0.046704408
3.07
1.75
0.285365073
-0.046476396
3.08
1.75
0.284901441
-0.046250234
3.09
1.76
0.284440062
-0.0460259
3.10
1.76
0.283980917
-0.045803374
3.11
1.76
0.283523989
-0.045582635
3.12
1.77
0.283069259
-0.045363663
3.13
1.77
0.282616709
-0.045146439
3.14
1.77
0.282166324
-0.044930943
3.15
1.77
0.281718085
-0.044717156
3.16
1.78
0.281271975
-0.044505059
3.17
1.78
0.280827978
-0.044294634
3.18
1.78
0.280386077
-0.044085861
3.19
1.79
0.279946255
-0.043878723
3.20
1.79
0.279508497
-0.043673203
3.21
1.79
0.279072786
-0.043469281
3.22
1.79
0.278639106
-0.043266942
3.23
1.80
0.278207442
-0.043066167
3.24
1.80
0.277777778
-0.042866941
3.25
1.80
0.277350098
-0.042669246
3.26
1.81
0.276924388
-0.042473066
3.27
1.81
0.276500632
-0.042278384
3.28
1.81
0.276078815
-0.042085185
3.29
1.81
0.275658923
-0.041893453
3.30
1.82
0.275240941
-0.041703173
3.31
1.82
0.274824855
-0.041514329
3.32
1.82
0.27441065
-0.041326905
3.33
1.82
0.273998312
-0.041140888
3.34
1.83
0.273587828
-0.040956262
3.35
1.83
0.273179182
-0.040773012
3.36
1.83
0.272772363
-0.040591125
3.37
1.84
0.272367355
-0.040410587
3.38
1.84
0.271964147
-0.040231383
3.39
1.84
0.271562723
-0.040053499
3.40
1.84
0.271163072
-0.039876922
3.41
1.85
0.270765181
-0.039701639
3.42
1.85
0.270369035
-0.039527637
3.43
1.85
0.269974624
-0.039354901
3.44
1.85
0.269581933
-0.03918342
3.45
1.86
0.269190951
-0.039013181
3.46
1.86
0.268801665
-0.038844171
3.47
1.86
0.268414064
-0.038676378
3.48
1.87
0.268028134
-0.038509789
3.49
1.87
0.267643864
-0.038344393
3.50
1.87
0.267261242
-0.038180177
3.51
1.87
0.266880256
-0.038017131
3.52
1.88
0.266500895
-0.037855241
3.53
1.88
0.266123148
-0.037694497
3.54
1.88
0.265747002
-0.037534887
3.55
1.88
0.265372446
-0.037376401
3.56
1.89
0.26499947
-0.037219027
3.57
1.89
0.264628062
-0.037062754
3.58
1.89
0.264258211
-0.036907571
3.59
1.89
0.263889907
-0.036753469
3.60
1.90
0.263523138
-0.036600436
3.61
1.90
0.263157895
-0.036448462
3.62
1.90
0.262794166
-0.036297537
3.63
1.91
0.262431941
-0.03614765
3.64
1.91
0.262071209
-0.035998792
3.65
1.91
0.261711961
-0.035850954
3.66
1.91
0.261354187
-0.035704124
3.67
1.92
0.260997875
-0.035558294
3.68
1.92
0.260643018
-0.035413453
3.69
1.92
0.260289603
-0.035269594
3.70
1.92
0.259937622
-0.035126706
3.71
1.93
0.259587066
-0.03498478
3.72
1.93
0.259237924
-0.034843807
3.73
1.93
0.258890187
-0.034703778
3.74
1.93
0.258543845
-0.034564685
3.75
1.94
0.25819889
-0.034426519
3.76
1.94
0.257855312
-0.03428927
3.77
1.94
0.257513101
-0.034152931
3.78
1.94
0.25717225
-0.034017493
3.79
1.95
0.256832748
-0.033882948
3.80
1.95
0.256494588
-0.033749288
3.81
1.95
0.25615776
-0.033616504
3.82
1.95
0.255822255
-0.033484588
3.83
1.96
0.255488065
-0.033353533
3.84
1.96
0.255155182
-0.033223331
3.85
1.96
0.254823596
-0.033093973
3.86
1.96
0.254493299
-0.032965453
3.87
1.97
0.254164284
-0.032837763
3.88
1.97
0.253836541
-0.032710894
3.89
1.97
0.253510063
-0.032584841
3.90
1.97
0.253184842
-0.032459595
3.91
1.98
0.252860869
-0.032335149
3.92
1.98
0.252538136
-0.032211497
3.93
1.98
0.252216636
-0.032088631
3.94
1.98
0.251896361
-0.031966543
3.95
1.99
0.251577303
-0.031845228
3.96
1.99
0.251259454
-0.031724679
3.97
1.99
0.250942807
-0.031604887
3.98
1.99
0.250627354
-0.031485848
3.99
2.00
0.250313087
-0.031367555
4.00
2.00
0.25
-0.03125
4.01
2.00
0.249688085
-0.031133178
4.02
2.00
0.249377334
-0.031017081
4.03
2.01
0.249067741
-0.030901705
4.04
2.01
0.248759298
-0.030787042
4.05
2.01
0.248451997
-0.030673086
4.06
2.01
0.248145833
-0.030559832
4.07
2.02
0.247840799
-0.030447273
4.08
2.02
0.247536886
-0.030335403
4.09
2.02
0.247234088
-0.030224216
4.10
2.02
0.246932399
-0.030113707
4.11
2.03
0.246631812
-0.03000387
4.12
2.03
0.24633232
-0.029894699
4.13
2.03
0.246033916
-0.029786188
4.14
2.03
0.245736594
-0.029678333
4.15
2.04
0.245440347
-0.029571126
4.16
2.04
0.245145169
-0.029464564
4.17
2.04
0.244851053
-0.02935864
4.18
2.04
0.244557994
-0.029253349
4.19
2.05
0.244265984
-0.029148685
4.20
2.05
0.243975018
-0.029044645
4.21
2.05
0.243685089
-0.028941222
4.22
2.05
0.243396192
-0.028838411
4.23
2.06
0.243108319
-0.028736208
4.24
2.06
0.242821466
-0.028634607
4.25
2.06
0.242535625
-0.028533603
4.26
2.06
0.242250792
-0.028433191
4.27
2.07
0.241966959
-0.028333368
4.28
2.07
0.241684122
-0.028234126
4.29
2.07
0.241402275
-0.028135463
4.30
2.07
0.241121411
-0.028037373
4.31
2.08
0.240841525
-0.027939852
4.32
2.08
0.240562612
-0.027842895
4.33
2.08
0.240284666
-0.027746497
4.34
2.08
0.24000768
-0.027650654
4.35
2.09
0.239731651
-0.027555362
4.36
2.09
0.239456571
-0.027460616
4.37
2.09
0.239182437
-0.027366412
4.38
2.09
0.238909241
-0.027272744
4.39
2.10
0.23863698
-0.02717961
4.40
2.10
0.238365647
-0.027087005
4.41
2.10
0.238095238
-0.026994925
4.42
2.10
0.237825747
-0.026903365
4.43
2.10
0.237557169
-0.026812322
4.44
2.11
0.237289499
-0.02672179
4.45
2.11
0.237022732
-0.026631768
4.46
2.11
0.236756862
-0.026542249
4.47
2.11
0.236491885
-0.026453231
4.48
2.12
0.236227796
-0.026364709
4.49
2.12
0.235964589
-0.02627668
4.50
2.12
0.23570226
-0.02618914
4.51
2.12
0.235440805
-0.026102085
4.52
2.13
0.235180217
-0.026015511
4.53
2.13
0.234920493
-0.025929414
4.54
2.13
0.234661627
-0.025843792
4.55
2.13
0.234403615
-0.025758639
4.56
2.14
0.234146453
-0.025673953
4.57
2.14
0.233890135
-0.02558973
4.58
2.14
0.233634657
-0.025505967
4.59
2.14
0.233380014
-0.025422659
4.60
2.14
0.233126202
-0.025339805
4.61
2.15
0.232873216
-0.025257399
4.62
2.15
0.232621053
-0.025175439
4.63
2.15
0.232369706
-0.025093921
4.64
2.15
0.232119173
-0.025012842
4.65
2.16
0.231869448
-0.024932199
4.66
2.16
0.231620527
-0.024851988
4.67
2.16
0.231372407
-0.024772206
4.68
2.16
0.231125082
-0.024692851
4.69
2.17
0.230878548
-0.024613918
4.70
2.17
0.230632802
-0.024535404
4.71
2.17
0.230387839
-0.024457308
4.72
2.17
0.230143654
-0.024379624
4.73
2.17
0.229900245
-0.024302351
4.74
2.18
0.229657606
-0.024225486
4.75
2.18
0.229415734
-0.024149025
4.76
2.18
0.229174624
-0.024072965
4.77
2.18
0.228934273
-0.023997303
4.78
2.19
0.228694677
-0.023922037
4.79
2.19
0.228455831
-0.023847164
4.80
2.19
0.228217732
-0.02377268
4.81
2.19
0.227980376
-0.023698584
4.82
2.20
0.227743759
-0.023624871
4.83
2.20
0.227507878
-0.02355154
4.84
2.20
0.227272727
-0.023478588
4.85
2.20
0.227038305
-0.023406011
4.86
2.20
0.226804606
-0.023333807
4.87
2.21
0.226571627
-0.023261974
4.88
2.21
0.226339365
-0.023190509
4.89
2.21
0.226107816
-0.023119409
4.90
2.21
0.225876976
-0.023048671
4.91
2.22
0.225646841
-0.022978293
4.92
2.22
0.225417409
-0.022908273
4.93
2.22
0.225188675
-0.022838608
4.94
2.22
0.224960635
-0.022769295
4.95
2.22
0.224733287
-0.022700332
4.96
2.23
0.224506628
-0.022631716
4.97
2.23
0.224280652
-0.022563446
4.98
2.23
0.224055357
-0.022495518
4.99
2.23
0.223830741
-0.02242793
5.00
2.24
0.223606798
-0.02236068
5.01
2.24
0.223383526
-0.022293765
5.02
2.24
0.223160921
-0.022227183
5.03
2.24
0.222938981
-0.022160933
5.04
2.24
0.222717702
-0.02209501
5.05
2.25
0.22249708
-0.022029414
5.06
2.25
0.222277112
-0.021964142
5.07
2.25
0.222057796
-0.021899191
5.08
2.25
0.221839127
-0.02183456
5.09
2.26
0.221621104
-0.021770246
5.10
2.26
0.221403721
-0.021706247
5.11
2.26
0.221186978
-0.021642561
5.12
2.26
0.220970869
-0.021579186
5.13
2.26
0.220755393
-0.02151612
5.14
2.27
0.220540546
-0.02145336
5.15
2.27
0.220326325
-0.021390905
5.16
2.27
0.220112727
-0.021328753
5.17
2.27
0.219899749
-0.0212669
5.18
2.28
0.219687388
-0.021205346
5.19
2.28
0.219475641
-0.021144089
5.20
2.28
0.219264505
-0.021083125
5.21
2.28
0.219053977
-0.021022455
5.22
2.28
0.218844055
-0.020962074
5.23
2.29
0.218634735
-0.020901982
5.24
2.29
0.218426014
-0.020842177
5.25
2.29
0.21821789
-0.020782656
5.26
2.29
0.21801036
-0.020723418
5.27
2.30
0.217803421
-0.020664461
5.28
2.30
0.21759707
-0.020605783
5.29
2.30
0.217391304
-0.020547382
5.30
2.30
0.217186121
-0.020489257
5.31
2.30
0.216981518
-0.020431405
5.32
2.31
0.216777492
-0.020373824
5.33
2.31
0.216574041
-0.020316514
5.34
2.31
0.216371161
-0.020259472
5.35
2.31
0.216168851
-0.020202696
5.36
2.32
0.215967106
-0.020146185
5.37
2.32
0.215765926
-0.020089937
5.38
2.32
0.215565307
-0.02003395
5.39
2.32
0.215365246
-0.019978223
5.40
2.32
0.215165741
-0.019922754
5.41
2.33
0.21496679
-0.019867541
5.42
2.33
0.21476839
-0.019812582
5.43
2.33
0.214570538
-0.019757876
5.44
2.33
0.214373231
-0.019703422
5.45
2.33
0.214176468
-0.019649217
5.46
2.34
0.213980246
-0.019595261
5.47
2.34
0.213784562
-0.01954155
5.48
2.34
0.213589414
-0.019488085
5.49
2.34
0.2133948
-0.019434863
5.50
2.35
0.213200716
-0.019381883
5.51
2.35
0.213007161
-0.019329144
5.52
2.35
0.212814133
-0.019276642
5.53
2.35
0.212621628
-0.019224379
5.54
2.35
0.212429644
-0.019172351
5.55
2.36
0.21223818
-0.019120557
5.56
2.36
0.212047232
-0.019068996
5.57
2.36
0.211856799
-0.019017666
5.58
2.36
0.211666878
-0.018966566
5.59
2.36
0.211477467
-0.018915695
5.60
2.37
0.211288564
-0.01886505
5.61
2.37
0.211100165
-0.018814632
5.62
2.37
0.21091227
-0.018764437
5.63
2.37
0.210724876
-0.018714465
5.64
2.37
0.21053798
-0.018664715
5.65
2.38
0.210351581
-0.018615184
5.66
2.38
0.210165676
-0.018565872
5.67
2.38
0.209980263
-0.018516778
5.68
2.38
0.20979534
-0.0184679
5.69
2.39
0.209610904
-0.018419236
5.70
2.39
0.209426954
-0.018370785
5.71
2.39
0.209243488
-0.018322547
5.72
2.39
0.209060503
-0.018274519
5.73
2.39
0.208877997
-0.018226701
5.74
2.40
0.208695968
-0.018179091
5.75
2.40
0.208514414
-0.018131688
5.76
2.40
0.208333333
-0.018084491
5.77
2.40
0.208152724
-0.018037498
5.78
2.40
0.207972583
-0.017990708
5.79
2.41
0.207792909
-0.01794412
5.80
2.41
0.2076137
-0.017897733
5.81
2.41
0.207434953
-0.017851545
5.82
2.41
0.207256668
-0.017805556
5.83
2.41
0.207078842
-0.017759763
5.84
2.42
0.206901472
-0.017714167
5.85
2.42
0.206724558
-0.017668766
5.86
2.42
0.206548096
-0.017623558
5.87
2.42
0.206372086
-0.017578542
5.88
2.42
0.206196525
-0.017533718
5.89
2.43
0.206021411
-0.017489084
5.90
2.43
0.205846742
-0.017444639
5.91
2.43
0.205672517
-0.017400382
5.92
2.43
0.205498734
-0.017356312
5.93
2.44
0.205325391
-0.017312428
5.94
2.44
0.205152485
-0.017268728
5.95
2.44
0.204980015
-0.017225211
5.96
2.44
0.20480798
-0.017181878
5.97
2.44
0.204636377
-0.017138725
5.98
2.45
0.204465205
-0.017095753
5.99
2.45
0.204294462
-0.01705296
6.00
2.45
0.204124145
-0.017010345
6.01
2.45
0.203954254
-0.016967908
6.02
2.45
0.203784786
-0.016925647
6.03
2.46
0.203615741
-0.016883561
6.04
2.46
0.203447115
-0.016841649
6.05
2.46
0.203278907
-0.01679991
6.06
2.46
0.203111116
-0.016758343
6.07
2.46
0.20294374
-0.016716947
6.08
2.47
0.202776776
-0.016675722
6.09
2.47
0.202610225
-0.016634665
6.10
2.47
0.202444083
-0.016593777
6.11
2.47
0.202278349
-0.016553056
6.12
2.47
0.202113021
-0.016512502
6.13
2.48
0.201948098
-0.016472112
6.14
2.48
0.201783578
-0.016431887
6.15
2.48
0.20161946
-0.016391826
6.16
2.48
0.201455741
-0.016351927
6.17
2.48
0.201292421
-0.01631219
6.18
2.49
0.201129497
-0.016272613
6.19
2.49
0.200966968
-0.016233196
6.20
2.49
0.200804832
-0.016193938
6.21
2.49
0.200643088
-0.016154838
6.22
2.49
0.200481735
-0.016115895
6.23
2.50
0.20032077
-0.016077108
6.24
2.50
0.200160192
-0.016038477
6.25
2.50
0.2
-0.016
6.26
2.50
0.199840192
-0.015961677
6.27
2.50
0.199680766
-0.015923506
6.28
2.51
0.199521721
-0.015885487
6.29
2.51
0.199363056
-0.01584762
6.30
2.51
0.199204768
-0.015809902
6.31
2.51
0.199046857
-0.015772334
6.32
2.51
0.198889321
-0.015734915
6.33
2.52
0.198732158
-0.015697643
6.34
2.52
0.198575368
-0.015660518
6.35
2.52
0.198418948
-0.015623539
6.36
2.52
0.198262896
-0.015586706
6.37
2.52
0.198107213
-0.015550017
6.38
2.53
0.197951896
-0.015513471
6.39
2.53
0.197796943
-0.015477069
6.40
2.53
0.197642354
-0.015440809
6.41
2.53
0.197488126
-0.01540469
6.42
2.53
0.197334259
-0.015368712
6.43
2.54
0.197180752
-0.015332873
6.44
2.54
0.197027602
-0.015297174
6.45
2.54
0.196874808
-0.015261613
6.46
2.54
0.196722369
-0.01522619
6.47
2.54
0.196570283
-0.015190903
6.48
2.55
0.19641855
-0.015155752
6.49
2.55
0.196267168
-0.015120737
6.50
2.55
0.196116135
-0.015085857
6.51
2.55
0.19596545
-0.01505111
6.52
2.55
0.195815112
-0.015016496
6.53
2.56
0.19566512
-0.014982015
6.54
2.56
0.195515472
-0.014947666
6.55
2.56
0.195366166
-0.014913448
6.56
2.56
0.195217202
-0.01487936
6.57
2.56
0.195068579
-0.014845402
6.58
2.57
0.194920294
-0.014811572
6.59
2.57
0.194772347
-0.014777872
6.60
2.57
0.194624736
-0.014744298
6.61
2.57
0.19447746
-0.014710852
6.62
2.57
0.194330519
-0.014677532
6.63
2.57
0.194183909
-0.014644337
6.64
2.58
0.194037631
-0.014611267
6.65
2.58
0.193891684
-0.014578322
6.66
2.58
0.193746065
-0.0145455
6.67
2.58
0.193600773
-0.014512802
6.68
2.58
0.193455808
-0.014480225
6.69
2.59
0.193311168
-0.01444777
6.70
2.59
0.193166852
-0.014415437
6.71
2.59
0.193022859
-0.014383223
6.72
2.59
0.192879187
-0.01435113
6.73
2.59
0.192735836
-0.014319156
6.74
2.60
0.192592804
-0.0142873
6.75
2.60
0.19245009
-0.014255562
6.76
2.60
0.192307692
-0.014223942
6.77
2.60
0.192165611
-0.014192438
6.78
2.60
0.192023843
-0.01416105
6.79
2.61
0.191882389
-0.014129778
6.80
2.61
0.191741247
-0.014098621
6.81
2.61
0.191600416
-0.014067578
6.82
2.61
0.191459895
-0.014036649
6.83
2.61
0.191319683
-0.014005833
6.84
2.62
0.191179778
-0.01397513
6.85
2.62
0.19104018
-0.013944539
6.86
2.62
0.190900887
-0.013914059
6.87
2.62
0.190761898
-0.01388369
6.88
2.62
0.190623213
-0.013853431
6.89
2.62
0.190484829
-0.013823282
6.90
2.63
0.190346747
-0.013793243
6.91
2.63
0.190208964
-0.013763311
6.92
2.63
0.19007148
-0.013733488
6.93
2.63
0.189934294
-0.013703773
6.94
2.63
0.189797405
-0.013674165
6.95
2.64
0.18966081
-0.013644663
6.96
2.64
0.189524511
-0.013615267
6.97
2.64
0.189388505
-0.013585976
6.98
2.64
0.189252791
-0.01355679
6.99
2.64
0.189117369
-0.013527709
7.00
2.65
0.188982237
-0.013498731
7.01
2.65
0.188847394
-0.013469857
7.02
2.65
0.188712839
-0.013441085
7.03
2.65
0.188578572
-0.013412416
7.04
2.65
0.18844459
-0.013383849
7.05
2.66
0.188310894
-0.013355383
7.06
2.66
0.188177482
-0.013327017
7.07
2.66
0.188044354
-0.013298752
7.08
2.66
0.187911507
-0.013270587
7.09
2.66
0.187778942
-0.013242521
7.10
2.66
0.187646656
-0.013214553
7.11
2.67
0.18751465
-0.013186684
7.12
2.67
0.187382922
-0.013158913
7.13
2.67
0.187251472
-0.013131239
7.14
2.67
0.187120297
-0.013103662
7.15
2.67
0.186989398
-0.013076182
7.16
2.68
0.186858773
-0.013048797
7.17
2.68
0.186728422
-0.013021508
7.18
2.68
0.186598343
-0.012994314
7.19
2.68
0.186468535
-0.012967214
7.20
2.68
0.186338998
-0.012940208
7.21
2.69
0.186209731
-0.012913296
7.22
2.69
0.186080732
-0.012886477
7.23
2.69
0.185952001
-0.012859751
7.24
2.69
0.185823537
-0.012833117
7.25
2.69
0.185695338
-0.012806575
7.26
2.69
0.185567405
-0.012780124
7.27
2.70
0.185439735
-0.012753764
7.28
2.70
0.185312329
-0.012727495
7.29
2.70
0.185185185
-0.012701316
7.30
2.70
0.185058303
-0.012675226
7.31
2.70
0.18493168
-0.012649226
7.32
2.71
0.184805318
-0.012623314
7.33
2.71
0.184679214
-0.012597491
7.34
2.71
0.184553368
-0.012571755
7.35
2.71
0.184427778
-0.012546107
7.36
2.71
0.184302445
-0.012520547
7.37
2.71
0.184177367
-0.012495072
7.38
2.72
0.184052543
-0.012469685
7.39
2.72
0.183927973
-0.012444382
7.40
2.72
0.183803656
-0.012419166
7.41
2.72
0.18367959
-0.012394034
7.42
2.72
0.183555775
-0.012368988
7.43
2.73
0.18343221
-0.012344025
7.44
2.73
0.183308894
-0.012319146
7.45
2.73
0.183185826
-0.012294351
7.46
2.73
0.183063006
-0.012269639
7.47
2.73
0.182940433
-0.012245009
7.48
2.73
0.182818106
-0.012220462
7.49
2.74
0.182696024
-0.012195996
7.50
2.74
0.182574186
-0.012171612
7.51
2.74
0.182452591
-0.01214731
7.52
2.74
0.182331239
-0.012123088
7.53
2.74
0.182210129
-0.012098946
7.54
2.75
0.18208926
-0.012074885
7.55
2.75
0.181968631
-0.012050903
7.56
2.75
0.181848242
-0.012027
7.57
2.75
0.181728091
-0.012003176
7.58
2.75
0.181608178
-0.011979431
7.59
2.75
0.181488502
-0.011955764
7.60
2.76
0.181369063
-0.011932175
7.61
2.76
0.181249858
-0.011908663
7.62
2.76
0.181130889
-0.011885229
7.63
2.76
0.181012154
-0.011861871
7.64
2.76
0.180893651
-0.01183859
7.65
2.77
0.180775382
-0.011815384
7.66
2.77
0.180657343
-0.011792255
7.67
2.77
0.180539536
-0.011769201
7.68
2.77
0.180421959
-0.011746221
7.69
2.77
0.180304611
-0.011723317
7.70
2.77
0.180187493
-0.011700487
7.71
2.78
0.180070602
-0.01167773
7.72
2.78
0.179953938
-0.011655048
7.73
2.78
0.1798375
-0.011632439
7.74
2.78
0.179721289
-0.011609902
7.75
2.78
0.179605302
-0.011587439
7.76
2.79
0.17948954
-0.011565048
7.77
2.79
0.179374001
-0.011542728
7.78
2.79
0.179258685
-0.011520481
7.79
2.79
0.179143591
-0.011498305
7.80
2.79
0.179028719
-0.0114762
7.81
2.79
0.178914067
-0.011454166
7.82
2.80
0.178799635
-0.011432202
7.83
2.80
0.178685422
-0.011410308
7.84
2.80
0.178571429
-0.011388484
7.85
2.80
0.178457653
-0.011366729
7.86
2.80
0.178344094
-0.011345044
7.87
2.81
0.178230751
-0.011323428
7.88
2.81
0.178117625
-0.01130188
7.89
2.81
0.178004714
-0.0112804
7.90
2.81
0.177892017
-0.011258988
7.91
2.81
0.177779534
-0.011237644
7.92
2.81
0.177667264
-0.011216368
7.93
2.82
0.177555206
-0.011195158
7.94
2.82
0.17744336
-0.011174015
7.95
2.82
0.177331726
-0.011152939
7.96
2.82
0.177220301
-0.011131928
7.97
2.82
0.177109087
-0.011110984
7.98
2.82
0.176998081
-0.011090105
7.99
2.83
0.176887284
-0.011069292
8.00
2.83
0.176776695
-0.011048543
8.01
2.83
0.176666313
-0.01102786
8.02
2.83
0.176556138
-0.011007241
8.03
2.83
0.176446168
-0.010986685
8.04
2.84
0.176336404
-0.010966194
8.05
2.84
0.176226844
-0.010945767
8.06
2.84
0.176117488
-0.010925403
8.07
2.84
0.176008336
-0.010905101
8.08
2.84
0.175899386
-0.010884863
8.09
2.84
0.175790638
-0.010864687
8.10
2.85
0.175682092
-0.010844574
8.11
2.85
0.175573747
-0.010824522
8.12
2.85
0.175465602
-0.010804532
8.13
2.85
0.175357656
-0.010784604
8.14
2.85
0.175249909
-0.010764736
8.15
2.85
0.175142361
-0.01074493
8.16
2.86
0.175035011
-0.010725184
8.17
2.86
0.174927857
-0.010705499
8.18
2.86
0.1748209
-0.010685874
8.19
2.86
0.174714139
-0.010666309
8.20
2.86
0.174607574
-0.010646803
8.21
2.87
0.174501203
-0.010627357
8.22
2.87
0.174395027
-0.01060797
8.23
2.87
0.174289044
-0.010588642
8.24
2.87
0.174183254
-0.010569372
8.25
2.87
0.174077656
-0.010550161
8.26
2.87
0.17397225
-0.010531008
8.27
2.88
0.173867036
-0.010511913
8.28
2.88
0.173762012
-0.010492875
8.29
2.88
0.173657178
-0.010473895
8.30
2.88
0.173552534
-0.010454972
8.31
2.88
0.173448078
-0.010436106
8.32
2.88
0.173343811
-0.010417296
8.33
2.89
0.173239732
-0.010398543
8.34
2.89
0.17313584
-0.010379847
8.35
2.89
0.173032135
-0.010361206
8.36
2.89
0.172928616
-0.010342621
8.37
2.89
0.172825282
-0.010324091
8.38
2.89
0.172722134
-0.010305617
8.39
2.90
0.17261917
-0.010287197
8.40
2.90
0.17251639
-0.010268833
8.41
2.90
0.172413793
-0.010250523
8.42
2.90
0.172311379
-0.010232267
8.43
2.90
0.172209148
-0.010214066
8.44
2.91
0.172107098
-0.010195918
8.45
2.91
0.172005229
-0.010177824
8.46
2.91
0.171903541
-0.010159784
8.47
2.91
0.171802033
-0.010141797
8.48
2.91
0.171700705
-0.010123862
8.49
2.91
0.171599556
-0.010105981
8.50
2.92
0.171498585
-0.010088152
8.51
2.92
0.171397793
-0.010070376
8.52
2.92
0.171297177
-0.010052651
8.53
2.92
0.171196739
-0.010034979
8.54
2.92
0.171096478
-0.010017358
8.55
2.92
0.170996392
-0.009999789
8.56
2.93
0.170896482
-0.009982271
8.57
2.93
0.170796746
-0.009964804
8.58
2.93
0.170697185
-0.009947388
8.59
2.93
0.170597798
-0.009930023
8.60
2.93
0.170498585
-0.009912708
8.61
2.93
0.170399544
-0.009895444
8.62
2.94
0.170300676
-0.009878229
8.63
2.94
0.170201979
-0.009861065
8.64
2.94
0.170103454
-0.00984395
8.65
2.94
0.1700051
-0.009826884
8.66
2.94
0.169906917
-0.009809868
8.67
2.94
0.169808903
-0.009792901
8.68
2.95
0.169711058
-0.009775983
8.69
2.95
0.169613383
-0.009759113
8.70
2.95
0.169515876
-0.009742292
8.71
2.95
0.169418537
-0.009725519
8.72
2.95
0.169321365
-0.009708794
8.73
2.95
0.169224361
-0.009692117
8.74
2.96
0.169127523
-0.009675488
8.75
2.96
0.169030851
-0.009658906
8.76
2.96
0.168934345
-0.009642371
8.77
2.96
0.168838003
-0.009625884
8.78
2.96
0.168741827
-0.009609443
8.79
2.96
0.168645814
-0.00959305
8.80
2.97
0.168549966
-0.009576703
8.81
2.97
0.16845428
-0.009560402
8.82
2.97
0.168358757
-0.009544147
8.83
2.97
0.168263397
-0.009527939
8.84
2.97
0.168168198
-0.009511776
8.85
2.97
0.168073161
-0.009495659
8.86
2.98
0.167978285
-0.009479587
8.87
2.98
0.167883569
-0.009463561
8.88
2.98
0.167789014
-0.00944758
8.89
2.98
0.167694618
-0.009431643
8.90
2.98
0.167600381
-0.009415752
8.91
2.98
0.167506303
-0.009399905
8.92
2.99
0.167412383
-0.009384102
8.93
2.99
0.16731862
-0.009368344
8.94
2.99
0.167225016
-0.00935263
8.95
2.99
0.167131568
-0.009336959
8.96
2.99
0.167038276
-0.009321332
8.97
2.99
0.166945141
-0.009305749
8.98
3.00
0.166852161
-0.009290209
8.99
3.00
0.166759336
-0.009274713
9.00
3.00
0.166666667
-0.009259259
9.01
3.00
0.166574151
-0.009243849
9.02
3.00
0.16648179
-0.009228481
9.03
3.00
0.166389581
-0.009213155
9.04
3.01
0.166297526
-0.009197872
9.05
3.01
0.166205624
-0.009182631
9.06
3.01
0.166113874
-0.009167432
9.07
3.01
0.166022275
-0.009152275
9.08
3.01
0.165930828
-0.00913716
9.09
3.01
0.165839532
-0.009122086
9.10
3.02
0.165748386
-0.009107054
9.11
3.02
0.16565739
-0.009092063
9.12
3.02
0.165566545
-0.009077113
9.13
3.02
0.165475848
-0.009062204
9.14
3.02
0.1653853
-0.009047336
9.15
3.02
0.165294901
-0.009032508
9.16
3.03
0.16520465
-0.009017721
9.17
3.03
0.165114547
-0.009002974
9.18
3.03
0.16502459
-0.008988267
9.19
3.03
0.164934781
-0.008973601
9.20
3.03
0.164845118
-0.008958974
9.21
3.03
0.164755602
-0.008944387
9.22
3.04
0.16466623
-0.008929839
9.23
3.04
0.164577005
-0.008915331
9.24
3.04
0.164487924
-0.008900862
9.25
3.04
0.164398987
-0.008886432
9.26
3.04
0.164310195
-0.008872041
9.27
3.04
0.164221546
-0.008857689
9.28
3.05
0.164133041
-0.008843375
9.29
3.05
0.164044679
-0.0088291
9.30
3.05
0.163956459
-0.008814863
9.31
3.05
0.163868381
-0.008800665
9.32
3.05
0.163780446
-0.008786505
9.33
3.05
0.163692651
-0.008772382
9.34
3.06
0.163604998
-0.008758298
9.35
3.06
0.163517485
-0.008744251
9.36
3.06
0.163430113
-0.008730241
9.37
3.06
0.16334288
-0.008716269
9.38
3.06
0.163255787
-0.008702334
9.39
3.06
0.163168833
-0.008688436
9.40
3.07
0.163082018
-0.008674575
9.41
3.07
0.162995342
-0.008660751
9.42
3.07
0.162908803
-0.008646964
9.43
3.07
0.162822402
-0.008633213
9.44
3.07
0.162736139
-0.008619499
9.45
3.07
0.162650012
-0.008605821
9.46
3.08
0.162564022
-0.008592179
9.47
3.08
0.162478168
-0.008578573
9.48
3.08
0.162392451
-0.008565003
9.49
3.08
0.162306868
-0.008551468
9.50
3.08
0.162221421
-0.00853797
9.51
3.08
0.162136109
-0.008524506
9.52
3.09
0.162050931
-0.008511078
9.53
3.09
0.161965887
-0.008497686
9.54
3.09
0.161880977
-0.008484328
9.55
3.09
0.1617962
-0.008471005
9.56
3.09
0.161711557
-0.008457717
9.57
3.09
0.161627046
-0.008444464
9.58
3.10
0.161542667
-0.008431246
9.59
3.10
0.161458421
-0.008418062
9.60
3.10
0.161374306
-0.008404912
9.61
3.10
0.161290323
-0.008391796
9.62
3.10
0.16120647
-0.008378715
9.63
3.10
0.161122748
-0.008365667
9.64
3.10
0.161039157
-0.008352653
9.65
3.11
0.160955695
-0.008339673
9.66
3.11
0.160872363
-0.008326727
9.67
3.11
0.16078916
-0.008313814
9.68
3.11
0.160706087
-0.008300934
9.69
3.11
0.160623142
-0.008288088
9.70
3.11
0.160540325
-0.008275274
9.71
3.12
0.160457636
-0.008262494
9.72
3.12
0.160375075
-0.008249747
9.73
3.12
0.160292641
-0.008237032
9.74
3.12
0.160210334
-0.00822435
9.75
3.12
0.160128154
-0.0082117
9.76
3.12
0.1600461
-0.008199083
9.77
3.13
0.159964172
-0.008186498
9.78
3.13
0.15988237
-0.008173945
9.79
3.13
0.159800693
-0.008161425
9.80
3.13
0.159719141
-0.008148936
9.81
3.13
0.159637714
-0.008136479
9.82
3.13
0.159556412
-0.008124054
9.83
3.14
0.159475233
-0.00811166
9.84
3.14
0.159394178
-0.008099298
9.85
3.14
0.159313247
-0.008086967
9.86
3.14
0.159232439
-0.008074667
9.87
3.14
0.159151754
-0.008062399
9.88
3.14
0.159071191
-0.008050161
9.89
3.14
0.15899075
-0.008037955
9.90
3.15
0.158910432
-0.008025779
9.91
3.15
0.158830234
-0.008013634
9.92
3.15
0.158750159
-0.00800152
9.93
3.15
0.158670204
-0.007989436
9.94
3.15
0.15859037
-0.007977383
9.95
3.15
0.158510656
-0.00796536
9.96
3.16
0.158431063
-0.007953367
9.97
3.16
0.158351589
-0.007941404
9.98
3.16
0.158272234
-0.007929471
9.99
3.16
0.158192999
-0.007917568
10.00
3.16
0.158113883
-0.007905694
10.01
3.16
0.158034885
-0.00789385
10.02
3.17
0.157956006
-0.007882036
10.03
3.17
0.157877244
-0.007870251
10.04
3.17
0.157798601
-0.007858496
10.05
3.17
0.157720074
-0.00784677
10.06
3.17
0.157641665
-0.007835073
10.07
3.17
0.157563373
-0.007823405
10.08
3.17
0.157485197
-0.007811766
10.09
3.18
0.157407138
-0.007800155
10.10
3.18
0.157329194
-0.007788574
10.11
3.18
0.157251366
-0.007777021
10.12
3.18
0.157173653
-0.007765497
10.13
3.18
0.157096056
-0.007754001
10.14
3.18
0.157018573
-0.007742533
10.15
3.19
0.156941205
-0.007731094
10.16
3.19
0.156863951
-0.007719683
10.17
3.19
0.156786811
-0.007708299
10.18
3.19
0.156709785
-0.007696944
10.19
3.19
0.156632872
-0.007685617
10.20
3.19
0.156556073
-0.007674317
10.21
3.20
0.156479386
-0.007663045
10.22
3.20
0.156402812
-0.007651801
10.23
3.20
0.15632635
-0.007640584
10.24
3.20
0.15625
-0.007629395
10.25
3.20
0.156173762
-0.007618232
10.26
3.20
0.156097635
-0.007607097
10.27
3.20
0.15602162
-0.007595989
10.28
3.21
0.155945715
-0.007584908
10.29
3.21
0.155869922
-0.007573854
10.30
3.21
0.155794238
-0.007562827
10.31
3.21
0.155718665
-0.007551827
10.32
3.21
0.155643202
-0.007540853
10.33
3.21
0.155567848
-0.007529906
10.34
3.22
0.155492603
-0.007518985
10.35
3.22
0.155417468
-0.00750809
10.36
3.22
0.155342442
-0.007497222
10.37
3.22
0.155267524
-0.00748638
10.38
3.22
0.155192714
-0.007475564
10.39
3.22
0.155118012
-0.007464774
10.40
3.22
0.155043418
-0.00745401
10.41
3.23
0.154968932
-0.007443272
10.42
3.23
0.154894553
-0.00743256
10.43
3.23
0.154820281
-0.007421873
10.44
3.23
0.154746115
-0.007411212
10.45
3.23
0.154672056
-0.007400577
10.46
3.23
0.154598104
-0.007389967
10.47
3.24
0.154524257
-0.007379382
10.48
3.24
0.154450516
-0.007368822
10.49
3.24
0.15437688
-0.007358288
10.50
3.24
0.15430335
-0.007347779
10.51
3.24
0.154229925
-0.007337294
10.52
3.24
0.154156604
-0.007326835
10.53
3.24
0.154083388
-0.0073164
10.54
3.25
0.154010276
-0.00730599
10.55
3.25
0.153937268
-0.007295605
10.56
3.25
0.153864364
-0.007285244
10.57
3.25
0.153791563
-0.007274908
10.58
3.25
0.153718865
-0.007264597
10.59
3.25
0.153646271
-0.007254309
10.60
3.26
0.153573779
-0.007244046
10.61
3.26
0.15350139
-0.007233807
10.62
3.26
0.153429103
-0.007223592
10.63
3.26
0.153356918
-0.007213402
10.64
3.26
0.153284835
-0.007203235
10.65
3.26
0.153212853
-0.007193092
10.66
3.26
0.153140973
-0.007182972
10.67
3.27
0.153069194
-0.007172877
10.68
3.27
0.152997515
-0.007162805
10.69
3.27
0.152925938
-0.007152757
10.70
3.27
0.15285446
-0.007142732
10.71
3.27
0.152783083
-0.00713273
10.72
3.27
0.152711805
-0.007122752
10.73
3.28
0.152640628
-0.007112797
10.74
3.28
0.152569549
-0.007102865
10.75
3.28
0.15249857
-0.007092957
10.76
3.28
0.15242769
-0.007083071
10.77
3.28
0.152356909
-0.007073208
10.78
3.28
0.152286226
-0.007063369
10.79
3.28
0.152215641
-0.007053552
10.80
3.29
0.152145155
-0.007043757
10.81
3.29
0.152074766
-0.007033985
10.82
3.29
0.152004475
-0.007024236
10.83
3.29
0.151934281
-0.00701451
10.84
3.29
0.151864185
-0.007004806
10.85
3.29
0.151794185
-0.006995124
10.86
3.30
0.151724282
-0.006985464
10.87
3.30
0.151654476
-0.006975827
10.88
3.30
0.151584766
-0.006966212
10.89
3.30
0.151515152
-0.006956619
10.90
3.30
0.151445633
-0.006947047
10.91
3.30
0.15137621
-0.006937498
10.92
3.30
0.151306883
-0.006927971
10.93
3.31
0.151237651
-0.006918465
10.94
3.31
0.151168514
-0.006908981
10.95
3.31
0.151099471
-0.006899519
10.96
3.31
0.151030523
-0.006890079
10.97
3.31
0.15096167
-0.00688066
10.98
3.31
0.15089291
-0.006871262
10.99
3.32
0.150824244
-0.006861886
11.00
3.32
0.150755672
-0.006852531
11.01
3.32
0.150687194
-0.006843197
11.02
3.32
0.150618808
-0.006833884
11.03
3.32
0.150550516
-0.006824593
11.04
3.32
0.150482316
-0.006815322
11.05
3.32
0.150414209
-0.006806073
11.06
3.33
0.150346195
-0.006796844
11.07
3.33
0.150278272
-0.006787637
11.08
3.33
0.150210442
-0.00677845
11.09
3.33
0.150142703
-0.006769283
11.10
3.33
0.150075056
-0.006760138
11.11
3.33
0.150007501
-0.006751013
11.12
3.33
0.149940036
-0.006741908
11.13
3.34
0.149872662
-0.006732824
11.14
3.34
0.149805379
-0.00672376
11.15
3.34
0.149738187
-0.006714717
11.16
3.34
0.149671085
-0.006705694
11.17
3.34
0.149604073
-0.006696691
11.18
3.34
0.149537151
-0.006687708
11.19
3.35
0.149470319
-0.006678745
11.20
3.35
0.149403576
-0.006669803
11.21
3.35
0.149336923
-0.00666088
11.22
3.35
0.149270359
-0.006651977
11.23
3.35
0.149203883
-0.006643094
11.24
3.35
0.149137497
-0.00663423
11.25
3.35
0.149071198
-0.006625387
11.26
3.36
0.149004989
-0.006616563
11.27
3.36
0.148938867
-0.006607758
11.28
3.36
0.148872834
-0.006598973
11.29
3.36
0.148806888
-0.006590208
11.30
3.36
0.148741029
-0.006581461
11.31
3.36
0.148675258
-0.006572735
11.32
3.36
0.148609575
-0.006564027
11.33
3.37
0.148543978
-0.006555339
11.34
3.37
0.148478468
-0.00654667
11.35
3.37
0.148413044
-0.00653802
11.36
3.37
0.148347707
-0.006529389
11.37
3.37
0.148282456
-0.006520776
11.38
3.37
0.148217292
-0.006512183
11.39
3.37
0.148152213
-0.006503609
11.40
3.38
0.148087219
-0.006495053
11.41
3.38
0.148022312
-0.006486517
11.42
3.38
0.147957489
-0.006477999
11.43
3.38
0.147892752
-0.006469499
11.44
3.38
0.147828099
-0.006461018
11.45
3.38
0.147763531
-0.006452556
11.46
3.39
0.147699048
-0.006444112
11.47
3.39
0.147634649
-0.006435687
11.48
3.39
0.147570334
-0.006427279
11.49
3.39
0.147506103
-0.00641889
11.50
3.39
0.147441956
-0.00641052
11.51
3.39
0.147377893
-0.006402167
11.52
3.39
0.147313913
-0.006393833
11.53
3.40
0.147250016
-0.006385517
11.54
3.40
0.147186202
-0.006377218
11.55
3.40
0.147122472
-0.006368938
11.56
3.40
0.147058824
-0.006360676
11.57
3.40
0.146995258
-0.006352431
11.58
3.40
0.146931775
-0.006344204
11.59
3.40
0.146868374
-0.006335995
11.60
3.41
0.146805055
-0.006327804
11.61
3.41
0.146741818
-0.00631963
11.62
3.41
0.146678662
-0.006311474
11.63
3.41
0.146615588
-0.006303336
11.64
3.41
0.146552595
-0.006295215
11.65
3.41
0.146489684
-0.006287111
11.66
3.41
0.146426853
-0.006279025
11.67
3.42
0.146364103
-0.006270956
11.68
3.42
0.146301434
-0.006262904
11.69
3.42
0.146238845
-0.006254869
11.70
3.42
0.146176337
-0.006246852
11.71
3.42
0.146113908
-0.006238852
11.72
3.42
0.146051559
-0.006230869
11.73
3.42
0.145989291
-0.006222902
11.74
3.43
0.145927101
-0.006214953
11.75
3.43
0.145864991
-0.006207021
11.76
3.43
0.145802961
-0.006199105
11.77
3.43
0.145741009
-0.006191207
11.78
3.43
0.145679137
-0.006183325
11.79
3.43
0.145617343
-0.00617546
11.80
3.44
0.145555627
-0.006167611
11.81
3.44
0.14549399
-0.006159779
11.82
3.44
0.145432432
-0.006151964
11.83
3.44
0.145370951
-0.006144165
11.84
3.44
0.145309548
-0.006136383
11.85
3.44
0.145248223
-0.006128617
11.86
3.44
0.145186976
-0.006120867
11.87
3.45
0.145125806
-0.006113134
11.88
3.45
0.145064713
-0.006105417
11.89
3.45
0.145003698
-0.006097716
11.90
3.45
0.144942759
-0.006090032
11.91
3.45
0.144881897
-0.006082363
11.92
3.45
0.144821112
-0.006074711
11.93
3.45
0.144760403
-0.006067075
11.94
3.46
0.14469977
-0.006059454
11.95
3.46
0.144639214
-0.00605185
11.96
3.46
0.144578733
-0.006044261
11.97
3.46
0.144518328
-0.006036689
11.98
3.46
0.144457999
-0.006029132
11.99
3.46
0.144397746
-0.006021591
12.00
3.46
0.144337567
-0.006014065
12.01
3.47
0.144277464
-0.006006556
12.02
3.47
0.144217436
-0.005999061
12.03
3.47
0.144157483
-0.005991583
12.04
3.47
0.144097604
-0.00598412
12.05
3.47
0.1440378
-0.005976672
12.06
3.47
0.143978071
-0.00596924
12.07
3.47
0.143918416
-0.005961823
12.08
3.48
0.143858834
-0.005954422
12.09
3.48
0.143799327
-0.005947036
12.10
3.48
0.143739894
-0.005939665
12.11
3.48
0.143680534
-0.005932309
12.12
3.48
0.143621247
-0.005924969
12.13
3.48
0.143562034
-0.005917644
12.14
3.48
0.143502894
-0.005910333
12.15
3.49
0.143443828
-0.005903038
12.16
3.49
0.143384834
-0.005895758
12.17
3.49
0.143325912
-0.005888493
12.18
3.49
0.143267064
-0.005881242
12.19
3.49
0.143208288
-0.005874007
12.20
3.49
0.143149584
-0.005866786
12.21
3.49
0.143090952
-0.00585958
12.22
3.50
0.143032392
-0.005852389
12.23
3.50
0.142973904
-0.005845213
12.24
3.50
0.142915488
-0.005838051
12.25
3.50
0.142857143
-0.005830904
12.26
3.50
0.142798869
-0.005823771
12.27
3.50
0.142740667
-0.005816653
12.28
3.50
0.142682536
-0.00580955
12.29
3.51
0.142624476
-0.00580246
12.30
3.51
0.142566487
-0.005795386
12.31
3.51
0.142508569
-0.005788325
12.32
3.51
0.142450721
-0.005781279
12.33
3.51
0.142392943
-0.005774247
12.34
3.51
0.142335236
-0.00576723
12.35
3.51
0.142277598
-0.005760227
12.36
3.52
0.142220031
-0.005753237
12.37
3.52
0.142162534
-0.005746262
12.38
3.52
0.142105106
-0.005739302
12.39
3.52
0.142047747
-0.005732355
12.40
3.52
0.141990459
-0.005725422
12.41
3.52
0.141933239
-0.005718503
12.42
3.52
0.141876088
-0.005711598
12.43
3.53
0.141819007
-0.005704707
12.44
3.53
0.141761994
-0.005697829
12.45
3.53
0.14170505
-0.005690966
12.46
3.53
0.141648175
-0.005684116
12.47
3.53
0.141591368
-0.00567728
12.48
3.53
0.141534629
-0.005670458
12.49
3.53
0.141477959
-0.005663649
12.50
3.54
0.141421356
-0.005656854
12.51
3.54
0.141364822
-0.005650073
12.52
3.54
0.141308355
-0.005643305
12.53
3.54
0.141251955
-0.00563655
12.54
3.54
0.141195624
-0.00562981
12.55
3.54
0.141139359
-0.005623082
12.56
3.54
0.141083162
-0.005616368
12.57
3.55
0.141027032
-0.005609667
12.58
3.55
0.140970969
-0.00560298
12.59
3.55
0.140914972
-0.005596305
12.60
3.55
0.140859042
-0.005589645
12.61
3.55
0.140803179
-0.005582997
12.62
3.55
0.140747382
-0.005576362
12.63
3.55
0.140691652
-0.005569741
12.64
3.56
0.140635988
-0.005563132
12.65
3.56
0.140580389
-0.005556537
12.66
3.56
0.140524857
-0.005549955
12.67
3.56
0.14046939
-0.005543386
12.68
3.56
0.140413989
-0.005536829
12.69
3.56
0.140358654
-0.005530286
12.70
3.56
0.140303383
-0.005523755
12.71
3.57
0.140248178
-0.005517238
12.72
3.57
0.140193039
-0.005510733
12.73
3.57
0.140137964
-0.005504241
12.74
3.57
0.140082954
-0.005497761
12.75
3.57
0.140028008
-0.005491294
12.76
3.57
0.139973128
-0.00548484
12.77
3.57
0.139918312
-0.005478399
12.78
3.57
0.13986356
-0.00547197
12.79
3.58
0.139808872
-0.005465554
12.80
3.58
0.139754249
-0.00545915
12.81
3.58
0.139699689
-0.005452759
12.82
3.58
0.139645193
-0.00544638
12.83
3.58
0.139590761
-0.005440014
12.84
3.58
0.139536393
-0.00543366
12.85
3.58
0.139482088
-0.005427319
12.86
3.59
0.139427847
-0.005420989
12.87
3.59
0.139373668
-0.005414672
12.88
3.59
0.139319553
-0.005408368
12.89
3.59
0.139265501
-0.005402075
12.90
3.59
0.139211512
-0.005395795
12.91
3.59
0.139157585
-0.005389527
12.92
3.59
0.139103721
-0.005383271
12.93
3.60
0.13904992
-0.005377027
12.94
3.60
0.13899618
-0.005370795
12.95
3.60
0.138942504
-0.005364575
12.96
3.60
0.138888889
-0.005358368
12.97
3.60
0.138835336
-0.005352172
12.98
3.60
0.138781845
-0.005345988
12.99
3.60
0.138728416
-0.005339816
13.00
3.61
0.138675049
-0.005333656
13.01
3.61
0.138621743
-0.005327507
13.02
3.61
0.138568499
-0.005321371
13.03
3.61
0.138515316
-0.005315246
13.04
3.61
0.138462194
-0.005309133
13.05
3.61
0.138409133
-0.005303032
13.06
3.61
0.138356133
-0.005296942
13.07
3.62
0.138303194
-0.005290864
13.08
3.62
0.138250316
-0.005284798
13.09
3.62
0.138197498
-0.005278743
13.10
3.62
0.138144741
-0.0052727
13.11
3.62
0.138092044
-0.005266668
13.12
3.62
0.138039408
-0.005260648
13.13
3.62
0.137986831
-0.005254639
13.14
3.62
0.137934315
-0.005248642
13.15
3.63
0.137881858
-0.005242656
13.16
3.63
0.137829462
-0.005236682
13.17
3.63
0.137777125
-0.005230718
13.18
3.63
0.137724847
-0.005224767
13.19
3.63
0.137672629
-0.005218826
13.20
3.63
0.137620471
-0.005212897
13.21
3.63
0.137568371
-0.005206978
13.22
3.64
0.137516331
-0.005201072
13.23
3.64
0.13746435
-0.005195176
13.24
3.64
0.137412427
-0.005189291
13.25
3.64
0.137360564
-0.005183418
13.26
3.64
0.137308759
-0.005177555
13.27
3.64
0.137257013
-0.005171704
13.28
3.64
0.137205325
-0.005165863
13.29
3.65
0.137153696
-0.005160034
13.30
3.65
0.137102124
-0.005154215
13.31
3.65
0.137050611
-0.005148408
13.32
3.65
0.136999156
-0.005142611
13.33
3.65
0.136947759
-0.005136825
13.34
3.65
0.13689642
-0.00513105
13.35
3.65
0.136845138
-0.005125286
13.36
3.66
0.136793914
-0.005119533
13.37
3.66
0.136742747
-0.00511379
13.38
3.66
0.136691638
-0.005108058
13.39
3.66
0.136640586
-0.005102337
13.40
3.66
0.136589591
-0.005096627
13.41
3.66
0.136538653
-0.005090927
13.42
3.66
0.136487773
-0.005085237
13.43
3.66
0.136436949
-0.005079559
13.44
3.67
0.136386181
-0.005073891
13.45
3.67
0.136335471
-0.005068233
13.46
3.67
0.136284817
-0.005062586
13.47
3.67
0.136234219
-0.005056949
13.48
3.67
0.136183678
-0.005051323
13.49
3.67
0.136133193
-0.005045708
13.50
3.67
0.136082763
-0.005040102
13.51
3.68
0.13603239
-0.005034507
13.52
3.68
0.135982073
-0.005028923
13.53
3.68
0.135931812
-0.005023349
13.54
3.68
0.135881606
-0.005017785
13.55
3.68
0.135831456
-0.005012231
13.56
3.68
0.135781362
-0.005006687
13.57
3.68
0.135731322
-0.005001154
13.58
3.69
0.135681339
-0.004995631
13.59
3.69
0.13563141
-0.004990118
13.60
3.69
0.135581536
-0.004984615
13.61
3.69
0.135531717
-0.004979123
13.62
3.69
0.135481954
-0.00497364
13.63
3.69
0.135432245
-0.004968167
13.64
3.69
0.13538259
-0.004962705
13.65
3.69
0.13533299
-0.004957252
13.66
3.70
0.135283445
-0.00495181
13.67
3.70
0.135233954
-0.004946377
13.68
3.70
0.135184518
-0.004940955
13.69
3.70
0.135135135
-0.004935542
13.70
3.70
0.135085807
-0.004930139
13.71
3.70
0.135036532
-0.004924746
13.72
3.70
0.134987312
-0.004919363
13.73
3.71
0.134938145
-0.004913989
13.74
3.71
0.134889032
-0.004908626
13.75
3.71
0.134839972
-0.004903272
13.76
3.71
0.134790967
-0.004897928
13.77
3.71
0.134742014
-0.004892593
13.78
3.71
0.134693115
-0.004887268
13.79
3.71
0.134644269
-0.004881953
13.80
3.71
0.134595476
-0.004876648
13.81
3.72
0.134546736
-0.004871352
13.82
3.72
0.134498048
-0.004866065
13.83
3.72
0.134449414
-0.004860789
13.84
3.72
0.134400833
-0.004855521
13.85
3.72
0.134352304
-0.004850264
13.86
3.72
0.134303827
-0.004845015
13.87
3.72
0.134255403
-0.004839777
13.88
3.73
0.134207032
-0.004834547
13.89
3.73
0.134158712
-0.004829327
13.90
3.73
0.134110445
-0.004824117
13.91
3.73
0.13406223
-0.004818916
13.92
3.73
0.134014067
-0.004813724
13.93
3.73
0.133965956
-0.004808541
13.94
3.73
0.133917896
-0.004803368
13.95
3.73
0.133869888
-0.004798204
13.96
3.74
0.133821932
-0.004793049
13.97
3.74
0.133774027
-0.004787904
13.98
3.74
0.133726174
-0.004782767
13.99
3.74
0.133678372
-0.00477764
14.00
3.74
0.133630621
-0.004772522
14.01
3.74
0.133582921
-0.004767413
14.02
3.74
0.133535273
-0.004762314
14.03
3.75
0.133487675
-0.004757223
14.04
3.75
0.133440128
-0.004752141
14.05
3.75
0.133392632
-0.004747069
14.06
3.75
0.133345187
-0.004742005
14.07
3.75
0.133297792
-0.004736951
14.08
3.75
0.133250448
-0.004731905
14.09
3.75
0.133203154
-0.004726868
14.10
3.75
0.13315591
-0.004721841
14.11
3.76
0.133108717
-0.004716822
14.12
3.76
0.133061574
-0.004711812
14.13
3.76
0.133014481
-0.004706811
14.14
3.76
0.132967438
-0.004701819
14.15
3.76
0.132920444
-0.004696835
14.16
3.76
0.132873501
-0.004691861
14.17
3.76
0.132826607
-0.004686895
14.18
3.77
0.132779763
-0.004681938
14.19
3.77
0.132732968
-0.00467699
14.20
3.77
0.132686223
-0.00467205
14.21
3.77
0.132639527
-0.004667119
14.22
3.77
0.132592881
-0.004662197
14.23
3.77
0.132546283
-0.004657283
14.24
3.77
0.132499735
-0.004652378
14.25
3.77
0.132453236
-0.004647482
14.26
3.78
0.132406785
-0.004642594
14.27
3.78
0.132360384
-0.004637715
14.28
3.78
0.132314031
-0.004632844
14.29
3.78
0.132267727
-0.004627982
14.30
3.78
0.132221471
-0.004623128
14.31
3.78
0.132175264
-0.004618283
14.32
3.78
0.132129106
-0.004613446
14.33
3.79
0.132082995
-0.004608618
14.34
3.79
0.132036933
-0.004603798
14.35
3.79
0.131990919
-0.004598987
14.36
3.79
0.131944953
-0.004594184
14.37
3.79
0.131899036
-0.004589389
14.38
3.79
0.131853166
-0.004584602
14.39
3.79
0.131807344
-0.004579824
14.40
3.79
0.131761569
-0.004575054
14.41
3.80
0.131715842
-0.004570293
14.42
3.80
0.131670163
-0.00456554
14.43
3.80
0.131624532
-0.004560795
14.44
3.80
0.131578947
-0.004556058
14.45
3.80
0.13153341
-0.004551329
14.46
3.80
0.131487921
-0.004546609
14.47
3.80
0.131442478
-0.004541896
14.48
3.81
0.131397083
-0.004537192
14.49
3.81
0.131351734
-0.004532496
14.50
3.81
0.131306433
-0.004527808
14.51
3.81
0.131261178
-0.004523128
14.52
3.81
0.13121597
-0.004518456
14.53
3.81
0.131170809
-0.004513792
14.54
3.81
0.131125694
-0.004509137
14.55
3.81
0.131080626
-0.004504489
14.56
3.82
0.131035605
-0.004499849
14.57
3.82
0.130990629
-0.004495217
14.58
3.82
0.1309457
-0.004490593
14.59
3.82
0.130900817
-0.004485977
14.60
3.82
0.130855981
-0.004481369
14.61
3.82
0.13081119
-0.004476769
14.62
3.82
0.130766445
-0.004472177
14.63
3.82
0.130721746
-0.004467592
14.64
3.83
0.130677093
-0.004463015
14.65
3.83
0.130632486
-0.004458447
14.66
3.83
0.130587924
-0.004453886
14.67
3.83
0.130543408
-0.004449332
14.68
3.83
0.130498938
-0.004444787
14.69
3.83
0.130454513
-0.004440249
14.70
3.83
0.130410133
-0.004435719
14.71
3.84
0.130365798
-0.004431196
14.72
3.84
0.130321509
-0.004426682
14.73
3.84
0.130277265
-0.004422175
14.74
3.84
0.130233065
-0.004417675
14.75
3.84
0.130188911
-0.004413183
14.76
3.84
0.130144802
-0.004408699
14.77
3.84
0.130100737
-0.004404223
14.78
3.84
0.130056717
-0.004399754
14.79
3.85
0.130012742
-0.004395292
14.80
3.85
0.129968811
-0.004390838
14.81
3.85
0.129924925
-0.004386392
14.82
3.85
0.129881083
-0.004381953
14.83
3.85
0.129837286
-0.004377521
14.84
3.85
0.129793533
-0.004373097
14.85
3.85
0.129749824
-0.004368681
14.86
3.85
0.129706159
-0.004364272
14.87
3.86
0.129662539
-0.00435987
14.88
3.86
0.129618962
-0.004355476
14.89
3.86
0.129575429
-0.004351089
14.90
3.86
0.12953194
-0.004346709
14.91
3.86
0.129488495
-0.004342337
14.92
3.86
0.129445093
-0.004337972
14.93
3.86
0.129401735
-0.004333615
14.94
3.87
0.129358421
-0.004329264
14.95
3.87
0.12931515
-0.004324921
14.96
3.87
0.129271922
-0.004320586
14.97
3.87
0.129228738
-0.004316257
14.98
3.87
0.129185597
-0.004311936
14.99
3.87
0.1291425
-0.004307622
15.00
3.87
0.129099445
-0.004303315
15.01
3.87
0.129056433
-0.004299015
15.02
3.88
0.129013465
-0.004294723
15.03
3.88
0.128970539
-0.004290437
15.04
3.88
0.128927656
-0.004286159
15.05
3.88
0.128884816
-0.004281888
15.06
3.88
0.128842018
-0.004277623
15.07
3.88
0.128799263
-0.004273366
15.08
3.88
0.128756551
-0.004269116
15.09
3.88
0.128713881
-0.004264873
15.10
3.89
0.128671253
-0.004260638
15.11
3.89
0.128628668
-0.004256409
15.12
3.89
0.128586125
-0.004252187
15.13
3.89
0.128543624
-0.004247972
15.14
3.89
0.128501166
-0.004243764
15.15
3.89
0.128458749
-0.004239563
15.16
3.89
0.128416374
-0.004235369
15.17
3.89
0.128374041
-0.004231181
15.18
3.90
0.128331751
-0.004227001
15.19
3.90
0.128289501
-0.004222828
15.20
3.90
0.128247294
-0.004218661
15.21
3.90
0.128205128
-0.004214501
15.22
3.90
0.128163004
-0.004210348
15.23
3.90
0.128120921
-0.004206202
15.24
3.90
0.12807888
-0.004202063
15.25
3.91
0.12803688
-0.00419793
15.26
3.91
0.127994921
-0.004193805
15.27
3.91
0.127953004
-0.004189686
15.28
3.91
0.127911128
-0.004185574
15.29
3.91
0.127869292
-0.004181468
15.30
3.91
0.127827498
-0.004177369
15.31
3.91
0.127785745
-0.004173277
15.32
3.91
0.127744033
-0.004169192
15.33
3.92
0.127702361
-0.004165113
15.34
3.92
0.12766073
-0.004161041
15.35
3.92
0.12761914
-0.004156975
15.36
3.92
0.127577591
-0.004152916
15.37
3.92
0.127536082
-0.004148864
15.38
3.92
0.127494613
-0.004144818
15.39
3.92
0.127453185
-0.004140779
15.40
3.92
0.127411798
-0.004136747
15.41
3.93
0.127370451
-0.004132721
15.42
3.93
0.127329143
-0.004128701
15.43
3.93
0.127287876
-0.004124688
15.44
3.93
0.12724665
-0.004120682
15.45
3.93
0.127205463
-0.004116682
15.46
3.93
0.127164316
-0.004112688
15.47
3.93
0.127123209
-0.004108701
15.48
3.93
0.127082142
-0.00410472
15.49
3.94
0.127041115
-0.004100746
15.50
3.94
0.127000127
-0.004096778
Kernel density thresholds identification
Now, for the remaining part of the main methodology of this paper, using the kernel polynomial density function presented above, a density value for each cell id
of the grid is estimated, allowing the generation of an approximately continuous kernel density distribution, as displayed in Figure 6, for three different clusters across four CMAs. In Figure 6, KDE values appear on the horizontal axis and the frequency count in thousands (unweighted number of cells from the full CMA grid) is on the vertical axis. The higher the KDE value, the greater the density of the cell (i.e., the larger the number of job locations available in the cell's neighborhood). As expected, the KDE density distribution is highly skewed. Interestingly enough, the empirical kernel density distribution below, is similar in shape with the theoretical first derivative function of the dynamic between spatial density and spatial distance between points, described in the previous page above (Figure 5). For all CMAs and clusters, the cell counts drop significantly when the KDE exceeds a value of 1 (dashed vertical line on Figure 6). Important drops are also available for even smaller KDE values, but we ignore those since the KDE values are getting near a density of zero. Given the pattern shown in Figure 6, a threshold of 1 was used as a pre-processing step to filter out the unnecessary segment of the grid cell distribution. However, the analysis was further refined with a second set of thresholds, as described below.
Description for Figure 6
Figure 6
Distributions of KDE values in different CMAs and industry clusters
Table summary The information is grouped by "KDE Value" (appearing as row headers), , calculated using (appearing as column headers).
"KDE Value"
"Montreal 313233 Cluster KDE Value Count"
"Montreal 4445 Cluster KDE Value Count"
"Montreal 72 Cluster KDE Value Count"
"Toronto 313233 Cluster KDE Value Count"
"Toronto 4445 Cluster KDE Value Count"
"Toronto 72 Cluster KDE Value Count"
"Winnipeg 313233 Cluster KDE Value Count"
"Winnipeg 4445 Cluster KDE Value Count"
"Winnipeg 72 Cluster KDE Value Count"
"Vancouver 313233 Cluster KDE Value Count"
"Vancouver 4445 Cluster KDE Value Count"
"Vancouver 72 Cluster KDE Value Count"
Source : authors’ computations from the BR database.
0
156,751
293,240
211,432
200,598
410,092
340,608
36,703
80,943
55,391
103,995
184,021
158,367
0.1
102,083
180,251
125,950
122,674
251,485
207,568
24,247
53,142
35,162
62,169
107,399
92,179
0.2
76,101
133,458
91,478
87,194
179,063
146,812
15,721
34,600
23,106
45,196
77,621
65,104
0.3
52,667
88,629
65,788
68,469
140,665
114,955
13,505
29,228
19,045
35,366
59,701
50,716
0.4
58,651
100,673
65,608
59,851
122,506
100,546
10,704
23,012
15,432
31,833
53,742
44,755
0.5
42,187
67,581
49,531
50,421
103,085
86,004
10,516
22,279
14,401
26,560
44,424
37,865
0.6
43,509
71,378
48,530
51,195
103,298
83,601
8,888
19,397
12,694
27,580
45,478
37,072
0.7
42,184
68,125
46,495
41,566
84,359
73,288
8,935
19,172
12,255
21,283
35,849
31,072
0.8
37,865
59,281
41,639
45,903
92,502
74,127
8,217
17,506
11,331
24,119
39,290
31,703
0.9
39,152
61,164
41,211
38,956
78,832
67,396
7,499
16,337
10,461
20,618
33,796
28,750
1
31,429
48,460
34,012
31,396
63,407
57,034
6,277
13,279
8,991
16,434
27,470
24,713
1.1
21,840
34,211
25,533
22,395
46,648
45,542
4,588
9,532
6,633
11,949
20,383
19,775
1.2
20,035
31,162
24,088
20,489
42,922
42,154
4,190
8,856
6,271
11,078
18,779
18,340
1.3
18,989
29,143
22,441
19,069
40,337
39,172
4,000
8,128
5,699
10,179
17,294
17,164
1.4
17,776
27,154
21,136
18,241
37,767
37,512
3,779
7,657
5,422
9,523
16,241
16,128
1.5
16,948
25,462
19,857
16,916
35,251
35,390
3,467
7,375
5,125
8,998
15,047
15,080
1.6
15,673
23,988
18,710
15,946
32,610
32,666
3,429
6,853
4,730
8,540
14,377
14,435
1.7
14,885
22,129
17,909
14,804
31,204
31,454
3,142
6,328
4,518
7,835
13,420
13,781
1.8
13,965
20,830
16,666
13,984
28,911
29,342
3,060
6,049
4,326
7,342
12,384
12,707
1.9
12,878
19,143
15,619
12,727
26,451
27,291
2,737
5,537
4,118
6,831
11,181
11,862
2
11,721
17,017
14,300
11,464
23,603
24,885
2,672
4,884
3,676
6,087
10,072
11,095
2.1
11,052
15,690
13,497
10,994
22,138
23,801
2,498
4,489
3,494
5,776
9,578
10,331
2.2
10,417
14,872
12,968
10,278
20,763
22,350
2,276
4,422
3,335
5,577
9,023
9,998
2.3
9,824
14,370
12,174
9,751
19,984
21,602
2,213
4,200
3,137
5,298
8,614
9,378
2.4
9,336
13,555
11,677
9,486
18,812
20,509
2,046
4,075
3,095
5,115
8,037
8,987
2.5
8,997
12,819
11,230
9,016
17,925
19,366
1,975
3,832
2,915
4,933
7,694
8,481
2.6
8,607
12,200
10,593
8,949
17,325
18,690
1,896
3,597
2,796
4,620
7,187
8,155
2.7
8,216
11,541
10,065
8,250
16,521
17,831
1,823
3,451
2,685
4,549
6,871
7,659
2.8
7,747
11,036
9,567
7,926
15,454
16,748
1,827
3,299
2,704
4,163
6,703
7,392
2.9
7,464
10,274
9,097
7,678
14,665
16,092
1,554
3,022
2,537
4,094
6,349
6,961
3
7,122
9,936
8,971
7,273
13,676
15,277
1,576
2,932
2,373
3,940
6,002
6,846
3.1
7,000
9,446
8,291
7,136
13,148
14,436
1,519
2,967
2,437
3,773
5,590
6,491
3.2
6,599
9,162
8,292
6,757
12,911
13,896
1,495
2,681
2,230
3,512
5,436
6,467
3.3
6,371
8,684
7,803
6,596
12,107
13,632
1,427
2,575
2,259
3,406
5,235
6,018
3.4
5,999
8,257
7,518
6,372
11,760
12,887
1,401
2,564
2,112
3,276
5,011
5,851
3.5
5,929
8,046
7,065
6,092
11,334
12,539
1,306
2,423
2,121
3,205
4,829
5,747
3.6
5,699
7,689
6,979
5,905
10,856
11,866
1,278
2,396
2,014
3,026
4,767
5,541
3.7
5,678
7,479
6,599
5,729
10,520
11,288
1,251
2,173
1,960
2,925
4,470
5,333
3.8
5,393
7,100
6,410
5,422
9,934
10,974
1,176
2,145
1,888
2,908
4,458
5,196
3.9
5,310
6,980
6,092
5,218
9,753
10,576
1,170
2,005
1,868
2,660
4,154
4,988
4
5,042
6,680
5,970
5,171
9,300
10,176
1,182
1,918
1,697
2,807
4,084
4,874
4.1
4,906
6,511
5,731
4,990
9,042
10,127
1,124
1,821
1,645
2,557
3,953
4,745
4.2
4,858
6,317
5,287
4,885
8,739
9,604
1,091
1,774
1,625
2,458
3,819
4,613
4.3
4,631
6,099
5,294
4,719
8,480
9,197
999
1,699
1,607
2,426
3,626
4,426
4.4
4,478
5,981
5,145
4,629
8,175
9,046
1,026
1,698
1,646
2,353
3,531
4,220
4.5
4,429
5,737
5,134
4,633
8,053
8,690
982
1,636
1,496
2,319
3,477
4,158
4.6
4,449
5,470
4,893
4,347
7,519
8,501
989
1,550
1,426
2,109
3,321
3,975
4.7
4,220
5,422
4,699
4,199
7,426
8,178
968
1,524
1,435
2,111
3,262
3,995
4.8
4,007
5,096
4,527
4,164
7,185
7,786
904
1,456
1,434
2,084
3,188
3,907
4.9
3,998
4,894
4,310
4,082
6,954
7,549
877
1,367
1,321
2,005
2,998
3,696
5
3,823
4,906
4,318
3,886
6,790
7,398
833
1,310
1,333
1,881
2,894
3,643
Using the density distributions displayed in Figure 6, a set of KDE threshold values was tested for each combination of CMA and industry sector by computing the number of DBs, employees, and establishments retained in the clusters at each KDE threshold level. The objective was to iteratively augment the threshold value and remove low-density DBs, and then stop before a significant drop in total employment within the cluster was observed. Ultimately, the KDE threshold values that retained a minimum of about 80% of total CMA employees were applied for most location/industry combinations, which essentially filtered out many DBs populated with small businesses. The KDE threshold values, ranging for the most part between 1 and 3, are shown in Table 2 for each combination of CMA and cluster.
As we can see in Table 2, a large proportion of thresholds preserve the initial value of 1 (Figure 6) and doesn’t need an adjustment. This is true especially for the industry clusters defined by Delgado et al. (2014), with one exception. This is an indication that the initial threshold is useful and robust. This result also suggests that the KDE model in this research succeeds to understand the notion of small business neighborhood in Canada, without being ingested with any formal notion of a 3-category classification label system, that is, small, medium, and large business size, beforehand. Indeed, the BR distribution of spatial data point is continuous and not multinomial. However, 7 out of 24 thresholds were set to a higher KDE cut-off value of 3. Further investigation of these clusters suggests that higher values were more suitable for representing clusters in CMAs with a higher proportion of small businesses, leading to a more spatially dispersed distribution of employment. In contrast, the threshold value for Hospitality and Tourism in Winnipeg was set at 0.1, indicating a very low concentration of small businesses in this CMA (Table 2). Finally, it is worth noting that the initial setting of thresholds (around value 1) is purely statistics-based and relies only on the shape of the KDE statistical distributions, while the second set of thresholds (Table 2) is driven by economic considerations and the potential for generating meaningful economic indicators from the cluster boundaries. Hence, this composite cut-off setting benefits from a multi-dimensional perspective and is more likely to be stable over time and space. The next three sections finalize the remaining steps of the methodology.
Table 2
KDE thresholds applied to CMA/industry cluster combinations
Table summary The information is grouped by Industry cluster (appearing as row headers), Montreal, Toronto, Winnipeg and Vancouver, calculated using KDE threshold value units of measure (appearing as column headers).
Industry cluster
Montreal
Toronto
Winnipeg
Vancouver
KDE threshold value
Source: authors’ computations from the BR database.
Manufacturing Sector
3
1
1
1
Retail Trade Sector
3
3
1
3
Accommodations and Food Services Sector
3
3
1
3
Distribution and Electronic Commerce (cluster 10)
1
1
1
1
Financial Services (cluster 16)
1
1
1
1
Hospitality and Tourism (cluster 22)
1
1
0.1
1
Generalized KDE
By definition, the KDE results identify areas within a CMA that contain the majority of a specific industry type. Although the filtered results provide an indication of where a specific industry is dominant within a study area, they are not uniform, and small gaps or holes can appear within identified high-density areas. To smooth out these results, a generalization step is applied to the filtered KDE outputs.
Description for Map 1
Maps showing the difference between kernel density estimation and generalized kernel density estimation. First map shows an example of kernel density estimation. Red color represents grid cells combined into a polygon feature. An info box is located at the top right corner of the map. This includes the name of the method displayed in the map, the short description of the method displayed in the map, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83). Second map shows application of generalized kernel density estimation to the first map. Red color represents buffering of combined grid cells by 50 meters then de-buffering. An info box is located at the top right corner of the map. This includes the name of the method displayed in the map, the short description of the method displayed in the map, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Sources: Statistics Canada, 2021 Census – Dissemination blocks boundary file, 2021 Census – Population centres boundary file, and authors’ computations.
The generalization step begins with the union of all KDE result grid cells into large polygon features. Once the KDE result cells are combined, the polygons are generalized by applying a technique of buffering the union results by 50 meters and then de-buffering by the same amount. This process removes small gaps and holes in the polygons and produces a cleaner output product for conflation. An example of this process is shown in Map 1 above.
Conflation of KDE results to DBs
The results obtained through KDE and the selection of density thresholds provided an indication of where concentrations of industry clusters exist within each CMA. However, these results were generated at the grid output cell level, and thus lacked association with established boundary files used by Statistics Canada. By re-associating the results with DB boundaries, a more comprehensive analysis of industry clusters can be performed, which is not limited to merely identifying the presence or absence of an industry.
Conflating the KDE results was accomplished by intersecting DB centroids of a CMA with the generalized KDE results. All DBs with centroids that intersected the KDE results were retained to represent industry clusters. In other words, the grid cell overlapping with the DB centroid defines the representative density of the DB. An example of the output from this process (Map 2), illustrates the DB polygons associated with an industry cluster through centroid intersection.
It is worth noting that a weighted DB centroid, based on the BR spatial location of establishment within the DB, could have been used instead of a standard geometric centroid. This is the case, for example, for the Spatial Access Measures Open Database (SAM), engineered at the Data Exploration and Integration Lab (DEIL) of the Center for Special Business Projects (CSBP) of Statistics Canada. However, the methodology of our paper makes use of a uniform randomization process for job location data points within the DB, which, rejects the pertinence of preserving the BR establishment locations within the DB, and validates the use of a simple geometric centroid. As explained earlier in this paper, the expected event
of the randomization process RP is a spatially uniform allocation of data points within the DB where each spatial spot gets the same number of data points. Consequently, there is no point in prioritizing a weighted centroid where job locations are likely to be located. Indeed, even though neighborhood accuracy matters, it also must be partially compromised because the objective remains to smooth the discrete spatial distribution of job locations and facilitate the contribution done by the kernel density estimation.
Description for Map 2
Maps showing DB polygons associated with an industry cluster using generalized KDE. First map shows an example of generalized kernel density estimation. Red color represents generalized KDE grid output cell overlayed DB boundaries. An info box is located at the top right corner of the map. This includes the name of the method displayed in the map, the short description of the method displayed in the map, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83). Second map shows DB polygons associated with an industry cluster. Red color represents generalized KDE conflated to DB. An info box is located at the top right corner of the map. This includes the name of the method displayed in the map, the short description of the method displayed in the map, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Clustering and filtering of conflated results
The results conflated to DB spatial boundaries illustrate which DBs are associated with an industry. However, due to the spatial distribution of some industries, and the irregular shapes of DB polygons, the results of this conflation can generate small groupings of DBs that may be associated with only a few industry establishments. To ensure that the final results focus on the main concentrations of businesses and preserve confidentiality in further analysis, the outputs of the conflation process were grouped into clusters of connected DB polygons, and summary statistics for each polygon were computed.
Clustering of the conflated DB polygons was performed by combining all edge touching DBs. Polygons that only touch at a corner were not considered part of the cluster as shown with the clustering example in Map 3. This rule was implemented to limit the creation of large sprawling industry clusters that have locations scattered across a CMA. With DB clusters identified, counts of employees, establishments, and DB polygons were tabulated for each cluster to assess whether the cluster should be retained.
Description for Map 3
Maps showing industry clustering of conflated DB polygons. First map shows an example of conflated DBs. Red color represents connected conflated DBs. An info box is located at the top right corner of the map. This includes the name of the method displayed in the map, the short description of the method displayed in the map, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83). Second map shows clustering of DBs by combining all edge touching DBs into a cluster. Red represents DBs part of cluster 1, yellow represents DBs part of cluster 2, and purple represents DBs part of cluster 3. An info box is located at the top right corner of the map. This includes the name of the method displayed in the map, the short description of the method displayed in the map, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Using the tabulated cluster results from the previous step, clusters that present a confidentiality risk due to being dominated by a single establishment are filtered out. This was done by removing any clusters that had too few establishments (less than 5) or where too many employees (+80%) by a single business establishment in a cluster. This process highlights the idea that a small cluster isolated from its own CMA largest clusters is not necessarily acknowledged as a confidentiality concern if both confidentiality conditions are met.
Recapitulation of steps from the perspective of a grid cell
Before moving to the next section of this paper, and to stay consistent with our notation, we recapitulate the recent few steps from the perspective of a grid cell of interest
. is initially matched with an estimated density numerical value
included in the non-negative real numbers
(Figure 6). If this value
is below its CMA and industry sector threshold (either a threshold value of 1 or 3 based on table 2), then
is excluded from the rest of the methodology. If
is equal to or above its respective threshold, then
is preserved for the next step. Also, at this point, the
representation of
doesn’t matter anymore, and
takes an arbitrary value
equal to the same arbitrary value of all other grid cell
included in the process of the CMA so far. For the generalization step, if
with value
is included in the process so far, then the generalization also preserves the existence of
and cannot exclude it. That is, the buffering and de-buffering step smooths cluster frontiers but cannot suppress existing cells of a frontier. However, grid cells that have been excluded so far can now become included as the generalization step is about removing small gaps and holes in the cluster shape. For the conflation to DBs,
is either overlapping with the centroid of its respective DB or not. If not overlapping, then
has no purpose anymore in the process. If overlapping, then the purpose of
is to accept its respective DB into the DB-level representation of the industry cluster heat mapping. Note that we now notate the DB of
as
because this explanation is from the perspective of
and the DB of
is included in the process so far because of the existence of . At this point, the role of grid cell
is over, and the rest of the decisions are cluster-based, and not cell-based or DB-based. For the final filtering step, the double confidentiality condition will either preserve or suppress the industry cluster where
and grid cell
are so far included. If the cluster is suppressed, then the contribution of grid cell , no matter how significant in the previous steps, now becomes zero because is no longer part of a cluster. If the cluster is preserved, then the grid cell
preserves its contribution because
is the cell whose role was to add within a cluster part of the final heat map.
Non-technical summary of the methodology
This methodology section concludes with a non-technical summary of the steps for the generation of the industry clusters. We decompose the approach into 4 main steps: 1) data, 2) kernel, 3) thresholds, and 4) post kernel processing.
Data:
Extract establishment level data from internal Statistics Canada BR database
Define location of each unique job part of an establishment (Randomize spatially all job locations of a DB within its own DB if the establishment is located within the DB)
Define weight of each job (in our case, a weight of 1 for each job)
Kernel:
Define grid cell shape and dimension for kernel
Define functional form for kernel
Define bandwidth length for kernel
Define geometric centroid of each cell of the grid within CMA
Compute a unique density value for each cell of the grid within CMA using distance between cell geometric centroid, and job location
Thresholds:
Calculate cut-offs based on statistical distribution of KDE for each CMA and industry
Calculate cut-offs based on retention ratio approach for each CMA and industry
Prioritize cut-offs based on retention ratio approach if different from cut-offs based on statistical distribution of KDE
Post kernel processing:
Generalize KDE by smoothing cluster boundaries and filling cluster inner gaps
Conflate KDE results to DBs to get a final DB representation of the clusters
Filter out clusters if not meeting the double condition confidentiality requirement
Results
Figure 7 presents an overview of the descriptive statistics for the cluster results. Number of DBs, establishments and employees per cluster are available for each of the 4 CMA study areas. The x-axis is for the count of DBs, establishments and employees, respectively. The y-axis is the frequency or the count of distinct clusters. Statistical distributions are similar in shape with the highly right-skewed KDE distribution of Figure 6. That is, similar to a power law distribution, which is in line with (Gabaix, 1999) that city size is power law distributed. Furthermore, for all 3 histograms of Figure 7, the frequency is more conservative for Montreal (blue highlighted histograms), which happens to be the CMA with the smallest and lowest range of DB superficies of Figure 1 and one of the lowest KDE distributions of Figure 6 with manufacturing NAICS 31-32-33. Based on these empirical observations, the number of clusters per CMA seems to be correlated with the DB population structure of the CMA. This is an intuitive result. That is, CMAs with a high proportion of large DBs would contribute to generating sparser segments of high-density grid output cells, and consequently produce a larger number of distinct spatial clusters. Alternatively, CMAs largely dominated with small squared DBs would tend to aggregate several clusters into a single large cluster during the generalization, buffering-de-buffering, and conflation process, due to the proximity of DBs to each other. This raises some interesting research questions, and we take the time to describe two here. 1) Does the DB configuration of a CMA explain a significant proportion of the population of clusters? 2) Are there other more important factors affecting the way the clusters and shapes emerge from the spatial data?
Description for Figure 7
Figure 7
Distribution of number of DBs, establishments, and employees for each CMA study area
Table summary The information is grouped by "DB Per Cluster (Range)" (appearing as row headers), , calculated using (appearing as column headers).
"DB Per Cluster (Range)"
"Montreal CMA Cluster Count"
"Toronto CMA Cluster Count"
"Winnipeg CMA Cluster Count"
"Vancouver CMA Cluster Count"
Source : authors’ computations from the BR database.
0-4
26
75
2
41
5-9
94
149
7
53
10-14
94
92
9
49
15-19
45
55
9
36
20-24
42
45
4
26
25-29
25
38
5
11
30-34
20
37
4
6
35-39
15
16
1
5
40-44
6
15
2
4
45-49
12
5
2
2
50-54
6
3
2
8
55-59
8
11
2
4
60-64
2
11
1
2
65-69
3
5
1
3
70-74
3
6
1
1
75-79
3
2
NA
1
80-84
2
4
NA
4
85-89
1
5
NA
3
90-94
NA
4
2
2
95-99
2
4
NA
1
100-104
NA
1
NA
NA
1050-1099
3
3
NA
2
1100-1149
3
3
NA
1
1150-1199
1
4
NA
3
1200-1249
2
3
NA
2
1250-1299
1
3
NA
1
1300-1349
2
2
NA
1
1350-1399
1
3
NA
1
1400-1449
2
2
1
1
1450-1499
NA
1
NA
1
1500-1549
1
3
NA
NA
"Employees Per Cluster (Range)"
0-49
25
98
9
14
50-99
74
89
14
52
100-149
52
59
6
30
150-199
30
47
2
23
200-249
26
41
3
18
250-299
23
24
2
16
300-349
19
18
2
8
350-399
17
19
2
7
400-449
13
16
1
4
450-499
9
8
2
7
500-549
4
12
1
5
550-599
8
12
1
4
600-649
8
7
NA
2
650-699
4
3
NA
6
700-749
7
8
NA
5
750-799
6
7
NA
6
800-849
10
9
1
5
850-899
7
2
NA
1
900-949
5
6
NA
4
950-999
3
3
1
3
1000-1049
3
6
NA
2
1050-1099
3
3
NA
2
1100-1149
3
3
NA
1
1150-1199
1
4
NA
3
1200-1249
2
3
NA
2
1250-1299
1
3
NA
1
1300-1349
2
2
NA
1
1350-1399
1
3
NA
1
1400-1449
2
2
1
1
1450-1499
NA
1
NA
1
1500-1549
1
3
NA
NA
"Establishments Per Cluster (Range)"
5-9
186
249
34
104
10-14
61
107
2
51
15-19
42
49
6
24
20-24
28
30
3
7
25-29
21
22
1
8
30-34
17
22
NA
8
35-39
12
15
3
5
40-44
9
10
NA
8
45-49
7
7
2
3
50-54
2
11
1
5
55-59
1
8
NA
3
60-64
5
6
1
5
65-69
3
4
NA
7
70-74
3
4
NA
2
75-79
1
6
1
2
80-84
1
5
NA
3
85-89
1
5
NA
2
90-94
1
3
NA
NA
95-99
2
2
NA
1
100-104
NA
2
NA
2
105-109
NA
1
NA
1
115-119
1
1
NA
NA
120-124
1
1
NA
1
135-139
NA
2
NA
1
140-144
1
2
NA
NA
145-149
NA
2
NA
NA
Key results of the spatial clusters for Montreal, Toronto, Winnipeg, and Vancouver along with six specifications of industry clusters are summarized in Table 3, while the mapping of each cluster is reported in the Appendix. In each map, the clusters for the reference industry sector are highlighted in red.
As expected, spatial distribution of clusters varies substantially among industries, with the Retail Trade and Accommodations and Food Services sectors generally covering the largest shares of DBs in each metropolitan area. This is also reflected in the count of DBs that belong to each cluster (Table 3). For all the CMAs, Retail Trade and Accommodations and Food Services sectors account for the largest number of DBs. Similarly, the DB retention ratio, that is the percentage of DBs that were included in the corresponding cluster in each CMA, ranges from 36% to 55%, and is generally higher (in the vicinity or above 40%) for Retail Trade and Accommodations and Food Services sectors (Table 3).
The clusters capture most of the establishments and, even more significantly, the majority of employment in the corresponding industries (Table 3). As described in the previous sections, this was a key criterion for determining the density value threshold for the inclusion of DBs in the cluster. The establishment retention ratio represents the share of establishments within the industry cluster DBs relative to the total count of establishments in that industry within the CMA; similarly, the employment ratio represents the share of employment generated by businesses within the industry cluster areas relative to the total employment in that industry within the CMA.
Table 3
Results: cluster summary statistics, 2023
Table summary The information is grouped by CMA Cluster (appearing as row headers), DB, Establishment , DB, Employment and Establishment, calculated using Count and retention ratio units of measure (appearing as column headers).
CMA Cluster
DB
Establishment
DB
Employment
Establishment
Count
retention ratio
Source: authors’ computations from the BR database.
Montréal
Manufacturing Sector
3,775
4,136
26.1
89.7
70.9
Retail Trade Sector
9,861
10,431
38.9
92
76.7
Accommodations and Food Services Sector
7,888
6,806
36.5
85.4
77.6
Distribution and Electronic Commerce (cluster 10)
5,660
5,406
34.4
94.9
74.6
Financial Services (cluster 16)
2,902
1,512
26.2
88.5
55.4
Hospitality and Tourism (cluster 22)
2,472
860
23.6
60.3
47.1
Toronto
Manufacturing Sector
5,732
7,594
34.8
94.4
81.6
Retail Trade Sector
9,716
15,075
35.5
97.1
76
Accommodations and Food Services Sector
9,946
10,085
38.5
89.2
77.2
Distribution and Electronic Commerce (cluster 10)
7,045
10,291
36.2
97.8
81.9
Financial Services (cluster 16)
5,972
4,856
35.6
95.7
72
Hospitality and Tourism (cluster 22)
2,842
1,329
25.3
71.2
50.9
Winnipeg
Manufacturing Sector
1,054
616
44.2
91.9
78.5
Retail Trade Sector
2,700
2,130
55.3
94.6
87.2
Accommodations and Food Services Sector
2,217
1,313
54.9
91.2
87.7
Distribution and Electronic Commerce (cluster 10)
855
850
33.1
93.2
73.5
Financial Services (cluster 16)
666
349
28.4
82.7
55.2
Hospitality and Tourism (cluster 22)
1,363
337
60.3
81.3
75.2
Vancouver
Manufacturing Sector
2,432
3,236
35.9
90.8
79.1
Retail Trade Sector
3,981
7,260
36.7
81.3
77.5
Accommodations and Food Services Sector
4,686
5,624
43.7
91.6
83
Distribution and Electronic Commerce (cluster 10)
2,662
4,610
32
93.7
77.1
Financial Services (cluster 16)
2,166
2,130
32.3
91.2
68.2
Hospitality and Tourism (cluster 22)
1,709
1,041
31.1
79.8
60.1
Except for Hospitality and Tourism (cluster 22), all other industry clusters capture well over 80% of the employment of that industry within their respective industries in the CMAs of reference, with some clusters capturing 95% or more of the employment (Table 3). For instance, the manufacturing sector cluster contains 89.7% of total manufacturing employment in Montreal, 94.4% in Toronto, 91.9% in Winnipeg, and 90.8% in Vancouver.
Similar percentages are computed for counts of businesses that are within the cluster areas. Although these percentages are slightly smaller than those for employment, they remain around 80% for most clusters, except Hospitality and Tourism (cluster 22) and Financial Services (cluster 16). For example, the manufacturing sector clusters contain 70.9% in Montreal, 81.6% in Toronto, 78.5% in Winnipeg, and 79.1% in Vancouver of total manufacturing establishments. This partition of businesses, within the cluster versus those operating outside the cluster, could also be used to monitor trends in the evolution of metropolitan clusters.
It should be noted that, for some clusters, the count of establishments is significantly smaller than the count of DBs comprising the cluster. For instance, the Hospitality and Tourism (cluster 22) cluster in Montreal includes 860 establishments and 2,472 DBs. This result is driven by the buffering and conflation methods used in the analysis and the concentration of businesses in areas with small DBs. The methodological approach developed in this analysis is designed to provide a neighborhood level representation, as opposed to representation of individual DBs. Thus, businesses within DBs that are in proximity to each other, but still separated by other DBs, are clustered together including DBs that are surrounded by such businesses but do not contain any within their boundaries. The opposite situation is also possible. That is, the DB include establishments and job locations but is not part of the final cluster. This can happen, if the grid cell overlapping with our DB centroid is not dense enough and gets filtered out during the 2-stage threshold process explained earlier. Such DB centroid is usually located at the boundaries of a large cluster and get rejected by the process due to the bandwidth aggregating a too large proportion of DB with no establishment location. A quick analogy would be fitting a non-linear regression model over an assembly of data points. The fitted model is sometimes above or below the actual data points, but overall, contributes to provide a good and continuous approximation of the data phenomena, and provides an anonymized representation of the actual spatial distribution of establishments and job locations. The fitted curve enables the analyst to do research about trends and patterns without directly observing the confidential data.
It is worth noting that a higher retention ratio (closer to 100%) for a CMA is not equivalent to a higher level of economic performance of the CMA. Retention ratios are rather closer, in meaning, to the distribution of business sizes, where lower ratios imply larger proportion of small firms within the CMA. The distribution of business sizes is publicly available and not an exclusive information provided by this research but remains a source of validation for the results of this research.
The results obtained with the mapping of business clusters were visually validated by comparing them with land zoning maps of municipalities included in the CMAs of this study. This validation was particularly feasible for manufacturing clusters and industrial zones or industrial parks. For example, for the municipality of Toronto, the industrial zoning of land parcels aligns closely with the location of manufacturing clusters. The use of data from third party sources, such as OSM or Google Maps, was not included in the current analysis; however, it could be considered for future validation efforts.
A few more notes for the understanding of our retention ratio. The retention threshold is based on a simple marginal gain principle. That is, we consider it worthwhile to keep filtering out DBs if the removal is sufficiently larger than the removal of employment. In other words, neatness and visual clarity matters for cluster heat map representation, until it affects economic performance. For cluster 22 Montreal, the retention ratio is 60% of employment. Consequently, choosing 80% of employment was not a motivation because a significant proportion of DB kept dropping at the expense of a small enough proportion of job locations. However, going beyond 60% of employment would not only make the job drop considerably larger than the DB drop, but it would also reach a state where the marginal decrease of jobs is exponential, say, from 60% to 30%, which would be unacceptable. This exponential drop is intuitive and explainable considering the very nature of our statistical BR distributions. As discussed in the previous sections of this report, the BR data are highly skewed with Power Law related shapes. Therefore, the non-flat segment of the distribution includes a chaotic segment, where a very small variation of information leads to a very large drop. Clusters defined over a relatively smaller proportion of establishments are likely to be smaller in cluster size and number of clusters. Investigating cluster 22 for Montreal and Toronto in the annex of this report, the clusters are indeed smaller, more sparse, and less voluminous in number.
Table 4
Percentage of co-location between 2 types of industry cluster for each CMA study area, 2023
Table summary The information is grouped by CMA Cluster (appearing as row headers), % of Co-Location Among Industry Clusters, Manufacturing , Retail Trade and Accommodations and Food , calculated using percent units of measure (appearing as column headers).
CMACluster
% of Co-Location Among Industry Clusters
Manufacturing
Retail Trade
Accommodations and Food
percent
Source: authors’ computations from the BR database.
Montréal
Manufacturing
-
37.3
20.5
Retail Trade
30.6
-
43.6
Accommodations and Food Services
24.6
63.6
-
Toronto
Manufacturing
-
50.1
37.0
Retail Trade
47.5
-
52.7
Accommodations and Food Services
39.9
59.9
-
Winnipeg
Manufacturing
-
47.5
37.4
Retail Trade
24.6
-
46.0
Accommodations and Food Services
25.4
60.4
-
Vancouver
Manufacturing
-
39.5
30.8
Retail Trade
46.6
-
61.1
Accommodations and Food Services
30.3
51.0
-
Finally, a simple co-location analysis of business clusters was performed by overlapping boundary files of two clusters and computing the percentage of area of one cluster (row) that is also included in another cluster (column). Table 4 shows the results for an assessment of co-location of Manufacturing, Retail Trade and Accommodations and Food services. A visualization of this co-location between two sets of clusters is also reported in Map 4.
Two key elements emerge from these examples. First, the percentage of co-location varies between approximately 20% and 60% of the cluster area, depending on the specific clusters. However, the strength of co-location between business clusters provides meaningful economic insights. In particular, across all the CMAs analyzed, the overlap between retail trade, and accommodations and food services is more pronounced compared to either of these two industries and manufacturing. This observation aligns with the general perception that these types of services tend to co-locate within the same geographic area within a metropolitan area.
Second, the type of co-location provides insights on potential differences that the same industry cluster may present across the metropolitan area. As an example, Map 4 shows that the Accommodations and Food Services sectors are concentrated in neighborhoods that also have a high concentration of retail trade, and separately, a high concentration of manufacturing sector neighbourhood. It is likely that these areas of overlap indicate Accommodations and Food service cater to different clienteles and present different types of sector linkages or dependencies.
Finally, it should be noted that the concentration of certain types of businesses, like manufacturing, in areas that are outside of what would appear to be their designated municipal zoning areas, such as industrial zoning, may indicate the concentration of specific business functions in commercial areas. For instance, the clustering of manufacturing businesses in downtown Toronto, in areas designated as commercial zoning, suggests that the manufacturing businesses in this area may be related to headquarters or office functions, as opposed to production establishments.
Description for Map 4
Maps showing co-location of industry clusters in Toronto CMA. First map shows manufacturing, and accommodation and food sectors. Red represents manufacturing sector, blue represents accommodation and food sector, and purple represents the two overlapping sectors. An info box is located at the top right corner of the map. This includes the name of the CMA, the name of the clusters, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83). Second map shows retail trade, and accommodation and food sectors. Green represents retail trade sector, blue represents accommodation and food sector, and purple represents the two overlapping sectors. An info box is located at the top right corner of the map. This includes the name of the CMA, the name of the clusters, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Directions for further research, analysis, and applications
The methodology presented in this paper can be further developed and applied to generate insights on business conditions at the local level. The most immediate development involves scaling the work to all metropolitan areas of Canada. Further extension to medium-sized agglomerations (Census Agglomerations) and rural and small-town areas should also be considered.
Increasing geographic coverage can be done in parallel with further refinement of the NAICS groupings that define each cluster. To develop the methodology, the focus of this paper was on simple 2-digit NAICS codes or pre-existing NAICS groupings, with specific reference to the work of Delgado et al. (2014) on U.S. industry clusters. Although there is validity in continuing with this approach, the use of BR microdata provides the opportunity to implement alternative aggregations of NAICS codes and different digit levels. Custom-specific groupings could be considered. For instance, some of the existing literature has focused on artistic and cultural clusters and neighborhood vitality; hence, custom specific definitions of these types of clusters could be considered for implementation.
Further analysis should consider profiling business cluster performance with additional BR variables. In addition to the 6-digit NAICS code, number of employees, and latitude and longitude position for each business establishment, the BR provides several additional fields of interest, which could be integrated into the analysis for profiling and indexing of business performances. In particular, the use of business revenue, total expense, total assets, and date of birth of the business could be explored to generate aggregate spatial indices. Engineering financial performance indicator ratios from the BR could also be possible. There are several of them; 1) profitability ratio ((total revenue minus total expenses) / total assets), 2) liquidity ratio (current total assets / current total debt), 3) tangibility ratio (total fixed assets / total assets), 4) total sales growth rate. It is important to emphasize that BR data comes with technical challenges. The extremely high level of granularity, both in terms of geography and industry group, would limit the type of information that could be extracted. However, the use of indices (e.g., heat maps) or categorical groups (e.g., simple classifications into high, medium, and low values) could mitigate these issues, while providing valuable insights on performance and trends of business clusters at the neighborhood level. This classification would involve clustering highly skewed statistical BR distributions where most of the observations are located within a small range of low values.
On the other hand, the BR remains a very important dataset for Statistics Canada and Canadians. As noted on the Statistics Canada’s BR official web page, ‘’As a statistical register, it provides listings of units and related attributes required for survey sampling frames, data integration, stratification, and business demographic statistics. The BR is a major pillar of the agency's economic statistics programs, including the Census of Agriculture.’’
The BR is maintained on a continuous basis, with inclusion of new businesses and with business employment and financial information updated at regular intervals. New vintages of the BR, with business counts by employment size, are released twice a year; hence, some statistics of cluster maps could be updated at regular intervals.
Finally, the results generated with neighborhood-level clusters could be combined with other data sources, including both Statistics Canada’s data holdings and alternative data sources from external providers. As an example, cluster mapping at the neighborhood level could be overlaid with proximity measures to services and amenities and spatial access measures to understand how the presence of amenities or accessibility interacts with the clustering of specific businesses.Note Similarly, cluster boundaries can be overlaid with mobility flows or commodity flows at a similar geographic scale. This information may provide insights into the level of economic activity within each cluster. Specifically, the combination of cluster boundary files with data on mobility flows is an example of data integration that may yield significant business insights. In this case, the boundary files would be used as geofences to estimate inward-outward mobility, using mobile device data, or mobility through the road network of cluster areas. This information could be used to estimate or monitor the economic activities in the business clusters. Spatial econometric models would be needed to acknowledge the spatial dependencies across these datasets.
As a final note, we observed a similarity of results between Figures 1, 6, and 10, which seems to suggest that the number of clusters per CMA is correlated with the DB population structure of the CMA. Interesting research question arises: to what extend can the DB configuration of a CMA explain the way clusters emerge and shape, and what would be the remaining factor explaining the phenomena? Furthermore, the current BR year assessed in this research could be compared with previous BR years available in the Statistics Canada’s Database. This analysis would be beyond a simple co-location analysis and exploit the signal located within the space-time BR. Spatial regression and autocorrelation (Moran' I) between several versions of the BR, at neighborhood level, could robustly measure the magnitude of change in the shape of the clusters across years. More specifically, we could identify neighborhoods of a CMA where the current clustered industrial expansion is explained by the spatial historical foundation of other industry clusters. The same for industrial stagnation and contraction at the neighborhood level. If the comparison of several years of the BR is a difficult exercise, then a methodology would need to be formalized following the steps of Statistics Canada methodologies for gross flow statistics and the approaches of sampling across time for the creation of unbiased and consistent statistics.
Also, the current kernel bandwidth acknowledges the DB distribution and the median DB dimension of the CMA but doesn’t consider the dispersion of job location. Therefore, a direct improvement would be a bandwidth including the 2 dimensions of information into the calculation. That is, a composite bandwidth weighting 2 components, our approach and the traditional Silverman’s approach.
Finally, granular spatial heat mapping engineering, like time series analysis, needs to be validated for the presence of random walks and spurious correlation. For time series, the correlation of 2 series with non-stationary statistical moment, such as the average and variance, can lead to correlations that appear strong but are unlikely to persist after simple transformation like a first time-step lag differentiation. The same applies for spatial data. Two granular spatial heat maps can be non-stationary, get a very strong spatial correlation, or in some cases, a very strong overlap or co-location of clusters, and not preserve the strong correlation after simple spatial transformation, such as spatial lag differentiation. Modern methods to identify robust spatial unit-root and spurious correlation can be found in Muller and Watson, 2023, and Hassan, 2012.
Conclusions
This paper presents an approach to developing business clusters at the neighborhood level using BR establishment-level data for four CMAs in Canada. A new approach to define the kernel bandwidth is detailed since the traditional Silverman’s rule bandwidth method fails, in the case of our applications, to directly recognize the configuration of the DB structure within the cities of interest. The generated business clusters fill a data gap on business cluster analysis at a highly granular level of geography, providing a framework that can support decision-making and policy at the local level. It also offers a comparative framework for analyzing trends across metropolitan areas in Canada.
The existing literature suggests that business location choices extend beyond the selection of metropolitan areas; neighborhood characteristics are also relevant determinants of location choices. As a result, business trends can vary greatly across neighborhoods within the same metropolitan areas. In turn, the clustering of businesses in a neighborhood may influence the overall economic prosperity and quality of life in that neighborhood, contributing to either expanding or reducing spatial disparities across the metropolitan area.
The methodology proposed in this analysis draws from existing literature on the identification of central districts. In simple terms, it transforms the discrete and fragmented spatial distribution of job locations (based on the geolocation of establishments) into a relatively smoother or more continuous representation using a fine-level grid to support spatial KDE, and subsequently conflates the results at the DB level. The data for this analysis comes from business records from Statistics Canada’s BR. The model introduced a new bandwidth approach to acknowledge more efficiently the configuration of DB within a CMA and be more robust to outliers present in the highly skewed distribution of BR’s job count per establishment. The results generated through this process are filtered to remove single DBs or small aggregations of DBs that would not meet basic confidentiality thresholds for business data. The results show that the proposed method is effective in capturing the large majority, if not nearly the totality, of employment in the respective industries, while simultaneously filtering out a large proportion of DB where the density of job locations is relatively weaker and less interesting. Moreover, while the focus of this analysis has been on the delineation of the boundaries (hence, the cluster areas are monochrome in the appendix maps), future analytical applications can highlight different trends with the cluster areas by using, for instance, heat maps.
As the demand for a more granular spatial information on business continues to grow, the use of business clusters at neighborhood level, as spatial framework of reference, can support the work of local economic agents and policy stakeholders. Neighborhood-level analysis can support local business associations interested in understanding and monitoring local businesses within a specific neighborhood. In addition, local cluster boundaries could be integrated with other measures of spatial distribution at the DB-level, such as population density and proximity to amenities, to generate more comprehensive analysis of local conditions and development opportunities.
References
Aitkin, M. (2022). Introduction to Statistical Modelling and Inference. The multinomial distribution. 1st Edition. Chapman and Hall/CRC.
Bekar, C., & Lipsey, R. G. (2001). Clusters and economic policy. Canadian Journal of Policy Research, 3(1), 62–70.
Bartle, R. G, & Sherbert, D. R. (2001). Introduction to real analysis. Wiley; 4th edition.
Bathelt, H., & Li, P.-F. (2014). Global cluster networks: Foreign direct investment flows from Canada to China. Journal of Economic Geography, 14(1), 45–71. https://doi.org/10.1093/jeg/lbt020
Campaniaris, C., Hayes, S., Jeffrey, M., & Murray, R. (2011). The applicability of cluster theory to Canada's small and medium-sized apparel companies. Journal of Fashion Marketing and Management: An International Journal, 15(1), 8–26.
Christensen, P., McIntyre, N., & Pikholz, L. (2002). Bridging community and economic development: A strategy for using industry clusters to link neighbourhoods to the regional economy. Shorebank Enterprise Group.
Delgado, M., Porter, M. E., & Stern, S. (2014). Defining clusters of related industries (NBER Working Paper No. 20375). National Bureau of Economic Research. https://www.nber.org/papers/w20375
Gabaix, X. (1999). Zipf's law for cities: An explanation. The Quarterly Journal of Economics, 114(3), 739–766.
Gabaix, X., & Ioannides, Y. (2004). The evolution of city size distributions. In J. V. Henderson & J. F. Thisse (Eds.), Handbook of regional and urban economics (Vol. 4, pp. 2341–2378). Elsevier.
Gabaix, X., Lasry, J., Lions, P., & Moll, B. (2016). The dynamics of inequality. Econometrica, 84(6), 2071–2111.
Grimaldi, R, P,. (2003). Discret and combinatorial mathematics. 5th edition.
Grodach, C., Currid-Halkett, E., Foster, N., & Murdoch III, J. (2014). The location patterns of artistic clusters: A metro- and neighborhood-level analysis. Urban Studies, 51(13), 2822–2843.
Hassan, A. (2012). Stationarity and unit roots in spatial auto-regressive models (Doctoral dissertation, Universidad Nacional de Colombia). Retrieved from https://repositorio.unal.edu.co/bitstream/handle/unal/9882/Stationarity_and_unit_roots_in_spatial_autoregressive_models.pdf?sequence=1&isAllowed=y
Patel, P. C. (2024). Nurturing neighborhoods, cultivating local businesses: The effects of amenities-to-infrastructure spending on new business licenses in Chicago’s wards. Journal of Business Venturing Insights, 21, e00467.
E. Parzen. On the estimation of a probability density and the mode. Ann. Math. Statist., 33:1965–1976, 1962.
Rao, J.N.K., Scott, A.J. (1981). The analysis of categorical data from complex sample surveys: chi-squared tests for goodness of fit and independence in two-way tables. Journal of American statistical association, Vol. 76, No. 374 (Jun., 1981), pp. 221-230
M. Rosenblatt. A central limit theorem and a strong mixing condition. Proc. Nat. Acad. Sci. USA, 42:43–47, 1956.
Sergerie, F., Chastko, K., Saunders, D., & Charbonneau, P. (2021). Defining Canada's downtown neighbourhoods: 2016 boundaries. Statistics Canada, Issue 2021001. Retrieved from https://www150.statcan.gc.ca/n1/pub/91f0015m/91f0015m2021001-eng.pdf
Severini, A. T. (2005). Elements of distribution theory. Cambridge University Press. https://www.cambridge.org/core/books/elements-of-distribution-theory/1091612241C8BD626046E5F1E47D30A3
Silverman, B. W. (1986). Density estimation for statistics and data analysis. Monographs on Statistics and Applied Probability. Available online at: https://ned.ipac.caltech.edu/level5/March02/Silverman/paper.pdf.
Spencer, G. M., Vinodrai, T., Gertler, M. S., & Wolfe, D. A. (2010). Do clusters make a difference? Defining and assessing their economic performance. Regional Studies, 44(6), 697–715.
Taubenböck, H., Klotz, M., Wurm, J., Schmieder, B., Wooster, M., Esch, T., & Dech, S. (2013). Delineation of central business districts in mega city regions using remotely sensed data. Remote Sensing of Environment, 136, 386–401. https://doi.org/10.1016/j.rse.2013.05.019
Toronto Region Board of Trade (2021). Regional centres: A business district analysis. Report. Retrieved from https://bot.com/Resources/Resource-Library/Regional-Centres-A-Business-District-Analysis
Wang, B., & Wen, B. (2021). The spatial distribution of businesses and neighborhoods: What industries match or mismatch what neighborhoods? Habitat International, 117, 102440.
Wheeler, C. H. (2006). Businesses don't just choose a city, they choose a specific neighborhood. Bridges. Federal Reserve Bank of St. Louis. Retrieved from from https://fraser.stlouisfed.org/title/bridges-federal-reserve-bank-st-louis-6272/bridges-600970
Wolfe, D. A., & Gertler, M. S. (2004). Clusters from the inside and out: Local dynamics and global linkages. Urban Studies, 41(5-6), 1071–1093.
Zhang, X. (2019). Building effective clusters and industrial parks. In Célestin Monga & Justin Yifu Lin (Eds.), The Oxford handbook of structural transformation. Oxford Handbooks. https://doi.org/10.1093/oxfordhb/9780198793847.013.13
Appendix 1: proof for the non-necessity of using probabilities of combinations for the comparative analysis of the binomial process
Acknowledging
as probability of head, the full expression is,
The expression is then equivalent to the following, if we substitute with the uniformity of ½,
Simplifying with the probabilities’ subtractions and exponents’ subtractions, we get,
Now, due to the uniformity of the binomial process, we use the probability exponent’s additive property to get,
finally, cancelling out common terms on both side of the inequality, and acknowledging some equalities, we get,
which explains why our comparative analysis focuses on the number of combinations and not the full expression. ■
This proof makes use of a denominator equal to the sum of all probabilities, on both sides of the inequality. This denominator is technically redundant and equal to 1. However, the denominator remains useful to visualize probabilities. That is, if we cancel out all probability terms
due to the uniformity of the random process, then only combinatorial terms remain at the numerator and denominator. The denominator becomes the sum of all possible combinations and the numerator becomes the number of combinations of interest. This provides a probability in itself.
Appendix 2: Justification for applying within DB instead of within DA random allocation of jobs during pre-kernel processing
Sergerie et al, 2021, applies a uniform randomization allocation of employment within the DA boundaries. This strategy is excellent considering the number of jobs within a DA is large and covers a large set of 2-digit NAICS. However, for the case of our applications, we treat a single 2-digit NAICS at a time for the generation of the clusters. Consequently, the number of jobs involved per DA is relatively more limited. The normality approximation theory documented in this paper is conditional to a large number of jobs. That is, if the geographical surface dedicated to jobs allocation is large compared to the number of jobs itself, then reaching out to an accurate normal distribution becomes difficult. This annex explains in detail the reason for processing within DB instead of within DA random allocation of jobs during pre-kernel processing. To do so, we decompose the global random allocation of the DA with the perspective of the several DBs of the same DA. In the next pages of this annex, we are going to present equations y*, y** and y*** before presenting our final explanation.
Let’s define space set
as the symmetric neighborhood around the centroid of the DB, the DB itself and the DA from which the DB belongs to, respectively. A is smaller than its DB and proportional in radius to the size of its DB, respectively. Let’s define supplementary space set
and
as the superficies subtracting the previous smaller superficies. That is,
and . Let’s define
as the finite number of spatial spots in available for random allocation of jobs. Let’s also define
and . Let’s assume that , which is by definition, always the case. Let’s define a share of normalized probability such as
and , which define the level of smoothness or fragmentation. The 3 shares should not be proportional to the relative volume of the space sets but rather proportional to their relative importance. Let’s define uniformity on . That is,
and
are the equal probabilities for all spatial spots located in , respectively. It is trivial to show that the sum of all probabilities is 1 in the sense that the space of probabilities measures is well define. That is,
. If the sum is 1 for the 3
terms, then the sum of the last expression remains 1, no matter the size of . An application allowing a limited amount of spread beyond the DB would be , , and . For room convenience matter, we label the probabilities such as, , , and . Also, , and . Such a random process
is not necessarily uniform anymore and still converges in distribution to the multivariate normal distribution. That is, , if , and in a finite context,
, if
is large. Note is the total job count of the DB of interest (spatial set A') and
is the new -dimension expected vector, that is, , is a row vector and , and , and is the identity matrix of dimension . That is, is a diagonal matrix whose diagonal elements are the items of vector . In other words,
And
For its part, is the new variance-covariance matrix, that is, the full uncertainty quantification can be represented as =
is a finite block matrix whose block diagonal includes 3 sub-matrices, one for , respectively. The dimensions of the 3 sub-matrices are , , and , respectively. Consequently, variance-covariance matrix
is symmetric and has dimension . is also a generalization of , presented in this paper. That is,
reduces to
and is equal to the first two of the 3 sub-matrices in the block diagonal of , if
is non existent, and , and
. Matrix terms below the block diagonal are omitted for clarity matter and also to make the 3 sub-matrices of the block diagonal more obvious. A matrix term below the block diagonal is equal to the term of the inverse coordinate in the above block diagonal.
is slightly more sophisticated than
due to the 3-level spatial structure involving . The new multivariate normal distribution of our converging multinomial random process can be expressedNote in a finite context in the following manner,
What is more? Let’s define ADA as the set of All DB of the dissemination Area. The material presented so far in this annex is about the random allocation of a population of v unique jobs of a DB in its ADA. However, the same should be done for all other DBs of the same ADA because each DB of the ADA will randomly allocate its own v jobs across the DBs of the same ADA. Consequently, we now must consider the joint multivariate normal distribution instead of a single multivariate normal distribution. Fundamentally, it is about a simple concatenation exercise of the material presented previously in this annex. To do so, it is sufficient to keep the existing notation of this annex and simply introduce the index and make it range from , where
is the number of total unique
in the ADA where jobs data points are available. From previous material in this annex, we know that each random process, one for each
, converges in distribution to the multivariate normal distribution. That is,
In a finite context, this is equivalent to say that,
It is also equivalent to say that,
Each of the multivariate normal distribution above refer to a distinct space of probabilities measures. That is, the probability vector and the set of employees (jobs) is distinct. However, the set of spatial spots is the same. Now, making use of the product operator, and assuming the independence of the several distributions, the multivariate normal distributions can be expressed in a single equation and generate the joint multivariate normal distribution , that is,
Looking at the third equality, we can notice that the term
doesn’t make use of the
index. That is, it is used repetitively for each multivariate normal equation because the value doesn’t change across . Looking at the fourth equality (the last equality), the summation is justified from the exponent property under the common term and is used to diminish the length of the full expression. It is worth noting that, the joint distribution
reaches its maximum density if each vector
is equal to their mean vector , respectively. The same way around, reaches its minimum density if each vector
allocates the totality of their resources to the element of the probability vector
where the probability is at its lowest possible, respectively. Furthermore, each input vector
of the joint distribution
is limited to its total number of available jobs , respectively. Consequently, the multivariate normal distributions
are fully independent of each other, in the sense that, the random allocation of jobs of one
within the ADA doesn’t provide any information about the way another
will allocate its own resources in the same ADA. In other words, we can write the following variance-covariance matrix
to summarize the setting. In other words,
Where,
is the variance-covariance matrix of the joint multivariate normal of vector of vectors , that is,
, if
is large, respectively.
Each element of the matrix
has been introduced previously. is only a way to re-package existing information to make ideas clearer. The dimension of is
Note . must be seen as a nested block matrix. That is, is a diagonal block matrix, where each element on the diagonal is itself a block matrix of dimension
representing inter-connections within a single , and each element off-diagonal is a block matrix of dimension
populated entirely by zeros representing the absence of inter-connection between a pair of distinct part of the same ADA and with available data. The latter is nothing more than the matrix representation of the full independence of the several multivariate normal distributions justifying the use of the product operator in the previous joint multivariate normal equation.
The material presented over the previous page of this annex propose a joint distribution of
distinct multivariate normal where each one is dedicated to a specific of the same ADA. The input variable is a vector of vectors
and require the analyst to input
distinct value. This number of inputs to provide is large and a burden for the analyst. For this reason, we re-package the same original ‘’joint’’ idea into an alternative model where the input variable is of dimension
only. This alternative model make use of a single multivariate normal, however the mean vector is a sum of all
mean vectors and the variance-covariance matrix is a sum of all variance-covariance matrices of the multivariate normal distributions. Formally, in a finite context, we have,
If
is large, respectively.
Where,
is a row vector of dimension , where the sum of terms of
is the scalar , that is,
And where,
And where,
Under independence of
level variables, and where the new multivariate normal distribution can be expressed asNote ,
differ significantly from presented earlier in this annex. The input variable
specify the allocation of unique jobs within a single input vector of dimension , instead of
distinct vectors of dimension
each.
is approximated by a multivariate normal distribution since the exact distribution would be a convolution of several multinomials not available in closed form.
Finally, we can now describe the randomization process of Sergerie et al, 2021 as a special case, where for all DB of a same DA with data available, we have,
.
The last condition guarantees that all unique jobs of the DA are randomly and uniformly allocated within the DA. Also, the main randomization method of this paper is a special case if,
for all DB of the ADA with available data. (doesn’t apply to y***).
Now that y*, y** and y*** have been introduced in detail, we can present our justification of why a within DB random allocation of jobs is preferred for the practical case of our application.
With y**, under a large geospatial surface such as a DA, there is 2 items to assume; 1) each multivariate normal distribution
is assumed to be well specified as a normal distribution, respectively. More precisely, for each multivariate distribution, the expected rate of occurrence λ of the polygonal space interval of interest, namely, a single spatial spot or a single pixel in the DA, is more than moderate. That is, expected outcomes vectors
are respectively made of
distinct elements equal or greater than 10. In other words, we assume here that a large enough number of jobs is present in each DB of the ADA where data are available and the number of spatial spots available in the DA is not large enough to lower down the ratio of the average number of jobs per spatial spot, below 10. Otherwise, a multivariate Poisson distribution (a multivariate right skewed normal getting squeezed toward its multivariate origin coordinate) would be more suitable for some of the distributions composing the joint distribution and the shape of the joint normal could be altered. 2) On the same topic, some probability vector could potentially be made of probability terms close to zero due to the large number of categories. This could challenge the good approximation toward the normal distribution. However, we count on the availability of a large per DB where data is available to compensate for the issue of sparcity and bring back the required quality of normal approximation. The intuition here, is similar to a dataset made of several variable, where the number of observations is large enough to handle the many dimensions of information.
Now, for the cases of our own application, if each DB of the ADA with available data allocates jobs data points from a single 2-digit NAICS within the DA, then it is unlikely that a minimum of 10 jobs will succeed to cover each spatial spot of the DA. Therefore, it is unlikely to produce accurate marginal normal approximations. It is also likely that the accuracy of the joint normal will be altered, even under uniformity of probabilities for the full DA. Consequently, the expected outcome won’t be a uniform allocation of jobs within the DA. It will be a sparce allocation not stable over random realizations. The property of smoothness will be lost. A larger kernel bandwidth will be needed to compensate for the instability and scarcity, and local accuracy of clusters will be lost. For this first reason, a random allocation within DB seems to be preferable.
A solution could be y***. With y***, it becomes more likely that a category will get a larger number of jobs randomly allocated because it is about the sum of counts from the DBs of the ADA with available data. However, for the cases of our own application, it is still unlikely that a single 2-digit NAICS will succeed to reach 10 jobs per spatial spot. For this second reason, a random allocation within DB seems to be preferable again.
A final alternative solution could be to aggregate uniformly the spatial spots of the DA and reduce the number of categories in such a way that the total number of jobs available in the DA divided per the number of aggregated categories is now equal to 10. However, for the cases of our own application, under a single 2-digit NAICS, the number of aggregated categories is likely to be not so large. A normal approximation will be reliable, and the expected outcome will be a uniform and stable allocation of jobs within the DA. However, again, a larger kernel bandwidth will be needed to compensate for the lack of high resolution and local accuracy of clusters will be lost. For this third reason, a random allocation within DB seems to be preferable again. An accurate normal approximation is not always easy to reach out to. However, in the case of our applications, a random allocation within DB remains our best option. Documenting the normal approximation in this paper remain essential, because if you have a deviation, the best way to understand this deviation is to understand what you are deviating from in the first place. ∎
Appendix 3: cluster mapping results
Description for Map 5
Map showing the CMA of Montréal. Dark red color represents Dissemination blocks part of the Manufacturing sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 6
Map showing the CMA of Montréal. Dark red color represents Dissemination blocks part of the Retail Trade sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 7
Map showing the CMA of Montréal. Dark red color represents Dissemination blocks part of the Accommodations and Food Services sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 8
Map showing the CMA of Montréal. Dark red color represents Dissemination blocks part of the Distribution and Electronic Commerce cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 9
Map showing the CMA of Montréal. Dark red color represents Dissemination blocks part of the Financial Services sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 10
Map showing the CMA of Montréal. Dark red color represents Dissemination blocks part of the Hospitality and Tourism sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 11
Map showing the CMA of Toronto. Dark red color represents Dissemination blocks part of the Manufacturing sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 12
Map showing the CMA of Toronto. Dark red color represents Dissemination blocks part of the Retail Trade sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 13
Map showing the CMA of Toronto. Dark red color represents Dissemination blocks part of the Accommodations and Food Services sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 14
Map showing the CMA of Toronto. Dark red color represents Dissemination blocks part of the Distribution and Electronic Commerce sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 15
Map showing the CMA of Toronto. Dark red color represents Dissemination blocks part of the Financial Services sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 16
Map showing the CMA of Toronto. Dark red color represents Dissemination blocks part of the Hospitality and Tourism sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 17
Map showing the CMA of Winnipeg. Dark red color represents Dissemination blocks part of the Manufacturing sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 18
Map showing the CMA of Winnipeg. Dark red color represents Dissemination blocks part of the Retail Trade sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 19
Map showing the CMA of Winnipeg. Dark red color represents Dissemination blocks part of the Accommodations and Food Services sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 20
Map showing the CMA of Winnipeg. Dark red color represents Dissemination blocks part of the Distribution and Electronic Commerce sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 21
Map showing the CMA of Winnipeg. Dark red color represents Dissemination blocks part of the Financial Services sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 22
Map showing the CMA of Winnipeg. Dark red color represents Dissemination blocks part of the Hospitality and Tourism sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 23
Map showing the CMA of Vancouver. Dark red color represents Dissemination blocks part of the Manufacturing sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 24
Map showing the CMA of Vancouver. Dark red color represents Dissemination blocks part of the Retail Trade sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the definition of a geographical spatial distance and a spatial projection definition.
Description for Map 25
Map showing the CMA of Vancouver. Dark red color represents Dissemination blocks part of the Accommodations and Food Services sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 26
Map showing the CMA of Vancouver. Dark red color represents Dissemination blocks part of the Distribution and Electronic Commerce sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 27
Map showing the CMA of Vancouver. Dark red color represents Dissemination blocks part of the Hospitality and Tourism sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
Description for Map 28
Map showing the CMA of Vancouver. Dark red color represents Dissemination blocks part of the Hospitality and Tourism sector cluster. An info box is available at the top right corner of the map. This includes the name of the CMA, the name of the cluster, a legend for the length of a geographical spatial distance in the map, and a spatial projection definition, which correspond to Lambert Conformal Conic Projection Datum: North American 1983 (NAD83).
More information
ISSN: 0000-0000
Note of appreciation
Canada owes the success of its statistical system to a long-standing partnership between Statistics Canada, the citizens of Canada, its businesses, governments and other institutions. Accurate and timely statistical information could not be produced without their continued co-operation and goodwill.
Standards of service to the public
Statistics Canada is committed to serving its clients in a prompt, reliable and courteous manner. To this end, the Agency has developed standards of service which its employees observe in serving its clients.
Copyright
Published by authority of the Minister responsible for Statistics Canada.