- Open Access
Monitoring county-level chlamydia incidence in Texas, 2004 – 2005: application of empirical Bayesian smoothing and Exploratory Spatial Data Analysis (ESDA) methods
© Owusu-Edusei and Owens; licensee BioMed Central Ltd. 2009
- Received: 02 January 2009
- Accepted: 26 February 2009
- Published: 26 February 2009
Chlamydia continues to be the most prevalent disease in the United States. Effective spatial monitoring of chlamydia incidence is important for successful implementation of control and prevention programs. The objective of this study is to apply Bayesian smoothing and exploratory spatial data analysis (ESDA) methods to monitor Texas county-level chlamydia incidence rates by examining spatiotemporal patterns. We used county-level data on chlamydia incidence (for all ages, gender and races) from the National Electronic Telecommunications System for Surveillance (NETSS) for 2004 and 2005.
Bayesian-smoothed chlamydia incidence rates were spatially dependent both in levels and in relative changes. Erath county had significantly (p < 0.05) higher smoothed rates (> 300 cases per 100,000 residents) than its contiguous neighbors (195 or less) in both years. Gaines county experienced the highest relative increase in smoothed rates (173% – 139 to 379). The relative change in smoothed chlamydia rates in Newton county was significantly (p < 0.05) higher than its contiguous neighbors.
Bayesian smoothing and ESDA methods can assist programs in using chlamydia surveillance data to identify outliers, as well as relevant changes in chlamydia incidence in specific geographic units. Secondly, it may also indirectly help in assessing existing differences and changes in chlamydia surveillance systems over time.
- Spatial Dependence
- Local Outlier
- Chlamydia Screening
- Exploratory Spatial Data Analysis
- Monte Carlo Randomization
Chlamydia is the most prevalent reportable disease in the United States with an estimated 2.8 million cases each year [1, 2]. Untreated chlamydial infections in women have been associated with more serious reproductive complications such as pelvic inflammatory disease (PID), ectopic pregnancy, tubal infertility, and chronic pelvic pain [3–6]. In men, chlamydia has been associated with urethritis and other complications such as epididymitis and acute proctitis [7–9]. Thus, it is a public health problem that has attracted public attention, albeit not as much as would be desired.
Several previous studies have recommended that the design and implementation of effective interventions to control or prevent sexually transmitted diseases (STDs) should be grounded on a good understanding of the existing and emerging spatiotemporal patterns because STDs are characterized by geographic patterns [10–16]. An emerging approach to achieving this end is the application of Exploratory Spatial Data Analysis (ESDA) methods which draws from the field of spatial statistics . At the state-level, ESDA methods can be used by state health officials to monitor spatial and temporal variations in rates using counties as spatial units. ESDA can also assist in identifying and monitoring hot spots ("problem counties") that may not be obvious otherwise. These methods can aid health officials to design more location-specific prevention programs that take into account global and local spatial influences. It is also valuable to be able to assess and develop surveillance systems that can immediately and effectively pick up warning signs of increases in any particular STD. The ideas and motivation for the application of these methods to STD were drawn from pioneering works in the area of ESDA by Luc Anselin and others on juvenile crime and cancer rates, among others [18–21].
The primary objective of this study was to use ESDA methods to identify and monitor Bayesian-smoothed chlamydia incidence rates using county-level data from the state of Texas. Our choice of counties as the unit of analysis was based on availability of data. Finer spatial units (cities or census tracts) may provide more location-specific information that can inform the design and implementation stages of existing or future interventions. Majority of chlamydia cases are asymptomatic prompting recommendations for routine screening for young women by individuals and organizations [22–30]. In view of this, differences in the incidence rates may be the result of differences in existing surveillance systems. Thus, indirectly, ESDA may help to identify disparities in chlamydia surveillance systems.
Data used in this study was obtained from the National Electronic Telecommunications System for Surveillance (NETSS) which is maintained by the Centers for Disease Control and Prevention (CDC). We used the overall incidence rates (per 100,000 residents, for all race, sex and age groups) for each county provided by the surveillance system.
Spatial relationship concept
We used the standardized 1st- order Queen Neighbors (all counties that share a border with the referent county) as the criteria for identifying neighbors. Spatial relationship through out this study was executed by the use of a spatial weight matrix.
Empirical Bayesian smoothing
Summary statistics of Bayesian-smoothed chlamydia incidence rates (n = 254)
Measuring spatial dependence
We used Moran's I , a statistical test for global spatial autocorrelation (dependence) in group-level data to identify departures from spatial randomness revealing the existing spatial patterns, such as clusters. The hypotheses for this test are:
Null: Smoothed rates in different regions are spatially independent,
Alternative: Smoothed rates are not spatially independent.
If the resulting value was positive, then there was spatial autocorrelation – nearby areas had similar rates, indicating global spatial clustering. Conversely, if the value was negative, then nearby areas were dissimilar. A value close to zero indicates random spatial units. We used Anselin's Local Moran test which was an extension of Moran's test to identify local spatial autocorrelation . This test was used to identify local outliers by comparing counties to their contiguous counties – how different the rates were for any spatial unit (county in this case) from its immediate neighbors.
After computing the appropriate statistic from the smoothed rates, a Monte Carlo Randomization (MCR) procedure was used to recalculate the statistic from the randomized data observations to generate a reference distribution using 999 permutations. The p-values were computed by comparing the observed statistic to the distribution generated by the MCR process. We used Simes Correction  to adjust the p-values to account for the lack of independence in the statistics computed. We used GeoDa (version 0.9.5-i) software application by Luc Anselin, 2004.
The basic steps in ESDA used in this study were:
1. Mapped the Bayesian-smoothed rates in classes for preliminary visual analysis – identified outliers and spatial association.
2. Conducted a statistical test to confirm or reject spatial dependence, and
3. Computed local Moran statistics to map local outliers.
We also analyzed trends by computing relative changes (from 2004 to 2005) for each county to identify global and local outliers using the same steps outlined above.
Bayesian-smoothed chlamydia incidence rates among the ten counties with the highest incidence in Texas, 2004 and 2005
Relative change 2004–2005
A visual inspection of Figure 1 (panels a and b) gives an indication of spatial association – counties with similar color shades had the tendency to be near each other. We carried out the formal statistical test using Moran's I. The test statistics were 0.19 (p < 0.01) and 0.21 (p < 0.01) for 2004 and 2005, respectively. Thus we rejected the null of spatial independence and concluded that there was sufficient evidence of spatial dependence.
The local Moran significance maps (Figure 1, panels d and e) indicated that in both years, smoothed chlamydia incidence rates for one county (Erath county) was significantly (p < 0.05) higher than its contiguous neighbors. Erath county had smoothed rates of more than 300, while its neighbors had 195 or less in the two years examined. For 2004, the smoothed rates for Runnels county was significantly (p < 0.05) lower than its immediate neighbors (see Figure 1, panel d). The rest of the rates for the other counties were either similar (high-high or low-low) to its neighbors or not statistically different from them.
To examine temporal patterns, we computed the relative changes (percent changes) in the smoothed rates from 2004 to 2005. Figure 1, panel c shows a graduated scheme map of the relative change in chlamydia rates. Our formal test also indicated evidence of spatial dependence in the relative changes that occurred from 2004 to 2005 (0.14, p < 0.01).
The highest ten values of relative change in the smoothed rates are also presented in Table 2. Gaines county had the highest relative change. The chlamydia rates for Gaines county increased from 139 to 379 (cases per 100, 000 residents); a 173 percent increase. Except for Gaines county all the other nine counties had z-values within 3 standard deviations of the mean.
Local Moran significance maps
Figure 1, panel f shows local Moran significance maps used to identify changes in smoothed rates for counties that were significantly higher or lower than their neighbors. Percent change in Newton county (15% increase) was significantly (p < 0.05) higher than its contiguous counties (Orange, 54% decrease; Jasper, 6% increase; Sabine, 37% decrease). The relative changes in smoothed rates for Andrews, Glasscock and Mills counties were significantly (p < 0.05) lower than their contiguous counties. The rest were either similar to their immediate neighbors, or not significantly different.
We applied empirical Bayesian smoothing and ESDA methods in GIS to study the most commonly reported STD in the United States – chlamydia. County-level chlamydia incidence rates for 2004 and 2005 in the state of Texas were used to characterize spatiotemporal patterns. Based on data from the National Electronic Telecommunications System for Surveillance (NETSS), our results indicated that empirical Bayesian-smoothed chlamydia rates were spatially dependent for the two years examined. Furthermore, the relative changes that occurred between 2004 and 2005 were also spatially dependent. Six counties (Bell, Falls, Potter, Taylor, Kleberg, and Lubbock) were among the highest ten counties for the two years examined. Also, Erath county had significantly (p < 0.05) higher smoothed rates than its immediate neighbors for the two years examined. The highest relative increase in chlamydia rate was in Gaines county, which experienced over 170% increase in smoothed rates. However, the increase in chlamydia rate in Gaines county was not significantly different from its contiguous neighbors. Relative change in chlamydia rates in Newton county was significantly (p < 0.05) higher than its contiguous counties. The counties identified suggest that they should be considered as the targets for further appraisal. Thus, more detailed examination of the data is required for these counties.
Surveillance data are not perfect. However, for chlamydia, the problem may be more pronounced. Most of chlamydia cases are asymptomatic so the data on incidence may largely be based on adherence to screening recommendations by individuals and organizations [22–30] that vary from county to county and from group to group. Nonetheless, the strength in this type of analyses is that it has the potential to prompt health officials to investigate the data further and subsequently help identify the disparities in the existing screening patterns. For instance, it is possible that Erath county has a relatively better chlamydia screening program compared to its neighbors, or that Gaines county substantially increased screening from 2004 to 2005.
One of the limitations of the Bayesian empirical smoothing method used in this study is that it can potentially overestimate the test statistic for spatial dependence. However, the conclusion reached (i.e., existence of spatial dependence) in this study was the same as in previous studies in which different smoothing methods and tests were used. Previous studies found spatial dependence in chlamydia rates using formal statistical tests on county-level data for 2000 from Texas  and 2000–2002 data on census blocks from Richmond, Virginia . Our review of the literature did not provide any information on the application of a formal test for spatial dependence on the relative changes of chlamydia rates. Thus more studies should be conducted to study the existence of spatial dependence of the relative changes in chlamydia rates.
By default, focusing on any particular spatial jurisdiction, such as a state, precludes one from studying the effects of spatial association with contiguous counties in neighboring states. One limitation of this study was that analyses of border counties did not include spatial effects from the counties in the bordering state. As an example, Newton County which was on the border with Louisiana would have to be examined closely to understand the pattern discovered in this study. Therefore, there was the need to further investigate the spatial relationships that existed for outliers that were situated on the border of the state, as their rates may be the result of interaction with counties in the states bordering them that may have been ignored.
The methods used in this study can be applied to any state, county or city, and for any age group within the chosen spatial/geographic unit. However, where data is available, smaller geographic units are preferable in such analyses. ESDA is one of the methods available for monitoring diseases. There are also other smoothing methods available in the literature, but there are no reports on a comprehensive objective assessment of all the available methods. Therefore, further research is needed in this area. The use of two or more methods on the same dataset may enhance validity if the final results are robust. Additionally, as shown in this study, simple mapping for any geographic units of interest by classes and by associated changes overtime, followed by critical inspection may help to correctly describe the existing spatiotemporal patterns. Identifying and describing the patterns can guide the design and implementation stages of interventions/programs, as well as indirectly help to evaluate existing chlamydia surveillance systems.
The authors thank Charlotte K. Kent, PhD, for constructive comments and suggestions and critical review of the final draft.
The findings and conclusions in this manuscript are those of the authors and do not necessarily represent the views of the Centers for Disease Control and Prevention (CDC). Mention of company names or products does not imply endorsement by CDC.
- Centers for Disease Control and Prevention: Sexually transmitted disease surveillance, 2006. 2007, Atlanta, GA: Centers for Disease Control and PreventionGoogle Scholar
- Weinstock H, Berman S, Cates W: Sexually transmitted diseases among American youth: Incidence and prevalence estimates, 2000. Perspectives on Sexual and Reproductive Health. 2004, 36 (1): 6-10. 10.1363/3600604.PubMedView ArticleGoogle Scholar
- Bakken IJ, Skjeldestad FE, Nordbo SA: Chlamydia trachomatis infections increase the risk for ectopic pregnancy: A population-based, nested case-control study. Sexually Transmitted Diseases. 2007, 34 (3): 166-169. 10.1097/01.olq.0000230428.06837.f7.PubMedView ArticleGoogle Scholar
- Hillis SD, Joesoef R, Marchbanks PA, Wasserheit JN, Cates W, Westrom L: Delayed care of pelvic inflammatory disease as risk factor for impaired fertility. Am J Obstet Gynecol. 1993, 168 (5): 1503-1509.PubMedView ArticleGoogle Scholar
- Scholes D, Stergachis A, Heidrich FE, Andrilla H, Holmes KK, Stamm WE: Prevention of pelvic inflammatory disease by screening for cervical chlamydial infection. New England Journal of Medicine. 1996, 334 (21): 1362-1366. 10.1056/NEJM199605233342103.PubMedView ArticleGoogle Scholar
- Westrom L, Joesoef R, Reynolds G, Hadgu A, Thompson SE: Pelvic inflammatory disease and fertility – a cohort study of 1,844 women with laparoscopically verified disease and 657 control women with normal laparoscopic results. Sexually Transmitted Diseases. 1992, 19 (4): 185-192.PubMedView ArticleGoogle Scholar
- Geisler WM, Krieger JN: Epididymitis. Sexually Transmitted Diseases. Edited by: Holmes KK, Sparling PF, Stamm WE, Piot P, Wasserheit JN, Corey L, Cohen MS, Watts HD. 2008, New York: McGraw Hill, 1127-1146.Google Scholar
- Martin DH: Urethritis in Males. Sexually Transmitted Diseases. Edited by: Holmes KK, Sparling PF, Stamm WE, Piot P, Wasserheit JN, Corey L, Cohen MS, Watts HD. 2008, New York: McGraw Hill, 1107-1126.Google Scholar
- Stamm WE: Chlamydia trachomatis infections in the adult. Sexually Transmitted Diseases. Edited by: Holmes KK, Sparling PF, Stamm WE, Piot P, Wasserheit JN, Corey L, Cohen MS, Watts HD. 2008, New York: McGraw Hill, 575-593.Google Scholar
- Bernstein KT, Curriero FC, Jennings JM, Olthoff G, Erbelding EJ, Zenilman J: Defining core gonorrhea transmission utilizing spatial data. American Journal of Epidemiology. 2004, 160 (1): 51-58. 10.1093/aje/kwh178.PubMedView ArticleGoogle Scholar
- Bush KR, Henderson EA, Dunn J, Read RR, Singh A: Mapping the core: Chlamydia and gonorrhea infections in Calgary, Alberta. Sexually Transmitted Diseases. 2008, 35 (3): 291-297.PubMedView ArticleGoogle Scholar
- Ellen JM, Brown BA, Chung SE, Potterat JJ, Muth SQ, Valente TW, Padian NS: Impact of sexual networks on risk for gonorrhea and chlamydia among low-income urban African American adolescents. Journal of Pediatrics. 2005, 146 (4): 518-522. 10.1016/j.jpeds.2004.11.023.PubMedView ArticleGoogle Scholar
- Jennings JM, Curriero FC, Celentano D, Ellen JM: Geographic identification of high gonorrhea transmission areas in Baltimore, Maryland. American Journal of Epidemiology. 2005, 161 (1): 73-80. 10.1093/aje/kwi012.PubMedView ArticleGoogle Scholar
- Law DCG, Serre ML, Christakos G, Leone PA, Miller WC: Spatial analysis and mapping of sexually transmitted diseases to optimise intervention and prevention strategies. Sexually Transmitted Infections. 2004, 80 (4): 294-299. 10.1136/sti.2003.006700.PubMedPubMed CentralView ArticleGoogle Scholar
- Niccolai LM, Stephens N, Jenkins H, Richardson W, Muth SQ, Rothenberg R: Early syphilis among men in Connecticut: Epidemiologic and spatial patterns. Sexually Transmitted Diseases. 2007, 34 (3): 183-187. 10.1097/01.olq.0000233708.27225.90.PubMedView ArticleGoogle Scholar
- Wylie JL, Cabral T, Jolly AM: Identification of networks of sexually transmitted infection: A molecular, geographic, and social network analysis. Journal of Infectious Diseases. 2005, 191 (6): 899-906. 10.1086/427661.PubMedView ArticleGoogle Scholar
- Anselin L: Exploratory spatial data analysis & geographic information systems. New Tools for Spatial Analysis. Edited by: Painho M. 1994, Luxemburg: Eurostat, 45-54.Google Scholar
- Anselin L: Local indicators of spatial association – LISA. Geographical Analysis. 1995, 27 (2): 93-115.View ArticleGoogle Scholar
- Anselin L: Exploratory spatial data analysis in geocomputational environment. GeoComputation, A Primer. Edited by: Longley P, Brooks S, Macmillan B, McDonnell R. 1998, New York: Wiley, 77-94.Google Scholar
- Anselin L: Interactive techniques & exploratory spatial data analysis. Geographic information systems: Principles, techniques, management & applications. Edited by: Longley PA, Goodchild MF, Maguire DJ, Rhind DW. 1999, New York: Wiley, 251-264.Google Scholar
- Anselin L, Cohen J, Cook D, Gorr W, Tita G: Spatial analyses of crime. Criminal Justice. Edited by: Duffee D. 2000, Washington D. C.: NIJ, 213-262.Google Scholar
- Centers for Disease Control and Prevention: Sexually transmitted diseases treatment guidelines 2002. Centers for Disease Control and Prevention. MMWR Recomm Rep. 2002, 51 (RR-6): 1-78.Google Scholar
- Hollblad-Fadiman K, Goldman SM: American college of preventive medicine practice policy statement – Screening for Chlamydia trachomatis. American Journal of Preventive Medicine. 2003, 24 (3): 287-292. 10.1016/S0749-3797(02)00636-0.PubMedView ArticleGoogle Scholar
- Lin KW, Ramsey L: Screening for chlamydial infection. Am Fam Physician. 2008, 78 (12): 1349-1350.PubMedGoogle Scholar
- Berg AO, Allan JD, Frame PS, Homer CJ, Lieu TA, Mulrow CD, Orleans CT, Peipert JF, Pender NJ, Sox HC: Screening for chlamydial infection – Recommendations and rationale. American Journal of Preventive Medicine. 2001, 20 (3): 90-94. 10.1016/S0749-3797(01)00254-9.View ArticleGoogle Scholar
- Centers for Disease Control and Prevention: Pelvic inflammatory disease: guidelines for prevention and management. MMWR Recomm Rep. 1991, 40 (RR-5): 1-25.Google Scholar
- Committee on Practice and Ambulatory Medicine: Recommendations for preventive pediatric health care. Pediatrics. 1995, 96 (2 Pt 1): 373-374.Google Scholar
- Elster AB, Kuznets NJ: Guidelines for adolescent preventive services. 1994, Baltimore, MD: American Medical AssociationGoogle Scholar
- Hillis SD, Wasserheit JN: Screening for chlamydia – A key to the prevention of pelvic inflammatory disease. New England Journal of Medicine. 1996, 334 (21): 1399-1401. 10.1056/NEJM199605233342111.PubMedView ArticleGoogle Scholar
- U.S. Preventive Services Task Force: Screening for chlamydial infection: U.S. Preventive Services Task Force recommendation statement. Annals of Internal Medicine. 2007, 147 (2): 128-134.View ArticleGoogle Scholar
- Clayton DG, Kaldor J: Empirical Bayesian estimates of age-standardized relative risks for use in disease mapping. Biometrics. 1987, 43: 671-682. 10.2307/2532003.PubMedView ArticleGoogle Scholar
- Waller LA, Gotway CA: Applied Spatial Statistics for Public Health Data. 2004, Hoboken, NJ: John Wiler & SonsView ArticleGoogle Scholar
- Moran PAP: Notes on continuous stochastic phenomena. Biometrika. 1950, 37: 17-23.PubMedView ArticleGoogle Scholar
- Simes RJ: An improved Bonferroni procedure for multiple tests of significance. Biometrika. 1986, 73: 751-754. 10.1093/biomet/73.3.751.View ArticleGoogle Scholar
- Owusu-Edusei K, Chesson HW: Association between County-Level Racial Composition and Reported Cases of Chlamydia and Gonorrhea: Application of Spatial Econometric Models. 2008 National STD Prevention Conference. Chicago, IL. 2008Google Scholar
- Delcher PC, Stover J: Geographic and Spatial Regression Analysis of Sexually Transmitted Diseases in Richmond, Virginia. 2006 National STD Prevention Conference. Jacksonville, FL. 2006Google Scholar
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.