Evaluation of sliding baseline methods for spatial estimation for cluster detection in the biosurveillance system
 Jian Xing^{1}Email author,
 Howard Burkom^{2},
 Linda Moniz^{2},
 James Edgerton^{3},
 Michael Leuze^{4} and
 Jerome Tokars^{1}
https://doi.org/10.1186/1476072X845
© Xing et al; licensee BioMed Central Ltd. 2009
Received: 21 April 2009
Accepted: 17 July 2009
Published: 17 July 2009
Abstract
Background
The Centers for Disease Control and Prevention's (CDC's) BioSense system provides nearreal time situational awareness for public health monitoring through analysis of electronic health data. Determination of anomalous spatial and temporal disease clusters is a crucial part of the daily disease monitoring task. Our study focused on finding useful anomalies at manageable alert rates according to available BioSense data history.
Methods
The study dataset included more than 3 years of daily counts of military outpatient clinic visits for respiratory and rash syndrome groupings. We applied four spatial estimation methods in implementations of spacetime scan statistics crosschecked in Matlab and C. We compared the utility of these methods according to the resultant background cluster rate (a false alarm surrogate) and sensitivity to injected cluster signals. The comparison runs used a spatial resolution based on the facility zip code in the patient record and a finer resolution based on the residence zip code.
Results
Simple estimation methods that account for dayofweek (DOW) data patterns yielded a clear advantage both in background cluster rate and in signal sensitivity. A 28day baseline gave the most robust results for this estimation; the preferred baseline is long enough to remove daily fluctuations but short enough to reflect recent disease trends and data representation. Background cluster rates were lower for the rash syndrome counts than for the respiratory counts, likely because of seasonality and the large scale of the respiratory counts.
Conclusion
The spatial estimation method should be chosen according to characteristics of the selected data streams. In this dataset with strong dayofweek effects, the overall best detection performance was achieved using subregion averages over a 28day baseline stratified by weekday or weekend/holiday behavior. Changing the estimation method for particular scenarios involving different spatial resolution or other syndromes can yield further improvement.
Background
The problem of detecting anomalous case clusters
Epidemiologists responsible for routine public health surveillance seek early indications of possible disease outbreaks. A growing collection of data types has become available to aid this surveillance, but with the analysis of these data comes the need to reduce the search possibilities. With cluster detection methodology, analysis of data sources can help by indicating the likely location and extent of a potential outbreak. Spatial and spatiotemporal scan statistics have the additional advantages of controlling the number of alarms resulting from multiple testing and accounting for usual clustering patterns, as found in the data [1].
Some natural clustering is expected, as in urban areas with high population density. If cluster detection methods do not account for this heterogeneity, customarily high concentrations will be signaled as anomalous. Only departures from the usual distribution should be flagged. In the implementation of the widely used SaTScan software [2], a detection statistic is applied to all candidate clusters, and the cluster with the highest statistical value is tested for significance. The statistic is a function of observed and expected values both inside and outside the candidate cluster. Thus, an accurate estimate of the expected spatial distribution is necessary for determination of relevant clusters with high sensitivity at reasonable false alarm rates. This requirement is consistent whether the detection approach uses scan statistics or another methodology.
A natural choice for the expected spatial distribution is the census population by subregion of the geographic area under surveillance. In the United States, for example, population counts stratified by 3digit zip code, 5digit zip code, or census tract are freely downloadable [3]. This choice is appropriate for some applications of scan statistics to chronic disease surveillance where the data give true population rates relative to census totals and longterm effects are of interest. One problem with using these counts is that they represent the population measured at the most recent census, taken only every 10 years, so they do not include more recent population shifts. A more serious problem is that some biosurveillance data sources are not representative of the entire population, so they might produce expectations that are biased relative to the data of interest. This bias typically occurs because the data represent only one segment of the population that is not evenly dispersed over the region. Examples are data covering only military personnel and their dependents, as in this study, participants in a health maintenance organization or other medical plan, or customers of a pharmacy chain. Organizations furnishing such data typically do not or cannot furnish details of their spatial coverage. Therefore, system developers have used past data to form spatial expectations.
If substantial (typically multiyear) quality historical data are available, detailed modeling of counts at the subregion level could be attempted to calculate expected distributed counts. Kleinman et al [4] obtained significantly reduced false alarm rates with this approach. However, health monitors typically need to apply data with little data history, and longer data history is likely to be unrepresentative of current spatial patterns because of changes in population behavior, in participation of data providers, and in data collection systems. This article focuses on databased methods of computing expected count distributions assuming that only 1–8 weeks of recent data are available.
For clarity, we distinguish two separate problems in applying data with spatial information to detection of disease clusters. The first problem is how to aggregate the data into subregions. Ozonoff et al and Olson et al. [5, 6] demonstrated that the level of aggregation can strongly affect the power of scan statistics for cluster detection. The second problem is that given the total number of cases occurring in a monitored region and the chosen partition of subregions, what is the expected distribution of cases among these subregions if no public health event of interest is underway? Reasonable agreement with this distribution might be seen as the null hypothesis that is rejected if a case cluster is found at a level of significance above a preset threshold. Practical data considerations often drive the answer to the first problem. Privacy laws typically preclude the use of patient addresses or other exact location fields, so zip codes, census tracts, and locations of clinics specified in data records have been used. The second problem is the subject of this article.
Study dataset
The choice of an effective method for estimating the case distribution is a function of the monitored data source. Even within the same data source, for example, estimation of the current distribution for a broad, seasonal syndromic filtering of clinical records is likely different from the same estimation for a rare subsyndrome occurring evenly throughout the year. For the current study, we used data from BioSense, a US national biosurveillance system operated by the Centers for Disease Control and Prevention (CDC) [7]. This system uses computer databases for automated disease surveillance. The source of data for the study was the BioSense collection of US Department of Defense (DoD) military outpatient clinic data records. Since 2004, these data records have been sent daily to CDC through the TriCare Management Activity, which manages clinical patient data for all US military treatment facilities. The patient base includes personnel on active duty, retirees, and dependents. Patient record fields provide key information including the visit date, age, gender, residence and treatment clinic zip codes, and for medical classification, a code from the International Classification of Diseases, Ninth Revision (commonly called ICD9), used by hospital and physicians' office personnel to report billing information to insurance companies [8].
One obvious difference between the geographic patient distribution for the study dataset and the general census distribution results from the large fraction of nonTexas residents among these patients who live near the military bases, with the result that zip codes near those bases are overrepresented relative to census figures. However, for reasons of both patient privacy and military operational security, the spatial distribution of people eligible for treatment at base clinics is unavailable for analysis. Because of commuting by Texas residents, it cannot be assumed from Figure 1 that nearly all patient residences are close to treatment facilities. Many commuters are eligible for treatment at military clinics throughout the continental United States, and in western states such as Texas, commuting distances of over 50 miles are common. For these reasons, databased spatial estimation methods are required. The Methods section discusses the data preparation and the estimation methods tested for cluster detection.
Methods
Treatment of data record spatial information

