# SIR Model for Dengue Disease with Effect of Dengue Vaccination

Authors: Pratchaya Chanprasopchai; I. Ming Tang; Puntani Pongsumpun
Journal: Computational and Mathematical Methods in Medicine
 (2018)

---

## Abstract

The dengue disease is caused by dengue virus, and there is no specific treatment. The medical care by experienced physicians and nurses will save life and will lower the mortality rate. A dengue vaccine to control the disease is available in Thailand since late 2016. A mathematical model would be an important way to analyze the effects of the vaccination on the transmission of the disease. We have formulated an SIR (susceptible-infected-recovered) model of the transmission of the disease which includes the effect of vaccination and used standard dynamical modelling methods to analyze the effects. The equilibrium states and their stabilities are investigated. The trajectories of the numerical solutions plotted into the 2D planes and 3D spaces are presented. The main contribution is determining the role of dengue vaccination in the model. From the analysis, we find that there is a significant reduction in the total hospitalization time needed to treat the illness.


---

## Body

## 1. Introduction

Dengue disease is a mosquito-borne viral infection caused by 4 serotypes of dengue virus, DEN-1, DEN-2, DEN-3, and DEN-4. Dengue disease is widely spread in tropical and subtropical regions of the world. Dengue virus is transmitted to human by the bite of the female mosquito of the speciesAedes aegypti and Aedes albopictus [1]. An estimated 3.9 billion people in 128 countries are at risk to this disease. The countries at danger to infection by the dengue viruses around the world are shown in Figure 1 [2].Figure 1
Distribution of global dengue risk [2].Thailand is located in the tropical region where dengue virus is widely circulating. Dengue is spreading nationwide in Thailand including the Bangkok metropolitan area. Thailand is in special danger since three of the four species of the dengue virus have been found in Thailand, and both of theAedes vector species are present. The Bureau of Epidemiology, Ministry of Public Health, has reported dengue cases in all provinces in 2016, a total of 63,931 cases with 64 deaths [3]. At the present time, there is no special treatment for dengue disease, but early detection and the appropriated medical care will decrease the fatality rates. A dengue vaccine would be another way to reduce the fatality rates. WHO reported the first dengue vaccine, called as Dengvaxia (CYD-TDV). It was registered in several countries in late 2015 and early 2016. It was recommended for use only in high dengue disease burden countries such as Thailand [1]. Dengue vaccine against four strains of the dengue virus was first launched in Thailand in late 2016. The vaccine would be suitable for use in individuals between 9 and 45 years of age living in endemic areas. Since the reported incidence of dengue peaks in the rainy season between June and September, the vaccination should be done in advance of the peak period in order for the immunity to develop.There were many mathematical models for describing and analyzing the behaviors of dengue disease. Esteva and Vargus [4] proposed an SIR (susceptible-infected-recovery) model to describe the transmission of dengue disease with constant human and vector populations while Chanprasopchai et al. [5] proposed a SEIR (susceptible-exposed-infected-recovered) model for Thailand to determine the effect of the rainfall on the spread of dengue in the Thailand model. The transmission of dengue disease is assumed to depend on the nature of the rainfall in different countries. The stability of the solution of the model was then analyzed. Numerical results taking into account the rainfall were obtained, and they were seen to correspond to the analytical results. Using standard dynamical analysis techniques, Chanprasopchai and Pongsumpun [6] established relations between the different variables in an SIR model of the dengue transmission model in which the biting rate of mosquito became as factor. Pongsumpun and Tang [7] analyzed the transmission of dengue hemorrhagic fever in an SIR model which included an age structure in human population.Recently, Shim [8] studied the recently approved dengue vaccination program in the Philippines and showed that with appropriated pricing of dengue vaccination, reduction of the burden of the dengue disease in the Philippines and a significant potential to confer the excellent value were possible. Recker et al. [9] reported that the availability of epidemiological and clinical data from the trials of vaccine provided a great opportunity for formulating mathematical models in which the vaccine efficacy depends on the serotype, age, host immune status, and severity. Mathematical modelling becomes a valuable tool in the policy-making process to estimate what the consequences of any decisions taken could be. Aguiar et at. [10] studied the impact of the tetravalent dengue vaccine (Dengvaxia™) in the dengue endemic countries. They divided the human population into four age groups (below 9 years, between 9 and 45 years and naturally infected with dengue fever, between 9 and 45 years and vaccinated, and older than 45 years). The reasons for dividing the humans into these four groups are the differences in the observed response of patients of different ages to the vaccine. Sridhar et al. [11] reported that the effect of the vaccine was different for people in these four groups. Aguiar et al. [10] also took into account the presence of two strains of the dengue virus in their model even though the vaccine was made with attenuated viruses of the four serotypes of the dengue virus (DEN-1, DEN-2, DEN-3, and DEN-4). In their model, there were thirty-one (31) population categories with differential equations defining the time evolution of each category. In their model, however, the transmission of the dengue virus was taken to be due to a contact between a susceptible person and an infected person even though dengue fever is abovirus disease where the transmission occurs when a susceptible vector (mosquito) bites an infected human and then inturn the infected mosquito bites a susceptible human. Sriprom et al. [12] have studied the transmission of DF when there are two dengue viruses present and no age structure in the human population and where the transmission of disease is via the mosquitoes. The transmission of a single strain of the virus in an age-structured human population has been studied by Pongsumpun and Tang [7].In this study, we propose an SIR mathematical model to analyze the behaviors of the transmission of dengue fever when vaccination with the Dengvaxia has been carried out. To be able to conduct a simple analytical analysis of the results, only one strain of the DF virus was taken into account and no age structure in the human populations was assumed. Otherwise the mathematics would become intractable. Only inferences based on the numerical simulations of the solutions to the multiple differential equations (31 for the model used in [10]) describing the time evolution of the many population categories can be drawn. The impact of the vaccination on the transmission in our model is expressed through the parameter α (the efficacy) appearing in (1). We note that this parameter appears in most of the analytical expressions for the various population categories and the basic reproduction rate R0. The standard analysis method using the Routh–Hurwitz criteria is applied to investigate the system stability in which the dynamical transmission model of dengue disease, equilibrium state, stability, numerical simulation, results, and conclusion are presented.


