Abstract
Temporal, spatial and spacetime scan statistics are commonly used to detect and evaluate the statistical significance of temporal and/or geographical disease clusters, without any prior assumptions on the location, time period or size of those clusters. Scan statistics are mostly used for count data, such as disease incidence or mortality. Sometimes there is an interest in looking for clusters with respect to a continuous variable, such as lead levels in children or low birth weight. For such continuous data, we present a scan statistic where the likelihood is calculated using the the normal probability model. It may also be used for other distributions, while still maintaining the correct alpha level. In an application of the new method, we look for geographical clusters of low birth weight in New York City.
Background
Spatial and spacetime scan statistics [14] have become popular methods in disease surveillance for the detection of disease clusters, and they are also used in many other fields. In most applications to date, the interest has been in count data such as disease incidence, mortality or prevalence, for which a Poisson or Bernoulli distribution is used to model the random nature of the counts. For example, in papers published in 2008, Chen et al. [5] studied cervical cancer mortality in the United States; Osei and Duker [6] studied cholera prevalence in Ghana; Oeltmann et al. [7] looked at multidrugresistant tuberculosis prevalence in Thailand; Mohebbi at al. [8] studied gastrointestinal cancer incidence in Iran; RubinskyElefant et al. [9] looked at human toxocariasis prevalence in Brazil; Frossling et al. [10] evaluated the Neospora caninum distribution in dairy cattle in Sweden; Heres et al. [11] studied madcow disease in the Netherlands; and Reinhardt et al. [12] developed a system for prospective meningococcal disease incidence surveillance in Germany.
It is also of interest to detect spatial clusters of individuals or locations with high or low values of some continuous data attribute. Gay et al. [13] developed a spatial hazard model which they applied to detect geographical clusters of dietary cows with a high somatic cell score, which is a continuous marker for udder inflamation. Stoica et al. [14] has proposed a cluster detection method based on a number of random disks that jointly cover the cluster pattern in a marked point process. Huang [15] and Cook et al. [16] have developed spatial scan statistics for survival type data with censoring. The former applied the method to prostate cancer survival while the latter used their method for the time from birth until to asthma, allergic rhinitis or exczema. Other continuous data, such as birth weight [17] or blood lead levels, may be better modeled using a normal distribution, sometimes after a suitable transformation.
In this paper we develop a scan statistic for continuous data that is based on the normal probability model. Under the null hypothesis, all observations come from the same distribution. Under the alternative hypothesis, there is one cluster location where the observations have either a larger or smaller mean than outside that cluster. A key feature of the method is that the statistical inference is still valid even if the true distribution is not normal, assuring that the correct alpha level is maintained. This is accomplished by evaluating the statistical significance of clusters through a permutation based Monte Carlo hypothesis testing procedure. The new method is applied to birth weight data from New York City. A simulation study is performed to evaluate the power for different types of clusters.
The application and simulation results presented in this paper are concerned with twodimensional spatial data, using a circular variable size scanning window. The new method is equally applicable to purely temporal and spatiotemporal data [1820], to be used for daily prospective disease surveillance to look for suddenly emerging clusters. In addition to circles, it may also be used with an elliptic scanning window [21], or with any collection of nonparametric shapes [2,2225].
The normal model has been incorporated into the freely available SaTScan software http://www.satscan.org webcite for spatial and sdpacetime scan statistics, so it is easy to use. While it requires the use of computer intensive Monte Carlo simulations, computing times are very reasonable, unless the data set is huge.
A Spatial Scan Statistic for Normal Data
Observations and Locations
The data consists of a number of continuous observations, such as birth weight, with values x_{i}, i = 1,...,N. Each observation is at a spatial location s, s = 1,...,S, with spatial latitude and longitude coordinates lat(s) and long(s). Each location has one or more observations, so that S ≤ N.
For each location s, define the sum of the observed values as x_{s }= ∑_{i∈s }x_{i }and the number of observations in the location as n_{s }The sum of all the observed values are X = ∑_{i}x_{i}.
Scanning Window
The circular spatial scan statistic is defined through a large number of overlapping circles [18]. For each circle z, a log likelihood ratio LLR(z) is calculated, and the test statistic is defined as the maximum LLR over all circles. The scanning window will depend on the application, but it is typical to define the window as all circles centered on an observation and with a radius varying continuously from zero up to some upper limit. To ensure that both small and large clusters can be found, the upper limit is often defined so that the circle contains at most 50 percent of all observations. It is never set above that number though, since a circular cluster with high values covering for example 80 percent of all observations is more appropriatly interpreted as a spatially disconnected 'cluster' with low values covering the 20 percent of observations that are located outside the circle, since it is those 20 percent that differ from the majority of observations. The maximum cluster size can also be defined using specific units of distance (e.g., 10 km). Circles with only one observation are ignored. Let n_{z }= ∑_{s∈z}n_{s }be the number of observations in circle z, and let x_{z }= ∑_{s∈z}x_{s }be the sum of the observed values in circle z.
Likelihood Calculations
Under the null hypothesis, the maximum likelihood estimates of the mean and variance are μ = X/N and respectively. The likelihood under the null hypothesis is then
and the log likelihood is
Under the alternative hypothesis, we first calculate the maximum likelihood estimators that are specific to each circle z, which is μ_{z }= x_{z}/n_{z }for the mean inside the circle and λ_{z }= (X  x_{z})/(N  n_{z}) for the mean outside the circle. The maximum likelihood estimate for the common variance is
The log likelihood for circle z is
This simplifies to
As the test statistic we use the maximum likelihood ratio
or more conveniently, but equivalently, the maximum log likelihood ratio
Only the last term depends on z, so from this formula it can be seen that the most likely cluster selected is the one that minimizes the variance under the alternative hypothesis, which is intuitive.
Randomization
The statistical significance of the most likely cluster is evaluated using Monte Carlo hypothesis testing [26]. Rather than generating random data from the normal distribution, a large set of random data sets are created by randomly permuting the observed values x_{i }and their corresponding locations s. That is, the analysis is conditioned on the collection of continuous observations that were observed, as well as on the locations at which they were observed, which are considered nonrandom. By doing the randomization this way, the correct alpha level will be maintained even if the observations do not truly come from a normal distribution. Note that it is the individual observations that are permuted, so two different observations in the same location will end up in two different locations in most of the random data sets.
For each random data set, the log likelihood lnL(z) is calculated for each circle. The most likely cluster is then found and its log likelihood ratio is noted. If the log likelihood ratio from the real data set is among the 5 percent highest of all the data set, then the most likely cluster from the real data set is statistically significant at the 0.05 alpha level. More specifically, if there are M random data sets, then the pvalue of the most likely cluster is R/(M + 1), where R is the rank of the log likelihood ratio from the real data set in comparison with all data sets. In order to obtain nice pvalues with a finite number of decimals, M should be chosen as for example 999, 4999 or 99999.
Note that these Monte Carlo based pvalues are exact in the sense that under the null hypothesis, the probability of observing a pvalue less than or equal to p is exactly p [26]. This is true irrespective of the number of random data sets M, but a higher M will provide higher statistical power.
If the random simulated data had instead been generated from a normal distribution with prespecified mean and variance, rather than through permutation, then one would test the null hypothesis that the observations come from exactly that normal distribution. We would then reject the null for many reasons other than the existance of spatial clusters. For example, the null may be rejected because the mean values are higher than specified uniformly throughout the whole study region.
Scanning for High or Low Values
As defined above, the normal scan statistic will search for clusters with exceptionally high values as well as clusters with exceptionally low values. Sometimes it makes more sense to only search for clusters with high values. The former is easily accomplished by adding an indicator function I(μ_{z }> λ_{z}) to the likelihood that is calculated under the alternative hypothesis. If one is only interested in cluster with low values, the indicator function is instead I(μ_{z }< λ_{z}).
Software
The normal scan statistic has been incorporated into the freely available SaTScan™ software package, version 7.0 http://www.satscan.org webcite. It can be used for temporal, spatial and/or spatiotemporal data. The spatial version may be applied using a circular or elliptic window in two dimensions or a spheric window in three or more dimensions. The spacetime version uses a cylindrical scanning window with either a circular or elliptic base. It is also possible for the user to define his/her own nonEuclidian neighborhood metric. The circle centroids can be identical to the collection of coordinates of the observations, or they may be specified by the user.
Detection of Low Birth Weight and Early Gestation Clusters in New York City
The New York City Department of Health and Mental Hygiene (NYCDOH) calculates infant mortality rates by neighborhood and reports them in its annual Summary of Vital Statistics [27]. Though the infant mortality rate has fallen dramatically over the past 20 years (from 13.4 per 1000 in 1988 to 5.4 per 1000 in 2007) neighborhood variation can be quite high and this has attracted much attention from public health officials, the press, and researchers [2830]. One approach to understanding the spatial pattern of infant mortality is to investigate the spatial patterns of known risk factors such as low birth weight, early gestation, or congenital conditions [28,17,32]. Attempts to identify clusters of low birth weight typically rely on dichotomized variables for birth weight (low: < 2500 grams; verylow: < 1500 grams) [17,31]. However, some researchers have noted a more complicated relationship among birth weight, gestation, and infant mortality with risk varying considerably within low birth weight and early gestation categories and modified by gender and other demographic characteristics [32].
This section examines spatial patterns of continuous measures of birth weight in New York City, using the spatial scan statistic with the normal probability model. Vital Records from NYCDOH were used to obtain data for all singleton births occurring in New York City in 2004. Births were georeferenced to the Mother's zip code. Births to mother's not residing in New York City and those with invalid zips were deleted. Birth weight was measured in grams. The normal spatial scan statistic was used to detect clusters of low birth weight, using a circular window shape and 50 percent of all birth as the maximum cluster size.
Two statistically significant geographical clusters of low birth weight were found (Table 1 and Figure 1). With a log likelihood ratio (LLR) of 125.8, the first one consists of 61 zipcode areas in eastern Brooklyn and southern Queens, where the birth weights were on average 60 grams less than the rest of the city (p < 0.001). The second cluster consists of 29 zipcode areas in northern Manhattan and southern Bronx, where the birth weights were on average 52 grams less (LLR = 62.7, p < 0.001). The two statistically significant clusters correspond closely to areas of increased risk for infant mortality and are highly correlated with clusters found using dichotomized variables. There was also a nonsignificant single zipcode cluster on Staten Island (60 grams less, LLR = 3.6, p = 0.90). Note that, while the weight difference is as large or larger in the State Island cluster, such a difference could easily be due to chance, due to the small number of births inside the cluster. Note also that since a circular scanning window is used, some low birth weight areas are just outside the BrooklynQueens cluster while some high birth weigth areas are just inside. The key thing to realize is that it is only the general area of the cluster that is detected, not its exact boundaries.
Figure 1. The geographical distribution of birth weight in New York City zip codes in 2004, with two statistically significant clusters found by the spatial scan statistic with the normal probability model.
Table 1. Geographical clusters of low birth weight in New York City in 2004.
For the City as a whole, the variance of the birth weights is 297250 and the standard deviation is 545.2. After accounting for the different means inside and outside of the most likely cluster, the variance is 296564 and the standard deviation is 544.6. These number are, by default, lower, but only marginally so. In fact, the most likely cluster only explains (297250  296564)/297250 = 0.23 percent of the total variance. This is not surprising for an outcome such as birth weight, since the natural variation is rather large.
The spatial pattern of birth weight may be largely driven by the spatial patterns of demographic and pregnancyspecific characteristics. If so, clusters are not surprising; simply reflecting the geographical distribution of known characteristics. As such, the public health utility of clustering in the raw data may be limited. In a substantive paper, we hope to reexamine the geographical clusters of birth weight adjusting for the mother's demographics, health status, and pregnancy characteristics that are known to correspond with low birth weight. This has two uses: first, it sharpens understanding of the relationship between demographic covariates and birth weight and second, it identifies areas with surprisingly low birth weight for investigation, which cannot be explained in terms of their underlying demographics.
Statistical Power and Spatial Precision
To evaluate the statistical power of the new method, we performed a simple simulation study. We simulated random normally distributed weights for infants born in New York City. The power to detect a cluster will depend on a number of factors, so data were generated using one standard baseline scenario and several variations: using different cluster locations within New York City, different cluster sizes, different sample size (total numbers of births), with different mean weights inside and outside the cluster, and with different variances inside and outside the true cluster.
In 2003, the average birth weight in New York City was around 3250 grams, with a standard deviation of approximately 600 grams. For a particular sample size, we fixed the total number of infants, and assigned them randomly to census tracts with the probability proportional to the census tract population size. All infants assigned into the same tract share the same latitude and longitude coordinates. This assignment was fixed and the same for all simulations with the same number of births.
We selected four cluster locations to have lower than average birth weight. These are shown in Figure 2. Cluster No. 1 is the baseline cluster, located in the center of Brooklyn. Cluster No. 2 is centered on The Rockaways, Queens, close to the Atlantic Ocean. Cluster No. 3 is located on Staten Island, far away from the rest of the City. Cluster No. 4 is split between southern Bronx and northern Queens. As the baseline, the maximum size of the clusters was defined to include 10 percent of all the births in the City. This means that the geographical size of the clusters varied, depending on the population density around the cluster centroid. For Cluster No. 1, we also evaluated clusters with 5 and 20 percent of the total number of births (dotted lines), while keeping the cluster centroid the same.
Figure 2. The location and size of the artificial low birth weight clusters used to evaluate the statistical power of the spatial scan statistic for normally distributed data.
Let Z denote the collection of census tracts in the true cluster and let Z^{c }denote the remaining New York City census tracts. For each of either 300, 600 (baseline), 900 or 1200 infants, we randomly simulate the birth weight from a normal distribution N(μ_{Z}, σ^{2}) if the infant was born inside the cluster and from N(, σ^{2}) if the infant was born outside the cluster. For all simulations, we set = 3250. We always choose μ_{Z }< , so that the simulated data has a cluster of low birth weight. For the baseline, we chose μ_{Z }to be ten percent less than , so that μ_{Z }= 3250  325 = 2925. We also evaluated clusters where μ_{Z }was 5, 8, 13, 15 and 20 percent less than . The variance σ^{2 }was the same inside and outside the cluster. For the baseline model we set the standard deviation to σ = 600, but we also evaluated σ = 300 and σ = 900. A complete list of all the evaluated cluster parameters are shown in the first five columns of Table 2.
Table 2. Estimated power, sensitivity and positive predictive value (PPV) when the normal scan statistic is used to detect different types of clusters, as described in the text.
For each cluster scenario, we simulated 1000 random data sets. The estimated power is calculated as the proportion of the 1000 random data sets for which the null hypothesis was rejected, expressed as a percentage.
Even when the null hypothesis is correctly rejected, the detected cluster is usually not exactly identical to the true cluster. The extant of the overlap, and hence, of the spatial accuracy of the detected cluster, can be evaluated using sensitivity and positive predicted value (PPV). The sensitivity is defined as the proportion of the infants in the true cluster that was included in the detected cluster. This obviously varies between the random data sets, and the estimated sensitivity is taken as the average over the 1000 random data sets. The positive predictive value is defined as the proportion of the infants in the detected cluster that are in the true cluster. Again, this is estimated by taking the average over the 1000 random data sets.
The results are presented in Table 2. The power is approximately the same for the four different cluster locations, which is a reflection of the fact that they are about the same population size. As expected, the power increase when the cluster size increase, when the sample size increase, when the mean weight difference increase and when the standard deviation decrease. Sensitivity and positive predictive value follow the same pattern. Note that the sensitivity is about the same as the positive predictive value. This means that we are about equally likely to leave out an infant that should be in the cluster as we are to include an infant that shouldn't be in the cluster. Note also that even when the power is 100, the sensitivity and positive predictive value are not. This means that while we can determine the general location of a cluster, there will almost always be uncertainty when it comes to the borders of the detected cluster.
Discussion
We have presented a scan statistic for continuous data. It is based on the normal distribution function, so if the data is truly normal, we have a likelihood ratio test. If the data follows some other distribution, it is no longer a likelihood ratio test, but it still maintains the correct alpha level. Hence, it can be used for a wide variety of continuous data, although we do not recommend it for exponential or other types of survival data, for which there are other scan statistics available [15,16].
The normal scan statistic performed well for the New York City birth weight data, finding two statistically significant clusters that corresponded to areas with high infant mortality.
The statistical power varies predictably with the type of cluster to be found. The same is true for sensitivity and the positive predictive value. One limitation of the simulation study is that we only evaluated the performance on data that were simulated from the normal distribution. While we know that the alpha level is correct for other distributions, we do not know about the power, sensitivity and positive predictive value.
As with most other scan statistics, the method is computer intensive, but not prohibitively so. The freely available SaTScan™ software http://www.satscan.org webcite is available to do the calculations in a purely temporal, purely spatial or spacetime setting, and when looking for clusters with either only high or only low values, or simultaneously for both. The normal probability model has been available in the SaTScan software since 2006, and the method has already been applied to study the epidemics of classical swine fever in Spain [33], the geographical differences in respondent and nonrespondents in epidemiological studies [34] and the geographical clustering of the time people spend walking and bicycling in Los Angeles and San Diego [35]. The method can also be used in other fields outside of medicine and public health. For example, the variable of interest could be the amount of rainfall in various geographical locations in a country, pollution levels in a city, the height of plants on a field or the size of stars in a galaxy.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
MK obtained funding and developed the statistical methods. KK selected and performed the SaTScan analysis on the real data. MK and LH designed and LH performed the simulated power evaluations. MK, KK and LH wrote the first draft for different parts of the manuscript. All authors revised the manuscript and approved the final version.
Acknowledgements
This work was funded by the National Institute of Child Health and Development, National Institutes of Health, grant number R01HD048852.
References