First, we included all records that have either a facility zip code or a home zip code inside Texas. We also included records from four facilities outside of Texas but near the TexasOklahoma and TexasLouisiana borders because of the large numbers of Texas residents commuting across the border and using those facilities.

Second, we calculated for each record the distance between the home and facility zip codes.

For the residencebased analysis, we used the patient home zip code if this distance was less than 100 miles. Otherwise, we used the facility zip code.

If the home zip code field was empty, the facility zip code was used as the residence zip code. If the facility zip code was missing or outside Texas, the home zip code was used.
We used this reasonably simple rule because a more detailed assignment would require information that was available neither for the study nor for routine monitoring. In a normal situation most people would probably not travel more than a couple of hours to a facility unless there were no other choices.
Calculation of expected counts
Baselinemean approach

Unstratified mean , or C2based: The expected count for a subregion was computed as the mean count over all baseline days for that subregion. For the 7day baseline, this method gives the same expected count for each zip code as that used for the temporal Early Aberration Reporting System (EARS) C2 algorithm [10].

Weekday/weekend mean , or W2based: The expected count for a zip code is over the only weekday (weekend) baseline day if the current day is a weekday (weekend). Days occurring on calendar holidays are classified with weekends because in several BioSense datasets, holiday counts are similar to weekend counts.

