# A Stochastic Model for Malaria Transmission Dynamics.

1. IntroductionMalaria is an infectious disease caused by the Plasmodium parasite and transmitted between humans through bites of female Anopheles mosquitoes. Approximately half of the worlds population is at risk of malaria. It remains one of the most prevalent and lethal human infections throughout the world. An estimated 40% of the world's population lives in malaria endemic areas. Most cases and deaths occur in subSaharan Africa. It causes an estimated 300 to 500 million cases and 1.5 to 2.7 million deaths each year worldwide. Africa shares 80% of the cases and 90% of deaths [1]. According to the website of the World Health Organization [2] there were approximately 214 million new cases of malaria and 438,000 deaths worldwide in 2015. Most cases were reported in the African region.

Recently, the incidence of malaria has been rising due to drug resistance. Various control strategies have been taken to reduce malaria transmissions. Since the first mathematical model of malaria transmission was introduced by Ross [3], quite a number of mathematical models have been formulated to investigate the transmission dynamics of malaria. Xiao and Zou [4] used mathematical models to explore a natural concern of possible epidemics caused by multiple species of malaria parasites in one region. They found that epidemics involving both species in a single region are possible.

Li and others [5] considered fast and slow dynamics of malaria model with relapse and analyzed the global dynamics by using the geometric singular perturbation theory. They suggested that a treatment should be given to symptomatic patients completely and adequately rather than to asymptomatic patients. On the other hand, for the asymptomatic patients, their results strongly suggested that, to control and eradicate the malaria, it is very necessary for governments to control the relapse rate strictly. Relapse is when symptoms reappear after the parasites had been eliminated from blood but persist as dormant hypnozoites in liver cells [6]. This commonly occurs between 8 and 24 weeks and is commonly seen with P. vivax and P. ovale infections. Other papers also consider the influence of relapse in giving up smoking or quitting drinking; please see [7].

Chitnis et al. [8] and Li et al. [5] assumed that the recovered humans have some immunity to the disease and do not get clinically ill but they still harbour low levels of parasite in their blood streams and can pass the infection to mosquitoes. After some period of time, they lose their immunity and return to the susceptible class. Unfortunately, Li and others did not consider that the recovered humans will return to their infectious state because of incomplete treatment.

Stochasticity is fundamental to biological systems. In some situations the system can be treated as a large number of similar agents interacting in a homogeneously mixing environment, and so the dynamics are well-captured by deterministic ordinary differential equations. However, in many situations, the system can be driven by a small number of agents or strongly influenced by an environment fluctuating in space and time [9-12].

Stochastic models incorporate discrete movements of individuals between epidemiological classes and not average rates at which individuals move between classes [13-15]. In stochastic epidemic models, numbers in each class are integers and not continuously varying quantities [13]. A significant possibility is that the last infected individual can recover before the disease is transmitted and the infection can only reoccur if it is reintroduced from outside the population [16]. In contrast, most deterministic models have the flaw that infections can fall to very low levels--well below the point at which there is only one infected individual only to rise up later [17]. In addition, the variability introduced in stochastic models may result in dynamics that differ from the predictions made by deterministic models [16].

For a large population size and a large number of infectious individuals, the deterministic threshold [R.sub.0] >1 provides a good prediction of a disease outbreak. However, this prediction breaks down when the outbreak is initiated by a small number of infectious individuals. In this setting, Markov chain (MC) models with a discrete number of individuals are more realistic than deterministic models where the number of individuals is assumed to be continuous-valued [18].

Motivated by these works, in this paper, we propose a model which is an extension of the model formulated by Huo and Qui (2014), who assumed that the pseudorecovered humans can recover and return to the susceptible class or relapse and become infectious again. Using the extended model, we will formulate the basic reproductive number [R.sub.0] and use it to compare the disease dynamics of the deterministic and stochastic models in order to determine the effect of randomness in malaria transmission dynamics.

This paper is organized as follows; in Section 2, we present a malaria transmission deterministic model with relapse, which is an extension of the model in [6]. We compute the basic reproduction number, [R.sub.0], of the malaria transmission deterministic model using the next-generation matrix approach. The stochastic version of the deterministic model and its underlying assumptions necessary for model formulation are presented and discussed in Section 3. In this section, we also compute the stochastic threshold for disease extinction or invasion by applying the multitype Galton-Watson branching process. In Section 4, we show the relationship between reproductive number of the deterministic model and the thresholds for disease extinction of the stochastic version; we also illustrate our results using numerical simulations. We conclude with a discussion of the results in Section 5.

2. The Malaria Deterministic Model

The interaction of the Mosquito-host is shown in Figure 1.

2.1. Model Formulation. In this section, we introduce a deterministic model of malaria with relapse, an extended form of the model in [6]. The model proposes a more realistic mathematical model of malaria, in which it is assumed that the pseudorecovered humans interact with infected mosquitoes and may acquire more parasites causing reinfection or, due to incomplete treatment, the infection my reoccur (relapse) and return to infectious class. Also the pseudorecovered human may lose their immunity and return to the susceptible class. Given that humans might get repeatedly infected due to not acquiring complete immunity, then the population dynamics are assumed to be described by the SIRS model; hence we consider a deterministic compartmental model which divides the total human population size at time t, denoted by N(t), into susceptible individuals Sh(t) (those who are not currently harbouring the parasite but are liable to be infected), infectious individuals Ih(t) (those already infected and are able to transmit the disease to mosquitoes), and pseudorecovered individuals Rh(t) (those who are treated from the disease but with partial recovery and hence can transmit the disease to mosquitoes). Mosquitoes are assumed not to recover from the parasites so the mosquito population can be described by the SI model.

The structure of model is shown in Figure 2.

2.2. Variable and Parameter Description for the Model. The variables for the model are summarized in description of variables for the malaria transmission model in the Notations.

The parameters for the model are described as in description of parameters for the malaria transmission model in the Notations.

