This article has Open Peer Review reports available.
Comparisons between mild and severe cases of hand, foot and mouth disease in temporal trends: a comparative time series study from mainland China
- Xiong Xiao†1, 2View ORCID ID profile,
- Qiaohong Liao†3,
- Michael G. Kenward2,
- Yaming Zheng3,
- Jiao Huang3,
- Fei Yin1,
- Hongjie Yu1, 3Email author and
- Xiaosong Li1Email author
© The Author(s). 2016
Received: 11 August 2016
Accepted: 11 October 2016
Published: 21 October 2016
Over recent decades, hand, foot and mouth disease (HFMD) has emerged as a serious public health threat in the Asia-Pacific region because of its high rates of severe complications. Understanding the differences and similarities between mild and severe cases can be helpful in the control of HFMD. In this study, we compared the two types of HFMD cases in their temporal trends.
We retrieved the daily series of disease counts of mild and severe HFMD cases reported in mainland China in the period of 2009–2014. We applied a quasi-Poisson regression model to decompose each series into the long-term linear trend, periodic variations, and short-term fluctuations, and then we compared each component between two series separately.
A total of 11,101,860 clinical HFMD cases together with 115,596 severe cases were included into this analysis. We found a biennial increase of 24.46 % (95 % CI: 22.80–26.14 %) for the baseline of disease incidence of mild cases, whereas a biennial decrease of 8.80 % (95 % CI: 7.26–10.31 %) was seen for that of severe cases. The periodic variations of both two series could be characterized by a mixture of biennial, annual, semi-annual and eight-monthly cycles. However, compared to the mild cases, we found the severe cases vary more widely for the biennial and annual cycle, and started its annual epidemic earlier. We also found the short-term fluctuations between two series were still significantly correlated at the current day with a correlation coefficient of 0.46 (95 % CI: 0.43–0.49).
We found some noticeable differences and also similarities between the daily series of mild and severe HFMD cases at different time scales. Our findings can help us to deepen the understanding of the transmission of different types of HFMD cases, and also provide evidences for the planning of the associated disease control strategies.
KeywordsEnterovirus infection Hand, foot and mouth disease Case surveillance Comparative time series analysis Mainland China
Hand, foot and mouth disease (HFMD) is a common childhood infectious disease caused by viruses that belong to the Enterovirus group, mainly by Coxsackievirus A16 (CVA 16), Enterovirus 71 (EVA 71) and, more recently, Coxsackievirus A6 (CVA 6) . It is typically characterized by a febrile illness followed by rash on hands and feet, sometimes also with vesicular/ulcer in the mouth . In most cases, the disease is mild and self-limiting and requires no more treatment other than symptomatic relief . However, some patients, especially those infected with EVA 71, may develop severe complications involving the central nervous system (CNS) such as meningitis, encephalitis and acute flaccid paralysis [4, 5]. A few patients can even progress to fatal cardiopulmonary failure.
Following the first observation of the clinical syndrome of HFMD in 1957 in New Zealand , HFMD and associated severe CNS disease have been reported globally [7–13]. The outbreaks of HFMD reported before the mid-1990s are normally considered as “benign” with rare severe and fatal cases, or many severe cases have been initially misdiagnosed as some other diseases (like poliomyelitis or encephalitis) . However, over the last few decades, the Asia-Pacific region has experienced a series of large epidemics of HFMD accompanied by abnormal high rates of severe and fatal cases [14–20]. Since then HFMD has started to prompt huge public health concerns in the Asia-Pacific region given its threat to young lives and the potential for its emergence as a leading cause of enterovirus-related CNS disease after poliomyelitis .
Given the high prevalence and benign outcome of this common childhood illness, a crucial problem in HFMD prevention and control is how to characterize the differences between mild and severe HFMD cases, so that we can optimize the associated clinical management and public health response. Previous studies have focused principally on their differences in virology, demography, and clinical manifestations . Most of those studies attempt to distinguish severe from mild HFMD cases in early stage which can provide implications for the related clinical management. To the best of our knowledge, few studies have focused on their differences in temporal trends. The temporal trends have played a crucial role in the scientific research of infectious disease, which can be used to inform the mechanism of disease transmission, identify time-varying risk factors, predict disease outbreak, etc. Therefore, obtaining an understanding of the differences and similarities between mild and severe HFMD cases in temporal trends may provide some implications for the HFMD-related public health response. In this study, we conducted a comparative time series analysis based on the surveillance data in mainland China from 2009 to 2014 aiming to compare mild and severe HFMD cases in their temporal trends.
Data source and case definitions
We retrieved the case surveillance data of HFMD from the Chinese Centre for Disease Control and Prevention (China CDC). HFMD was made statutorily notifiable as a Class C infectious disease on 2 May 2008 in mainland China. All clinical cases of HFMD are reported to China CDC via the online reporting system within 24 h of diagnosis by a standardized form.
A clinical case of HFMD is defined as a patient with papular or vesicular rash on hands, feet, mouth or buttocks, with or without fever. A severe case is defined as a clinical case with any CNS complications (including aseptic meningitis, encephalitis, encephalomyelitis, acute flaccid paralysis, or autonomic nervous system dysregulation), or cardiopulmonary dysfunctions (pulmonary edema, pulmonary hemorrhage, or cardiorespiratory failure), or both. Otherwise, the patient is classified as a mild case . More details about case definitions can be found in the Chinese Guidelines for HFMD public health response .
We aggregated the daily counts of mild and severe HFMD cases reported between 1 January 2009 and 31 December 2014 based on the date of onset of symptom, respectively. Previous analysis has shown that data collected during the first year are less reliable than those from more recent years, mainly because of improvements in reporting and surveillance . Therefore, we excluded the data collected in 2008 and started our study period from 1 January 2009. Given that the proportion of severe cases is much less than that of mild cases, the huge difference in magnitude (or amplitude) will inevitably tend to mask their relationships in temporal trends. To address above issue, we rescaled both two series by scaling by the arithmetic average for each series separately. The rescaling makes the range of the two series comparable but preserve the temporal trends. We can interpret the values of rescaled time series as the relative daily increase or decrease to the overall mean.
Comparing the long-term linear trend
This figure can be interpreted as the accumulated relative increase in the baseline of the disease incidence over a biennial cycle after adjusting for the periodic variations. Then, we compared the biennial increase between two series and tested their equality (i.e. whether or not the difference between β il equals zero) by a Z test  suggested for comparing the slopes between two regressions.
Comparing the periodic variations
β ijs and β ijc are the estimated coefficients of sine and cosine functions for the cycle j, γ ij and φ ij are the estimated amplitude and phase time for the cycle j. For the overall periodic curve, we are unable to derive the parametric forms of its amplitude or phase time because it is a mixture of four types of cycles. Instead, we calculated some of its key characteristics by definitions based on the fitted periodic curve, including the start timing (the temporal position of the minimum value), the peak value (the difference between the peak and minimum value) and the peak timing (the temporal position of the peak value). It could be very cumbersome to derive the parametric distributions of the characteristics that we compared. Therefore, we made the statistical inferences based on the Monte Carol simulations . In simplest terms, we re-generated 1000 replicates of the four types of cycles as well as the overall fitted periodic curves by sampling from the model estimations of the coefficients of Fourier terms (i.e. [β ijs, β ijc ]) assuming a multivariate normal distribution. And then, we obtained the confidence intervals and implemented the statistical tests that we were interested in based on the distributions of the sampled replicates.
Comparing the short-term fluctuations
Unlike the long-term linear trend and the periodic variations, the short-term fluctuations cannot be expressed as some comparable figures. To compare the short-term fluctuations, we implemented the cross-correlation analysis  between two series to see whether they were still temporally related and to explore if there was any potential delay in their temporal trend. Also, we implemented a partial cross-correlation analysis  to adjust for the holiday effect (including the summer and winter vacations for school and national public holidays), day of week and the autocorrelations. The selection of the autoregressive terms was informed by the autocorrelation plot for each series.
We did all statistical analysis with R software (version 3.2.3) using the packages stats and TSA.
The comparisons of time series components between the series of mild and severe HFMD cases
Mild HFMD cases
Severe HFMD cases
Long-term linear trend
Biennial increase (%)
124.46 (122.80, 126.14)
91.20 (89.69, 92.74)
peak time (days)
peak time (days)
peak time (days)
peak time (days)
First year cycle of the overall periodic curve
Start time (days)
Major peak time (days)
Major peak value
Minor peak time (days)
Minor peak value
Second year cycle of the overall periodic curve
Start time (days)
Major peak time (days)
Major peak value
Minor peak time (days)
Minor peak value
As mentioned earlier, the periodic variations of both two series consist of the same four types of cycles (Fig 2e–h). For both two series, the annual cycle made the largest contributions (measured with the size of peak value) to the periodic variations, followed by the semi-annual cycle, biennial cycle and 8-monthy cycle (Table 1). The annual peak value for mild and severe cases were 1.78 (95 % CI: 1.73–1.84) and 3.08 (95 % CI: 2.97–3.19), respectively. We estimated the overall periodic curves by combining the above four types of cycles collectively (Fig 2c). For both two series, the overall periodic curves cycled every 2 years with a high-low pattern (i.e. a relatively low-level epidemic in the first year followed by a high-level epidemic in the next year, denote as first and second year cycle in the following), which was mainly driven by the biennial cycle. Besides, we found two peaks in the annual epidemics of mild cases, including a major peak in early summer and a minor peak in autumn, which were mainly driven by the annual and semi-annual cycle. Whereas, only the major peak was seen in the series of severe cases with a peak value of 2.10 (95 % CI: 2.00–2.19) and 4.20 (95 % CI: 4.07–4.34) for the first and second year cycle respectively, which was nearly twice that of mild cases (Table 1). The discrepancy could be explained by the relatively great peak size of annual and biennial cycle in the series of severe cases, which can conceal the peak of semi-annual cycle. With regards to the peak timing, we found the series of severe cases started the annual epidemic earlier than the series of mild cases with a leading time of 9 (95 % CI: 7–11) and 14 (95 % CI: 11–16) days for the first and second year cycle respectively. However, they reached their major annual peak nearly at the same time (Table 1).
In this study, we have characterized the temporal trends of HFMD for both the mild and severe HFMD cases. We have found sustained epidemics of HFMD in our study period of 2009–2014. In comparison with reports from other countries or districts in the Asia-Pacific region, the epidemics of HFMD in mainland China exhibited a more regular temporal pattern. In Sarawak (Malaysia)  and Japan , epidemics of HFMD occurred every 3 to 4 years, while in Taiwan , Hong Kong  and Singapore  the epidemic pattern was more irregular with the inter-epidemic interval varying from 0 to 3 years. The relatively constant growth of population in China  could partly explain the regular behaviour of the epidemics in mainland China. Besides, the related vaccines were available until recently after two inactivated monovalent EV-A71 vaccines had been licensed in China in November 2015 . The lack of specific treatments or vaccines could be another reason for the sustained epidemics of HFMD.
Normally, we would expect the baseline of disease incidence (i.e. long-term trend) for a specific infection decreased over time because of the accumulation of immune individuals, just as we saw in the series of severe cases. However, an upward trend was observed in the series of mild cases. Unlike the severe cases which are overwhelmingly caused by EV-A71 , the mild symptoms of HFMD can be caused by a wide range of enterovirus. Even though there is evidence of cross-protection between different serotypes of enterovirus, the cross-protection is believed to be temporal . Therefore, the replacement of circulating serotypes of Enterovirus in different years might be one reason for the continued growth of mild cases. Besides, it is worth noting that the number of mild cases is more likely biased by the under reporting compared to that of severe cases because its symptom is benign. According to a retrospective study in Jiangsu province (China) in 2009, only 11.2 % of HFMD cases were reported to the National Infectious Disease Information Management System . Therefore, we cannot rule out the possibility that the upward trend in the series of mild cases may be due just to an improvement in reporting.
By comparing the periodic variations between two series, most of the severe cases were reported around early summer with only one dominant peak every year, whereas the mild cases were more evenly distributed all year around with two annual peaks. Our findings imply that the specific condition in early summer may be more favourable to the transmission of severe cases than any other time of year. The more dominant major peak of severe cases in early summer can lead to a stronger depletion of the pool of susceptible individuals than that of mild cases, which can then partly explain the absence of the minor peak in the following autumn. Consequently, because there is no minor peak in the autumn, we expect a more sufficient replenishment of susceptible individuals for the severe cases in the second half of the year. This is probably why we see an even stronger and earlier major peak in the next year in the series of severe cases. So far, the exact causes of HFMD related severe complications are still unclear . Apart from the infection of EV-A71, previous epidemiological studies suggest that some host factors may also increase the risk of severe complication, such as male gender, young age, etc. . Therefore, the specific condition in early summer may have some impacts on the virulence of EV-A71 or the immunity of hosts. However, further studies are needed to verify our hypothesis.
Even after removing the long-term linear trend and periodic variations, the short-term fluctuations of mild and severe cases have been found to be correlated with each other, especially at the current day. As mentioned earlier, the long-term linear trend and periodic variation of HFMD are mainly driven by some long-term and stable factors, such as the seasonal change of climate and immunity level, the growth of population, the improvement of reporting, etc. However, for the short-term fluctuations of HFMD incidences, it is more likely they are driven by some transient risk factors, for instance, the daily change of temperature and humidity [40, 41]. Our findings suggest that such kind of transient risk factors may have similar impacts on both the mild and severe HFMD cases.
Few previous studies have compared the differences in temporal trends between mild and severe HFMD cases. That is probably because it is very hard to obtain the temporal trends of severe cases given the relatively rare occurrence of severe complications, especially for a small area. We only found one study in Taiwan  which examined the temporal associations between the weekly series of mild and severe HFMD cases. It suggested a significantly correlation (r = 0.553) between two types of cases occurring in the same week. Their results are consistent with our findings in the short-term fluctuations. However, the Taiwan study only applied a simple correlation analysis between two series without distinguishing the short-term and long-term temporal trends. As mentioned earlier, the temporal trends at different time scale can be driven by different factors. Therefore, we believe our study can provide a more explicit understanding of HFMD temporal trends.
In this analysis, we did not implement the conventional methods of time series decompositions which was mainly based on the nonparametric smooth function . Instead, we decomposed the time series by the application of a quasi-Poisson regression model. We did this for two main reasons. First, our study interest is different from the conventional purpose of time series decompositions which aim to prediction. In our study, we are more interested in the comparison of temporal trends between two series. By using the regression model, we can extract the parametric forms of temporal trends which allow us to make statistical inferences between two series. Second, the regression form model is more widely used in the field of epidemiological research. Therefore, the methods we used can make our study to be better understood by the potential readers.
There are several limitations in this study. First, we did not conduct any subgroup analyse either by regions or subgroups. As some literature has shown [23, 43–45], the temporal trends of HFMD can vary in different places and subpopulations. Therefore, our conclusions may be not applicable in some specific settings. However, subgroup analyses do raise issues of data sparsity with such a rare occurrence of severe cases. Second, we are unable to analyse the temporal trends of laboratory data in detail to further help with the interpretation of our findings. The laboratory data of HFMD in China were collected on a monthly basis, and so could not provide information at the short-term time scale. Also, the collection of laboratory specimens was not based on the random sampling, and the proportion of laboratory-confirmed cases was very small, hence severely limiting its values in our study . Third, 6 years is still a relatively short period for the time series analysis. As the accumulation of case surveillance data, a longer series would be likely to provide a better and more robust understanding of the temporal trends of HFMD.
In conclusion, we found some noticeable differences between the daily series of mild and severe HFMD cases in their long-term linear trend and periodic variations, but we also found the short-term fluctuations of these two series were still significant correlated. The different temporal trends between mild and severe HFMD cases can provide us some insights into their potentially different mechanisms of transmission. Our findings can be used to generate hypotheses for future studies to gain a better understanding of the transmission of HFMD, in particular for different types of clinical cases. Besides, the differences between mild and severe cases in temporal trends should also be considered when planning the disease control strategies to optimize the cost-effective. For instance, the more intense major peak in early summer in the series of severe cases may require additional attention when preparing to response the HFMD epidemic.
The authors thank the Chinese Centre for Disease Control and Prevention for providing data.
This work was supported by the National Science Fund for Distinguished Young Scholars (No. 81525023) (HY), the National Natural Science Foundation of China (No. 81473031) (HY), the National Natural Science Foundation of China (No. 81402766) (FY). Xiong Xiao was supported by the LSHTM-CMB scholarship.
Availability of data and materials
The data that support the findings of this study are available from the Chinese Centre for Disease Control and Prevention but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available.
HY and XL conceptualized and designed this experiment. XX and MK implemented the data analysis. XX, QL and MK interpreted the results and drafted the manuscript. YZ, JH and FY retrieved the dataset and did the preliminary analysis. MK, HY and XL revised 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
Hand, foot and mouth disease belongs to the class C notifiable disease in China. The National Health and Family Planning Commission of China decided that the collection of individual data for all notifiable diseases, according to the national surveillance protocol was part of an ongoing public health response and was thus exempt from institutional review board assessment.
The China CDC granted permission for the authors to use the data. The China CDC has deleted the personal identifiers before the co-authors of this article were given access to the surveillance data for the purposes of research. Data were also aggregated by day before analysis.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Stewart CL, Chu EY, Introcaso CE, Schaffer A, James WD. Coxsackievirus A6–induced hand-foot-mouth disease. JAMA Derm. 2013;149(12):1419–21.View ArticleGoogle Scholar
- World Health Organization. A guide to clinical management and public health response for hand, foot and mouth disease (HFMD). Manila: World Health Organization; 2011. http://www.wpro.who.int/emerging_diseases/documents/HFMDGuidance/en/index.html.Google Scholar
- Frydenberg A, Starr M. Hand, foot and mouth disease. Aust Fam Physician. 2003;32(8):594–5.PubMedGoogle Scholar
- Sabanathan S, Tan le V, Thwaites L, Wills B, Qui PT, Rogier van Doorn H. Enterovirus 71 related severe hand, foot and mouth disease outbreaks in South-East Asia: current situation and ongoing challenges. J Epidemiol Community Health. 2014;68(6):500–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Wong S, Yip C, Lau S, Yuen K. Human enterovirus 71 and hand, foot and mouth disease. Epidemiol Infect. 2010;138(08):1071–89.View ArticlePubMedGoogle Scholar
- Duff M. Hand-foot-and-mouth syndrome in humans: coxackie A10 infections in New Zealand. Br Med J. 1968;2(5606):661–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Blomberg J, Lycke E, Ahlfors K, Johnsson T, Wolontis S, Von Zeipel G. New enterovirus type associated with epidemic of aseptic meningitis and/or hand, foot, and mouth disease. Lancet. 1974;304(7872):112.View ArticleGoogle Scholar
- Kennett M, Birch C, Lewis FA, Yung A, Locarnini S, Gust I. Enterovirus type 71 infection in Melbourne. Bull World Health Organ. 1974;51(6):609.PubMedPubMed CentralGoogle Scholar
- Deibel R, Gross LL, Collins DN. Isolation of a new enterovirus. Exp Biol Med. 1975;148(1):203–7.View ArticleGoogle Scholar
- Shindarov L, Chumakov M, Voroshilova M, Bojinov S, Vasilenko S, Iordanov I, Kirov I, Kamenov E, Leshchinskaya E, Mitov G. Epidemiological, clinical, and pathomorphological characteristics of epidemic poliomyelitis-like disease caused by enterovirus 71. J Hyg Epidemiol Microbiol Immunol. 1978;23(3):284–95.Google Scholar
- Tagaya I, Takayama R, Hagiwara A. A large-scale epidemic of hand, foot and mouth disease associated with enterovirus 71 infection in Japan in 1978. Jpn J Med Sci Biol. 1981;34(3):191–6.View ArticlePubMedGoogle Scholar
- Nagy G, Takátsy S, Kukán E, Mihály I, Dömök I. Virological diagnosis of enterovirus type 71 infections: experiences gained during an epidemic of acute CNS diseases in Hungary in 1978. Arch Virol. 1982;71(3):217–27.View ArticlePubMedGoogle Scholar
- Samuda G, Chang W, Yeung C, Tang P. Monoplegia caused by Enterovirus 71: an outbreak in Hong Kong. Pediatr Infect Dis J. 1987;6(2):206.View ArticlePubMedGoogle Scholar
- Ho M, Chen E-R, Hsu K-H, Twu S-J, Chen K-T, Tsai S-F, Wang J-R, Shih S-R. An epidemic of enterovirus 71 infection in Taiwan. N Engl J Med. 1999;341(13):929–35.View ArticlePubMedGoogle Scholar
- Chan L, Parashar UD, Lye M, Ong F, Zaki SR, Alexander JP, Ho K, Han LL, Pallansch MA, Suleiman AB. Deaths of children during an outbreak of hand, foot, and mouth disease in Sarawak, Malaysia: clinical and pathological characteristics of the disease. Clin Infect Dis. 2000;31(3):678–83.View ArticlePubMedGoogle Scholar
- Chan KP, Goh KT, Chong CY, Teo ES, Lau G, Ling AE. Epidemic hand, foot and mouth disease caused by human enterovirus 71, Singapore. Emerg Infect Dis. 2003;9(1):78–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Van Tu P, Thao NTT, Perera D, Truong KH, Tien NTK, Thuong TC, How OM, Cardosa MJ, McMinn PC. Epidemiologic and virologic investigation of hand, foot, and mouth disease, southern Vietnam, 2005. Emerg Infect Dis. 2007;13(11):1733–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Ma E, Chan KC, Cheng P, Wong C, Chuang SK. The enterovirus 71 epidemic in 2008—public health implications for Hong Kong. Int J Infect Dis. 2010;14(9):e775–80.View ArticlePubMedGoogle Scholar
- Zhang Y, Zhu Z, Yang W, Ren J, Tan X, Wang Y, Mao N, Xu S, Zhu S, Cui A. An emerging recombinant human enterovirus 71 responsible for the 2008 outbreak of hand foot and mouth disease in Fuyang city of China. Virol J. 2010;7(1):94.View ArticlePubMedPubMed CentralGoogle Scholar
- Biswas T. Enterovirus 71 causes hand, foot and mouth disease outbreak in Cambodia. Natl Med J India. 2011;25(5):316.Google Scholar
- Regional Emerging Diseases Intervention Center. Forum on hand foot and mouth disease (HFMD) in Asia-Pacific region: epidemiological, laboratory, clinical, and public health aspects. Singapore: Regional Emerging Diseases Intervention Center; 2009. http://www.aseanplus3fetn.net/doc/hfmd_singapore_redi.pdf.Google Scholar
- Fang Y, Wang S, Zhang L, Guo Z, Huang Z, Tu C, Zhu B-P. Risk factors of severe hand, foot and mouth disease: a meta-analysis. Scand J Infect Dis. 2014;46(7):515–22.View ArticlePubMedGoogle Scholar
- Xing W, Liao Q, Viboud C, Zhang J, Sun J, Wu JT, Chang Z, Liu F, Fang VJ, Zheng Y. Hand, foot, and mouth disease in China, 2008–12: an epidemiological study. Lancet Infect Dis. 2014;14(4):308–18.View ArticlePubMedPubMed CentralGoogle Scholar
- Ministry of Health of China. Guideline for HFMD public health response. Beijing: Ministry of Health of China; 2009.Google Scholar
- Dagum EB. Time series modelling and decomposition. Statistica. 2010;70(4):433–57.Google Scholar
- Shumway RH, Stoffer DS. Time series analysis and its applications: with R examples. 3rd ed. New York: Springer Science & Business Media; 2010.Google Scholar
- Paternoster R, Brame R, Mazerolle P, Piquero A. Using the correct statistical test for the equality of regression coefficients. Criminology. 1998;36(4):859–66.View ArticleGoogle Scholar
- Naumova E, Jagai J, Matyas B, DeMaria A, MacNeill I, Griffiths J. Seasonality in six enterically transmitted diseases and ambient temperature. Epidemiol Infect. 2007;135(02):281–92.View ArticlePubMedGoogle Scholar
- Buckland ST. Monte Carlo confidence intervals. Biometrics. 1984;40(3):811–817.Google Scholar
- Chatfield C. The analysis of time series: an introduction. 6th ed. London: Chapman & Hall /CRC; 2013.Google Scholar
- Stark E, Drori R, Abeles M. Partial cross-correlation analysis resolves ambiguity in the encoding of multiple movement features. J Neurophysiol. 2006;95(3):1966–75.View ArticlePubMedGoogle Scholar
- Podin Y, Gias EL, Ong F, Leong Y-W, Yee S-F, Yusof MA, Perera D, Teo B, Wee T-Y, Yao S-C. Sentinel surveillance for human enterovirus 71 in Sarawak, Malaysia: lessons from the first 7 years. BMC Public Health. 2006;6(1):180.View ArticlePubMedPubMed CentralGoogle Scholar
- Infectious Agents Surveillance Report. Hand, foot and mouth disease in Japan, 2002–11. IASR. 2012;33:55–6.Google Scholar
- Chen S-C, Chang H-L, Yan T-R, Cheng Y-T, Chen K-T. An eight-year study of epidemiologic features of enterovirus 71 infection in Taiwan. AmJTrop Med Hyg. 2007;77(1):188–91.Google Scholar
- Ang LW, Koh B, Chan KP, Chua LT, James L, Goh KT. Epidemiology and control of hand, foot and mouth disease in Singapore. Ann Acad Med Singapore. 2009;38:106–12.PubMedGoogle Scholar
- The World Bank. Birth rate, crude (per 1,000 people) in China. 2015. http://data.worldbank.org/indicator/SP.DYN.CBRT.IN?page=2. Accessed 21 Apr 2016.
- China Food and Drug Administration. Announcement on licensed drugs approved by China Food and Drug Administration (No.4 in 2016). http://www.sfda.gov.cn/WS01/CL0087/142000.html. Accessed 01 May 2016.
- Takahashi S, Liao Q, Van Boeckel TP, Xing W, Sun J, Hsiao VY, Metcalf CJE, Chang Z, Liu F, Zhang J. Hand, foot, and mouth disease in China: modeling epidemic dynamics of enterovirus serotypes and implications for vaccination. PLoS Med. 2016;13(2):e1001958.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhu F-C, Liang Z-L, Meng F-Y, Zeng Y, Mao Q-Y, Chu K, Song X-F, Yao X, Li J-X, Ji H. Retrospective study of the incidence of HFMD and seroepidemiology of antibodies against EV71 and CoxA16 in prenatal women and their infants. PLoS One. 2012;7(5):e37206.View ArticlePubMedPubMed CentralGoogle Scholar
- Hii YL, Rocklov J, Ng N. Short term effects of weather on hand, foot and mouth disease. PLoS One. 2011;6(2):e16796.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang Y, Deng T, Yu S, Gu J, Huang C, Xiao G, Hao Y. 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:134.View ArticlePubMedPubMed CentralGoogle Scholar
- Chan T-C, Hwang J-S, Chen R-H, King C-C, Chiang P-H. Spatio-temporal analysis on enterovirus cases through integrated surveillance in Taiwan. BMC Public Health. 2014;14(1):1.View ArticleGoogle Scholar
- Lee C-CD, Tang J-H, Hwang J-S, Shigematsu M, Chan T-C. Effect of Meteorological and Geographical Factors on the Epidemics of Hand, Foot, and Mouth Disease in Island-Type Territory, East Asia. Biomed. Res. Int. 2014; doi:10.1155/2015/805039.
- 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. Glob Health Action. 2014;7:24664.View ArticlePubMedGoogle Scholar
- Zhu L, Yuan Z, Wang X, Li J, Wang L, Liu Y, Xue F, Liu Y. The impact of ambient temperature on childhood HFMD incidence in inland and coastal area: a Two-city study in Shandong province, China. Int J Environ Res Public Health. 2015;12(8):8691–704.View ArticlePubMedPubMed CentralGoogle Scholar