▪ Individual weekday mean , or dayofweek (DOW) based: The expected count for a zip code is more than only baseline days on the same day of the week as the current day. Days occurring on calendar holidays are classified with Sundays.
We used baseline lengths of 7, 28, and 56 days (1, 4, and 8 weeks) with a 2day buffer between baseline and test day. As in the EARS methods [10], the purpose of the buffer is to reduce the effect of any current outbreak on the baseline data. We excluded the 7day baseline from W2based and DOWbased methods because the 7day baseline is too short for weekends and holidays. For each test day, these averages were calculated and conditioned on the daily case total to obtain expected counts for the scan statistics (see Appendix).
The expected count for each subregion is calculated as a fraction of the total number of cases in the test period, set to 1 day in this study. It is critical to correctly estimate this fraction to the correct recognition of spatial aberrations. In the traditional scan statistics application of SaTScan, this fraction is the population of the subregion divided by the total population of the surveillance region. In the current data environment, where the denominator population is unknown, we estimate this fraction as the number of subregion baseline cases divided by the sum of baseline cases in the entire region. The only difference in this estimate for the C2, W2, and DOWbased prediction methods is which days are included in the baseline.
 a)
The test day is a weekday.
 b)
There are 50 cases in the entire region for the test day and 200 data subregions.
 c)
There are a total of 1,000 cases on weekdays in the baseline.
 d)
Subregion j has 75 cases on weekdays in the baseline.
If one or more of the cases above are true, then the expected number of cases is [(75 + 1)/(1000 + 200)] * 50 = 3.167. See Appendix 1 for the formal statement of this estimation method.
(1) SpaceTime permutation approach
Figure 4(b) shows the spacetime permutation approach. In this approach we get estimated counts by conditioning on spacetime marginal (SPM) totals for both geocoordinates and days as in the SaTScan spacetime permutation option [11]. In other words, expectations are based on keeping the total counts constant for both individual days for all subregions, and individual subregions for all days. For each randomization trial, the set of dates for all records is kept fixed and the set of subregions shuffled to keep these marginal totals constant. For a purely spatial scan statistic, this method is equivalent to the unstratified mean method with the baseline including the test period. The baseline might also be restricted as in the W2based and DOWbased methods. As noted above, this approach guarantees that the baseline sum will be positive for any subregion with positive count in the test period, but the disadvantage is that including the test day biases the baseline distribution.
Software implementation
We implemented the algorithms described in the previous sections into C and MatLab computer programs both for research purposes and for possible use in BioSense. Development of codes giving identical random draws in both software environments allowed perfect agreement of results despite the diverse coding methods. These methods were also developed to facilitate comparison of the spatial estimation methods, to allow inspection of all variables at all stages, and to incorporate two of the following features: (1) computation of the statistical significance of maximum likelihood clusters using the extreme value (or Gumbel) distribution [12], and (2) for the capability to inject realistic signals on the background data for detection performance testing. The C program runs in PC and/or Unix platforms and is callable from SAS code. This code does not have all of the features of SaTScan. Readers interested in obtaining the code should contact the corresponding author.
Cluster search space
Our implementation searches cylinders with circular bases. The cylinder height is the number of days in the current test period, so that for a 7day cylinder, a candidate cluster might have from 1 to 7 days of data. The analysis that follows focuses on comparing spatial estimation methods and is restricted to 1day cylinders. For the spatial search, any point grid of possible cluster centers could be used, and our study used the set of centroids of zip codes in the data. For analyses based on facility and presumed residence zip codes, the numbers of zip codes were 32 and 1,233, respectively. We defined the radius of a candidate cluster as the maximum distance from the center to another centroid in the cluster, and we allowed radii as large as 100 km for the residence zip code runs and 300 km for the clinic zip runs, based on possible care seeking behavior.
Test statistic
This statistic gives a measure of the likelihood of the observed counts inside and outside the cluster given the expected distribution and the total number of observed cases. Candidate clusters within the search space are ranked in descending order by their GLR values. The next step is to decide which if any of these clusters should be flagged as statistically significant.
Significance testing
The LLR values used as the test statistic do not comply with a known distribution in general, so thresholds for significance are usually determined empirically. The empirical approach is to randomly generate a large set of simulated data distributions under the null hypothesis of the expected spatial distribution for the test period. For each simulated distribution k, apply the same search procedure used for the observed data to obtain the maximum test statistics max(LLR(k)) over all candidate clusters. Derive the LLR threshold for significance from the set of max(LLR(k)) and apply it to decide if one or more of the clusters found using the observed data are significant. Specific steps include:
Step 1
from the set of subregion expected values E(j) derived by the estimation method of choice, and J is the number of subregions. Let N be the total number of cases observed on a day in which data are analyzed. For each of these cases, we choose a subregion with a multinomial random draw using {p(j), j = 1,,, J}. A random number x between 0 and 1 is drawn using the uniform distribution, and subregion j is chosen if x is between ∑^{ j1 }p(j) and ∑^{ j }p(j), so that the likelihood of choosing region j is proportional to p(j). We thus obtain a random distribution corresponding to the expected ratios. For seeking spacetime clusters in test periods of more than one day, we apply this same procedure for each test day.
For the spacetime permutation method, we generated trial distributions by implementing the procedure in [11], keeping the dates of all baseline cases fixed and shuffling their subregions. As noted above, the baseline includes the test period in this method, so this shuffling gives a new test distribution. For efficiency, we achieved this randomization using the FisherYates Shuffle, commonly called the Knuth Shuffle [13].
Step 2
We applied the search for the maximum test statistic for each random trial distribution exactly as for the observed data and collected the set of trial LLR maxima.
Step 3
The pvalue is then tested for significance against a threshold chosen to control the background cluster rate. We used a pvalue threshold of 0.0027, for a nominal daily background rate of one cluster per year, in the results below. Following [12], this significance method has two advantages over simply ranking LLR* among the trial maxima. First, fewer trial distributions are required to obtain a useful pvalue estimate, and second, the Gumbel pvalues are on a continuous scale so that strong significance can be inferred without a huge number of trials.
Detection performance assessment
Reliable quantitative measurement of the sensitivity and specificity of a cluster detection method requires multiple signals distributed in time over the subregions of interest. Documentation of authentic outbreaks sufficient to provide enough detail for this purpose is very rare. Previous power studies [15, 16] aimed at chronic disease applications focused on increased risk in a few subregions. However, infectious disease surveillance is concerned with detection of transient events in which cases are variably distributed in the region of interest. For a controlled, quantifiable detection performance capability, we therefore added to the software an injection component that allows the user the capability to: (1) specify the injected signal in detail, including the number of cases to inject, the epidemic curve of these injected cases by day, and the geometry of distribution of these cases; (2) produce the injected cases for each day using random draws and add them to the background data; (3) run the scan statistics for a set of days including the days of injected cases to determine whether clusters are detected at critical significance levels for the injected signals; and, (4) conduct sets of repeated trials of these tests at different start dates to allow computation of detection sensitivity and timeliness. The stochastic generation of signals and changing data background over time provides realistic variability for these trials.

