Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Detecting Spatial-Temporal Clusters of HFMD from 2007 to 2011 in Shandong Province, China

  • Yunxia Liu ,

    Contributed equally to this work with: Yunxia Liu, Xianjun Wang

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Shandong University, Jinan, China

  • Xianjun Wang ,

    Contributed equally to this work with: Yunxia Liu, Xianjun Wang

    Affiliation Shandong Center for Disease Control and Prevention, Jinan, China

  • Yanxun Liu,

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Shandong University, Jinan, China

  • Dapeng Sun,

    Affiliation Shandong Center for Disease Control and Prevention, Jinan, China

  • Shujun Ding,

    Affiliation Shandong Center for Disease Control and Prevention, Jinan, China

  • Bingbing Zhang,

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Shandong University, Jinan, China

  • Zhaohui Du,

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Shandong University, Jinan, China

  • Fuzhong Xue

    xuefzh@sdu.edu.cn

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Shandong University, Jinan, China

Abstract

Background

Hand, foot, and mouth disease (HFMD) has caused major public health concerns worldwide, and has become one of the leading causes of children death. China is the most serious epidemic area with a total of 3,419,149 reported cases just from 2008 to 2010, and its different geographic areas might have different spatial epidemiology characteristics at different spatial-temporal scale levels. We conducted spatial and spatial-temporal epidemiology analysis to HFMD at county level in Shandong Province, China.

Methods

Based on the China National Disease Surveillance Reporting and Management System, the spatial-temporal database of HFMD from 2007 to 2011 was built. The global autocorrelation statistic (Moran’s I) was first used to detect the spatial autocorrelation of HFMD cases in each year. Purely Spatial scan statistics combined with Space-time scan statistic were used to detect epidemic clusters.

Results

The annual average incidence rate was 93.70 per 100,000 in Shandong Province. Most HFMD cases (93.94%) were aged within 0–5 years old with an average male-to-female sex ratio 1.71, and the incidence seasonal peak was between April and July. The dominant pathogen was EV71 (47.35%), and CoxA16 (26.59%). HFMD had positive spatial autocorrelation at medium spatial scale level (county level) with higher Moran’s I from 0.31 to 0.62 (P<0.001). Seven spatial-temporal clusters were detected from 2007 to 2011 in the landscape of the whole Shandong, with EV71 or CoxA16 as the dominant pathogen for most hotspots areas.

Conclusions

The spatial-temporal clusters of HFMD wandered around the whole Shandong Province during 2007 to 2011, with EV71 or CoxA16 as the dominant pathogen. These findings suggested that a real-time spatial-temporal surveillance system should be established for identifying high incidence region and conducting prevention to HFMD timely.

Introduction

Hand, foot, and mouth disease (HFMD) is a common communicable disease which usually affects children, particularly those less than 5 years old. It is characterized by a distinct clinical presentation of fever, or vesicular exanthema on their hands, feet, mouths, or buttocks. HFMD can be caused by numerous human enteroviruses (EV) [1], with Coxsackievirus A16 (CoxA16) and Enterovirus 71(EV71) the major causative agents [2][6]. The transmission of HFMD occurs from person to person through direct contact with saliva, faeces, vesicular fluid or respiratory droplets of an infected person and indirectly by contaminated articles. At present, there are still no available effective vaccines or drugs against HFMD human use.

Recently, numerous large-scale outbreaks of HFMD in East and Southeast Asia [7][11], especially deaths caused by EV71 [12], [13], have caused major public health concerns worldwide. In particular, the epidemic situation of HFMD in China is quite serious. It has become one of the leading causes of child death and a public health priority in China [14]. There were 3419,149 HFMD cases and 1,384 fatal cases during 2008–2010 from the report of China’s Health Ministry.