2.3. Equations of the Model. Assuming that the disease transmits in a closed system which translates into the simplifying assumption of a constant population size with the birth rate equal to death rate, that is, [lambda] = [mu] (see [21]) so that N(t) = N, hence, the above assumptions lead to the following system of differential equations which describe the interaction between mosquitoes and humans:

[mathematical expression not reproducible] (1)

where [S.sub.h], [I.sub.h], [R.sub.h], [S.sub.m], and [I.sub.m] represent the number of susceptible humans, infectious humans, recovered humans, susceptible mosquitoes, and infectious mosquitoes, respectively.

2.4. Computation of [R.sub.0] Using the Next-Generation Matrix Approach. The basic reproduction number, [R.sub.0], is defined as the secondary infections produced by one infective agent that is introduced into an entirely susceptible population at the disease-free equilibrium. The next-generation matrix approach is frequently used to compute [R.sub.0]. The original nonlinear system of ODEs including these compartments can be written as [partial derivative] [X.sub.i]/[partial derivative]t = F - V, where F = ([F.sub.i]) and V = ([V.sub.i]) represent new infections and transfer between compartments, respectively [22-24].

The Jacobian matrices of F(X) and V(X) at the disease-free equilibrium [E.sub.0] are, respectively,

[mathematical expression not reproducible] (2)

The matrix J = F -V is the Jacobian matrix evaluated at the DFE.

[mathematical expression not reproducible] (3)

The inverse matrix of V is

[mathematical expression not reproducible] (4)

The matrix F[V.sup.-1] is called the next-generation matrix. The (i,k) entry of F[V.sup.-1] indicates the expected number of new infections in compartment i produced by the infected individual originally introduced into compartment k. The model reproduction number, [R.sub.0], which is defined as the spectral radius of F[V.sup.-1] and denoted by p(F [V.sup.-1]), is given by

[mathematical expression not reproducible] (5)

Here [R.sub.0] is associated with disease transmission by infected humans as well as the infection of susceptible humans by infected mosquitoes.

Simplifying (5) we find an equivalent form, which is defined as the product of the transmission from mosquito to human and from human to mosquito as follows:

[mathematical expression not reproducible] (6)

In (6), [R.sub.0] is associated with disease transmission by infected mosquitoes as well as infection of susceptible mosquitoes by infected humans. The term [beta]/[[eta].sub.1] represents the number of infected humans generated by infectious mosquito in its life span. The term [[alpha].sub.1]/([gamma] + [[mu].sub.1]) represents the number of infected mosquitoes generated by infectious human during the infectious period of the individual while the term [[alpha].sub.2]/([[rho].sub.1] + [[rho].sub.2] + [[mu].sub.1]) represents the number of infected mosquitoes generated by pseudorecovered human during his/her infectious period. Susceptible humans acquire infection following effective contacts with infected mosquitoes. Susceptible mosquitoes acquire malaria infection from infected humans in two ways, namely, by infectious human or pseudorecovered humans.

Also, from (6), we can rewrite [R.sup.0] as infection between mosquitoes and infectious humans or infection between mosquitoes and pseudorecovered humans.

[R.sup.0] = [square root of [R.sub.01] + [R.sub.02]], (7)

where

[R.sub.01] -([beta]/[[eta].sub.1]) ([[alpha].sub.1]/([gamma] + [[mu].sub.1]), (8)

[R.sub.02] -([beta]/[[eta].sub.1]) ([[alpha].sub.2]/([[rho].sub.1] + [[rho].sub.2] + [[mu].sub.1]) (9)

[R.sub.01] in (8) represents the product of individuals generated by infectious human and infected mosquito, respectively, during their infectious time.

[R.sub.02] in (9) represents the product of individuals generated by pseudorecovered human and infected mosquito, respectively, during their infectious time.

For model (1), the disease dies out if [R.sub.0] < 1 and the disease persists if [R.sub.0] > 1. Hence [R.sub.0] < 1, iff [R.sub.01] < 1 and [R.sub.02] < 1.

3. Malaria Stochastic Epidemic Model

For the mosquito-human malaria dynamics, since time is continuous and the disease states are discrete, we derive the stochastic version of the deterministic model (1) using a continuous-time discrete state Galton-Watson branching process (CTDSGWbp), which is a type of stochastic process. The malaria CTDSGWbp model is a time-homogeneous process with the Markov property. The model takes into account random effects of individual birth and death processes, that is, demographic variability. A stochastic process is defined by the probabilities with which different events happen in a small time interval At. In our model there are two possible events (production and death/removal) for each population. The corresponding rates in the deterministic model are replaced in the stochastic version by the probabilities that any of these events occur in a small time interval of length [DELTA]t [25,26].

3.1. Model Formulation. Let time be continuous, t [member of] [0; [infinity]), and let [S.sub.h](t), [I.sub.h](t), [R.sub.h](t), [S.sub.m](t), and [I.sub.m](t) represent discrete random variables for the number of susceptible humans, infectious humans, recovered humans, susceptible mosquitoes and infectious mosquitoes, respectively, with finite space.

[S.sub.h] (t), [I.sub.h] (*), [R.sub.h] (*), [S.sub.m] (t), [I.sub.m] (t) [member of] {0,1, 2, 3, ..., H}, (10)

where H is positive and represents the maximum size of the populations.

If a disease emerges from one infectious group with [R.sub.0] >1 and if i infective agents are introduced into a wholly susceptible population, then the probability of a major disease outbreak is approximated by 1- (1/[R.sub.0]).sup.i] while the probability of disease extinctionis approximately (1/[R.sub.0]).sup.i] [13]. However, this result does not hold if the infection emanates from multiple infectious groups [27]. For multiple infectious groups, the stochastic thresholds depend on two factors, namely, the number of individuals in each group and the probability of disease extinction for each group. Further, the persistence of an infection into a wholly susceptible population is not guaranteed by having [R.sub.0] greater than one.

For CTDSGWbp models, the transition from one state to a new state may occur at any time t. If the process begins in a state G(0), after a random time period [tau], it transits to a new state G([tau]). The process remains in state G([tau]) for a random time t after which it moves through to the new state G(p), with p = [tau] + t [27]. This process continues throughout the model. The state transitions and rates for the stochastic model are presented in Table 1.