For each trial, a signal start date is selected, and we generate a stochastic epidemic curve using random draws from a lognormal distribution [17, 18].

Parameters of this distribution are supplied by the user to control the median peak day and time spread of the injected cases.

These random draws are tabulated to give the total number of injected cases on each day. For each day with injected cases, these cases must be assigned to subregions, or zip codes in this study.

Three geometries have been implemented to reflect hypotheses of spread of an aerosol pathogen. The implemented geometries include radial, wedgeshaped, and hourglassshaped spread at rates that the user may specify.

For the results given here, we used the radial spread, in which the dropoff of cases from the centre zip code is proportional to e^{d/k}, for a dispersion factor k that determines the radial decay rate.
In the current study, total signals consisted of either 50 or 100 cases added to the original respiratory counts for either home or clinic zip codes. For epicurve generation, we used lognormal parameters ζ = 1.3, σ = .4. These parameters gave regionwide time series epidemic curves of 8–13 days in duration. The radial spatial distribution was calculated using distance from the signal centre in kilometres, and for our Texas DoD dataset, we used a dispersion factor of k = 18 for a plausible spread of cases. We modified this spatial distribution procedure to avoid producing obvious signals by injecting large numbers of counts in zip codes where there were few or no baseline counts. For this reason, the maximum allowed signal for each zip code for each day was at least 2, but no more than the maximum number of background cases in that zip code for the case day of week. We ran the distribution algorithm imposing this criterion recursively until all cases were distributed. For the repeated trials to measure detection sensitivity, we chose a start day for the first signal, and in each subsequent trial, the start date was advanced by 8 days to sample the background counts at different seasons, avoid dayofweek bias, and reduce correlation among successive trials. We generated these trials over approximately 3 years of background data to obtain 127 trials. For comparison of detection capability using residence zip codes and facility zip codes, we first produced injected signals at the resolution of the 1,233 residence zip codes and ran the scan statistics at that resolution. We then distributed these cases to the 32 facilities using distributions based on the distance to each facility, for comparison runs at the coarser resolution.
For the sequential runs without injects to assess background cluster rates, the update procedure is depicted in Figure 4. Our primary purpose was to compare spatial estimation methods, so the test period for all runs was a single day with no lookback, the purely spatial case. For dates 1/1/2004 to 12/31/2006, we started the first run at the end of the 2day buffer after the first baseline period. The maximum spatial cluster sizes were 300 km and 50 km for DoD facility zip scanning and patient zip scanning, respectively.
Results
Crosschecking results of the C and Matlab software described above for implementing scan statistics, we compared the four estimation methods described in the Methods section using both facility and patient residence zip codes with 7, 28, and 56day sliding baselines. The measures of effectiveness were compared using the background cluster rate and sensitivity to simulated clusters at chosen significance thresholds.
Results of background cluster rate comparisons (a false alarm rate surrogate)
We applied the four estimation methods to scan all test dates and all coordinates in this study dataset. Because they are to be used for daily surveillance, and disease outbreaks are expected to be rare, a primary concern for using these methods is the actual cluster determination rate, or how often the method will indicate a need for investigation. If alarms occur too often, then the alarms corresponding to true health events might be ignored.
We ran the purely spatial cluster detection algorithm for the 3 years of data with no correction for temporal multiple testing to get the relative cluster determination rates of the four methods. In the absence of information about any true outbreaks during the data period, we could not refer to significant clusters as true or false alarms, so as a surrogate for a false alarm rate, we measured the background alarm rate as the number of significant clusters per unit time for the pvalue 0.003.
Number of clusters found per 100 weekdays/weekends, DoD Texas, 2004–2006, (pvalue = 0.003, subregion = 2).
Spatial Resolution  Syndrome  Weekday/Weekend  Method/Baseline Duration (days)  