Naus J: Clustering of random points in two dimensions.
Biometrika 1965, 52:263267. Publisher Full Text

Kulldorff M: A spatial scan statistic.
Communications in Statistics: Theory and Methods 1997, 26:14811496. Publisher Full Text

Glaz J, Naus J, Wallenstein S: Scan Statistics. Springer; 2001.

Chen J, Roth RE, Naito AT, Lengerich EJ, MacEachren AM: Geovisual analytics to enhance spatial scan statistic interpretation: An analysis of US cervical cancer mortality.
International Journal of Health Geographics 2008, 7:57. PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Osei FB, Duker AA: Spatial dependency of V. cholera prevalence on open space refuse dumps in Kumasi, Ghana: A spatial statistical modeling.
International Journal of Health Geographics 2008, 7:62. PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Oeltmann JE, Varma JK, Ortega L, Liu Y, O'Rourke T, Cano M, Harrington T, Toney S, Jones W, Karuchit S, Diem L, Rienthong D, Tappero JW, Ijaz K, Maloney S: Multidrugresistant tuberculosis outbreak among USbound Hmong refugees, Thailand, 2005.
Emerging Infectious Diseases 2008, 14:17151721. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Mohebbi M, Mahmoodi M, Wolfe R, Nourijelyani K, Mohammad K, Zeraati1 H, Fotouhi A: Geographical spread of gastrointestinal tract cancer incidence in the Caspian Sea region of Iran: Spatial analysis of cancer registry data.
BMC Cancer 2008, 8:137. PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