3.2. The Branching Process Approximation. We use branching process to analyze the malaria dynamics near the diseasefree equilibrium. Since infectious human, pseudorecovered humans, and infectious mosquitoes are the only sources of infection, we apply the multitype branching process in the three variables [I.sub.h](t), [R.sub.h](t), and [I.sub.M](t). The susceptible humans and mosquitoes are assumed to be at the disease-free state [28]. We use the multitype Galton-Watson branching process (GWbp) to determine disease invasion and extinction probabilities. More review on the GWbp branching theory process can be accessed through [18, 27, 28]. We now define the offspring pgfs for the three variables, where each offspring pgf has a general form

[mathematical expression not reproducible] (11)

where [P.sub.i]([k.sub.1], [k.sub.2], ..., [k.sub.n]) = prob([X.sub.i1] = [k.sub.1], [X.sub.i2] = [k.sub.2], ..., [X.sub.in] = [k.sub.n]) is the probability that one infected individual of type i gives birth to [k.sub.j] individuals of type j; see [29].

For one malaria infectious human, there are three possible events: infection of a mosquito, recovery of the infectious human, or death of the infectious human. The offspring pgf for infectious humans define the probabilities associated with the "birth" of secondary infectious mosquito or the "death" of the initial infectious human, given that the process started with only one infectious human; that is [I.sub.h](0) = 1, [R.sub.h](0) = 0, and [I.sub.M](0) = 0.

The offspring pgf for [I.sub.h] is given by

[mathematical expression not reproducible] (12)

where [P.sub.1] ([k.sub.1],[k.sub.2], [k.sub.3]) = prob ([X.sub.11] = [k.sub.1], [X.sub.12] = [k.sub.2], [X.sub.13] = [k.sub.3]) is the probability that one infectious human through infection produces another infectious human [k.sub.1] or a pseudorecovered human k2 or an infectious mosquito [k.sub.3].

Similarly, the offspring pgf for [R.sub.h] is given by

[mathematical expression not reproducible]. (13)

where [P.sub.2] ([k.sub.1],[k.sub.2], [k.sub.3]) = prob ([X.sub.21] = [k.sub.1], [X.sub.22] = [k.sub.2], [X.sub.3] = [k.sub.3]) is

is the probability that one pseudorecovered human through infection produces an infectious human [k.sub.1] or another pseudo-recovered human [k.sub.2] or an infectious mosquito [k.sub.3].

Lastly, the offspring pgf for [I.sub.m] is given by

[mathematical expression not reproducible] (14)

where [P.sub.3] ([k.sub.1],[k.sub.2], [k.sub.3]) = prob ([X.sub.31] = [k.sub.1], [X.sub.32] = [k.sub.2], [X.sub.33] = [k.sub.3]) is is the probability that one infectious mosquito through infection produces an infections human k1 or a pseudo recovered human [k.sub.2] or another infectious mosquito [k.sub.3].

The power to which [s.sub.j] is raised is the number of infectious individuals generated from one infectious individual. If an individual recovers or dies, then no new infections are generated, hence ([s.sup.0.sub.j]).

The offspring pgfs for [I.sub.h], [R.sub.h], and [I.sub.m] are used to calculate the expected number of offsprings produced by a single infectious human or by pseudorecovered human or infectious mosquito. They are also used to calculate the probability of disease extinction.

The specific offspring pgfs for [I.sub.h], [R.sub.h], and [I.sub.m] are defined using the rates in description of parameters for the malaria transmission model in the Notations, when the initial susceptible populations are near disease-free equilibrium, [S.sub.h](0) [approximately equal to] N and [S.sub.m](0) [approximately equal to] M.

From (12), the offspring pgf for infectious human, given [I.sub.h](0) = 1, [R.sub.h](0) = 0, and [I.sub.M](0) = 0, is given by

[mathematical expression not reproducible] (15)

For [G.sub.1], one infectious human dies or is treated with probability ([[mu].sub.i] + [gamma])/([[alpha].sub.1] + [gamma] + [[mu].sub.1]); this means an infectious human dies before infecting a susceptible mosquito. The term [gamma]/([[alpha].sub.1] + [gamma] + [[mu].sub.1]) represents the probability that the infectious human is treated and moves to the pseudorecovered class, this results in "[X.sub.11] = 0, [X.sub.12] = 01, and [X.sub.13] = 0, though there is movement of an infectious human to pseudorecovered state due to partial treatment (this is not new offspring)." The infectious human infects a mosquito with probability [[alpha].sub.1]/([[alpha].sub.1] + [gamma] + [[mu].sub.1]). This means an infectious human infects a susceptible mosquito and remains infectious, which results in "[X.sub.11] = 1, [X.sub.12] = 0, and [X.sub.13] = 1." Note the term [s.sub.1] [s.sub.3] in (15) means one infectious human generates one infectious mosquito ([s.sub.3] raised to power one) and remains infectious ([s.sub.1] raised to power one).

For one pseudorecovered human, there are four events: infection of a mosquito, relapse to infected class, successful treatment of the pseudorecovered human, or death of the recovered host. Similarly, from (13), the offspring pgf for recovered humans given [I.sub.h](0) = 0, [R.sub.h](0) = 1, and [I.sub.m](0) = 0 is given by

[mathematical expression not reproducible]. (16)