C2  W2  DOW  STM  
7  28  56  28  56  28  56  7  28  56  
Facility Zip  Rash  Weekday  2.0  1.9  2.0  0.9  1.7  1.9  1.5  2.2  1.8  1.8 
Weekend  1.3  0.6  0.6  1.0  1.0  1.0  1.0  1.0  1.0  1.0  
Respiratory  Weekday  21.6  21.4  25.3  13.9  18.5  18.9  19.9  17.0  18.9  23.4  
Weekend  21.3  23.9  22.0  9.2  8.0  7.6  6.1  22.3  23.6  23.6  
Residence Zip  Rash  Weekday  4.3  1.8  1.4  1.9  1.5  5.4  2.7  1.9  1.3  1.5 
Weekend  2.2  1.6  1.6  1.6  0.6  3.2  1.6  2.2  2.2  1.9  
Respiratory  Weekday  20.1  13.0  15.3  13.3  13.7  23.7  18.9  10.9  11.1  14.5  
Weekend  44.3  36.6  35.0  19.7  20.4  19.4  15.9  44.6  34.7  34.4 
▪ From Table 1, in general, the rash syndrome yield lower cluster rates than the respiratory syndrome, probably because of the seasonality and high mean and overdispersion of the respiratory counts.
▪ For the rash syndrome, the W2 estimation method with the 28day baseline has a clear advantage over the unstratified averaging of the C2 method at the facility level, but not at the residence level, at which the dayofweek effect is less clear. The plotted levels of cluster alerts suggest that the full respiratory syndrome is too noisy for this dataset in the sense that its high mean and overdispersed variance cause much more significant but random clustering than would be expected for a sparser outcome variable, or especially for a dataset in which spatial distribution is similar to that of census data. Replacement by a more specific subsyndrome such as influenzalike illness (ILI) should be considered; even though ILI counts are also seasonal, counts are much smaller, and their spatial distribution is likely to be fairly stable for most of the year.
▪ For the respiratory syndrome, management of dayofweek (DOW) effects is a primary concern, and the W2 and DOW methods clearly yield fewer cluster alerts than methods that ignore these effects. See Figures 5(b) and 6(b).
▪ The longest baseline of 56 days gives reduced background clustering only in certain circumstances, such as when the DOW method is used at the residence zip level. The finer temporal DOW stratification, combined with the finer spatial stratification, appears to require a longer baseline. For coarser aggregation, the longer baseline can actually increase the number of cluster alerts, likely because of seasonal effects.
▪ Extending the temporal stratification from the W2 to the DOW method never yields a clear overall advantage in reduced alerting.
▪ The W2 method with a 28day baseline is a safe choice for controlling the number of cluster alerts overall and gives as clear advantage in some situations. However, for a sparse outcome variable such as the rash visit counts and a fine spatial subdivision of records, a longer baseline might reduce the alerting.
Sensitivity comparison with signal injection
As described in the Methods section, we randomly generated 127 visit count signals with a regionwide lognormal epidemic curve and a radically decaying spatial distribution. Separate signal sets were generated for totals of 50 and 100 outbreak cases. For each set, we added these signals to the authentic respiratory syndrome background data for a series of detection trials.
Discussion
The above results substantially support the idea that accurate spatial estimation methods are essential for robust detection of anomalous spatial clusters at controlled background cluster rates. Among the SaTScan user community, when an expected spatial distribution of cases is not available in the form of a population file with reliable estimates, a common method for deriving subregion expected counts is the conditioning on marginal distributions for both dates and subregions, as in the spacetime permutation scan statistic [11]. For the purely spatial scan statistic, this conditioning is equivalent to a special case of the moving baseline methods described above, with the test period data included in the baseline and with no buffer between baseline and test cases. The advantage of this method is that there is never a problem of assessing the likelihood of cases in a subregion that has no baseline cases. The disadvantage is that the baseline is biased by the test period cases. Effects of this disadvantage can become significant if the baseline case counts are sparse. The estimation methods in the current paper provide an alternative whereby, in SaTScan terms, the Poisson method could be used and the population file updated with subregion estimates from a moving baseline separated from the test period by a buffer. In this alternative, if case records occur in the test period from subregions that have no baseline cases, expected values are computed without bias by the test cases, as discussed in the Appendix. Given the current baseline interval, the problem is then to estimate an expected spatial distribution that accounts for key factors such as seasonality and dayofweek patterns. The above study compares several estimation methods for two syndrome groupings and two levels of spatial resolution of the Texas DoD data.
 a)