## 2. Materials and Methods

In our SIR model, the population is divided into 2 populations, a human and a vector population. The human population consists of three epidemiological states: susceptible humans (SH¯), infected humans (IH¯), and recovered humans (RH¯), while the vector population has two epidemiological states: susceptible vector (SV¯) and infected vector (IV¯). Mosquito has no recovery state since the mosquito dies before it can recover from the disease. The susceptible mosquito state is unimmune and uninfected, while in the infected state, it is infected with dengue virus and can transmit the virus. The recovery state in the human population is a person who has recovered from an infection by the dengue virus. We assume that the human and vector populations are constant. The dynamical transmission in human and mosquito population with effect of vaccination is shown in Figure 2.Figure 2
Dynamical transmission in humans and mosquitoes with the effect of vaccination incorporated.Note. SH¯t = number of susceptible human population who are unvaccinated at time t; IH¯t = number of infected human population who are unvaccinated at time t; RH¯t = number of recovered human population who are unvaccinated at time t; SHV¯t = number of susceptible human population who have been vaccinated at time t; IHV¯t = number of infected human population who have been vaccinated at time t; SV¯t = number of susceptible vector at any time t; IV¯t = number of infected vector at any time t; dH, dV = death rates of human and vector populations; NH, NV = total human and vector populations; βH, βV = transmission rate of dengue virus from vector to human and human to vector; bH = birth rate of human population; A = constant recruitment rate of vector population; α = vaccine efficacy; p = fraction of newborns vaccinated.The transmission model of dengue disease with effect of vaccination can be described by the following differential equations:(1)dSH¯dt=1−pbHNH−βHSH¯IV¯−dHSH¯,dIH¯dt=βHSH¯IV¯−γIH¯−dHIH¯,dRH¯dt=γIH¯+IHV¯−dHRH¯,dSHV¯dt=pbHNH−1−αβHSHV¯IV¯−dHSHV¯,dIHV¯dt=1−αβHVSH¯IV¯−γIHV¯−dHIHV¯,dSV¯dt=A−βVSV¯IH¯+IHV¯−dVSV¯,dIV¯dt=βVSV¯IH¯+IHV¯−dVIV¯.The effect of the vaccination is incorporated into our model by the presence of the efficacy coefficientα. Setting it to 1, there will be no infected mosquitoes to transmit the virus and so there will be no dengue fever epidemic. This is however not the way the vaccination works. When a person is vaccinated with a vaccine made from live-attenuated virus, weak attenuated viruses are introduced into the human where the body immune system produces the agents which will provide immunity to further infection of the virus. Invariably, some people will not have an immune response able to the kill the weak virus, and weak virus will develop into a strong virus able to infect the person, producing an infectious person. The effect of the vaccination is to create a new pathway for the transmission of the disease.The total human and vector populations are assumed to be governed by the following conditions:(2)SH¯+IH¯+RH¯+SHV¯+IHV¯=NH,SV¯+IV¯=NV.If total human and vector populations are constants, then the rate of change for total human and vector populations is 0. As the result, we will have the following equations:(3)dSH¯dt+dIH¯dt+dRH¯dt+dSHV¯dt+dIHV¯dt=0,dSV¯dt+dIV¯dt=0,NV=AμV,bH=dH.Normalizing the equations by introducing the following normalized variables:(4)SH=SH¯NH,IH=IH¯NH,RH=RH¯NH,SHV=SHV¯NH,IHV=IHV¯NH,SV=SV¯NV,IV=IV¯NV.Introducing these normalized variables into (1), we get the new set of equations of the following states:(5)dSHdt=1−pbH−βHSHIVNV−dHSH,dIHdt=βHSHIVNV−γIH−dHIH,dSHdt=pbH−1−αβHSHVIVNV−dHSHV,dIHVdt=1−αβHSHVIVNV−γIHV−dHIHV,dIVdt=βVSVIH+IHVNH−dVIV.The equilibrium states are obtained by setting the right-hand side of (5) to be 0. By doing this, we obtain an expression for something known as the basic production number R0. This number is defined as(6)R0=ε1−−2+αε2dV+NH1−pαdH+ε3NVβHβVε12α2ε22dV2+2αε2dVNH1+p−2+αdH+−1+2pε3NVβHβV+NH2−1+pαdH+ε3NVβH2βV2.WhenR0 ≤ 1, the equilibrium state will be the disease-free state E1 defined as(7)E1t=SH=1−p,IH=0,SHV=p,IHV=0,IV=0,and when R0 > 1, the equilibrium state is the endemic state defined as(8)E2t=SH∗t,IH∗t,SHV∗t,IHV∗t,IV∗t,where(9)SH∗t=ε1αε2dV+ε3ε4+ε52αdHε6,IH∗t=ε1−αε2dV+1+α−2pαε4−ε52αε2ε6,SHV∗t=ε1−αε2dV+ε3ε4+ε52ε3αdHε6,IHV∗t=ε1αε2dV+ε3−1+2pαε4−ε52ε3αε2ε6,IV∗t=ε1−−2+αε2dV+ε3ε4−ε52ε3ε2dVNV2βH2,with (10)ε1=dHNVβH,ε2=γ+dH,ε3=−1+α,ε4=NHNVSVβHβV,ε5=dH2NV2βH2α2γ+dH2dV2+2−1+2p−1+ααγ+dHdVNHNVSVβHβV+−1+α2NH2NV2SV2βH2βV2,ε6=NHNV2SVβH2βV.The equilibrium states are local asymptotically stable if all the eigenvalues have negative real parts. The eigenvalues (λ) are obtained by solving the eigenvalue matrix equation(11)detJ−λI=0,where J is the Jacobian matrix of each equilibrium point, λ is the eigenvalue, and I is the identity matrix.The Jacobian matrix of system (6) is as follows:(12)J=−βHIVNV−dH000−βHSHNVβHIVNV−γ−dH00βHSHNV00−1−αβHIVNV−dH0−1−αβHSHVNV001−αβHIVNV−γ−dH1−αβHSHVNV0βVSV0βVSV−dV.Constructing the Jacobian matrix from (5) and evaluating it at the two equilibrium points, we obtain the eigenvalue equation(13)−λ−γ−dHλ4+e1λ3+e2λ2+e3λ1+e4=0,for the disease-free state E1 and the eigenvalue equation(14)λ5+e1λ4+e2λ3+e3λ2+e4λ1+e5=0,for the endemic state E2.The eigenvalues of disease-free equilibrium state will have negative real parts when the coefficients of (13) have values satisfying the Routh–Hurwitz criteria(15)e1>0,e3>0,e4>0,e1e2e3>e32+e12e4.The eigenvalues of the endemic equilibrium state will have negative real parts when the coefficients of (14) have values which satisfy a different Routh–Hurwitz criterion(16)e1>0,e2>0,e3>0,e4>0,e5>0,e1e2e3−e32−e12e4>0,e1e4−e5e1e2e3−e32−e12e4−e5e1e2−e32−e1e52>0.