RubinskyElefant G, SilvaNunes M, Malafronte RS, Muniz PT, Ferreira MU: Human toxocariasis in rural Brazilian Amazonia: Seroprevalence, risk factors, and spatial distribution.
American Journal of Tropical Medicine and Hygiene 2008, 79:9398. PubMed Abstract  Publisher Full Text

Frossling J, Nodtvedt A, Lindberg A, Björkman C: Spatial analysis of Neospora caninum distribution in dairy cattle from Sweden.

Heres L, Brus DJ, Hagenaars TJ: Spatial analysis of BSE cases in the Netherlands.
BMC Veterinary Research 2008, 4:21. PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Reinhardt M, Elias J, Albert J, Frosch M, Harmsen D, Vogel U: EpiScanGIS: An online geographic surveillance system for meningococcal disease.
International Journal of Health Geographics 2008, 7:33. PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Gay E, Senoussi R, Barnouin J: A spatial hazard model for cluster detection on continuous indicators of disease: Application to somatic cell score.
Veterinary Research 2007, 38:585596. PubMed Abstract  Publisher Full Text

Stoica RS, Gay E, Kretzschmar A: Cluster pattern detection in spatial data based on Monte Carlo inference.
Biometrical Journal 2007, 49:505519. PubMed Abstract  Publisher Full Text