The complexity of spatial estimation should be appropriate for the time series structure and reflect what is known about the data. Simple methods are preferable unless the data structure is rich enough to derive advantage from more complex methods such as finer stratification (e.g., the DOW method) or the regression modeling of [4]. More complex methods have a greater dependence on data quality and involve expense for development of new data sources and convergence issues (e.g., for nonlinear regression).
 b)
Choosing the baseline interval length involves a tradeoff between stability and capturing recent data behavior. In the charted results for the flat average (C2 based) method, 7day baseline estimates produced consistently higher background cluster rates than the 28day estimates. Lengthening the baseline further gives more stability, but can actually degrade the clustering performance if the time series mean changes either because of unmodeled seasonal effects or because of altered data provider participation or other common changes in the data acquisition process. Lengthening the baseline from 28 days to 56 days yielded an advantage for certain situations, such as when the DOW stratification was used, but in some situations produced higher background cluster rates.
 c)
In some health monitoring systems, a single alerting method with fixed parameters might be required for operational reasons. Given such a requirement, the weekday or weekend/holiday stratified averages (W2based) with a 28day baseline provide a safe and sometimes the best choice among the methods tested for the DoD outpatient visit data. The unstratified averages might yield marginally better results if dayofweek effects are absent, but could fare much worse if such effects are present. Moreover, dayofweek effects are not always apparent from visual inspection.
 d)