## 3. Numerical Results

The transmission of dengue disease in this study is based on the SIR model with vaccination. The nonzero values ofα and p are the parameters pertaining to the vaccination program. The numerical simulations were done using the following values of parameters: dH = 1/(65∗365) per day corresponding to a life expectancy of 65 years for the Thai people and dV = 1/12 corresponding to a life expectancy of 12 days of mosquito population. For the disease-free equilibrium state, the parameter values were A = 1,000; NH = 1,000; γH = 1/3; βH = 0.000012; βV = 0.000012; p = 0.8; and α = 0.8, while the parameters value of the endemic equilibrium state were A = 500; NH = 500; γH = 0.03; βH = 0.000045; βV = 0.000045; p = 0.75; and α = 0.75. These numerical values in the first set gave R0 < 1, while the values in the second set gave R0 > 1. The trajectories of the numerical simulations for disease-free and endemic states of SH, IH, SHV, IHV, and IV are shown in Figures 3 and 4, respectively. The trajectories of the numerical simulation for disease-free and endemic states plotted in the 2D planes (SH, IH), (SH, SHV), (SH, IHV), (SH, IV), (IH, SHV), (IH, IHV), (IH, IV), (SHV, IHV), (SHV, IV), and (IHV, IV) planes are shown in Figures 5 and 6, respectively. The trajectories of the numerical solutions for disease-free and endemic states plotted in the 3D spaces (SH, IH, SHV), (SH, IH, IHV), (SH, IH, IV), (SH, SHV, IHV), (SH, SHV, IV), and (SH, IHV, IV) spaces are shown in Figures 7 and 8, respectively.Figure 3
The trajectory of (a)SH, (b) IH, (c) SHV, (d) IHV, and (e) IV towards the disease-free equilibrium state (E1).
(a)
(b)
(c)
(d)
(e)Figure 4
The trajectory of (a)SH, (b) IH, (c) SHV, (d) IHV, and (e) IV towards the endemic equilibrium state (E2).
(a)
(b)
(c)
(d)
(e)Figure 5
The trajectories of dengue disease for disease-free equilibrium projected onto (a) (SH, IH), (b) (SH, SHV), (c) (SH, IHV), (d) (SH, IV), (e) (IH, SHV), (f) (IH, IHV), (g) (IH, IV), (h) (SHV, IHV), (i) (SHV, IV), and (j) (IHV, IV) planes.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)Figure 6
The trajectories of dengue disease for endemic equilibrium projected onto (a) (SH, IH), (b) (SH, SHV), (c) (SH, IHV), (d) (SH, IV), (e) (IH, SHV), (f) (IH, IHV), (g) (IH, IV), (h) (SHV, IHV), (i) (SHV, IV), and (j) (IHV, IV) planes.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)Figure 7
The trajectories of dengue disease for disease-free equilibrium projected onto (a) (SH, IH, SHV), (b) (SH, IH, IHV), (c) (SH, IH, IV), (d) (SH, SHV, IHV), (e) (SH, SHV, IV), and (f) (SH, IHV, IV) spaces.
(a)
(b)
(c)
(d)
(e)
(f)Figure 8
The trajectories of dengue disease for endemic equilibrium projected onto (a) (SH, IH, SHV), (b) (SH, IH, IHV), (c) (SH, IH, IV), (d) (SH, SHV, IHV), (e) (SH, SHV, IV), and (f) (SH, IHV, IV) spaces.
(a)
(b)
(c)
(d)
(e)
(f)