Huang L, Kulldorff M, Gregorio D: A spatial scan statistic for survival data.
Biometrics 2007, 63:109118. PubMed Abstract  Publisher Full Text

Cook A, Gold DR, Li Y: Spatial cluster detection for censored outcome data.
Biometrics 2007, 63:540549. PubMed Abstract  Publisher Full Text

Ozdenerol E, Williams BL, Kang SY, Magsumbol MS: Comparison of spatial scan statistic and spatial filtering in estimating low birth weight clusters.
International Journal of Health Geographics 2005, 4:19. PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Kulldorff M, Athas W, Feuer E, Miller B, Key C: Evaluating cluster alarms: A spacetime scan statistics and brain cancer in Los Alamos.
American Journal of Public Health 1998, 88:13771380. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Kulldorff M: Prospective timeperiodic geographical disease surveillance using a scan statistic.
Journal of the Royal Statistical Society Series A 2001, 164:6172.

Kulldorff M, Heffernan R, Hartman J, Assunção R, Mostashari F: A spacetime permutation scan statistic for the early detection of disease outbreaks.
PLoS Medicine 2005, 2:e59. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Kulldorff M, Huang L, Pickle L, Duczmal L: An elliptic spatial scan statistic.
Statistics in Medicine 2006, 25:39293943. PubMed Abstract  Publisher Full Text