The results suggest that more robust clustering performance is possible in a system whereby the estimation method and its parameters could be chosen according to the data. For example, one could use dayofweek stratification time series with rich counts like the respiratory syndrome and large subregions, as in the facilitylevel runs. A longer baseline could be advantageous for more rare syndrome groupings such as rash, especially if seasonal effects are absent.
The experience of this study also suggested that the pvalue threshold for cluster significance should be determined by test runs using authentic data for different outcome variables. Finally, clustering runs on background data could help evaluate the utility of syndrome groupings. The yearround high mean and variance of the respiratory syndrome count data lead to a high background cluster rate, as measured by the sequential background runs. This observation suggests a more selective record filtering, such as a subsyndrome for influenzalike illness. The runs using the sparser rash syndrome counts support this approach as did additional runs with other syndrome groupings.
Conclusion
This study was conducted as part of an effort to bridge the gap between surveillancerelated research and the need for robust and sensitive routine health monitoring. The BioSense data environment has grown to include various types of data from an increasing number of hospitals and other clinical data providers. Statistical behaviors such as seasonal trends, dayofweek effects, and age distributions vary among these sources, and appropriate syndromic classifications are still under study. Next steps are to seek appropriate spatial estimation methods for the new BioSense clinical data sources and record groupings and to standardize the method selection criteria. Regardless of the cluster detection method of choice in the future, reliable spatial background estimation will remain essential for robust surveillance.
Appendix
Computing expected values for each subregion by the modified baseline method:
Declarations
Acknowledgements
We acknowledge Drs Kenneth Cox and Julie Pavlin, and other US Department of Defense personnel who have graciously shared data with the BioSense program for several years.
Authors’ Affiliations
References
 Kulldorff M, Athas W, Feuer E, Miller B, Key C: Evaluating cluster alarms: a spacetime scan statistic and brain cancer in Los Alamos. Am J Public Health. 1998, 88: 13771380. 10.2105/AJPH.88.9.1377.PubMedPubMed CentralView ArticleGoogle Scholar
 SaTScan v7.0: software for the spatial and spacetime scan statistics. 2007,http://www.satscan.org/
 US Census Bureau downloadable population data.http://factfinder.census.gov/home/saff/
 Kleinman KP, Abrams AM, Kulldorff M, Platt R: A modeladjusted spacetime scan statistic with an application to syndromic surveillance. Epidemiol Infect. 2005, 133: 409419. 10.1017/S0950268804003528.PubMedPubMed CentralView ArticleGoogle Scholar
 Ozonoff A, Jeffery C, Manjourides J, ForsbergWhite L, Pagano M: Effect of spatial resolution on cluster detection: a simulation study. Int J Health Geogr. 2007, 6: 5210.1186/1476072X652.PubMedPubMed CentralView ArticleGoogle Scholar
 Olson KL, Grannis SJ, Mandl KD: Privacy protection versus cluster detection in spatial epidemiology. Am J Public Health. 2006, 96 (11): 20022008. 10.2105/AJPH.2005.069526.PubMedPubMed CentralView ArticleGoogle Scholar
 CDC BioSense Web site.http://www.cdc.gov/BioSense/
 Hart AC, Hopkins CA, Eds: 2003 ICD9CM expert for hospitals. 2003, Salt Lake City, UT: St. Anthony Publishing, 6Google Scholar
 Syndrome groups definitions in CDC's BioSense surveillance system.http://www.bt.cdc.gov/surveillance/syndromedef/index.asp
 Hutwagner L, Thompson W, Seeman GM, Treadwell T: The bioterrorism preparedness and response early aberration reporting system (EARS). J Urban Health. 2003, 80 (2): i89i96.PubMedPubMed CentralGoogle Scholar
 Kulldorff M, Heffernan R, Hartman J, Assunção R, Mostashari F: A spacetime permutation scan statistic for disease outbreak detection. PLoS Med. 2005, 2 (3): e5910.1371/journal.pmed.0020059.PubMedPubMed CentralView ArticleGoogle Scholar
 Abrams A, Kulldorff M, Kleinman K: Empirical/asymptotic pvalues for Monte Carlobased hypothesis testing: an application to cluster detection using the scan statistic. Proceedings of the Fourth Conference on Extreme Value Analysis and Statistical Models and their Applications: August 15. 2005, ; GothenburgGoogle Scholar
 Knuth DE: The Art of Computer Programming volume 2: Seminumerical algorithms. 1969, Reading, MA: AddisonWesley, 124125.Google Scholar
 Evans M, Hastings N, Peacock B: Statistical distributions. 1993, New York: John Wiley and Sons, 682Google Scholar
 Kulldorff M, Tango T, Park P: Power comparisons for disease clustering tests. computational statistics and data analysis. Int J Health Geogr. 2003, 42: 665684.Google Scholar
 Song C, Kulldorff M: Power evaluation of disease clustering tests. Int J Health Geogr. 2003, 2: 910.1186/1476072X29.PubMedPubMed CentralView ArticleGoogle Scholar
 Evans M, Hastings N, Peacock B: Statistical distributions. 1993, New York: John Wiley and Sons, 102105. 2Google Scholar
 Brunk HD: An Introduction to Mathematical Statistics. 1965, Waltham, MA: Blaisdell Publishing Co, Chapter 14:Google Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.