Time Series Analyses of Hand, Foot and Mouth Disease Integrating Weather Variables

• Published in: PloS one Vol. 10; no. 3; p. e0117296
• Main Authors: Song, Yuanbin, Wang, Fan, Wang, Bin, Tao, Shaohua, Zhang, Huiping, Liu, Sai, Ramirez, Oscar, Zeng, Qiyi
• Format: Journal Article
• Language: English
• Published: United States Public Library of Science 02.03.2015; Public Library of Science (PLoS)
• Subjects: Adolescent; Child; Child, Preschool; China - epidemiology; Climate; Climate change; Climate models; Climatic analysis; Correlation; Correlation analysis; Coxsackievirus infections; Data processing; Disease control; Early warning systems; Enterovirus; Enterovirus A, Human - genetics; Enterovirus A, Human - isolation & purification; Epidemiology; Female; Foot and mouth disease; Hand, Foot and Mouth Disease - epidemiology; Hand, Foot and Mouth Disease - virology; Hand-foot-and-mouth disease; Health aspects; Hospitalization; Hospitals; Humans; Humidity; Impact prediction; Incidence; Independent variables; Infant; Infections; Laboratories; Male; Mathematical models; Models, Theoretical; Outbreaks; Pediatrics; Performance prediction; Precipitation (Meteorology); Prediction models; Rain; Rainfall; Repressors; Reverse Transcriptase Polymerase Chain Reaction; RNA, Viral - analysis; Serotypes; Serotyping; Time series; Variables; Viral infections; Weather; China; United States > US
• DOI: 10.1371/journal.pone.0117296
• ISSN: 1932-6203

The past decade witnessed an increment in the incidence of hand foot mouth disease (HFMD) in the Pacific Asian region; specifically, in Guangzhou China. This emphasized the requirement of an early warning system designed to allow the medical community to better prepare for outbreaks and thus minimize the number of fatalities. Samples from 1,556 inpatients (hospitalized) and 11,004 outpatients (non-admitted) diagnosed with HFMD were collected in this study from January 2009 to October 2013. Seasonal Autoregressive Integrated Moving Average (SARIMA) model was applied to establish high predictive model for inpatients and outpatient as well as three viral serotypes (EV71, Pan-EV and CA16). To integrate climate variables in the data analyses, data from eight climate variables were simultaneously obtained during this period. Significant climate variable identified by correlation analyses was executed to improve time series modeling as external repressors. Among inpatients with HFMD, 248 (15.9%) were affected by EV71, 137 (8.8%) were affected by Pan-EV+, and 436 (28.0%) were affected by CA16. Optimal Univariate SARIMA model was identified: (2,0,3)(1,0,0)52 for inpatients, (0,1,0)(0,0,2)52 for outpatients as well as three serotypes (EV71, (1,0,1)(0,0,1)52; CA16, (1,0,1)(0,0,0)52; Pan-EV, (1,0,1)(0,0,0)52). Using climate as our independent variable, precipitation (PP) was first identified to be associated with inpatients (r = 0.211, P = 0.001), CA16-serotype (r = 0.171, P = 0.007) and outpatients (r = 0.214, P = 0.01) in partial correlation analyses, and was then shown a significant lag in cross-autocorrelation analyses. However, inclusion of PP [lag -3 week] as external repressor showed a moderate impact on the predictive performance of the SARIMA model described here-in. Climate patterns and HFMD incidences have been shown to be strongly correlated. The SARIMA model developed here can be a helpful tool in developing an early warning system for HFMD.