Patil GP, Taillie C: Geographic and network surveillance via scan statistics for critical area detection.
Statistical Science 2003, 18:457465. Publisher Full Text

Duczmal L, Assunçao R: A simulated annealing strategy for the detection of arbitrarily shaped spatial clusters.
Computational Statistics and Data Analysis 2004, 45:269286. Publisher Full Text

Tango T, Takahashi K: A flexibly shaped spatial scan statistic for detecting clusters.
International Journal of Health Geographics 2005, 4:11. PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Assunçao RM, Costa M, Tavares A, Ferreira S: Fast detection of arbitrarily shaped disease clusters.
Statistics in Medicine 2006, 25:723742. PubMed Abstract  Publisher Full Text

Dwass M: Modified randomization tests for nonparametric hypotheses.
Annals of Mathematical Statistics 1957, 28:181187. Publisher Full Text

Office of Vital Statistics: Summary of Vital Statistics 2004, the City of New York. New York. New York: New York City Department of Health; 2004.

Sohler NL, Arno PS, Chang CJ, Fang J, Schechter C: Income inequality and infant mortality in New York City.

Grady SC, McLafferty S: Disentangling the effects of residential segregation and neighborhood poverty on low birthweight for immigrant and nativeborn black women in New York City.
Urban Geography 2007, 28:377397. Publisher Full Text