Several outbreaks, mainly caused by CoxA16 and/or EV71, have been reported since 2007 in China, such as Linyi in Shandong (2007) [15], Fuyang in Anhui (2008) [16], Shanghai (2009) [17], Nanchang in Jiangxi (2010) [18], etc., suggesting that the incidence of HFMD might have variability in different regions and times. Therefore, a better understanding of the spatial-temporal distribution patterns of HFMD would help to identify areas and population at high risk, and then to formulate and take appropriate regional public health intervention strategy to prevent and control the outbreak. For example, one study found that the occurrence of HFMD in Ningbo, Zhejiang Province, had an apparent seasonal distribution with a peak in June [5]. Recently, Shandong Province has been suffering from serious HFMD epidemic, the reported cases were totally up to 448,251 and ranked top 5 among 31 provinces in China during 2007 to 2011. However, the distribution of HFMD at medium spatial scale level (county level in China) is still not clear. Thus, based on the reported data of HFMD from the China Information System for Disease Control and Prevention (CISDCP, http://www.cdpc.chinacdc.cn), we conducted the spatial and space-time scan statistics analysis in Shandong Province to explore the distribution characteristics and detect spatial and spatial-temporal clusters (hotspots) of HFMD cases.

Materials and Methods

Data Collection

Shandong Province, located between latitude 34°25′ and 38°23′ north, and longitude 114°36′ and 112°43′ east, is a coastal province in Eastern China with a population of approximately 98 million people (Figure 1). It includes 140 counties (subdistricts) belonging to 17 regions (municipal districts) with a total land area of 156,700 square kilometers.

thumbnail
Figure 1. The location of study area, Shandong Province in China.

https://doi.org/10.1371/journal.pone.0063447.g001

Data of HFMD in Shandong Province during 2007 to 2011 were derived from CISDCP, including the basic social-demographic characteristics of HFMD cases, and the pathogen type (CoxA16, EV71 and other EV) of some HFMD cases. The diagnosis was based on the clinical criteria from the HFMD Control and Prevention Guide published by the Chinese Ministry of Health [19]. The report cards of HFMD cases were filled out by professional doctors, and collected by trained reporter and then input into the CISDCP within 24 hours based on the P. R. China infectious disease prevention and cure statute. EV71, CA16 and other enteroviruses were tested by real-time PCR in 17 laboratories located in 17 cities in Shandong guided by Shandong CDC, the pathogen data was also uploaded to the China CISDCP. In this paper, we focused on the children HFMD cases (0 to 5 years) which accounted for about 94% of the total cases (described in the results section). The corresponding demographic data of each county between 2007 and 2011 were obtained from Shandong Statistical Yearbook. All collected data were geographically referenced based on 140 counties of Shandong Province, i.e., 140 spatial units for analysis.

Statistical Analysis

The frequencies of HFMD were summarized monthly and annually by geographic area (i.e., county). The 0 to 5 years incidence rates of HFMD were calculated by HFMD counts aged 0 to 5 years divided by the population aged 0 to 5 years, which were expressed as the number of cases per 100,000. And the incidence rates of all ages were calculated using the total number of HFMD cases and total population.

The autocorrelation statistic (Moran’s I) [20] was used to detect the global spatial autocorrelation of HFMD cases in the study area and disclose the spatial pattern of HFMD with Z score at county level. The significance of Moran’s I was assessed by employing Monte Carlo randomization. A statistically significant (Z score ≥1.96) estimate of I indicates that neighboring counties have a similar prevalence rate of HFMD and the cases are likely to cluster at county level. The software GeoDa™ 0.9.5-i was used to conduct the analysis [21].

Scan statistics [22] were used to determine the presence of high rates clusters of HFMD. (1) Spatial scan statistic, based on discrete Poisson model, was used to identify purely spatial clusters of HFMD cases by year. The purely spatial scan statistic imposes a circular (or elliptic) window which was in turn centered on each geographical area throughout the study region. The radius of the window varies continuously in size according to the population range of the area (from zero to some upper limit specified by user). (2) Space-time scan statistic, based on Space-Time Permutation model, was adopted to determine the presence of space-time clusters of HFMD cases during the study period. Space-Time Permutation scan statistic, which is developed for the early detection of disease outbreaks and is especially useful when we have only the number of cases and no population-at-risk information [23]. This approach was defined by a cylindrical window with a circular (or elliptic) geographic base and with height corresponding to time, which undergo dynamic changes in space and time to detect possible spatial-temporal clusters [24]. The space-time scan statistic automatically adjusted for both purely spatial and purely temporal variation. The base was defined exactly as for the purely spatial scan statistic, while the height changes according to the defined time interval (less than or equal to half the total study period) which reflects the time period of potential clusters. In two scan statistics mentioned above, the null hypothesis is that the rate (Poisson model), or the independence of cases in space and time (Space-Time Permutation model), was the same within and outside the scanning window. For each window, a likelihood ratio and the relative risk were calculated to test the hypothesis. And the P values for detected clusters were calculated by using Monte Carlo hypothesis testing to generate a number of random replications of the data set under the appropriate null hypothesis. The scanning window with the maximum likelihood constituted most likely cluster, other windows for which the likelihood value was statistically significant were defined as secondary clusters ranked according to their likelihood ratio test statistic. The statistical analyses were done by the free SaTScan™ v9.1.1 software.

In this study, we used 140 counties of Shandong Province as spatial units, 60 months from January 2007 to December 2011 as time unit. In order to scan for small to large clusters, the largest radius was set to 50% of the total population at risk, the largest height was set to 50% of the total study period. To ensure excellent statistical power and consider the computation times, 999 Monte Carlo replications were set, and clusters with statistical significance of P<0.05 were all reported, including secondary clusters that did not overlap with a previously reported cluster (i.e., they had no location IDs in common). Furthermore, we used ArcGIS v9.0 (ESRI, Redlands, CA, USA) to visualize the results of scan statistic analysis.

Results

Prevalence of HFMD

There were 448,251 HFMD cases reported in Shandong Province from 2007 to 2011, including 9,006 severe and 79 fatal cases. The annual average incidence rate of HFMD was 93.70 per 100,000 (ranged from 35.34 per 100,000 in 2008 to 148.82 per 100,000 in 2010).

Table 1 summarized the social-demographic characteristics of HFMD cases and the pathogen types of some cases from 2007 to 2011. It showed that 0–5 year old age group was majority of the victims in the outbreaks, which accounted for 93.94% (421,072 cases) of all reported cases over the study period, with the annual average incidence rate 1442.50 (ranged from 550.11 in 2008 to 2249.54 in 2009) per 100,000. Thus, we mainly focused on 0–5 year old age group in the following spatial analysis.

thumbnail
Table 1. Social-demographic characteristics of HFMD cases and the pathogen types of some cases in Shandong Province, 2007–2011.

https://doi.org/10.1371/journal.pone.0063447.t001

Of 421,488 HFMD cases, 266,185 were boys and 155,303 were girls, with an average male-to-female sex ratio 1.71 (1.73 in 2007, 1.82 in 2008, 1.68 in 2009, 1.73 in 2010, and 1.69 in 2011). In China, all 0–5 year old age children spent their daytime at home, kindergarten or school. Table 1 displayed that most of HFMD cases were preschoolers (74.52% scattered children and 25.35% nursery children), the rest (0.13%) were students. From 2008, pathogen detection for HFMD had been conducting using nucleic acid testing in Shandong Province. Among 13,532 genotyped cases between 2008 and 2011, CoxA16, EV71 and other EV accounted for 30.82%, 45.48% and 23.71% respectively. Clearly, the annual distribution of HFMD pathogens varied heavily (Table 1).

Figure 2 illustrated the monthly distribution of HFMD cases, which indicated that the occurrence of HFMD presented significant seasonality. It was obvious that the incidence peak appeared between April and July in the study years.

thumbnail
Figure 2. Monthly distribution of HFMD cases (0–5 years), 2007–2011.

https://doi.org/10.1371/journal.pone.0063447.g002

Spatial Autocorrelation of HFMD Cases

Figure 3 showed the annual incidence rate of 0–5 year old age per county accounting for the spatial variability of population size. It clearly indicated that the distribution of HFMD was heterogeneous at county level. Table 2 listed the results of the spatial autocorrelation test, which demonstrated that high global spatial autocorrelation of HFMD was detected at county level in Shandong Province within each epidemic year during 2007 to 2011 (Moran’s I >0.3, P<0.001).

thumbnail
Figure 3. Annual incidence rates of HFMD in 0–5 years old children (/100,000) per county in Shandong Province, 2007–2011.

https://doi.org/10.1371/journal.pone.0063447.g003

thumbnail
Table 2. The results of the spatial autocorrelation test on HFMD cases in Shandong Province, 2007–2011.

https://doi.org/10.1371/journal.pone.0063447.t002

Spatial Clusters of HFMD

Figure 4 showed the statistically significant spatial clusters (including the most likely cluster and several secondary clusters) for a high occurrence of HFMD identified by purely Spatial scan statistic based on discrete Poisson model. It was obvious that the locations and sizes of these clusters varied by year. The most likely clusters were in Qingdao area (2007, 2010), Jinan area (2010), Linyi area (2008, 2010), Liaocheng area (2009), Zibo, Dongying and Binzhou areas (2010, 2011), Yantai, Weifang, Taian, Weihai, Rizhao, Laiwu, Dezhou areas (2011). Table 3 summarized the size, relative risk, and P value, etc. of most likely clusters, indicating hotspots of HFMD wandered around in the whole Shandong from 2007 to 2011.

thumbnail
Figure 4. The detected purely spatial clusters of HFMD in Shandong Province, 2007–2011.

https://doi.org/10.1371/journal.pone.0063447.g004

thumbnail
Table 3. The most likely high risk clusters of HFMD cases detected using the purely spatial analysis.

https://doi.org/10.1371/journal.pone.0063447.t003

Space-time Clusters of HFMD

Using Space-Time Permutation model, we detected the statistically significant monthly spatial-temporal clusters for a high occurrence of HFMD during 2007 to 2011, shown in Table 4 and Figure 5. Seven spatial-temporal clusters were detected from 2007 to 2011, including 1 most likely cluster and 6 secondary clusters. The first (most likely) cluster was located in southwest Shandong from February 2009 to April 2010 (cluster 1, Radius = 97.01 km) with 36235 cases (RR = 1.83, P<0.0001). The dominant pathogen was EV71 (57.81%). Followed by cluster in north Shandong from May 2011 to October 2011 (cluster 2, Radius = 53.84 km) with 14257 cases (RR = 2.46, P<0.0001) and EV71 as the dominant pathogen (64.37%); cluster in south Shandong from May 2007 to May 2009 (cluster 3, Radius = 59.85 km) with 23147 cases (RR = 1.74, P<0.0001) and EV71 as dominant (68.42%); cluster in northeast Shandong from June to October 2007 (cluster 4, Radius = 145.51 km) with 10488 cases (RR = 2.23, P<0.0001) and no pathogen data; cluster in northwest Shandong from June 2010 to July 2010 (cluster 5, Radius = 82.40 km) with 12097 cases (RR = 1.68, P<0.0001) and CoxA16 as dominant (50.09%); cluster in middle west Shandong from September 2008 to January 2009 (cluster 6, Radius = 28.44 km) with 1234 cases (RR = 6.88, P<0.0001) and EV71 as dominant (83.33%); cluster in middle east Shandong from June 2010 to July 2010 (cluster 7, Radius = 49.09 km) with 10180 cases (RR = 1.70, P<0.0001) and combined EV71 (29.45%) with CoxA16 (26.71%) as dominant. The above spatial pattern of clusters indicated that the occurrence of HFMD was with spatial-temporal heterogeneity. In chronological order, the earliest high occurrence area was cluster 4, followed by cluster 3, cluster 6, cluster 1, cluster 5, cluster 7 and cluster 2. Although different HFMD pathogen existed in different cluster areas, EV71 was the dominant pathogen for most clusters.

thumbnail
Figure 5. The detected spatial-temporal clusters of HFMD in Shandong Province during 2007 to 2011, and the distribution of HFMD pathogens within these clusters.

Each pie shows the proportion of different HFMD pathogen for each spatial-temporal cluster, but no pathogen data in 2007.

https://doi.org/10.1371/journal.pone.0063447.g005

thumbnail
Table 4. The spatial-temporal high risk clusters of HFMD cases detected using Space-Time Permutation model, 2007–2011.

https://doi.org/10.1371/journal.pone.0063447.t004

Discussion

Shandong Province has been one of the most serious HFMD epidemic areas since 2007, with annual average incidence rate 93.70 per 100,000. Most HFMD cases (93.94%) were aged less than 5 years old of both scattered children (74.52%) and nursery children (25.35%), with an average male-to-female sex ratio 1.71, and the incidence single peak season was between April and July (see Table 1 and Figure 2). The dominant pathogen was EV71 (45.48%), followed by CoxA16 (30.82%), and other EVs (23.70%). Although these epidemiological characteristics was similar to the epidemic areas in other districts of China, some different epidemic features emerged in different districts [25]. For example, though with similar pathogen pattern to Shandong Province, there were double seasonal peaks with the highest occurrence between April and June and the second occurring in November in Jiangsu Province [26]. In Hong Kong, both warmer seasonal peak (May-July) and winter peak (October-December) were detected with EV71 as dominant pathogen, but the number of older children (>5 years) infected increased from 25.4% in 2001 to 33.0% in 2009 [27]. In Guangdong Province, HFMD incidence peaked in April/May and September/October, also with EV71 (22.4%) and CoxA16 (23%) as dominant pathogen [28]. These difference might be partly attributed to climatic, geographic, social factors, etc. [28][30].

In Shandong Province, HFMD had positive spatial autocorrelation at medium spatial scale level (county level) with higher Moran’s I from 0.31 to 0.62 (P<0.001), and the high incidence areas wandered around in the whole Shandong from 2007 to 2011 (see Table 2 and Figure 3). This was similar to the spatial autocorrelation patterns in the whole mainland China at province spatial scale level [31]. Further spatial autocorrelation analysis should be done to demonstrate its spatial epidemic behavior at much smaller spatial scale level, including town and village level.

The hotspots of HFMD wandered around in the whole Shandong from 2007 to 2011 from the results of Spatial scan analysis (see Figure 4). Furthermore, using spatial-temporal analysis, 7 spatial-temporal clusters were detected, including 1 most likely cluster and 6 secondary clusters with different pathogen proportion, indicating that the occurrence of HFMD was with spatial-temporal heterogeneity (see Figure 5). Therefore, HFMD propagates in a composite space-time domain rather than showing a purely spatial and purely temporal variation in Shandong. In addition, EV71 was the dominant pathogen for most clusters, though different HFMD pathogen existed in different hotspot areas (see Figure 5). Analogously, the spatial-temporal distribution of HFMD was also nonrandom in Jiangsu Province [26] and in the landscape of whole mainland China [32]. A study had shown that this spatial heterogeneity was associated with the monthly precipitation types of the region [33]. For establishing the appropriate regional prevention strategy and measure, further study need to be done on risk factors of HFMD.

In this study, Scan statistics were used to detect the spatial and spatial-temporal clusters. The method had been widely used to detect clusters of infectious disease and may serve as a useful adjunct to disease surveillance, particularly in areas of limited resources, e.g., malaria [34], dengue [35], Japanese encephalitis [36], tuberculosis [37] and HFMD [32], etc. However, the Scan statistics are based on the assumption of circular spatial scanning windows and space-time cylinders to detect clusters, while the actual shapes of clusters were not all like that. This limitation prompts us to explore better statistic model to conduct further spatial epidemiology analysis.

In summary, this study highlighted the spatial epidemiological characteristics of HFMD at medium spatial scale level from 2007 to 2011 in Shandong Province, China. It indicated that the spatial-temporal hotspots of HFMD wandered around in whole Shandong from 2007 to 2011 with EV71 or CoxA16 as the dominant pathogen. These findings suggested that a real-time spatial-temporal surveillance system should be established for identifying high incidence region and for conducting prevention to HFMD timely.

Author Contributions

Conceived and designed the experiments: YXL XJW FX. Performed the experiments: YXL YL BBZ FX. Analyzed the data: YXL BBZ ZHD. Contributed reagents/materials/analysis tools: XJW DPS SJD. Wrote the paper: YXL.

References

  1. 1. Gopalkrishna V, Patil PR, Patil GP, Chitambar SD (2012) Circulation of multiple enterovirus serotypes causing hand, foot and mouth disease in India. J Med Microbiol 61: 420–425.
  2. 2. Yang F, Zhang T, Hu Y, Wang X, Du J, et al. (2011) Survey of enterovirus infections from hand, foot and mouth disease outbreak in China, 2009. Virol J 8: 508.
  3. 3. Wei SH, Huang YP, Liu MC, Tsou TP, Lin HC, et al. (2011) An outbreak of coxsackievirus A6 hand, foot, and mouth disease associated with onychomadesis in Taiwan, 2010. BMC Infect Dis 11: 346.
  4. 4. Bracho MA, González-Candelas F, Valero A, Córdoba J, Salazar A (2011) Enterovirus coinfections and onychomadesis after hand, foot, and mouth disease, Spain, 2008. Emerg Infect Dis 17: 2223–2231.
  5. 5. Ni H, Yi B, Yin J, Fang T, He T, et al. (2012) Epidemiological and etiological characteristics of hand, foot, and mouth disease in Ningbo, China, 2008–2011. J Clin Virol 54: 342–348.
  6. 6. Gopalkrishna V, Patil PR, Patil GP, Chitambar SD (2012) Circulation of multiple enterovirus serotypes causing hand, foot and mouth disease in India. J Med Microbiol 61: 420–425.
  7. 7. Tu PV, Thao NT, Perera D, Huu TK, Tien NT, et al. (2007) Epidemiologic and Virologic investigation of hand, foot, and mouth disease, southern Vietnam, 2005. Emerg Infect Dis 13: 1733–1741.
  8. 8. Fujimoto T, Chikahira M, Yoshida S, Ebira H, Hasegawa A, et al. (2002) Outbreak of central nervous system disease associated with hand, foot, and mouth disease in Japan during the summer of 2000: detection and molecular epidemiology of enterovirus 71. Microbiol Immunol 46: 621–627.
  9. 9. Wang Y, Feng Z, Yang Y, Self S, Gao Y, et al. (2011) Hand, foot, and mouth disease in China: patterns of spread and transmissibility. Epidemiology 22: 781–792.
  10. 10. Ang LW, Phoon MC, Wu Y, Cutter J, James L, et al. (2011) The changing seroepidemiology of enterovirus 71 infection among children and adolescents in Singapore. BMC Infect Dis 11: 270.
  11. 11. Puenpa J, Theamboonlers A, Korkong S, Linsuwanon P, Thongmee C, et al. (2011) Molecular characterization and complete genome analysis of human enterovirus 71 and coxsackievirus A16 from children with hand, foot and mouth disease in Thailand during 2008–2011. Arch Virol 156: 2007–2013.
  12. 12. Chong CY, Chan KP, Shah VA, Ng WY, Lau G, et al. (2003) Hand, foot and mouth disease in Singapore: acomparison of fatal and non-fatal cases. Acta Paediatr 92: 1163–1169.
  13. 13. McMinn PC (2002) An overview of the evolution of enterovirus71 and its clinical and public health significance. FEMS Microbiol Rev 26: 91–107.
  14. 14. Chinese Center for Disease Control and Prevention (China CDC) (2010) National Incidence and death cases of notifiable class A or class B infectious diseases. (http://www.chinacdc.cn).
  15. 15. Zhang Y, Tan XJ, Wang HY, Yan DM, Zhu SL, et al. (2009) An outbreak of hand, foot, and mouth disease associated with subgenotype C4 of human enterovirus 71 in Shandong, China. J Clin Virol 44: 262–267.
  16. 16. Zhang Y, Zhu Z, Yang W, Ren J, Tan X, et al. (2010) An emerging recombinant human enterovirus 71 responsible for the 2008 outbreak of hand foot and mouth disease in Fuyang city of China. Virol J 7: 94.
  17. 17. Yan XF, Gao S, Xia JF, Ye R, Yu H, et al. (2012) Epidemic characteristics of hand, foot, and mouth disease in Shanghai from 2009 to 2010: Enterovirus 71 subgenotype C4 as the primary causative agent and a high incidence of mixed infections with coxsackievirus A16. Scand J Infect Dis 44: 297–305.
  18. 18. Liu MY, Liu W, Luo J, Liu Y, Zhu Y, et al. (2011) Characterization of an outbreak of hand, foot, and mouth disease in Nanchang, China in 2010. PLOS One 6: e25287.
  19. 19. Chinese Ministry of Health. Hand, Foot and Mouth Disease Control and Prevention Guide, 2008.
  20. 20. Pfeiffer D, Robinson T, Stevenson M, Stevens K, Rogers D, et al.. (2008) A spatial analysis in epidemiology. New York: Oxford university press. 35–37 p.
  21. 21. Anselin,L(2004). GeoDa0.95i Release Notes. Spatial Analysis Laboratory (SAL). Department of Agricultural and Consumer Economics, University of Illinois, Urbana-Champaign, IL.
  22. 22. Kulldorff M (1997) A spatial scan statistic. Communications in Statistics - Theory and Methods 26: 1481–1496.
  23. 23. Kulldorff M, Heffernan R, Hartman J, Assunção R, Mostashari F (2005) A space-time permutation scan statistic for disease outbreak detection. PLoS Med 2: e59.
  24. 24. Kulldorff M, Athas WF, Feuer EJ, Miller BA, Key CR (1998) Evaluating cluster alarms: A space-time scan statistic and brain cancer in Los Alamos, New Mexico. Am J Public Health 88: 1377–1380.
  25. 25. Chang ZR, Zhang J, Sun JL, Zhang WD, Wang ZJ (2011) Epidemiological features of hand, foot and mouth disease in China, 2008–2009. Zhonghua Liu Xing Bing Xue Za Zhi 32: 676–680.
  26. 26. Liu WD, Wu Y, Liang Q, Hu JL, Dai QG, et al. (2012) Epidemiological characteristics and temporal-spatial clustering analysis on hand-foot-mouth disease in Jiangsu province, 2009–2011. Zhonghua Liu Xing Bing Xue Za Zhi 33: 813–817.
  27. 27. Ma E, Lam T, Chan KC, Wong C, Chuang SK (2010) Changing epidemiology of hand, foot, and mouth disease in Hong Kong, 2001–2009. Jpn J Infect Dis 63: 422–426.
  28. 28. Zou XN, Zhang XZ, Wang B, Qiu YT (2012) Etiologic and epidemiologic analysis of hand, foot, and mouth disease in Guangzhou city: a review of 4,753 cases. Braz J Infect Dis 16: 457–465.
  29. 29. Cao LH, Ren M, Zhao PL, Ma JB, Sun SL, et al. (2011) A exploration and study of the relationships of hand-foot-mouth disease (HFMD) and the climate. Zhonghua Shi Yan He Lin Chuang Bing Du Xue Za Zhi 25: 227–229.
  30. 30. Onozuka D, Hashizume M (2011) The influence of temperature and humidity on the incidence of hand, foot, and mouth disease in Japan. Sci Total Environ 410–411: 119–125.
  31. 31. Xiao GX, Hu YH, Ma JQ, Hao YT, Wang XF, et al. (2012) Spatial clustering and changing trend of hand-foot-mouth disease during 2008–2011 in China. Zhonghua Liu Xing Bing Xue Za Zhi 33: 808–812.
  32. 32. Zhu Q, Hao Y, Ma J, Yu S, Wang Y (2011) Surveillance of hand, foot, and mouth disease in mainland China (2008–2009). Biomed Environ Sci 24: 349–356.
  33. 33. Wang JF, Guo YS, Christakos G, Yang WZ, Liao YL, et al. (2011) Hand, foot and mouth disease: spatiotemporal transmission and climate. Int J Health Geogr 10: 25.
  34. 34. Coleman M, Coleman M, Mabuza AM, Kok G, Coetzee M, et al. (2009) Using the SaTScan method to detect local malaria clusters for guiding malaria control programmes. Malar J 8: 68.
  35. 35. Morrison AC, Getis A, Santiago M, Rigau-Perez JG, Reiter P (1998) Exploratory space-time analysis of reported dengue cases during an outbreak in Florida, Puerto Rico, 1991–1992. Am J Trop MedHyg 58: 287–298.
  36. 36. Tang XY, Gan WY, Xu B, Chen CY, Chu M, et al. (2011) Study on the spatiotemporal trend of Japanese encephalitis in Guangxi, based on geographic information system and space-time permutation scan statistic. Zhonghua Liu Xing Bing Xue Za Zhi 32: 274–278.
  37. 37. Higgs BW, Mohtashemi M, Grinsdale J, Kawamura LM (2007) Early detection of tuberculosis outbreaks among the San Francisco homeless: trade-offs between spatial resolution and temporal scale. PLoS One 2: e1284.