For [G.sub.2], one pseudorecovered host dies or relapses or is fully treated with probability ([mathematical expression not reproducible] or infects a mosquito with probability [[alpha].sub.2]/([[alpha].sub.2] + [[rho].sub.1] + [[rho].sub.2] + [[mu].sub.1]). This means a pseudorecovered human infects a susceptible mosquito and remains infectious, which results in "[X.sub.21] = 0, [X.sub.22] = 1, and [X.sub.23] = 1." Note the term [s.sub.2][s.sub.3] in (16) means one recovered human generates one infectious mosquito ([s.sub.3] raised to power one) and remains infectious ([s.sub.2] raised to power one).

For one infectious mosquito, there are only two events: infection of a susceptible human or death of the mosquito. From (14), the offspring pgf for infected mosquito given [I.sub.h](0) = 0, [R.sub.h](0) = 0, and [I.sub.M](0) = 1 is given by

[mathematical expression not reproducible] (17)

For [G.sub.3], one infectious mosquito dies with probability [[eta].sub.1]/([beta] + [[eta].sub.1]) or infects a human with probability [beta]/([beta] + [[eta].sub.1]). This means an infectious mosquito infects a susceptible human and remains infectious, which results in "[X.sub.31] = 1, [X.sub.32] = 0, and [X.sub.3]3 = 1." Note the term [s.sub.1][s.sub.3] in (17) means one infectious mosquito generates one infectious human ([s.sub.1] raised to power one) and remains infectious ([s.sub.3] raised to power one).

3.3. The Relationship between [R.sub.0] and the Stochastic Threshold [S.sub.0]. The offspring pgfs, evaluated at (1,1,1), gives the expectation matrix with elements [m.sub.ji]. Below are the offspring pgfs evaluated at (1,1,1).

[mathematical expression not reproducible] (18)

The expectation matrix of the offspring pgfs, evaluated at (1,1,1), is given by

[mathematical expression not reproducible] (19)

The entries [m.sub.11], [m.sub.21], and [m.sub.31] represent the expected number of infectious humans, pseudorecovered humans, and infectious mosquitoes, respectively, produced by one infectious human. Similarly, the entries [m.sub.12], [m.sub.22], and [m.sub.32] represent the expected number of infectious humans, pseudo-recovered humans, and infectious mosquitoes, respectively, produced by one pseudorecovered human. Lastly, the entries [m.sub.13], [m.sub.2]3, and [m.sub.33] represent the expected number of infectious humans, pseudorecovered humans, and infectious mosquitoes, respectively, produced by one infectious mosquito.

The matrix C = W(M - I) is the Jacobian matrix evaluated for the stochastic version.

[mathematical expression not reproducible] (20)

where [mathematical expression not reproducible] is a diagonal matrix and 1 is the identity matrix.

From (3) and (20), we show that

[mathematical expression not reproducible] (21)

The spectral radius of matrix M obtained by finding the eigenvalues of matrix M is given by

[mathematical expression not reproducible] (22)

Taking the first expression of [S.sub.0] = [rho](M) in (24), we have

[mathematical expression not reproducible] (23)

This gives the probability of malaria transmission by either infectious human or by infectious mosquito. From (23),

[mathematical expression not reproducible] (24)

The probability of disease extinction is one if p(M) < 1. Hence from (24) we have

[mathematical expression not reproducible] (25)

Expanding and simplifying the inequality, we have

[beta][[alpha].sub.1] < [[eta].sub.1] ([gamma] + [[mu].sub.1]) (26)

which reduces to

[mathematical expression not reproducible] (27)

When the infection is between infectious humans and infectious mosquitoes, (27) is true. The result in (27) agrees with the deterministic reproduction number for disease elimination. Hence we conclude that the probability of disease elimination in the CTDSGWbp model is one iff

[mathematical expression not reproducible] (28)

3.4. Deriving Probability of Disease Extinction P0 Using Branching Process Approximation. To find the probability of extinction (no outbreak), we compute the fixed points of the system ([q.sub.1], [q.sub.2], [q.sub.3]) [member of] (0,1) of the offspring pgfs for the three infectious stages; that is, we solve [mathematical expression not reproducible], and [mathematical expression not reproducible]. The solutions of these systems are (1,1,1) and ([q.sub.1], [q.sub.2], [q.sub.3]); see [19]. Equating [G.sub.1]([s.sub.1], [s.sub.2], [s.sub.3]) in (15) to [q.sub.1], then letting [s.sub.1] = [q.sub.1] and [s.sub.3] = [q.sub.3], and solving for [q.sub.1], we have

[mathematical expression not reproducible] (29)

Equating [G.sub.3]([s.sub.1], [s.sub.2], [s.sub.3]) in (16) to [q.sub.2], letting [s.sub.2] = [q.sub.2] and [s.sub.3] = [q.sub.3], and solving for [q.sub.2], we have

[mathematical expression not reproducible] (30)

Equating [G.sub.3]([s.sub.1], [s.sub.2], [s.sub.3]) in (17) to [q.sub.3], then letting [s.sub.1] = [q.sub.1] and [s.sub.3] = [q.sub.3], and solving for [q.sub.3], we have

[mathematical expression not reproducible] (31)

The expression for [q.sub.1] in (29) has a biological interpretation. Beginning from one infectious human, there is no outbreak if the infectious human recovers or dies with probability ([gamma] + [[mu].sub.1])/([][alpha].sub.1] + [gamma] + [[mu].sub.1]) or if there is no successful transmission to a susceptible mosquito with probability ([[alpha].sub.1]/([[alpha].sub.1] + [gamma] + [[mu].sub.1]))(1/[R.sub.01]). This implies that if there is successful contact, then the probability of successful transmission from infectious human to susceptible mosquito is 1- 1/[R.sub.01].

The expression for q3 in (31) has a biological interpretation. Beginning from one infectious mosquito, there is no outbreak if the infectious mosquito dies with probability [[eta].sub.1]/([beta] + [[eta].sub.1]) or if there is no successful transmission to a susceptible human with probability ([beta]/([beta] + [[eta].sub.1]))(1/[R.sub.01]).

There are some other important relationships; if disease transmission is by infectious mosquitoes as well as infection of susceptible mosquitoes by infectious humans, then

[q.sub.1] * [q.sub.3] = 1/[R.sub.01] (32)

From (32), we see that, in both mosquito and human populations, the probability of no successful transmission from infectious human to susceptible mosquito and from infected mosquito to susceptible human is 1/ [R.sub.01].

If disease transmission is by infectious mosquitoes as well as infection of susceptible mosquitoes by pseudorecovered humans, then

[q.sub.2] x [q.sub.3] = 1/[R.sub.02] (33)