Grady SC: Racial disparities in low birthweight and the contribution of residential segregation: A multilevel analysis.
Social Science and Medicine 2006, 63:30133029. PubMed Abstract  Publisher Full Text

Rushton G, Krishnamurthy R, Krishnamurt R, Lolonis P, Song H: The spatial relationship between infant mortality and birth defect rates in a U.S. city.
Statistics in Medicine 1996, 15:19071919. PubMed Abstract  Publisher Full Text

Solis P, Pullman SG, Frisbie WP: Demographic models of birth outcomes and infant mortality: An alternative measurement approach.
Demography 2000, 37:489498. PubMed Abstract  Publisher Full Text

MartínezLópez B, Perez AM, SánchezVizcaíno JM: A stochastic model to quantify the risk of introduction of classical swine fever virus through import of domestic and wild boars.
Epidemiology and Infection 2009, 137:15051515. PubMed Abstract  Publisher Full Text

Shen M, Cozen W, Huang L, Colt J, De Roos AJ, Severson RK, Cerhan JR, Bernstein L, Morton LM, Pickle L, Ward MH: Census and geographic differences between respondents and nonrespondents in a casecontrol study of nonHodgkin lymphoma.

Huang L, Stinchcomb D, Pickle L, Dill J, Berrigan D: Identifying clusters of active transportation using spatial scan statistics.
American Journal of Preventive Medicine 2009, 37:157166. PubMed Abstract  Publisher Full Text  PubMed Central Full Text