## 4. Discussion and Conclusion

In this study, the dynamical transmission of dengue disease based on an SIR model where a dengue vaccination campaign in the human population has occurred is studied. Again, it is found that the model system has two equilibrium points, a disease-free and an endemic state. The occurrence of the two equilibrium states depend on whetherR0 < 1 and R0 > 1 where R0 is the basic reproduction number or number of secondary infection caused by an initial infection. The conditions for the stability of the disease-free and endemic equilibrium states were established.It should be noted that the population was not divided into age groups or the presences of a second, third, or fourth serotype of the dengue virus were taken into account. The different population categories include all peoples between the age of one and sixty-five. The time series solutions of the disease-free and endemic equilibrium states are presented in Figures3 and 4, respectively. In Figure 3, we see that the different population groups except for the recovered population decay monotonically to their equilibrium for R0 < 1. Figure 4 shows that the time evolutions of the population categories oscillate before reaching their equilibrium values when R0 > 1. Interestingly, the solutions for infected populations with and without vaccination exhibit oscillations in two time intervals with the oscillation between two intervals becoming weak. The oscillation of the infected population without vaccination in first interval is stronger than that in the second time interval. This is opposite to the behavior of the oscillation of the infected population with vaccination. The trajectories of the disease-free and endemic equilibrium states onto 2D planes are shown in Figures 5 (when R0 < 1) and 6 (when R0 > 1), while the trajectories of the disease-free and endemic equilibrium states projected onto 3D planes are shown in Figures 7 (when R0 < 1) and 8 (when R0 > 1), respectively. In both the 2D and 3D projections, the trajectories towards the equilibrium values are smooth when R0 < 1. However, when R0 > 1, the trajectories exhibit oscillatory behavior in both the 2D and 3D projections. The conditions for the stability of disease-free and endemic equilibrium states were established. The time series solution of disease-free and endemic equilibrium states are presented in Figures 3 and 4, respectively. The trajectories of disease-free and endemic equilibrium projected onto 2D planes are shown in Figures 5 and 6, while the trajectories of disease-free and endemic equilibrium projected onto 3D planes are shown in Figures 7 and 8, respectively.In order to analyze the effect of dengue vaccination, we have investigated both the disease-free and endemic equilibrium states using different values of parameters which would giveR0 < 1 and R0 > 1. These same set of numerical values were used for numerical simulation with and without the influence of dengue vaccination campaign. α = 0 and p = 0 were used in the simulation to get the trajectories in the case where there was no vaccine administered. The influence of dengue vaccination is seen in Figures 9 and 10. Figure 9 shows that the disease-free state is sooner when there are dengue virus vaccines administered than when there are not vaccines administered. This means that the hospitalization time can be reduced. Figure 10 shows the effects of the vaccine when the parameters are such that the endemic state is in the equilibrium state.Figure 9
The infected human without vaccination (a) and with vaccination (b) by comparison of time series to disease-free equilibrium point.
(a)
(b)Figure 10
Infected human without vaccination (a) and with vaccination (b) by comparison of time series to endemic equilibrium point.
(a)
(b)The presence of oscillations around the endemic equilibrium stateE2 means that the imaginary part of the eigenvalue is not zero. For the simulation shown in Figure 4, the imaginary part of the complex roots is approximately 0.000238428. This leads to an estimate of the period of the oscillations or Tperiod = 2π/ω, where ω = imaginary part of λ or 2π/0.000238428 ≈ 72.20 years. This value is the approximation to the period of the solutions [4].In any vaccination campaigns, one must take into account the difference in the efficacy of the vaccine. It may not be the same for all age groups. Since one is not sure about the safety of the vaccine to children, the vaccination has been recommended only for people between the ages of 9 and 45. As of now, the vaccination schedule consists of 3 injections of 0.5 mL administered at 6-month intervals, given on a 0/6/12 month schedule [13].The campaign in Thailand began in December 2016, and information on efficacy of the vaccine against the different serotypes and the difference in the efficacy for different age groups is being collected. Dengue disease in Thailand occurs in urban and suburban areas [14–17] with peak transmission rates during the rainy season [5, 18]. Seasonal and climate affect the dengue fluctuation [19–22]. At present, it is not recommended to give dengue vaccination to pregnant women and travelers or health-care workers at this time due to lack of sufficient data.Since the present model does not take into account the age structure of the human population, the presence of more than one serotype of the dengue virus, and a proper treatment of how the vaccine interacts with a susceptible human being after it is administered to a human, one cannot answer the present question now facing the Public Health Community: Is the tetravalent vaccine safe or dangerous to the communities in which a vaccination program being carried out? In a community in which no dengue epidemic exist, for example, no dengue virus of a particular serotype exist, Does the small possibility of introducing that virus into the community where there is no chance of the virus entering into the community through natural means worth the risk? Does the chance of the more virulent form of the dengue fever (dengue hemorrhagic fever DHF or dengue shock syndrome, DSS) occurring in community in which only one serotype of DV is circulating worth the risk? In Thailand, three of the serotypes are circulating, so the use of the vaccine does not impose additional risk. These questions cannot be answered on the basis of the results of our model. We have shown that the use of the vaccine is beneficial.



---

*Source: 9861572-2018-08-28.xml*