From (33), we see that, in both mosquito and human populations, the probability of no successful transmission from pseudorecovered humans to susceptible mosquito and from infected mosquito to pseudorecovered humans is 1/[R.sub.02].

To compute the probability of disease extinction and of an outbreak for our malaria model, we recall that, for multiple infectious groups, the stochastic thresholds depend on two factors, namely, the number of initial individuals in each group and the probability of disease extinction for each group. Using [q.sub.1], [q.sub.2], and [q.sub.3] in (29)--(31) and assuming initial individuals for infectious humans, pseudorecovered humans and infected mosquitoes are [I.sub.h](0) = [h.sub.0], [R.sub.h](0) = [r.sub.0], and [I.sub.m](0) = [m.sub.0], respectively. Then the probability of malaria clearance is given by

[mathematical expression not reproducible] (34)

4. Numerical Simulations

In this section, we illustrate numerically the disease dynamics of model (1) using parameter values in Table 2. The numerical simulations are done using Maple codes.

4.1. Effects of Model Parameters on [R.sub.0]. Using parameter values in Table 2, we identified how different input parameters affect the reproduction number [R.sub.0] of our model as shown in Figure 3.

From Figure 3, we see that increase in [gamma], [[rho].sub.1], [[mu].sub.1], and [q.sub.1] will decrease [R.sub.0]. Also decrease in [beta], [[alpha].sub.1], [[alpha].sub.2], and [[rho].sub.2] will decrease [R.sub.0].

