Research Article  Open Access
A HandFootandMouth Disease Model with Periodic Transmission Rate in Wenzhou, China
Abstract
We establish an SEIQRS epidemic model with periodic transmission rate to investigate the spread of seasonal HFMD in Wenzhou. The value of this study lies in two aspects. Mathematically, we show that the global dynamics of the HFMD model can be governed by its reproduction number ; if , the diseasefree equilibrium of the model is globally asymptotically stable, which means that the disease will vanish after some period of time; while if , the model has at least one positive periodic solution and is uniformly persistent, which indicates that HFMD becomes an endemic disease. Epidemiologically, based on the statistical data of HFMD in Wenzhou, we find that the HFMD becomes an endemic disease and will break out in Wenzhou. One of the most interesting findings is that, for controlling the HFMD spread, we must increase the quarantined rate or decrease the treatment cycle.
1. Introduction
Hand, foot, and mouth disease (HFMD) was first reported in New Zealand in 1957 and is caused by Coxsackievirus A16 (CVA16) and human enterovirus 71 (HEV71) and occasionally by Coxsackie virus A4A7, A9, A10, B1B3, and B5. HFMD is an acute viral illness that usually affects infants and children younger than 5 years old, and, however, it can sometimes occur in adults. HFMD is characterized by fever, intraoral vesicles and erosions, and papulovesicles that favor the palms and soles. HFMD is spread from person to person through nose and throat secretions (such as saliva, sputum, or nasal mucus), blister fluid, or stool of infected persons [1, 2].
HFMD is moderately contagious and usually not a serious illness among the infected population; however, recent outbreak of HFMD in countries such as China, Singapore, and Finland has brought the world attention to HFMD due to complications of death related cases [3ā10]. In 2008, there were 488,955 cases reported, with a morbidity of 37/100,000, mortality of 0.0095/100,000, and illdeath rate of 0.26/1000; while, in 2009, there were 1,155,525 cases reported [11]. Due to the severity of this disease, the Ministry of Health of the Peopleās Republic of China upgraded HFMD to a Class C communicable disease on May 2, 2008.
It is obvious that HFMD not only causes health problems but also has great social and economic impacts which are not easily quantifiable. So it is important to understand the spread dynamics of HFMD among the susceptible populations and to enable policy makers to take effective measures to curb the disease spread and reduce the adverse impact of the disease [12, 13].
There are several types of analytical models that are valuable to understand and predict transmission of HFMD. One is the statistical models which can help us find novel information concerning pathogen detection and some probable coinfection factors in HFMD and have been applied to understand HFMDās spatiotemporal transmission and discover the relationship between HFMD occurrence and climate [2ā11, 14, 15]. The other is compartmental differential equation model [12, 13, 16ā20]. Of them, Tiing and Labadin [16] were the first to propose a simple SIR epidemic model to predict the number of the infected and the duration of an outbreak when HFMD occurs and found that the disease spread quite rapidly and the parameter that may be able to be controlled would be the number of susceptible persons. Roy and Haider [12] and Roy [13] established a simple SEIR model to understand the dynamics of HFMD among young children and found that disease transmission depended more on the number of actively infected people in the population at the initial time and could be controlled by quarantine of more actively infected individuals. Wang and Xue [18] and Yang et al. [19] established SEIQR (susceptibleexposedinfectiousquarantinedrecovered) model, respectively, and estimated the basic reproduction number for the HFMD transmission in Mainland China and predicted that HFMD was an endemic disease in China.
After a susceptible individual is infected, he/she firstly enter the incubation period of HFMD, which period is about 3 to 7 days. After the incubation period, the infected will show some clinical symptoms, such as having a fever, poor appetite, malaise, and sore throat, and few people may develop dehydration, febrile seizures, encephalitis, meningitis, cardiomyopathy, and so forth. And the infected people will fully recover after 7 to 10 days [1]. The occurrence of HFMD has seasonal characteristics and is often associated with climate changes [11, 21]. Periodic changing in the birth rate [22] and seasonally changing in the contact rate [23] are often regarded as sources of periodicity.
In the sense of seasonally changing in the contact rate, Liu [17] and Ma et al. [20] established HFMD with periodic transmission rate, respectively. In [17], Liu constructed a periodic SEIQR epidemic model to simulate the dynamics of HFMD transmission and showed that quarantine in the children population had a positive impact on controlling the spread of HFMD. And Ma et al. [20] proposed a SEIIeQR HFMD model, and found that the recessive subpopulation played an important role in the spread of HFMD, and only quarantining the infected is not an effective measure in controlling the disease.
Wenzhou is a prefecturelevel city in southeastern Zhejiang, province in China. At the time of the 2010 Chinese census, 9,122,100 people lived in Wenzhou [24]. Since Wenzhou has a humid subtropical climate zone with an annual average 18.08Ā°C (64.5Ā°F), it is of particular public health significance to update molecular epidemiology of HFMD in Wenzhou. It is reported that, in 2012, there were 147,941 HFMD cases and 17 deaths in Zhejiang province, and it is the first of the āten legal infectious diseasesā; there were 41,438 HFMD cases and 4 deaths in Wenzhou [25]. The numbers of weekly reported HFMD cases in Wenzhou from March 5, 2010, to December 27, 2013, are listed in Table 1.

