Different effects of meteorological factors on hand, foot and mouth disease in various climates: a spatial panel data model analysis
BMC Infectious Diseases volume 16, Article number: 233 (2016)
Major outbreaks of hand, foot and mouth disease (HFMD) have been reported in China since 2008, posing a great threat to the health of children. Although many studies have examined the effect of meteorological variables on the incidence of HFMD, the results have been inconsistent. This study aimed to quantify the relationship between meteorological factors and HFMD occurrence in different climates of mainland China using spatial panel data models.
All statistical analyses were carried out according to different climate types. We firstly conducted a descriptive analysis to summarize the epidemic characteristics of HFMD from May 2008 to November 2012 and then detected the spatial autocorrelation of HFMD using a global autocorrelation statistic (Moran’s I) in each month. Finally, the association between HFMD incidence and meteorological factors was explored by spatial panel data models.
The 353 regions were divided into 4 groups according to climate (G1: subtropical monsoon climate; G2: temperate monsoon climate; G3: temperate continental climate; G4: plateau mountain climate). The Moran’s I values were significant with high correlations in most months of group G1 and G2 and some months of group G3 and G4. This suggested the existence of a high spatial autocorrelation with HFMD. Spatial panel data models were more appropriate to describe the data than fixed effect models. The results showed that HFMD incidences were significantly associated with average atmospheric pressure (AAP), average temperature (AT), average vapor pressure (AVP), average relative humidity (ARH), monthly precipitation (MP), average wind speed (AWS), monthly total sunshine hours (MSH), mean temperature difference (MTD), rain day (RD) and average temperature distance (ATD), but the effect of meteorological factors might differ in various climate types.
Spatial panel data models are useful and effective when longitudinal data are available and spatial autocorrelation exists. Our findings showed that meteorological factors were related to the occurrence of HFMD, which were also affected by climate type.
Hand, foot and mouth disease (HFMD) is a common infectious disease, which is mainly caused by the enteroviruses coxsackie A16 and enterovirus 71 . In most cases, the disease is mild and self-limiting, but sometimes serious neurological and cardiopulmonary complications may occur, particularly among those aged 5 years and younger . HFMD was listed as a notifiable Class-C communicable disease since May 2008 . In recent years, the trend of HFMD outbreaks has increased among children in China, which is regarded as an important public health problem [4, 5].
Many studies have been performed to analyze the association between meteorological factors and HFMD, but the results have been inconsistent. It has been found that temperature and relative humidity were positively associated with HFMD in most studies [6–11], whereas studies in Japan found that the number of days per week with an average temperature above 25°C was negatively associated with HFMD incidence . Wind velocity was found to be positively associated with HFMD in Ma  and Li’s studies , but in Huang’s publication , no relationship was found between wind speed and HFMD. Furthermore, Wang  found a negative association. Other possible risk factors (total sunshine, difference in temperature, atmospheric pressure, vapor pressure, etc.) were examined in a few studies and require further research.
Most previous studies have not considered temporal or spatial effects (Multiple Linear Regression Model /Negative Binomial Regression /Generalized Additive Poisson Model ), or merely focused on the spatial dimension (Geographically Weighted Regression Models (GWR) ) or on time dimension (time-series analysis ) approaches, which might cause loss of information by ignoring the heterogeneity in both time and space and result in different conclusions.
Compared to traditional methods on the basis of time-series or cross-sectional data alone, spatial panel data models can control for both spatial dependency and unknown heterogeneity [17, 18]. In this study, spatial panel data models were used to explore the relationship between meteorological variables and HFMD incidence in 353 regions of mainland China according to climate types.
Thirty-one provinces of mainland China were divided into 353 regions in this study. There were four main types of climate in mainland China. All the data were grouped by climate types, as shown in Fig. 1: subtropical monsoon climate (SMC, G1), temperate monsoon climate (TMC, G2), temperate continental climate (TCC, G3) and plateau mountain climate (PMC, G4).
Surveillance data of HFMD
Data of monthly reported HFMD cases in each region in mainland China from May 2008 to November 2012 were obtained from the National Center for Public Health Surveillance and Information Services, at the China Center for Disease Control and Prevention (China CDC). All cases were diagnosed according to the clinical criteria provided in a guidebook published by the National Health and Family Planning Commission of the People’s Republic of China in 2008 : with or without fever, a probable case of HFMD was defined as a patient with a papular or vesicular rash on or in the hands, feet, mouth, and/or buttocks. A confirmed case was defined as a probable case with laboratory evidence of enterovirus infection. The demographic data for each region were obtained from the National Bureau of Statistics of China . The ethical approval and the consent from each individual subject was not required because we used only aggregated data (number of cases for each county in months) but not any individualized data in this study.
Meteorological data from 328 monitoring stations that were nearest to the centers of 353 regions were obtained from the China Meteorological Data Sharing Service System (http://cdc.nmic.cn/home.do), which was publicly accessible. Monthly meteorological variables in this study included average atmospheric pressure (AAP), average temperature (AT), average vapor pressure (AVP), average relative humidity (ARH), monthly precipitation (MP), average wind speed (AWS), monthly total sunshine hours (MSH), mean temperature difference (MTD), monthly rainfall days (MRD), and average temperature distance (ATD).
Global spatial autocorrelation analysis
The autocorrelation statistic (global Moran’s I)  was used to detect the global spatial autocorrelation of reported HFMD incidences in the study area by climate types. The significance of Moran’s I was assessed by employing Monte Carlo randomization. A higher positive Moran’s I value with a statistically significant p-value (P < 0.05) indicates that the values of neighboring areas tend to cluster . ArcGIS 10.1 (ESRI Inc. Redlands, CA, USA) was used to perform the analysis.
Spatial panel data models
Spatial panel data typically refer to data containing continuous observations of a number of spatial units. Spatial panel models, which could address data with spatial dependence and also enable researchers to consider spatial and/or temporal heterogeneity, were used to examine the role of different meteorological factors in this study. Spatial panel data models are more informative and contain more variation and less collinearity between the variables compared with cross-sectional or time series models .
The basic form of the simple panel data model with a spatial and temporal specific effect is:
where i and t are indices for the cross-sectional dimension (spatial units) and time dimension (time periods), respectively; y it is the dependent variable at i and t; X it is the group of explanatory variables; β is the vector of regression coefficients that explains the relationship between X it and y it ; ε it is an independently and identically distributed error term with zero mean and variance σ2; μ i denotes a spatial specific effect and γ t represents temporal specific effects. The spatial and/or temporal specific effects may be treated as fixed effects or random effects. A random effect is appropriate if a certain number of individuals is randomly sampled from a large population. The fixed effect model is favored [23–25] when the regression analysis is applied to a precise set of regions. For this reason, because our data contained all regions of the study area, we established fixed effects panel data models that included spatial error autocorrelation or a spatially lagged dependent variable. We also compared the random effects specification against fixed effects specification by Housman’s specification test [26, 27], which suggested that fixed effect specification was more appropriate.
The simple panel data models with specific effects can be extended to the spatial lag (including spatially lagged dependent variables) and the spatial error model (including spatial error autocorrelation terms). The dependent variable in the spatial lag model depends on the dependent variable observed in neighboring units . The spatial lag model with spatial and temporal fixed effects could be specified as follows :
The spatial error model assumed that the effect of spatial correlation is a representation of the ignored variable. The model with spatial and temporal fixed effects could be specified as below :
where W is an N × N positive non-stochastic spatial weight matrix; ρ is the spatial autoregressive coefficient and λ is usually called the spatial autocorrelation coefficient.
Likelihood ratio (LR) tests were used to determine the spatial and/or time-period fixed effects in the extension of the model. The Lagrange multiplier (LM) and robust Lagrange multiplier (robust LM) test were used to determine whether the spatial lag model and/or the spatial error model are more appropriate than a simple panel data model to describe the data. The spatial error model is more appropriate if the LM error is more significant than the LM lag and the robust LM error test is significant and the robust LM lag is not significant and vice versa. R2 and log-likelihood were the commonly used effective criteria to evaluate the model [29–31]. Matlab R2014a (Mathworks Inc., Natick, MA, USA) was used to perform the analysis of the spatial panel data models, LR, and LM tests.
A total of 7,061,525 HFMD cases were reported in 353 regions of mainland China from May 2008 to November 2012. The number of cases of HFMD per month ranged from 0 to 11,730 in each region. The number of reported cases for each month is shown in Fig. 2, which indicates a potential seasonality of the incidence of HFMD, since more cases occurred in April-July.
Spatial autocorrelation of HFMD incidence
The spatial autocorrelation test for each month is shown in Fig. 2. The figure demonstrates that a high global spatial autocorrelation of HFMD was detected at the regional level in mainland China within each epidemic month during May 2008 to November 2012 in group G1 (Moran’s I from 0.046 to 0.848) and group G2 (Moran’s I from 0.020 to 0.550). In group G3 and group G4, the global Moran’s I value was only significant in some months. Because of the existence of a high spatial dependency on the occurrence of HFMD, a spatial autocorrelation should be included in the regression model and thus, the panel data models were more appropriate.
Spatial panel data models
Logarithmic transformation of the reported incidence was used due to the distribution of the reported disease incidence having a non-normal distribution. Because of some zeroes in the dependent variable, we assumed that there were 0.5 cases in the zones with no reported cases before calculating the incidence and the logarithmic transformation. Table 1 shows the basic descriptive statistics for meteorological variables and dependent variables (logarithmic transformation of the reported incidence) in each group.
Table 2 shows the results to determine which specific fixed effect and which type of spatial dependency term should be included in the model for each group. The LR tests showed that the extension of the model with both spatial and time-period fixed effects was more suitable (P < 0.001). LM and robust LM test results demonstrated that the spatial lag models were more appropriate than the spatial error models. Overall, the test results implied that the spatial lagged model with spatial and time-period fixed effects was more appropriate to process the data.
The results of the three models (fixed effects model, spatial lag panel model and spatial error panel model) for each group are shown in Table 3. It can be seen in Table 3 that the spatial lag panel model and the spatial error panel model were better than the classic fixed effects model, and the spatial lag model is more appropriate than the spatial error model at comparing the values of R2 and the log-likelihood in each group. Both the spatial autoregressive coefficient (ρ = 0.452, 0.615, 0.241 and 0.201 for G1, G2, G3, and G4, respectively) and the spatial autocorrelation coefficient (λ = 0.450, 0.616, 0.178 and 0.092 for G1, G2, G3, and G4, respectively) were positive and statistically significant, suggesting the HFMD incidence of a spatial unit correlates positively with the incidence of surrounding spatial units and unmeasured variables.
Different significant factors were found in different groups. The average atmospheric pressure (AAP), average vapor pressure (AVP), monthly total sunshine hours (MSH) and mean temperature difference (MTD) had a significant correlation with the HFMD incidence in the three models for regions, and average wind speed (AWS) was also significant in the spatial lag model for south/east China. Average temperature (AT), average vapor pressure (AVP), average relative humidity (ARH), rain day (RD) and average temperature distance (ATD) were risk factors in the spatial lag model in regions with a temperate monsoon climate (TMC). Average temperature (AT) has a significant positive correlation with the HFMD incidence, but monthly total sunshine hours (MSH) was found to have a significant negative correlation in both G3 and G4. Monthly precipitation (MP) and average wind speed (AWS) also had a significant negative association with the HFMD incidence in G4, although average vapor pressure (AVP) was positively significant. The collected variables and the spatial neighborhood effects could jointly explain approximately 80.0 % of the variation of the HFMD incidence in the first 3 groups (G1 [R2 = 0.803], G2 [R2 = 0.802] and G3 [R2 = 0.797]) and 57.0 % in group G4 (R2 = 0.570). We can also find that most coefficients of spatial panel models were smaller than the fixed effect model, which indicated that if the spatial autocorrelation is ignored, the effects of meteorological factors would be overestimated.
The relationship between HFMD and meteorological factors was quantified using spatial panel data models based on longitudinal data from 353 regions of mainland China according to climate type from May 2008 to November 2012.
According to Elhorst , spatial panel data models contain more variation and less co-linearity among variables and are more informative than purely cross-sectional models or time-series models. The model tests justified the spatial and temporal fixed effects, which indicate that spatial and temporal heterogeneity do affect the robustness of statistical models, and taking the panel structure and spatial autocorrelation terms in the model simultaneously result in a better fit for the data . Previous studies have confirmed the existence of spatial correlations in the incidence of HFMD in various regions, which is consistent with our study [32–35]. Both the spatial autoregressive coefficient in the spatial lag panel model and the spatial autocorrelation coefficient in the spatial error panel model are positive and statistically significant. This further indicates the importance and necessity to include the spatial correlation in exploring the risk factors of HFMD.
These results show that every meteorological factor might have an association with HFMD incidence, and the effect may differ with different climates. A positive relationship between temperature and HFMD incidence (in group G2, G3 and G4) as well as between relative humidity (in group G2) and HFMD incidence was found in our study, which was consistent with the finding of other recent studies [8–10, 13–16]. Average wind speed has a negative association with HFMD incidence (in group G1 and G4) in this study, which is inconsistent with Ma’s, Li’s and Yong Huang’s studies [8, 13, 14]. A negative effect was observed for sunshine at lag days 3–4 in Wu’s study , which is consistent with the results in groups G3 and G4 in our research. Deng’s study  showed a positive correlation between sunshine and HFMD, which was also consistent with our result in group G1. The differences in the results may be caused by different climate types, which need further research. AAP was found to be positively significant with HFMD in group G1 in this study, which was different from Li’s study in Guangzhou . In Wang’s research , vapor pressure had a negative correlation with HFMD, which is consistent with our findings in G2. We also found opposite results in G1 and G4. We also found that monthly precipitation (MP) and mean temperature difference (MTD) were negatively correlated with HFMD in G4 and G1, respectively. Both rain day (RD) and average temperature distance (ATD) had a positive relationship with HFMD in G2, but not in G1, G3 and G4.
There are some limitations in this study. Although the spatial fixed effect and temporal fixed effect were included in the spatial panel models, the models did not consider the effect of space-time interaction, and we assumed that the effect of meteorological variables on HFMD was consistent in all regions. However, the population and climate distribution vary greatly in the world’s third biggest country. It would be difficult to keep the spatial stationarity assumption in such a complicated area. The geographical weighted regression model (GWR) or random coefficients models would be more useful to explore the local effect of factors of interest in the future. We selected the month as the temporal scale; more precise results or lag effects of meteorological factors may have been attained if detailed data were obtained and used for the analysis. Another limitation is that other factors, such as social/economic status, which may affect the HFMD, were not quantified in this study.
Spatial panel data models are useful and effective when longitudinal data are available and spatial autocorrelation exists. This study provides quantitative evidence that the incidence of HFMD could be affected by average atmospheric pressure, average temperature, average vapor pressure, average relative humidity, monthly precipitation, average wind speed, monthly total sunshine hours, mean temperature difference, rain day and average temperature distance. The effects of meteorological factors may be different in different climate types. This can facilitate a better understanding of epidemic trends and a preparedness for HFMD prevention and control strategies.
average atmospheric pressure
average relative humidity
average temperature distance
average vapor pressure
average wind speed
- China CDC:
China Center for Disease Control and Prevention
geographically weighted regression models
hand, foot and mouth disease
monthly total sunshine hours
mean temperature difference
plateau mountain climate
- Robust LM:
robust Lagrange multiplier
subtropical monsoon climate
temperate continental climate
temperate monsoon climate
Yang F, Jin Q, Zhang T, Hu Y, Wang S, Sun X, et al. Survey of enterovirus infections from hand, foot and mouth disease outbreak in China, 2009. Virol J. 2011;8:508.
World Health Organization. A guide to clinical management and public health response for Hand, Foot and Mouth disease (HFMD). Geneva: WHO Press; 2011.
Chinese Center for Disease Control and Prevention (China CDC). National incidence and death cases of notifiable class A or class B infectious disease (2008, 2009, 2010, 2011, 2012). (http://www.chinacdc.cn). Accessed 03 Apr 2014.
Xing W, Liao Q, Viboud C, Zhang J, Sun J, Wu JT, et al. Hand, foot, and mouth disease in China, 2008–12: an epidemiological study. Lancet Infect Dis. 2014;14(4):308–18.
Wang C, Li X, Zhang Y, Xu Q, Huang F, Cao K, et al. Spatiotemporal cluster patterns of hand, foot, and mouth disease at the county level in mainland China, 2008–2012. PLoS One. 2016;11(1):e0147532.
Onozuka D, Hashizume M. The influence of temperature and humidity on the incidence of hand, foot, and mouth disease in Japan. Sci Total Environ. 2011;410–411:119–25.
Ma E, Lam T, Chan KC, Wong C, Chuang SK. Changing epidemiology of hand, foot, and mouth disease in Hong Kong, 2001–2009. Jpn J Infect Dis. 2010;63:422–6.
Huang Y, Deng T, Yu S, Gu J, Huang C, Xiao G, et al. Effect of meteorological variables on the incidence of hand, foot, and mouth disease in children: a time-series analysis in Guangzhou, China. BMC Infect Dis. 2013;13(1):134.
Deng T, Huang Y, Yu S, Gu J, Huang C, Xiao G, et al. Spatial-temporal clusters and risk factors of hand, foot, and mouth disease at the district level in Guangdong Province, China. 2013. e56943.
Wang H, Du Z, Wang X, Liu Y, Yuan Z, Liu Y, et al. Detecting the association between meteorological factors and hand, foot, and mouth disease using spatial panel data models. Int J Infect Dis. 2015;34:66–70.
Hii YL, Rocklov J, Ng N. Short term effects of weather on hand, foot and mouth disease. PLoS One. 2011;6(2):e16796.
Urashima M, Shindo N, Okabe N. Seasonal models of herpangina and hand-foot-mouth disease to simulate annual fluctuations in urban warming in Tokyo. Jpn J Infect Dis. 2003;56:48–53.
Ma E, Lam T, Wong C, Chuang SK. Is hand, foot and mouth disease associated with meteorological parameters? Epidemiol Infect. 2010;138(12):1779–88.
Li T, Yang Z, Liu X, Kang Y, Wang M. Hand-foot-and-mouth disease epidemiological status and relationship with meteorological variables in Guangzhou, southern China, 2008–2012. Rev Inst Med Trop Sao Paulo. 2014;56(6):533–9.
Wu H, Wang H, Wang Q, Xin Q, Lin H. The effect of meteorological factors on adolescent hand, foot, and mouth disease and associated effect modifiers. Global health action, 2014, 7. doi: http://dx.doi.org/10.3402/gha.v7.24664.
Hu M, Li Z, Wang J, Jia L, Liao Y, Lai S, et al. Determinants of the incidence of hand, foot and mouth disease in China using geographically weighted regression models. PLoS One. 2012;7(6):e38978.
Elhorst JP. Spatial panel data models. Handbook of applied spatial analysis: 377–407; 2010.
Lee L, Yu J. Some recent developments in spatial panel data models. Reg Sci Urban Econ. 2010;40(5):255–71.
The Ministry of Health of the People's Republic of China. Guide for the preparedness and control measures of hand, foot, and mouth disease in China. (2008 version) (in Chinese). (http://www.gov.cn/gzdt/2008-05/03/content_960347.htm).
National Bureau of Statistics of the People’s Republic of China. http://www.stats.gov.cn/tjsj/tjgb/rkpcgb/dfrkpcgb/. Accessed 14 Mar 2014.
Pfeiffer DU, Robinson TP, Stevenson M, Stevens KB, Rogers DJ, Clements AC. Spatial analysis in epidemiology. New York: Oxford University Press; 2008.
Chan TC, King CC. Surveillance and epidemiology of infectious diseases using spatial and temporal clustering methods. In: Castillo-Chavez C, Lober WB, Thurmond M, Zeng D, editors. Infectious disease informatics and biosurveillance: research, systems and case studies. U.S.A.: Springer; 2010. p. 208–34.
Wooldridge JM. Econometric analysis of cross section and panel data. 2nd ed. Boston: The MIT press; 2010.
Hsiao C. Analysis of panel data. 2nd ed. London: Cambridge University Press; 2002.
LeSage J, Pace RK. Introduction to spatial econometrics: CRC press. 2010.
Hausman JA. Specification tests in econometrics. Econometrica. 1978;1251–1271.
Baltagi BH, Song SH, Jung BC, Koh W. Testing for serial correlation, spatial autocorrelation and random effects using panel data. J Econ. 2007;140(1):5–51.
Elhorst JP. Specification and estimation of spatial panel data models. Int Reg Sci Rev. 2003;26:244–68.
Lee L, Yu J. Estimation of spatial autoregressive panel data models with fixed effects. J Econ. 2010;154:165–85.
Elhorst JP. Matlab software for spatial panels. Int Reg Sci Rev. 2012;15905357.
Kelejian HH, Prucha IR. Specification and estimation of spatial autoregressive models with autoregressive and heteroskedastic disturbances. J Econ. 2010;157:53–67.
Liu Y, Wang X, Liu Y, Sun D, Ding S, Zhang B, et al. Detecting spatial-temporal clusters of HFMD from 2007 to 2011 in Shandong Province, China. PLoS One. 2013;8(5):e63447.
Wang J, Cao Z, Zeng DD, Wang Q, Wang X, Qian H. Epidemiological analysis, detection, and comparison of space-time patterns of Beijing hand-foot-mouth disease (2008–2012). PLoS One. 2014;9(3):e92745.
Xie YH, Chongsuvivatwong V, Tang Z, McNeil EB, Tan Y. Spatio-temporal clustering of hand, foot, and mouth disease at the county level in Guangxi, China. PLoS One. 2014;9(2):e88065.
Xiao GX, Hu YH, Ma JQ, Hao YT, Wang XF, Zhang E. Spatial clustering and changing trend of hand-foot-mouth disease during 2008–2011 in China. Zhonghua liu xing bing xue za zhi. 2012;33:808–12.
We thank the local health departments, hospitals, and CDC China for assistance in coordinating the data. The views expressed herein are those of the authors and do not necessarily reflect Chinese CDC policies.
This work was supported by The National S&T Major Project (2012ZX10005009-003 and 2014ZX10004005-001), National 973 Project (2011CB505404), and National Natural Science Foundation of China (81302516). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The dataset(s) supporting the conclusions of this article is(are) available in the Data-center of China Public Health Science repository (hyperlink to dataset(s) in http://www.phsciencedata.cn/Share/en/index.jsp), the China Meteorological Data Sharing Service System (hyperlink to dataset(s) in http://cdc.nmic.cn/home.do) and the National Science & Technology Infrastructure of China, Data Sharing Infrastructure of Earth System Science (hyperlink to dataset(s) in http://www.geodata.cn), which are publicly accessible.
CW participated in the design of the study, collected the data, performed the statistical analysis and drafted the manuscript. KC participated in the design of the study, collected the disease data, participated in the statistical analysis and helped to draft the manuscript. QX, FH, LT, JG and QG collected meteorological data, contributed data analysis and helped to draft the manuscript. YZ, LF and XL participated in the design of the study, collected disease data and the base map of China, and helped to draft the manuscript. XG conceived of the study and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
About this article
Cite this article
Wang, C., Cao, K., Zhang, Y. et al. Different effects of meteorological factors on hand, foot and mouth disease in various climates: a spatial panel data model analysis. BMC Infect Dis 16, 233 (2016). https://doi.org/10.1186/s12879-016-1560-9
- Hand, foot and mouth disease
- Spatial panel data model
- Meteorological factors
- Climate type