From the Tornado plot, the infection of susceptible humans by infected mosquitoes (denoted by [beta]) is a major factor in the malaria transmission dynamics. Reducing [beta[ would reduce [R.sub.0] significantly hence reducing the possibility of disease outbreak. Vector control is the main way to prevent and reduce malaria transmission. If coverage of vector control interventions within a specific area is high enough, then a measure of protection will be conferred across the community. WHO recommends protection for all people at risk of malaria with effective malaria vector control which agrees with the information from Figure 3. Two forms of vector control insecticide-treated mosquito nets and indoor residual spraying are effective in a wide range of circumstances. Another factor which affects the malaria transmission significantly is the death rate of infected mosquitoes (denoted by [[eta]/sub.1]). Increasing [[eta]/sub.1] would decrease [R.sub.0] hence reducing the chances of disease outbreak.

4.1.1. Graphical Representation of Parameters Effects on [R.sub.0]. Figure 4 shows the effects of relapse and recover rates on the basic reproduction number. From the simulations, we find that [R.sub.0] is increasing with increase in relapse rate, while it is decreasing with increase in recovery rate. To control and eradicate the malaria epidemic, it is important and necessary for governments of endemic countries to decrease the relapse rate and increase the treatment and recovery rate.

4.2. Probability of Disease Extinction. Using parameter values in Table 2, we compute numerically the probability of disease extinction [P.sub.0] and of an outbreak 1- [P.sub.0] for our malaria model using different initial values for the infectious classes.

The probability of disease extinction is high if the disease emerges from infected humans. It is very low if the infection emerges from infected mosquitoes. However, as the initial number of infected humans grows largely, there is a high probability of disease outbreak as illustrated in Table 3. The probability of disease extinction is significantly low if the disease emerges from infected mosquitoes. Therefore, the disease dynamics for model system (1) at the beginning of the epidemics are being driven by initial number of infected mosquitoes.

A female mosquito will continue to bite and draw blood until her abdomen is full. If she is interrupted before she is full, she will fly to the next person. After feeding, the mosquito rests for two or three days before laying her eggs, and then it is ready to bite again. Since mosquitoes do not recover from the infection, then a single infected mosquito is capable of biting and infecting so many susceptible humans in their lifespan which may in turn infect so many uninfected mosquitoes there by reducing the probability of disease clearance and increasing the probability of a major disease outbreak. Moreover, mosquitoes are the reservoir host for the parasites that cause malaria and it takes a long time for them to be malaria-free, hence the high probability of malaria outbreak if the parasite is introduced by infected mosquito.

Table 3 depicts that, at the beginning of malaria outbreak, any policy or intervention to control the spread of malaria should focus on controlling the infected mosquito population as well as the infected humans. If more effort to control the disease is only focused on the infected humans, then it is very unlikely that malaria will be eliminated. This is an interesting insight from the stochastic threshold that could not be provided by the deterministic threshold.

4.3. Numerical Simulation of Malaria Model. Using parameter values in Table 2, we numerically simulate the behavior of model (1). Initial conditions are [S.sub.h](0) = 99, [I.sub.h](0) = 1, [R.sub.h](0) = 0, [S.sub.m](0) = 999, and [I.sub.m](0) = 1.

4.3.1. Numerical Simulation When [R.sub.0] < 1. From Figure 5, the analysis shows that when [R.sub.0] < 1 and [S.sub.0] < 1, then the probability of disease extinction is [P.sub.0] = 0.9476 [congruent to] 1, although this agrees with (30), which points out that the probability of disease elimination in the CTDSGWbp model is one iff

[mathematical expression not reproducible] (35)

There is still a small probability 1 - [P.sub.0] = 0.0524 of disease outbreak. [R.sub.0] has been widely used as a measure of disease dynamics to estimate the effectiveness of control measures and to inform on disease management policy. However, from the analysis in Figure 5, it is evident that R0 can be flawed and disease can persist with [R.sub.0] <1 depending on the kind of disease being modeled.

4.3.2. Numerical Simulation When [R.sub.0] > 1. Using parameter values in Table 2, we numerically simulate the behavior of model (1) when [R.sub.0] > 1. Initial conditions are [S.sub.h](0) = 99, [I.sub.h](0) = 1, [R.sub.h](0) = 0, [S.sub.m](0) = 999, and IJ0) = 1.

The analysis from Figure 6 suggests that when [R.sub.0] >1, there is still some probability of disease extinction ([P.sub.0] = 0.3713). Although the probability is low, there is still a chance to clear the disease. Therefore, from the analysis in Figure 6, disease can be eliminated with [R.sub.0] >1 and hence the R0 threshold should not be the only parameter to consider in quantifying the spread of a disease.

4.3.3. Effects of Relapse on Human Population. Figure 7 shows changing effects of relapse on the human populations.

To control and eradicate the malaria epidemic, it is important and necessary for governments of endemic countries to decrease the relapse rate as can be seen in Figure 7.

5. Discussions and Recommendations

In this study, we investigated the transmission dynamics of malaria using CTDSGWbp model. The disease dynamic extinction thresholds from the stochastic model were compared with the corresponding deterministic threshold. We derived the stochastic threshold for disease extinction [S.sub.0] and showed the relationship that exists between [R.sub.0] and [S.sub.0] in terms of disease extinction and outbreak in both deterministic and stochastic models.

Our analytical and numerical results showed that both deterministic and stochastic models predict disease extinction when [R.sub.0] <1 and [S.sub.0] <1. However, the predictions by these models are different when [R.sub.0] > 1. In this case, deterministic model predicts with certainty disease outbreak while the stochastic model has a probability of disease extinction at the beginning of an infection. Hence, with stochastic models, it is possible to attain a disease-free equilibrium even when [R.sub.0] >1. Also we noticed that initial conditions do not affect the deterministic threshold while the stochastic thresholds are affected. Thus, the dynamics of the stochastic model are highly dependent on the initial conditions and should not be ignored.

The probabilities of disease extinction for different initial sizes of infected humans and infected mosquitoes were approximated numerically. The results indicate that the probability of eliminating malaria is high if the disease emerges from infected human as opposed to when it emerges from infected mosquito at the beginning of the disease. The analysis has shown that any policy or intervention to control the spread of malaria at the beginning of an outbreak must focus not only on infected humans but also on controlling the infected mosquito population. Also our results strongly suggest that, to eradicate malaria, governments of endemic countries should increase the recovery rate and reduce the relapse rate.

In conclusion, to achieve the WHO vision 2030 strategic global targets, which include (i) reducing malaria case incidence by at least 90% by 2030 and (ii) reducing malaria mortality rates by at least 90% by 2030, the governments of endemic countries should embark on increasing the coverage of vector control interventions and reduce the malaria relapse rate as well as controlling the infected mosquito population.

For future work, the study can be extended by considering vertical transmission of the mosquitoes. Another extension of the study can be the inclusion of immigration and emigration of individuals in both populations in order to investigate the effect of movement on the disease transmission dynamics on the two thresholds. Also the study can be extended by considering the effect of climatic conditions, such as rainfall patterns, temperature, and humidity on the disease transmission dynamics.

Notations

Description of Variables for the Malaria Transmission Model

[S.sub.h]: Susceptible humans [I.sub.h] : Infectious humans [R.sub.h]: Pseudorecovered humans [S.sub.m]: Susceptible mosquito [I.sub.m]: Infectious mosquito

Description of Parameters for the Malaria Transmission Model

[lambda]: Natural birth rate of humans [mu]: Natural death rate of humans [[mu].sub.1]: Disease induced death rate of humans [beta]: Transmission rate from an infectious mosquito to a susceptible human [[rho].sub.1]: Recovery rate [[rho].sub.2]: Relapse rate [gamma]: Treatment rate [delta]: Natural birth rate of mosquitoes [eta]: Natural death rate of mosquitoes [[eta].sub.1]: Parasite induced death rate of mosquitoes [[alpha].sub.1]: Transmission rate from an infectious human to a susceptible mosquito [[alpha].sub.2]: Transmission rate from a pseudorecovered human to a susceptible mosquito N: The total size of human population M: The total size of mosquito population.

https://doi.org/10.1155/2018/2439520

Conflicts of Interest

All authors declared no conflicts of interest regarding the publication of this paper.

Authors' Contributions

Rachel Waema Mbogo produced the first draft of the manuscript. All authors critically reviewed the paper and approved the final version.

References

[1] WHO, World Malaria, World Health Organization, Geneva, Switzerland, 2012.

[2] WHO, World Malaria, World Health Organization, Geneva, Switzerland, 2015.

[3] R. Ross, "An application of the theory of probabilities to the study of a priori pathometry," in Proceedings of the Royal Society A, vol. 92, pp. 204-230,1916.

[4] Y. Xiao and X. Zou, "Can multiple malaria species co-persist?" SIAM Journal on Applied Mathematics, vol. 73, no. 1, pp. 351-373, 2013.

[5] J. Li, Y. Zhao, and S. Li, "Fast and slow dynamics of malaria model with relapse," Mathematical Biosciences, vol. 246, no. 1, pp. 94-104, 2013.

[6] H.-F. Huo and G.-M. Qiu, "Stability of a mathematical model of malaria transmission with relapse," Abstract and Applied Analysis, vol. 2014, Article ID 289349, 9 pages, 2014.

[7] H.-F. Huo and C.-C. Zhu, "Influence of relapse in a giving up smoking model," Abstract and Applied Analysis, vol. 2013, Article ID 525461, 12 pages, 2013.

[8] N. Chitnis, J. M. Cushing, and J. M. Hyman, "Bifurcation analysis of a mathematical model for malaria transmission," SIAM Journal on Applied Mathematics, vol. 67, no. 1, pp. 24-45, 2006.

[9] L. J. S. Allen, "Stochastic population and epidemic models:persistence and extinction," in Mathematical Biosciences Institute lecture series, Stochastic in Biological systems, Springer International Publishing, Switzerland, 2015.

[10] T. Britton, "Stochastic epidemic models: A survey," Mathematical Biosciences, vol. 225, no. 1, pp. 24-35, 2010.

[11] Y. Cai, Y. Kang, M. Banerjee, and W. Wang, "A stochastic epidemic model incorporating media coverage," Communications in Mathematical Sciences, vol. 14, no. 4, pp. 893-910, 2016.

[12] Y. Cai, Y. Kang, and W. Wang, "A stochastic SIRS epidemic model with nonlinear incidence rate," Applied Mathematics and Computation, vol. 305, pp. 221-240, 2017

[13] M. Maliyoni, F. Chirove, H.. Gaff, and K. S. Govinder, "A stochastic tick-borne disease model: exploring the probability of pathogen persistence," Bulletin of Mathematical Biology, vol. 79, no. 9, pp. 1999-2021, 2017

[14] M. S. Bartlett, "The relevance of stochastic models for large-scale epidemiological phenomena," Journal of Applied Statistics, vol. 13, no. 13, pp. 2-8, 1964.

[15] M. S. Bartlett, Stochastic population models, Methuen, London, UK, 1960.

[16] L. J. S. Allen, "An Introduction to stochastic Epidemic models," in Mathematical Epidemiology, Springer, Berlin, Germany, 2008.

[17] L. J. S. Allen and P. Van den Driessche, "The basic reproduction number in some discrete-time epidemic models," Journal of Difference Equations and Applications, vol. 14, pp. 11-27, 2008.

[18] L. J. Allen and P. van den Driessche, "Relations between deterministic and stochastic thresholds for disease extinction in continuous- and discrete-time infectious disease models," Mathematical Biosciences, vol. 243, no. 1, pp. 99-108, 2013.

[19] L. J. S. Allen, "A primer on stochastic epidemic models: Formulation, numerical simulation, and analysis," Infectious Disease Modelling, vol. 243, pp. 1-15, 2017

[20] C. Manore and M. Hyman, "Mathematical models for fighting zika virus," SIAM News, vol. 2016, pp. 1-5, 2016.

[21] S. Olaniyi, M. A. Lawal, and O. S. Obabiyi, "Stability and sensitivity analysis of a deterministic epidemiological model with pseudo-recovery," IAENG International Journal of Applied Mathematics, vol. 46, no. 2, pp. 1-8, 2016.

[22] L. Xue and C. Scoglio, "Network-level reproduction number and extinction threshold for vector-borne diseases," Mathematical Biosciences and Engineering, vol. 12, no. 3, pp. 565-584,2015.

[23] P. van den Driessche and J. Watmough, "Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission," Mathematical Biosciences, vol. 180, pp. 29-48, 2002.

[24] O. Diekmann, J. A. Heesterbeek, and J. A. Metz, "On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations," Journal of Mathematical Biology, vol. 28, no. 4, pp. 365-382, 1990.

[25] W. R. Mbogo, L. S. Luboobi, and J. W. Odhiambo, "Stochastic model for in-host HIV dynamics with therapeutic intervention," ISRN Biomathematics, vol. 2013, Article ID 103708, 11 pages, 2013.

[26] W. R. Mbogo, L. S. Luboobi, and J. W. Odhiambo, "Mathematical model for HIV and CD4+ cells dynamics in vivo," International Journal of Pure and Applied Mathematics, vol. 6, no. 2, pp. 83-103, 2013.

[27] L. J. S. Allen and G. E. Lahodny Jr., "Extinction thresholds in deterministic and stochastic epidemic models," Journal of Biological Dynamics, vol. 6, no. 2, pp. 590-611, 2012.

[28] J. Lahodny, R. Gautam, and R. Ivanek, "Estimating the probability of an extinction or major outbreak for an environmentally transmitted infectious disease," Journal ofBiological Dynamics, vol. 9, no. 1, pp. 128-155, 2015.

[29] L. J. S. Allen, "Branching processes," in Encyclopaedia of Theoretical Ecology, University of California Press, Calif, USA, 2012.

Rachel Waema Mbogo [ID], Livingstone S. Luboobi, John W. Odhiambo

Institute of Mathematical Sciences (IMS), Strathmore University, Box 59857 00200, Nairobi, Kenya

Correspondence should be addressed to Rachel Waema Mbogo; rmbogo@strathmore.edu

Received 14 September 2017; Accepted 26 December 2017; Published 11 February 2018

Academic Editor: Zhidong Teng

Caption: Figure 1: The mosquito-human transmission cycle. An infectious mosquito bites a susceptible, uninfected human and transmits the virus via saliva. Once the human becomes infectious (usually accompanied by symptoms), the human can transmit the pathogen to an uninfected mosquito via the blood the mosquito ingests. Source. Figure 1 is reproduced from Manore and Hyman (2016) [20] [under the Creative Commons Attribution License/public domain].

Caption: Figure 2: Schematic representation of the mosquito-human malaria transmission dynamics.

Caption: Figure 3: Tornado plot for parameters that influence [R.sub.0].

Caption: Figure 4: Effects of relapse and recovery rate on [R.sub.0].

Caption: Figure 5: Malaria dynamics when [R.sub.0] = 0.4564 and [S.sub.0] = 0.5661 and [P.sub.0] = 0.9476.

Caption: Figure 6: Malaria dynamics when [R.sub.0] = 2.4629 and [P.sub.0] = 0.3713.

Caption: Figure 7: Changing effects of relapse on the human populations.

Table 1: State transitions and rates for the CTDSGWbp malaria model. Event Population Population components at t t + [DELTA]t Birth ofhumans ([S.sub.h], ([S.sub.h] + [I.sub.h], 1,[I.sub.h], [R.sub.h], [R.sub.h], [S.sub.m], [S.sub.m], [I.sub.m]) [I.sub.m]) Death of susceptible ([S.sub.h], 4, ([S.sub.h]-1, humans [R.sub.h], [I.sub.h], [S.sub.m], [R.sub.h], [I.sub.m]) [S.sub.m], [I.sub.m]) Infection of humans ([S.sub.h], ([S.sub.h] [I.sub.h], 1,[I.sub.h] + [R.sub.h], 1, [R.sub.h], [S.sub.m], [S.sub.m], [I.sub.m]) [I.sub.m]) Recovery rate of ([S.sub.h], ([S.sub.h] + humans [I.sub.h], 1,[I.sub.h], [R.sub.h], [R.sub.h] -1, [S.sub.m], [S.sub.m],[I.sub.m]) [I.sub.m]) Relapse rate of ([S.sub.h], ([S.sub.h], humans [I.sub.h], [I.sub.h] +1, [R.sub.h], [R.sub.h] -1, [S.sub.m], [S.sub.m], [I.sub.m]) [I.sub.m]) Treatment rate ([S.sub.h], ([S.sub.h], [I.sub.h], [I.sub.h] -1 [R.sub.h], ,[R.sub.h] + 1, [S.sub.m], [S.sub.m], [I.sub.m]) [I.sub.m]) Death of infected ([S.sub.h], ([S.sub.h], humans [I.sub.h], Rh, [I.sub.h]-1, [I.sub.m]) [R.sub.h], [S.sub.m], [I.sub.m]) Death of recovered ([S.sub.h], ([S.sub.h], humans [I.sub.h], Rh, [I.sub.h], [S.sub.m], [R.sub.h]-1, [I.sub.m]) [S.sub.m], [I.sub.m]) Birth of mosquitoes ([S.sub.h], ([S.sub.h], [I.sub.h], Rh, Th[R.sub.h], [I.sub.m]) [S.sub.m] +1, [I.sub.m]) Death of susceptible ([S.sub.h], ([S.sub.h], mosquitoes [I.sub.h], Rh, [I.sub.h], [S.sub.m], [R.sub.h], [I.sub.m]) [S.sub.m][I.sub.m]) Infection of ([S.sub.h], ([S.sub.h], mosquitoes from [I.sub.h], Rh, [I.sub.h], humans [S.sub.m], [R.sub.h], [S.sub.m] [I.sub.m]) -[I.sub.m] + 1) Death of infected ([S.sub.h], ([S.sub.h], mosquitoes [I.sub.h], Rh, [I.sub.h], [S.sub.m], [R.sub.h], [I.sub.m]) [S.sub.m], [I.sub.m]-1 Event Transition probabilities Birth ofhumans [lambda] [DELTA]t Death of susceptible [mu][S.sub.h][DELTA]t humans Infection of humans [beta]([S.sub.h]/ N)[I.sub.m] [DELTA]t Recovery rate of [[rho].sub.1][R.sub.h] humans [DELTA]t Relapse rate of [[rho].sub.2] humans [R.sub.h][DELTA]t Treatment rate [gamma] [l.sub.h] [DELTA]t Death of infected [[mu].sub.1][I.sub.h] humans [DELTA]t Death of recovered [[mu].sub.1][R.sub.h] humans [DELTA]t Birth of mosquitoes [delta] [DELTA]t Death of susceptible [eta][S.sub.m][DELTA]t mosquitoes Infection of ([[alpha].sub.1][I.sub.h] mosquitoes from + [[alpha].sub.2] humans [R.sub.h])([S.sub.m]/ N)[DELTA]t Death of infected [[eta].sub.1][I.sub.m] mosquitoes [DELTA]T Table 2: Model parameter values. Parameter Description Units [lambda] Birth rate ofhumans Per day [mu] Death rate of susceptible humans Per day [beta] Infection rate ofhumans Per day [[mu].sub.1] Death rate of infected humans Per day [[rho].sub.1] Recovery rate ofhumans Per day [[rho].sub.2] Relapse rate ofhumans Per day [gamma] Treatment rate Per day [[mu].sub.1] Death rate of recovered humans Per day [delta] Birth rate of mosquitoes Per day [eta] Death rate of susceptible mosquitoes Per day [[alpha].sub.1] Infection rate from infectious human Per day [[alpha].sub.2] Infection rate from recovered human Per day [[eta].sub.1] Date rate of infected mosquitoes Per day Parameter Parameter Source values [lambda] 0.000039 [5] [mu] 0.000039 [5] [beta] 0.02 [19] [[mu].sub.1] 0.00039 [5] [[rho].sub.1] 0.01 [19] [[rho].sub.2] 0.002 [5] [gamma] 0.037 [5] [[mu].sub.1] 0.00034 [5] [delta] 0.143 [5] [eta] 0.143 [5] [[alpha].sub.1] 0.072 [5] [[alpha].sub.2] 0.0072 [5] [[eta].sub.1] 0.143 [19] Table 3: Probability of disease extinction [P.sub.0] and of an outbreak 1-[P.sub.0] for the malaria model. [h.sub.0] [r.sub.0] [m.sub.0] [P.sub.0 1-[P.sub.0] 1 0 0 0.6838 0.3162 0 1 0 0.9834 0.0166 0 0 1 0.1365 0.8635 1 1 0 0.6724 0.3276 1 1 1 0.0918 0.9082 10 0 0 0.0223 0.9777 0 10 0 0.8459 0.1541 0 0 10 0 1 10 10 0 0.0188 0.9812 10 10 10 0 1 [h.sub.0] [R.sub.01] [R.sub.02] [R.sub.0] 1 [R.sub.01] > 1 [R.sub.02] < 1 [R.sub.0] > 1 0 [R.sub.01] > 1 [R.sub.02] < 1 [R.sub.0] > 1 0 [R.sub.01] > 1 [R.sub.02] < 1 [R.sub.0] > 1 1 [R.sub.01] > 1 [R.sub.02] < 1 [R.sub.0] > 1 1 [R.sub.01] > 1 [R.sub.02] < 1 [R.sub.0] > 1 10 [R.sub.01] > 1 [R.sub.02] < 1 [R.sub.0] > 1 0 [R.sub.01] > 1 [R.sub.02] < 1 [R.sub.0] > 1 0 [R.sub.01] > 1 [R.sub.02] < 1 [R.sub.0] > 1 10 [R.sub.01] > 1 [R.sub.02] < 1 [R.sub.0] > 1 10 [R.sub.01] > 1 [R.sub.02] < 1 [R.sub.0] > 1

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Mbogo, Rachel Waema; Luboobi, Livingstone S.; Odhiambo, John W. |

Publication: | Journal of Applied Mathematics |

Article Type: | Report |

Date: | Jan 1, 2018 |

Words: | 8079 |

Previous Article: | Accessing the Power of Tests Based on Set-Indexed Partial Sums of Multivariate Regression Residuals. |

Next Article: | A Comparison of Algorithms for Finding an Efficient Theme Park Tour. |

Topics: |