The aim of this study is to use mathematical modeling to gain some insights into the transmission dynamics of HFMD in Wenzhou and to assess the aforementioned preventive strategies. In particular, we aim to answer the following questions through our analytic and numerical results of the HFMD model.(1)What is the disease dynamics of the HFMD in Wenzhou?(2)How can we control the HFMD spread in Wenzhou?
The paper is organized as follows. In Section 2, we derive a SEIQRS HFMD model. In Section 3, we give the main theoretical results of the HFMD model. In Section 4, based on the HFMD data of the Wenzhou Center for Disease Control and Prevention, we perform some simulation results of the model and sensitivity analysis. In Section 5, we give a brief discussion.
2. Model Derivation
In order to establish an HFMD model in Wenzhou, we classify the population into five compartments according to their states: susceptible, exposed, infected, quarantined, and recovered, which are denoted by , and , respectively. We denote the total population by , that is . And the nonautonomous differential equations for HFMD model are and the meaning and units of each variable and constant in model (1) are as follows:ā: the recruitment rate of susceptible (),ā: the per capita natural mortality rate (),ā: the periodic transmission rate coefficient, which is a continuous, positive periodic function,ā (): the āpsychologicalā effect, which can be interpreted as follows: when an infectious disease appears and spreads in a region, people deter from risky behavior or from taking precautionary measures in relation to the disease outbreak; human behavior change consequently leads to reduction in number of real susceptible individuals or effective contact rates [26],ā: the mean incubation period (week),ā: the quarantined rate (); this means that some of the infected people will be hospitalized for treatment, and, thus, these quarantined are isolated from other subpopulations,ā: the infectious recover rate, which means that the infectious individuals recover and return to recovered compartment from compartments (week),ā: the quarantined rate, which means that the quarantined individuals recover and return to recovered compartment from compartments (week),ā: the diseaseinduced mortality for the infective individuals (),ā: the diseaseinduced mortality for the quarantined individuals (),ā: the progression rate of the recovered individual .
3. Dynamics Analysis of the HFMD Model
3.1. Preliminaries
Let be the standard ordered dimensional Euclidean space with a norm . For , we denote , if ; , if ; and , if .
We briefly present some main results of Floquet theory. Consider the following linear periodic system: where is an vector and is an matrix of principal periodic , such that .
Let be a fundamental solution matrix of (2); that is, a nonsingular matrix each of the columns of which is a solution of equation such that , the identify matrix. Floquet theory show that the fundamental matrix has the form , where is a periodic matrix with initial values and is a constant matrix. The matrix is the socalled principle or monodromy or Floquet transition matrix of (2). Let be the spectral radius of . It then follows from [27, 28] that is a matrix with all entries positive for each . By the PerronFrobenius theorem [29], is the principal eigenvalue of in the sense that it is simple and admits an eigenvector . The following result is useful for our subsequent comparison arguments.
Lemma 1 (see [30, Lemma 2.1]). Let . Then there exists a positive, periodic function such that is a solution of (2).
3.2. Positively Invariant Sets
The study of the dynamics of model (1) requires the introduction of the following important sets:
And we can easily obtain the following theorem of the positively invariant sets.
Theorem 2 (positively invariant sets). and are positive invariant sets for model (1).
Proof. It shows that is positively invariant. From model (1), the total population satisfies the following equation: Hence, by integration, we check Hence, which implies that is positively invariant with respect to model (1).
3.3. The Basic Reproduction Number
It is easy to see that model (1) always has one diseasefree equilibrium . By the definition of BacaĆ«r and Guernaoui [31] and the general calculation procedure in Wang and Zhao [32], we have It follows that
Furthermore, is nonnegative, and is cooperative in the sense that the offdiagonal elements of are nonnegative.
It is obvious that and are the eigenvalues of and are negative. Thus is stable; namely, . We can see that the eigenvalues of are the diagonal elements and are negative. So is stable; namely, .
Let , be the evolution operator of the linear periodic system. Consider That is, for each , the matrix satisfies where is the identity matrix. Thus, the monodromy matrix of (9) equals .
Following the method in Wang and Zhao [32], we let be periodic in and the initial distribution of infectious individuals. So is the rate of new infections produced by the infected individuals who are introduced at time . When , gives the distribution of those infected individuals who are newly infected by and remain in the infected compartments at time . Naturally, is the distribution of accumulative new infections at time produced by all those infected individuals introduced at time previous to .
Let be the ordered Banach space of all periodic functions from to , which is equipped with the maximum norm and the positive cone . Then we can define a linear operator. Consider is called the next infection operator and the spectral radius of is defined as the basic reproduction number. Consider for the periodic epidemic model.
Lemma 3 (see [32]). The following statements are valid:(i) if and only if ;(ii) if and only if ;(iii) if and only if .Thus, the diseasefree equilibrium of model (1) is asymptotically stable if and unstable if .
Now we introduce the linear periodic system with parameter . Let , be the evolution operator of system (14) on . Clearly, . Moreover, Hence, we derive Following the general calculation procedure in Wang and Zhao [32], the basic reproduction number is the unique solution of .
3.4. Extinction of the Disease
Theorem 4. The diseasefree equilibrium of model (1) is globally asymptotically stable if and unstable if .
Proof. If , according to Lemma 3, we have . From model (1), for , we know that Consider the following comparison system: that is, By Lemma 1, there exists a positive periodic function such that is a solution of (18), where . Then we conclude that as , which implies that the zero solution of model (18) is globally asymptotically stable. Applying the comparison principle, we known that ,āāand as . It follows that ,āāand as . Therefore, the diseasefree equilibrium is globally asymptotically stable. This completes the proof.
Remark 5. From Theorem 4, one can know that if , then is globally asymptotically stable, which means that the disease will vanish after some period of time. Therefore, it is interesting to study the diseasefree equilibrium for controlling infectious disease.
3.5. Uniform Persistence of the Disease
Define and . Denote as the unique solution of model (1) with the initial value . Let be the PoincarĆ© map associated with model (1); that is, where is the period. By applying the fundamental existenceuniqueness theorem [33], we know that is the unique solution of model (1) with . From Theorem 2, we know that is dissipative point. We then introduce the following lemma.
Lemma 6. If , then there exists a such that when , for any , one has
Proof. We proceed by contradiction to prove that If not, then for some . Without loss of generality, we can assume that , for all . By the continuity of the solutions with respect to the initial values, for all , we obtain For any , let , where and , which is the greatest integer less than or equal to . Then we have It follows that , for . Then . Thus, for , we have We consider the following comparison system Set . Then, model (28) can be rewritten as If , we obtain by Lemma 3. Then we can choose small enough such that . Once again by Lemma 1, it follows that there exists a positive periodic such that is a solution of model (28), where , which implies that as . By applying the comparison principle, we have , and before as . This is a contraction.
By using Lemma 6, we can obtain the following theorem.
Theorem 7. If , model (1) has at least one positive periodic solution and is uniformly persistent.
Proof. We first prove that is uniformly persistent with respect to (). First of all, we show that and are positively invariant, for any . From the first equation of model (1), we derive that
So, we have
Similarly,
Thus, is positively invariant. Clearly, is relatively closed in . Set
We claim that
Note that
We only need to prove that
In fact, for any , in the case where , it is clear that , , for all . From the second equation of model (1), we have
In the case where , then , for any , . Therefore, for sufficiently small . That is to say, for any , . This implies that and is the only fixed point of and acyclic in . Moreover, Lemma 6 implies that is an isolated invariant set in and , where is the stable set of . By the acyclicity theorem on uniform persistence for maps [34], it follows that is uniformly persistent with respect to .
Furthermore, taking advantage of Theorem 1.3.6 in Zhao [34], has a fixed point. Consider
Then we see that , , , ,āāand . We further prove that , , and . Suppose not, if , from the first equation of model (1), we derive that
with , . Then we have
The periodicity of implies , for all . Similarly, and . Thus, and is a positive periodic solution of model (1). This completes the proof.
Remark 8. Theorem 7 tells us that when the basic reproduction number , model (1) admits at least one positive periodic solution, and the infectives are ultimately bounded below by some positive constant , which indicates that HFMD becomes an endemic disease.
Remark 9. From Theorems 4 and 7, we can conclude that is the threshold parameter for the extinction and the uniform persistence of the disease, which partially answer our proposed question in the introduction: āwhat is the disease dynamics of the HFMD in Wenzhou?ā.
4. Simulation Results via Statistical Data in Wenzhou
In this section, by using model (1), we simulate the reported data of HFMD in Wenzhou and carry out the sensitivity analysis based on the parameters.
Thanks to the insightful work of Liu [17], Ma et al. [20], and Yang et al. [19], we estimate the values of all parameters of model (1) in details:(i)the recruitment rate of susceptible : from the statistical data of Wenzhou population [24], one can know that the average birth number of Wenzhou population from 2005 to 2010 is 102729; then we adopt the weekly birth population number ;(ii)the natural mortality rate : in this paper, we take the average death rate of Wenzhou population as the natural mortality rate [24];(iii)the transmission rate : there are some epidemiological models using sinusoidal function () for the seasonal varying transmission rate [16, 17, 19, 20]; in this paper, we use the periodic function with period (here weeks) as the transmission rates; by using the leastsquare fitting of the numbers of HFMD , we can estimate the values of , , and ;(iv)the āpsychologicalā effect : in this paper, we adopt ; it would be likely to note that when , the transmission rate is the classical standard incidence rate;(v)the mean incubation period : as discussed in the introduction section, a person who is exposed to HFMD viruses will exhibit the symptoms after 3 to 7 days, for simplicity; in this paper, we assume that the average incubating time is 4 days, then the mean incubation period can be determined by ;(vi)the quarantined rate : some of the infected people will be hospitalized to isolate for treatment; the ratio of these quarantined in infectious is called the quarantined rate; based on the statistical data, we can determine ;(vii)the infected recover rate : we assume that the average recovery time of infectious is 8 days; then the recovery rate is ;(viii)the quarantined recover rate : we assume that the average recovery time of quarantined is 16 days; then the recovery rate is ;(ix)the diseaseinduced mortality and : we assume that the death comes from the quarantined; that is, ; and by using the leastsquare fitting of HFMD number , we can determine ;(x)the progression rate : by using the leastsquare fitting of HFMD number , we can determine .
And the number of the initial susceptible population at the end of 2010, , is obtained from the Wenzhou City Population and Family Planning Commission [24]. The numbers of the initial infected and quarantined can be obtained from the reported data of HFMD. Because the numbers of the initial exposed population and the recovered population cannot be obtained directly, we derive by and by and .
The numerical simulation results of the model (1) about the number of HFMD infectious are shown in Figures 2ā4.
In Figure 1, we show the spread dynamics of HFMD in Wenzhou from March 2010 to November 2013, and the parameters in model (1) undefined in the above are taken as In this case, ; the diseasefree equilibrium of model (1) is globally asymptotically stable, which means the HFMD disease will vanish after long time. The HFMD will die out, but is it always true?
(a)
(b)
For the sake of learning the spread dynamics of HFMD in Wenzhou further, we display the simulation results via the statistical data 20112012 (Figure 2) and 20122013 (Figure 4), respectively. In Figure 2, with parameters then , which means that HFMD becomes an endemic disease in Wenzhou. In this case, from Theorems 4 and 7, one can know that model (1) has at least one positive periodic solution which is uniformly persistent. Figures 3(a) and 3(b), with the same parameters as in Figure 2, confirm this conclusion.
In Figure 4, with parameters and, hence, , which means the HFMD disease will vanish after long time.
Comparing Figure 1 with Figures 2 and 4, one can see that with partial or the whole HFMD data in Wenzhou from 2010 to 2013, with different time periodic data, then by using leastsquare fitting method, we can obtain different periodic transmission rate , and, hence, obtain different , so we can obtain different disease spread dynamics and draw different conclusions. This difference is induced by different data and initial values.
For the sake of learning the effects of the quarantined rate () and infected recovery rate () on the controlling of disease spreading, we consider the relations between the reproduction number with and , respectively. Figure 5(a) shows that the larger quarantined rate is, the less is; and that is to say, quarantine has a positive impact on controlling the spread of disease. Figure 5(b) has the similar tendency with Figure 5(a); that is, the larger quarantined rate is, the less is. And ( is the cycle of treatment); hence, we can conclude that the less the cycle of treatment is, the less is.
(a)
(b)
5. Conclusions and Discussions
In this paper, we construct an SEIQRS epidemic model with periodic transmission rates to investigate the spread of seasonal HFMD in Wenzhou. We show that the global dynamics of the HFMD model (1) can be governed by its reproduction number : if ; the diseasefree equilibrium of model (1) is globally asymptotically stable, which means the disease will die out; while , model (1) has at least one positive periodic solution and is uniformly persistent, which indicates that HFMD becomes an endemic disease; that is, the HFMD will break out.
Epidemiologically, we partially provide answers to the two questions proposed in the introduction: what is the spread dynamics of the HFMD in Wenzhou? How can we control the HFMD spread in Wenzhou?
From Figures 1, 2, and 4, we can obtain different disease spread dynamics and draw different conclusions based on the different data and initial values. In fact, a mathematical model may help to explain a system, to study the effects of different components, and to make predictions about disease spreading behaviour. We need to point out that all these three figures are partly explaining and predicting the spread dynamics of the HFMD in Wenzhou. Besides, in the cases of Figures 1, 2, and 4, we obtain different basic reproduction numbers , which can mathematically help us to determine whether an infectious disease can spread or not. There are two values of less than 1, and close to 1. It is worthy to note that, in the HFMD model (1), we assume that the whole population of Wenzhou (more than 9 million people) is homogeneously mixing. If we try to take children under fiveā17 years old and their family as the susceptible population (about 1 million people), the value of would be larger. All in all, we can conclude that the HFMD in Wenzhou persists under current circumstances; that is, the HFMD becomes an endemic disease and will break out in Wenzhou.
From Figure 5, we find that the quarantined rate () and the infected recovery rate () are very beneficial to control the disease outbreak. That is to say, for controlling the HFMD spread in Wenzhou, we must increase the quarantined rate or decrease the treatment cycle.
In [20], based on the statistical data HFMD in Shandong province, China, Ma et al. indicated that the quarantined rate has little effect on controlling the disease spread. This is very different from our present results. It is reported that the HFMD incidence strongly depends on the climate [19]. We think this difference may be induced by different climate: Shandong province is in north China, and Wenzhou is in south China.
Furthermore, from Figures 1, 2, and 4, one can see that the fittings have some flaws because of oscillation peak in the data of summer 2011 and 2012. We think this phenomenon may represent an alarm for next year prevalence. In addition, there may be more suitable incidence rate to fit the spread dynamics of the HFMD in Wenzhou. This is desirable in future studies.
It is well known that no vaccine currently exists to protect individuals from infection with the viruses that cause HFMD, but such vaccines are being developed [35]. HFMD is highly contagious and is transmitted by oropharyngeal secretions such as saliva or nasal mucus, by direct contact or fecaloral transmission. Preventive measures include avoiding direct contact with infected individuals (including keeping infected children home from school), proper cleaning of shared utensils, disinfecting contaminated surfaces, and proper hand hygiene. These measures have been proved to be effective in decreasing the transmission of the viruses responsible for HFMD [1, 35].
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to thank JingAn Cui and Meng Ye (Beijing University of Civil Engineering and Architecture) and Maoxing Liu and Yangjun Ma (North University of China) for very helpful suggestions and comments which led to the improvement of the original paper. This research is supported by the National Science Foundation of China (61373005) and Zhejiang Provincial Natural Science Foundation (LY12A01014).
References
 Centers for Disease Control and Prevention, āAbout hand, foot, and mouth disease,ā http://www.cdc.gov/handfootmouth. View at: Google Scholar
 J. P. Lott, K. Liu, M. Landry et al., āAtypical handfootandmouth disease associated with Coxsackievirus A6 infection,ā Journal of the American Academy of Dermatology, vol. 69, no. 5, pp. 736ā741, 2013. View at: Publisher Site  Google Scholar
 L.X. Mao, B. Wu, W.X. Bao et al., āEpidemiology of hand, foot, and mouth disease and genotype characterization of Enterovirus 71 in Jiangsu, China,ā Journal of Clinical Virology, vol. 49, no. 2, pp. 100ā104, 2010. View at: Publisher Site  Google Scholar
 M. Y. Liu, W. Liu, J. Luo et al., āCharacterization of an outbreak of hand, foot, and mouth disease in nanchang, China in 2010,ā PLoS ONE, vol. 6, no. 9, Article ID e25287, 2011. View at: Publisher Site  Google Scholar
 Q. Zhu, Y. Hao, J. Ma, S. Yu, and Y. Wang, āSurveillance of hand, foot, and mouth disease in Mainland China (20082009),ā Biomedical and Environmental Sciences, vol. 24, no. 4, pp. 349ā356, 2011. View at: Publisher Site  Google Scholar
 M. Zeng, Y. F. Li, X. H. Wang et al., āEpidemiology of hand, foot, and mouth disease in children in Shanghai 2007–2010,ā Epidemiology and Infection, vol. 140, no. 6, pp. 1122ā1130, 2012. View at: Publisher Site  Google Scholar
 L. W. Ang, B. K. Koh, K. P. Chan, L. T. Chua, L. James, and K. T. Goh, āEpidemiology and control of hand, foot and mouth disease in Singapore, 2001–2007,ā Annals of the Academy of Medicine Singapore, vol. 38, no. 2, pp. 106ā112, 2009. View at: Google Scholar
 S. Blomqvist, P. Klemola, S. Kaijalainen et al., āCocirculation of coxsackieviruses A6 and A10 in hand, foot and mouth disease outbreak in Finland,ā Journal of Clinical Virology, vol. 48, no. 1, pp. 49ā54, 2010. View at: Publisher Site  Google Scholar
 Y. L. Hii, J. Rocklöv, and N. Ng, āShort term effects of weather on Hand, foot and mouth disease,ā PLoS ONE, vol. 6, no. 2, Article ID e16796, 2011. View at: Publisher Site  Google Scholar
 S. S. Y. Wong, C. C. Y. Yip, S. K. P. Lau, and K. Y. Yuen, āHuman enterovirus 71 and hand, foot and mouth disease,ā Epidemiology and Infection, vol. 138, no. 8, pp. 1071ā1089, 2010. View at: Publisher Site  Google Scholar
 J.F. Wang, Y.S. Guo, G. Christakos et al., āHand, foot and mouth disease: spatiotemporal transmission and climate,ā International Journal of Health Geographics, vol. 10, article 25, 2011. View at: Publisher Site  Google Scholar
 N. Roy and N. Haider, āCompartmental modeling of hand, foot and mouth infectious disease (HFMD),ā Research Journal of Applied Sciences, vol. 5, no. 3, pp. 177ā182, 2010. View at: Publisher Site  Google Scholar
 N. Roy, āMathematical modeling of handfootmouth disease: quarantine as a control measure,ā International Journal of Advanced Scientific Engineering and Technological Research, vol. 1, no. 2, pp. 34ā44, 2012. View at: Google Scholar
 E. Ma, T. Lam, K. C. Chan, C. Wong, and S. K. Chuang, āChanging epidemiology of hand, foot, and mouth disease in Hong Kong, 2001–2009,ā Japanese Journal of Infectious Diseases, vol. 63, no. 6, pp. 422ā426, 2010. View at: Google Scholar
 D. Zhu, X. Zhao, Y. Yao et al., āA new factor influencing pathogen detection by molecular assay in children with both mild and severe hand, foot, and mouth disease,ā Diagnostic Microbiology and Infectious Disease, vol. 76, no. 2, pp. 162ā167, 2013. View at: Publisher Site  Google Scholar
 F. C. S. Tiing and J. Labadin, āA simple deterministic model for the spread of hand, foot and mouth disease (HFMD) in Sarawak,ā in Proceedings of the 2nd Asia International Conference on Modelling and Simulation (AMS '08), pp. 947ā952, IEEE, May 2008. View at: Publisher Site  Google Scholar
 J. Liu, āThreshold dynamics for a HFMD epidemic model with periodic transmission rate,ā Nonlinear Dynamics, vol. 64, no. 12, pp. 89ā95, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 X. Wang and Y. Xue, āAnalysis of HFMD in China using a SEIQR model,ā Far East Journal of Mathematical Sciences, vol. 58, no. 1, pp. 83ā98, 2011. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 J.Y. Yang, Y. Chen, and F.Q. Zhang, āStability analysis and optimal control of a handfootmouth disease (HFMD) model,ā Journal of Applied Mathematics and Computing, vol. 41, no. 12, pp. 99ā117, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Ma, M. Liu, Q. Hou, and J. Zhao, āModelling seasonal HFMD with the recessive infection in Shandong, China,ā Mathematical Biosciences and Engineering, vol. 10, no. 4, pp. 1159ā1171, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Ni, B. Yi, J. Yin et al., āEpidemiological and etiological characteristics of hand, foot, and mouth disease in Ningbo, China, 2008–2011,ā Journal of Clinical Virology, vol. 54, no. 4, pp. 342ā348. View at: Publisher Site  Google Scholar
 J. Ma and Z. Ma, āEpidemic threshold conditions for seasonally forced SEIR models,ā Mathematical Biosciences and Engineering, vol. 3, no. 1, pp. 161ā172, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Dushoff, J. B. Plotkin, S. A. Levin, and D. J. D. Earn, āDynamical resonance can account for seasonality of influenza epidemics,ā Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 48, pp. 16915ā16916, 2004. View at: Publisher Site  Google Scholar
 Wenzhou City Population and Family Planning Commission, āA survey of wenzhou population,ā http://www.wzrkjs.gov.cn, (in Chinese). View at: Google Scholar
 Zhejiang Online Source for Disease Control and Prevention, āThe statistic results of “Ten legal infectious diseases” in Zhejiang province in 2012,ā http://www.cdc.zj.cn, (in Chinese). View at: Google Scholar
 D. Xiao and S. Ruan, āGlobal analysis of an epidemic model with nonmonotone incidence rate,ā Mathematical Biosciences, vol. 208, no. 2, pp. 419ā429, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Aronsson and R. B. Kellogg, āOn a differential equation arising from compartmental analysis,ā Mathematical Biosciences, vol. 38, no. 12, pp. 113ā122, 1978. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. W. Hirsch, āSystems of differential equations that are competitive or cooperative. II. Convergence almost everywhere,ā SIAM Journal on Mathematical Analysis, vol. 16, no. 3, pp. 423ā439, 1985. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. L. Smith and P. Waltman, The Theory of the Chemostat, vol. 13 of Cambridge Studies in Mathematical Biology, Cambridge University Press, Cambridge, UK, 1995. View at: Publisher Site  MathSciNet
 F. Zhang and X.Q. Zhao, āA periodic epidemic model in a patchy environment,ā Journal of Mathematical Analysis and Applications, vol. 325, no. 1, pp. 496ā516, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 N. Bacaër and S. Guernaoui, āThe epidemic threshold of vectorborne diseases with seasonality. The case of cutaneous leishmaniasis in Chichaoua, Morocco,ā Journal of Mathematical Biology, vol. 53, no. 3, pp. 421ā436, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 W. Wang and X.Q. Zhao, āThreshold dynamics for compartmental epidemic models in periodic environments,ā Journal of Dynamics and Differential Equations, vol. 20, no. 3, pp. 699ā717, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Perko, Differential Equations and Dynamical Systems, Springer, New York, NY, USA, 2000.
 X.Q. Zhao, Dynamical Systems in Population Biology, Springer, New York, NY, USA, 2003. View at: MathSciNet
 N. Sarma, āHand, foot, and mouth disease: current scenario and indian perspective,ā Indian Journal of Dermatology, Venereology, and Leprology, vol. 79, no. 2, pp. 165ā175, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Yeting Zhu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.