Abstract and Applied Analysis

Abstract and Applied Analysis / 2014 / Article

Research Article | Open Access

Volume 2014 |Article ID 234509 | 11 pages | https://doi.org/10.1155/2014/234509

A Hand-Foot-and-Mouth Disease Model with Periodic Transmission Rate in Wenzhou, China

Academic Editor: Kaifa Wang
Received11 Jan 2014
Revised18 Feb 2014
Accepted20 Feb 2014
Published25 Mar 2014

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 disease-free 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 A4-A7, A9, A10, B1-B3, 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 [310]. In 2008, there were 488,955 cases reported, with a morbidity of 37/100,000, mortality of 0.0095/100,000, and ill-death 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 [211, 14, 15]. The other is compartmental differential equation model [12, 13, 1620]. 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 (susceptible-exposed-infectious-quarantined-recovered) 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 prefecture-level 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.


WeekDataWeekData WeekData WeekData WeekData WeekData WeekData WeekData

1 368 26 331 51 46 76 366 101 164 126 469 151 238 176 749
2 380 27 320 52 47 77 371 102 193 127 413 152 189 177 794
3 457 28 475 53 136 78 317 103 234 128 413 153 145 178 708
4 615 29 415 54 153 79 291 104 226 129 467 154 141 179 600
5 581 30 339 55 197 80 441 105 264 130 429 155 131 180 518
6 694 31 352 56 366 81 549 106 456 131 439 156 118 181 505
7 924 32 201 57 322 82 395 107 478 132 676 157 223 182 523
8 1124 33 254 58 401 83 478 108 743 133 861 158 268 183 520
9 1200 34 226 59 460 84 410 109 876 134 793 159 351 184 446
10 1539 35 243 60 549 85 565 110 1035 135 765 160 455 185 510
11 1744 36 182 61 574 86 508 111 1403 136 645 161 517 186 714
12 1706 37 173 62 767 87 543 112 1660 137 491 162 596 187 917
13 1528 38 187 63 1165 88 538 113 1908 138 542 163 577 188 672
14 1329 39 214 64 1247 89 787 114 2240 139 600 164 517 189 589
15 1252 40 186 65 1168 90 689 115 1974 140 644 165 552 190 482
16 1447 41 166 66 1060 91 854 116 2163 141 579 166 546 191 369
17 1663 42 192 67 1388 92 785 117 1545 142 691 167 668 192 326
18 1346 43 190 68 1686 93 618 118 1611 143 777 168 874 193 366
19 1247 44 181 69 1532 94 476 119 1868 144 783 169 1068 194 491
20 1094 45 171 70 1323 95 414 120 1683 145 732 170 1266 195 402
21 794 46 95 71 819 96 365 121 1631 146 561 171 1001 196 439
22 673 47 66 72 743 97 308 122 1228 147 636 172 1037 197 371
23 584 48 55 73 636 98 206 123 758 148 441 173 1363 198 358
24 448 49 37 74 602 99 162 124 502 149 345 174 1051 199 382
25 420 50 42 75 478 100 161 125 519 150 225 175 995

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 disease-induced mortality for the infective individuals (),: the disease-induced 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 so-called 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 Perron-Frobenius 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 disease-free 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 off-diagonal 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 disease-free 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 disease-free 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 disease-free 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 disease-free 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 existence-uniqueness 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 least-square 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 disease-induced mortality and : we assume that the death comes from the quarantined; that is, ; and by using the least-square fitting of HFMD number , we can determine ;(x)the progression rate : by using the least-square 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 24.

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 disease-free 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?

For the sake of learning the spread dynamics of HFMD in Wenzhou further, we display the simulation results via the statistical data 2011-2012 (Figure 2) and 2012-2013 (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 least-square 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.

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 disease-free 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 fecal-oral 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 Jing-An 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

  1. Centers for Disease Control and Prevention, “About hand, foot, and mouth disease,” http://www.cdc.gov/hand-foot-mouth. View at: Google Scholar
  2. J. P. Lott, K. Liu, M. Landry et al., “Atypical hand-foot-and-mouth 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
  3. 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
  4. 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
  5. Q. Zhu, Y. Hao, J. Ma, S. Yu, and Y. Wang, “Surveillance of hand, foot, and mouth disease in Mainland China (2008-2009),” Biomedical and Environmental Sciences, vol. 24, no. 4, pp. 349–356, 2011. View at: Publisher Site | Google Scholar
  6. 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
  7. 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
  8. S. Blomqvist, P. Klemola, S. Kaijalainen et al., “Co-circulation 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
  9. 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
  10. 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
  11. 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
  12. 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
  13. N. Roy, “Mathematical modeling of hand-foot-mouth 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
  14. 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
  15. 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
  16. 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
  17. J. Liu, “Threshold dynamics for a HFMD epidemic model with periodic transmission rate,” Nonlinear Dynamics, vol. 64, no. 1-2, pp. 89–95, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  18. 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
  19. J.-Y. Yang, Y. Chen, and F.-Q. Zhang, “Stability analysis and optimal control of a hand-foot-mouth disease (HFMD) model,” Journal of Applied Mathematics and Computing, vol. 41, no. 1-2, pp. 99–117, 2013. View at: Publisher Site | Google Scholar | MathSciNet
  20. 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
  21. 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
  22. 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
  23. 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
  24. Wenzhou City Population and Family Planning Commission, “A survey of wenzhou population,” http://www.wzrkjs.gov.cn, (in Chinese). View at: Google Scholar
  25. 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
  26. 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
  27. G. Aronsson and R. B. Kellogg, “On a differential equation arising from compartmental analysis,” Mathematical Biosciences, vol. 38, no. 1-2, pp. 113–122, 1978. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  28. 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
  29. 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
  30. 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
  31. N. Bacaër and S. Guernaoui, “The epidemic threshold of vector-borne 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
  32. 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
  33. L. Perko, Differential Equations and Dynamical Systems, Springer, New York, NY, USA, 2000.
  34. X.-Q. Zhao, Dynamical Systems in Population Biology, Springer, New York, NY, USA, 2003. View at: MathSciNet
  35. 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 © 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.


More related articles

1826 Views | 967 Downloads | 13 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.