# Extinction of disease pathogenesis in infected population and its subsequent recovery: a stochastic approach.

1. IntroductionThe modeling of epidemic diseases is an age-old problem [1] and makes a sincere effort in understanding the development of mathematical models for epidemics from the 18th century to the present day. These models are shown to be of used in predicting and controlling the spread of infection. But in the recent context of epidemiological research, microbial pathogenesis reflects the interaction between two entities, host and pathogen, which is somewhat related to the predator-prey model [2, 3]. These classes of models are also relevant for host-parasite type of models. Host-pathogen models are mathematical prototypes pertaining to epidemiology, and are of immense importance in view of the emergence and reemergence of epidemiological diseases in the present day global scenario [4, 5].

Application of mathematical concepts and techniques to analyze host-pathogen interactions was done by various researchers [6-8]. A very recent study throws light on the antiviral drug treatment along with immune activator IL-2 and optimal control in disease pathogenesis by deterministic model formulation [9]. In this research paper, a conventional host-pathogen model has been considered including the recovery of the infected individuals to the healthy susceptible organisms. In this case, the host population is divided into two classes: susceptible (S), that is, healthy organisms, and infected individuals (I).Pathogens (V) cause infection to host population by transforming S to I. Over the last several years, many researchers focused their attention on the mathematical and biological aspects of host-pathogen interactions. Beltrami and Carroll [10] as well as Venturino, worked on the role of viral disease considering a three-species model of susceptible and infected phytoplankton as well as their predator [4]. Hethcote and Driessche [11] formulated SIS epidemiological models where delay has been incorporated corresponding to the infectious period and disease-related deaths. Another pioneering work was reported by [12] where the transport of coevolved host-pathogen systems into new environment leads to the evolution of altered levels of pathogen aggressiveness, if transmission rates are different in the new environment. With the growing research in the prey-predator and other prototypical systems, it is apparent that the pathogen or viral growth through replication influences the model dynamics. This has been emphasized by Bairagi et al. in a subsequent communication [13]. In another pioneering communication in recent years, Bairagi et al. carried a comparative study of the prey-predator model with several response functions [14].

Epidemiological modeling uses both deterministic and stochastic approaches in host-pathogen interactions [15]. Both model types have their respective advantages and weaknesses. Deterministic models may be considered as an approximation of a corresponding stochastic model. An important difference between the two is that the stochastic model deals with a finite population size, while the deterministic model deals only with proportion [16]. The deterministic version of the model can be derived as an approximation of the stochastic version using the state variables which are susceptible host (S), infected host (7), and pathogen (V) where the total population comprising the susceptible and infected host is N. Recurrence can be explained by the combined influence of epidemic and demographic forces. In stochastic models, infection will eventually become extinct, and time to extinction is an important measure of the persistence of the infection [17].

This paper specifies a stochastic approach of our earlier work [18] designed for some widespread infection in closed population. It is a well-known fact that beyond a threshold value, that is, basic reproduction ratio, which is determined by the parameters of the model, the deterministic model predicts that the proportion of infected individuals will approach a positive endemic level as time approaches infinity. However, the stochastic model predicts that the infection will become extinct. Generally, stochastic version of the model is compared with prey-predator interaction or more precisely the host-pathogen interaction, and thus the concept of the quasi-stationary state is an important aspect of this model [19]. It should be mentioned here that the stochastic logistic model can be interpreted as SIS model and that SIS model is used for infection that gives no immunity. In this case, it is assumed that an individual who recovers from an infection remains susceptible to further infection. The time to extinction from quasi-stationary distribution is a measure of the persistence of the infection. The stationary distribution is to be an additional importance since the expected time to extinction from quasi-stationary can be expressed in terms of this distribution. A major goal of the analysis of the stochastic model is therefore to derive an approximation of the quasi-stationary distribution. This derivation is based on a diffusion approximation of the stochastic discrete state model.

The model is analyzed in two different avenues, analytical and numerical. In our research paper, we have found an approximation for the marginal distribution of the infected population in quasi-stationary condition and its time to extinction. The time to extinction of the infected population depends on the total population, and the study has been carried out through stochastic approach so that infected population is to be extinct and increase the susceptible population. However, being an extension of our previous work [18], where deterministic model was used and optimal control theory applied to predict the decrease of infected population along with cost-effective combination of treatment, our present work emphasizes on the time to extinction of the infected host population or the disease by stochastic model formulation and its subsequent analysis and evaluation.

2. The Deterministic Model

We consider the three components of the basic three-dimensional deterministic host pathogen model [18] consisting of a host population, whose concentration is denoted by N ([N] = number of host per designated area) and a pathogen population inflicting infection in the host population whose concentration is denoted by V ([V] = number of pathogens per designated area). The following differential equations are formed at initial dynamics on the change of host pathogen interaction with time t:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (1)

Here in the presence of pathogenic infection, the host population is divided into two disjoint classes, susceptible host S and infected host I. In the ideal case of no pathogen, the growth of susceptible host population follows the logistic law [19] implying that this growth is entirely controlled by an intrinsic birth rate constant r([member of] [R.sub.+]) with a carrying capacity K([member of] [R.sub.+]). [gamma]([member of] [R.sub.+]) is the force of infection through contact with pathogens, and pathogens maximally infect [gamma] susceptible hosts per day. This infection rate is half maximal at susceptible host population density of [h.sub.[gamma]] host. [lambda]([member of] [R.sub.+]) is the intensity of infection by infected host, and [d.sub.I]([member of] [R.sub.+]) is the death rate constant of I. Rate of cell lysis (replication of pathogens) is [eta]([member of] [R.sub.+]), and the natural death rate of pathogens is denoted as [mu]([member of] [R.sub.+]). We assume that the infected hosts do not grow or reproduce, but they can recover from pathogenic infection and move to the susceptible host population. Such recovery would stem out from immunization or vaccination. We consider a recovery rate of infected host (I) to be denoted by [delta]([member of] [R.sub.+]).

In spirit of Bonhoefer et al., we employed a simplified system with two components, the susceptible and the infected host. It is assumed that at equilibrium point [??] = 0, so V can be eliminated by putting V = [eta][d.sub.I]([h.sub.[gamma]] + S) I/([mu]([h.sub.[gamma]] + S) + [gamma]S). With this choice of V model, (1) reduces to

ds/dt = rS(1 - S + I/K) - [lambda]SI - [alpha] SI/[beta] + S + [delta]I, dI/dt = [lambda]SI + [alpha] SI/[beta] + S - [d.sub.I]I - [delta]I, (2)

where [gamma][eta][d.sub.I]/([gamma] + [mu]) = [alpha] and [h.sub.[gamma]][mu]/([gamma] + [mu]) = [beta].

An alternative deterministic formulation of the reduced model (2), one assumes birth and death rate functions B(S) and D(S), respectively, of

B (S) = [b.sub.1]S - [b.sub.2](S + I) S, D (S) = [d.sub.1]S + [d.sub.2](S + I) S, (3)

where [b.sub.1], [b.sub.2] and [d.sub.1], [d.sub.2] > 0. Here, [b.sub.1] and [d.sub.1] are intrinsic rate, [b.sub.2] and [d.sub.2] are crowding coefficient [20], and thus intrinsic growth rate (r) and carrying capacity (K) are defined as follows: r = [b.sub.1] - [d.sub.1] and K = ([b.sub.1] - [d.sub.1])/([b.sub.2] + [d.sub.2]) [20],

dS/dt = (B(S) - D(S)) - [lambda]SI - [alpha] SI/[beta] + S + [delta]I, dI/dt = [lambda]SI + [alpha] SI/[beta] + S - ([d.sub.I] + [delta]) I. (4)

3. The Stochastic Model Formulation

There are two state variables, namely, the number of susceptible hosts S(t) and the number of infected hosts I(t) at time t. They jointly take values in the state space [S.sub.p] = {(m, n): m = 0, 1, 2, n = 0, 1, 2, ...}. The joint distribution of S(t) and I(t) at time t is denoted by

[p.sub.mn] (t) = P{S(t) = m, I(t) = n}. (5)

We use this notation even when m and/or n are negative, with the convention that [p.sub.mn](t) is equal to zero. The model is based on the following four basic events, that is, birth of a susceptible host, death of a susceptible host, infection of an uninfected host, and death or recovery of an infected host. The transition rates of the model are shown in Table 1.

4. Description of the Transition States

The total number of population N is increased by unity, if there is a birth of susceptible host (S) for a small time interval [DELTA]t. But to make the population to be the same, we should assume that there must be a natural death of susceptible host (S). These phenomena are captured through the first two rows of the transition matrix. On the other hand, if there is an infection in the susceptible host it can be balanced by an increase of an infected host, (I). The susceptible class is infected either by direct infection with the pathogen or by replication of viral generated within the infected population.

This occurs one at a time, and so the increases of the infected class are reflected by the rise of unity in the transition state. If there is an infected host, the natural death should be reflected through a natural birth of susceptible host. At the end, naturally, the recovery of infected host must be reinstated to the susceptible host.

5. Formulation of Kolmogorov's Forward Equation

We are supposing that in a time interval of infinitesimally little length ([DELTA]t), the probability of precisely one birth (or one death) is birth rate (or death rate) x ([DELTA]t) + intuitively more than one occasion (birth and/or death) in o([DELTA]t). We as well believe the prospect [p.sub.mn](t + [DELTA]t), where [DELTA]t [down arrow] 0. The Kolmogorov's forward equations for the representation can be written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (6)

where

[kappa](m, n) = [[lambda].sub.1](m, n) + [[mu].sub.1] (m, n) + [[mu].sub.2] (m, n) + [[gamma].sub.2] (m, n). (7)

Note that all proceedings consisting of more than one birth or more than one death are incorporated in the o([DELTA]t) expression

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (8)

6. The Quasi-Stationary Distribution

We obtain primarily a deferential equation that will be used afterward in this subsection and in the discussion of time to extinction in the subsequent subsection. Place n = 0 in Kolmogorov's Forward equations (8), and put in the explicit expressions for the evolution rates given above. In addition, commence

[p.sub..1](t) = [[infinity].summation over (m = 0)][p.sub.mn] = P{I(t) = n} (9)

to indicate the subsidiary allocation of the number of infected individuals at time t. By summing the forward equations intended for n = 0 over all m-values, we acquire

[p'.sub..0](t) = ([delta] + [d.sub.I])[p.sub..1](t). (10)

After that, we develop an organization of equations for the quasi-stationary distribution. The state possibility conditioned on not being engrossed is signified [q.sub.mn](t) and specified by

[q.sub.mn](t) = P{S(t) = m, I(t) = n|I(t) [not equal to] 0} = [p.sub.mn](t)/1 - [p.sub..0](t), m = 0, 1, 2, ..., n = 1, 2, ... (11)

We initiate [q.sub..n](t) = [[summation].sup.[infinity].sub.m = 0] [p.sub.mn](t) to stand for the marginal distribution for the number of infected individuals at time t, habituated on not having attained any condition in the engrossing position. By differentiating the expression for [q.sub.mn](t) in (11) and relating (10), we get hold of

[q'.sub.mn](t) = [[p'.sub.mn](t)/1 - [p.sub..0](t)] + ([delta] + [d.sub.I])[q.sub..1](t) [[p.sub.mn](t)/1 - [p.sub..0](t)]. (12)

By pertaining the forward equations for [q.sub.mn](t) in (8), we achieve the following scheme of differential equations for the conditional state probabilities [q.sub.mn](t):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (13)

The quasi-stationary distribution [q.sub.mn](t) is the stationary solution of this system of equations.

7. The Distribution of the Time to Extinction

Two initial distributions are principally fascinating. One is the quasi-stationary distribution, and another communicates to one infected individual. The distribution of the time to extinction [tau] can be unwavering if we can resolve the Kolmogorov forward equations (8) for [p.sub.mn](t). The reason is that the happening that [tau] surpasses t is equal to the event that the number of infected individuals at time t is positive. By allowing for the balancing events we attain

P{[tau] [less than or equal to] t} = P{I(t) = 0} = [p.sub..0] (t). (14)

Therefore, the cumulative distribution function of the time to extinction at time t equals the marginal possibility in which the number of infected individuals at time t equals 0. The distribution of the time to extinction is particularly straightforward when the opening distribution is equal to the quasi-stationary distribution. Let us indicate the time to extinction from quasi-stationarity by [[tau].sub.Q] explicitly. We demonstrate that [[tau].sub.Q] has an exponential distribution and that its predictable value is equal to

E([[tau].sub.Q]) = [1/[delta] + [d.sub.I]] [q.sub..1](t). (15)

To gain this consequence, we put [q'.sub.mn](t) = 0 in (12). Thus, we are guided to the initial value problems

[p'.sub.mn] = -([delta] + [d.sub.I])[q.sub..1][p.sub.mn], [p.sub.mn](0) = [q.sub.mn], m = 0, 1, 2,..., n = 1, 2, 3... (16)

with solutions

[p.sub.mn] = [q.sub.mn] exp (-([delta] + [d.sub.I])[q.sub..1]t), m = 0, 1, 2..., n = 1, 2, 3.... (17)

By adding these expressions of [p.sub.mn](t) over all m, we obtain

[p.sub..n] = [q.sub..n] exp (-([delta] + [d.sub.I])[q.sub..1]t), n = 1, 2, 3.... (18)

The differential equation for [p.sub..0] (t) in (10) can now be answered since the right-hand side of this equation is known from above. Remembering that we have the initial value [p.sub..0](0) = 0, we acquire [p.sub..0](t) = 1 - exp (-([delta] + [d.sub.I])[q.sub..1]t). This institutes the claim that [[tau].sub.Q] has an exponential distribution with expected value agreed by (15).

8. Diffusion Approximation and the Approximation of Quasi-Stationary Distribution

In this section, we have derived the diffusion approximation for the process formulated in Section 2. In order to approximate the quasi-stationary distribution, we consider the two-dimensional process which is represented by the set of differential equation (4)

dS/dt = (B(S) - D(S)) - [lambda]SI - [alpha] SI/[beta] + S + [delta]I, dI/dt = [lambda]SI + [alpha] SI/[beta] + S - [d.sub.I]I - [delta]I. (19)

The main result is the quasi-stationary distribution is approximated by a bivariate normal distribution, if N is sufficiently large.

The critical point of the rescaled deterministic model that corresponds to a pathogenic infection is denoted by x = [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where

[??] = [([d.sub.I] + [delta] - [alpha] -[lambda][beta]) + [square root of ([([d.sub.I] + [delta] - [alpha] - [lambda][beta]).sup.2] + 4[beta][lambda] ([delta] + [d.sub.I]))]/2[lambda]], (20)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (21)

The changes in the scaled state variables S and I during the time interval from t to t + [DELTA]t are denoted by [DELTA]S and [DELTA]I using Table 2: [DELTA]S = S(t + [DELTA]t) - S(t) and [DELTA]I = I(t + [DELTA]t) - I(t). From the hypotheses of the original process, we determine the mean and the covariance of the vector with components [DELTA] and [DELTA]I. We begin with the mean

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (22)

The Jacobian matrix of the vector b(x) with respect to x is denoted by B(x)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (23)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

Next, we determine the covariance matrix of the vector of changes in the state variables during the time interval (t, t + [DELTA]t)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (24)

The matrix M(x) is approximated by evaluating it at the critical point [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] corresponding to the pathogenic infection level

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (25)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and [A.sub.5] = -([delta] + [d.sub.I])[??].

For large N, the process [N.sup.1/2](x - [??]) is approximated by a stable bivariate Ornstein-Uhlenbeck process, with local drift matrix B([??]) and local covariance matrix M([??]). Then, the stationary distribution of the Ornstein-Uhlenbeck process is bivariate normal with mean 0 and variance [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] through the relationship

B([??])[summation] + [summation][B.sup.T] ([??]) = -M ([??]), (26)

where the superscript T is used to denote the transpose. After solving the above equation, we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (27)

Note that the parameter should also satisfy [[sigma].sub.1] > 0, [[sigma].sub.3] > 0. Thus, the diffusion approximation led to the conclusion that the marginal distribution of the infected population size in quasi-stationarity is approximately [q.sub..n]. To achieve consistency with the fact I [greater than or equal to] 0, the approximating normal distribution is modified by truncation at 1/2. Hence, we have the following approximation

[q.sub..n] [approximately equal to] [1/[square root of ([[sigma].sub.3]/N)]] [[PHI]((n - [??])/[square root of ([[sigma].sub.3]/N)])/(([??] - 0.5)/[square root of ([[sigma].sub.3]/N)])], (28)

where [PHI] and [phi] are, respectively, the standard normal c.d.f. and the standard normal p.d.f.

8.1. The Expected Time to Extinction. We find the expected time to extinction (E([[tau].sub.Q])) from quasi-stationary distribution. The expected time to extinction is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (29)

We see that the expected time to extinction is a function of population size N, which is a function when the population size is increasing (Figure 2).

9. Model Modification under Immune Host

In model (4), we consider the recovery class to be joined with the susceptible host, and it may increase the population size of uninfected host. It is to be noted that sometimes for a specific disease, the recovery from infected host may not come back to the susceptible class. For such a case, after the recovery from infected host, this specific recovery class becomes immunized and can be classified into a new category. Under this assumption the modified model can be written from model (4) by just ignoring the term [delta]I, in the growth equation of susceptible host. The algebraic manipulations for obtaining the quasi-stationary distribution and time to extinction for revised model are pretty similar with that of model (4).

10. Numerical Illustration

In the present study, a stochastic approach has been adopted to eradicate the infected population from a system and add on to the susceptible population. This is achieved by immunization or vaccination, and infected population recovers to susceptible population using the recovery rate [delta]. The stochastic model in this case predicts the extinction of infection, and the time to extinction is the important measure of the persistence of infection.

Figure 1 represents the normalized simulated marginal distribution profile of the infected host (I) in quasi-stationary state for N = 100, N = 200, and N = 300. It is observed that when the population size is N = 100, the quasi-stationary distribution is truncated and skewed, while when N is increased to N = 200 and subsequently to N = 300, the distribution displays higher kurtosis and enhanced symmetry. The population distribution of the infected class thus becomes much narrower for higher total strength of the susceptible and infected individuals. From Figure 2, we get the expected time to extinction of the infected host from quasi-stationary state when the total population (N) varies from N = 100 to N = 500. In this case, it is observed that the expected time to extinction increases with the increase in the total population N. When N is small, the number of infected population is also expected to be small, and at that time, the natural immune system can decelerate the rapid growth of the infected cells and add it to the susceptible class. This results in the lesser time to extinction at the initial stages of N. But when the population size is high, the number of infected class is also expected to rise, and the immune system cannot control the rapid growth of the infected host population, and this leads to the high extinction time as observed. Thus, it is evident from the numerical analysis that at the onset of the time to extinction, the infected individuals get transformed and move to the susceptible population.

11. Discussion and Conclusion

In this paper, we have presented a basic mathematical model of host-pathogen interaction using the stochastic approach based on the concept of quasi-stationary and the diffusion approximation result. The time to extinction has been predicted as a function of the total population size (N). The numerical simulation reveals that as the total population increases, the quasi-stationary distribution inclines to reduced skewness and narrower distribution of the infected class under consideration. The time to extinction of the infected class and its transformation to the susceptible class also vary with the total population size (N) and are found to exhibit a gradual rise with increasing value of the total population. Since the deterministic version is an approximation of the stochastic model, the estimation of the time to extinction is unlikely to be feasible with the deterministic version. So, we are inclined to conclude that the deterministic version of the model gives an approximation to the stochastic model.

The diffusion approximation led to the conclusion that the marginal distribution of the infected host population size in quasi-stationarity is approximately N([??], [square root of ([[sigma].sub.3]/N)]). To achieve consistency with the fact that I [greater than or equal to] 0, the approximating normal distribution is modified by truncating at 1/2. We conclude from the above result that as the total population size of the system becomes large (N [right arrow] [infinity]), the marginal distribution of S(t) is approximately N([??], [square root of ([[sigma].sub.1]/N)]), and the marginal distribution of I(t) is approximately N([??], [square root of ([[sigma].sub.3]/N)]). Further, the covariance is approximated by [[sigma].sub.3]. The numerical simulation confirms the analysis which is relevant in Figures 1 and 2. Further, the numerical simulation reveals that as the total population increases, the quasi-stationary distribution changes from positively skewed to symmetrical nature along with increases in the time to extinction.

The entire study has been carried out in a different angle to focus on the time to extinction of disease pathogenesis, and this reflects an advancement of our earlier reported work based on the decreasing infected population by deterministic modeling. Thus, we can conclude that stochastic modeling provides a more accurate prediction in finding out the expected time to extinction of infected population and hence disease pathogenesis.

http://dx.doi.org/ 10.1155/2013/381286

References

[1] J. Gani, "Modelling epidemic diseases," Australian and New Zealand Journal of Statistics, vol. 52, no. 3, pp. 321-329, 2010.

[2] H. I. Freedman, "A model of predator-prey dynamics as modified by the action of a parasite," Mathematical Biosciences, vol. 99, no 2, pp 143-155, 1990

[3] J. Chattopadhyay and N. Bairagi, "Pelicans at risk in Salton sea--an eco-epidemiological model," Ecological Modelling, vol.136, no. 2-3, pp. 103-112, 2001.

[4] E. Venturino, "Epidemics in predator-prey models: disease in the predators," IMA Journal of Mathematics Applied in Medicine and Biology, vol.19, no.3, pp.185-205, 2002.

[5] D. Mukherjee, "Persistence in a prey-predator system with disease in the prey," Journal of Biological Systems, vol. 11, no. 1, pp. 101-112, 2003.

[6] J. Chattopadhyay, A. Mukhopadhyay, and P. K. Roy, "Effect of viral infection on the generalized Gause model of predator-prey system," Journal of Biological Systems, vol.11, no.1, pp.19-26, 2003.

[7] I. Nasell, "Stochastic models of some endemic infections," Mathematical Biosciences, vol. 179, no. 1, pp. 1-19, 2002.

[8] P. K. Roy and B. Chattopadhyay, "Host pathogen interactions with recovery rate: a mathematical study," in Proceedings of the World Congress on Engineering (WCE '10), Lecture Notes in Engineering and Computer Science, pp. 521-526, London, UK, June-July 2010.

[9] A. N. Chatterjee and P. K. Roy, "Anti-viral drug treatment along with immune activator IL-2: a control-based mathematical approach for HIV infection," International Journal of Control, vol. 85, no. 2, pp. 220-237, 2012.

[10] E. Beltrami and T. O. Carroll, "Modelling the role of viral disease in recurrent phytoplankton blooms," Journal of Mathematical Biology, vol.32, no.8, pp.857-863, 1994.

[11] H. W. Hethcote and P. van den Driessche, "Two SIS epidemiologic models with delays," Journal of Mathematical Biology, vol. 40, no.1, pp.3-26, 2000.

[12] R. A. Ennos, "The introduction of lodgepole pine as a major forest crop in Sweden: implications for host-pathogen evolution," Forest Ecology and Management, vol. 141, no. 1-2, pp. 85-96, 2001.

[13] N. Bairagi, P. K. Roy, R. R. Sarkar, and J. Chattopadhyay, "Virus replication factor may be a controlling agent for obtaining disease-free system in a multi-species eco-epidemiological system," Journal of Biological Systems, vol.13, no.3, pp.245-259, 2005.

[14] N. Bairagi, P. K. Roy, and J. Chattopadhyay, "Role of infection on the stability of a predator-prey system with several response functions--a comparative study," Journal of Theoretical Biology, vol. 248, no. 1, pp. 10-25, 2007.

[15] T. Britton and P. Neal, "The time to extinction for a stochastic SIS-household-epidemic model," Journal of Mathematical Biology, vol.61, no.6, pp.763-779, 2010.

[16] I. Nasell, "A new look at the critical community size for childhood infections," Theoretical Population Biology, vol.67, no.3, pp.203-216, 2005.

[17] I. Nasell, "On the time to extinction in recurrent epidemics," Journal of the Royal Statistical Society B, vol. 61, part 2, pp. 309-330, 1999.

[18] P. K. Roy and J. Mondal, "Host pathogen interactions: insight of delay response recovery and optimal control in disease pathogenesis," Engineering Letters, vol. 18, no. 4, 2010.

[19] M. Carletti, "On the stability properties of a stochastic model for phage-bacteria interaction in open marine environment," Mathematical Biosciences, vol.175, no.2, pp.117-131, 2002.

[20] J. H. Matis, T. R. Kiffe, E. Renshaw, and J. Hassan, "A simple saddle point approximation for the equilibrium distribution of the stochastic logistic model of population growth," Ecological Modelling, vol.161, no.3, pp.239-248, 2003.

[21] J. H. Matis, T. R. Kiffe, and P. R. Parthasarathy, "On the cumulants of population size for the stochastic power law logistic model," Theoretical Population Biology, vol. 53, no. 1, pp. 16-29, 1998.

Priti Kumar Roy, (1) Jayanta Mondal, (2) Rupa Bhattacharyya, (3) Sabyasachi Bhattacharya, (4) and Tamas Szabados (5)

(1) Centre for Mathematical Biology and Ecology, Department of Mathematics, Jadavpur University, Kolkata, West Bengal 700032, India

(2) Department of Mathematics, Barasat College, Kolkata 700126, India

(3) Department of Chemistry, Narula Institute of Technology, Kolkata 700109, India

(4) Agricultural and Ecological Research Unit, Indian Statistical Institute, 203 B. T. Road, Kolkata 700108, India

(5) Department of Stochastics, Institute of Mathematics, Budapest University of Technology and Economics, Budapest 1521, Hungary

Correspondence should be addressed to Priti Kumar Roy; pritiju@gmail.com

Received 5 February 2013; Accepted 21 May 2013

Academic Editor: Xinyu Song

TABLE 1: Hypothesized transition rates for the stochastic version. Event Transition Transition rates Birth of a (m,n) [right [[lambda].sub.1] susceptible arrow] (m + 1, (m, n) = host n) [b.sub.1]m - [b.sub.2] (m + n)m + [delta]n Death of a (m,n) [right [[mu].sub.1] (m, susceptible arrow] (m -1, n) n) = [d.sub.1]m host + [d.sub.2] (m + n)m Infection of an (m,n) [right [[gamma].sub.2] uninfected host arrow] (m -1, n (m, n) = + 1) ([lambda] + ([alpha]/ ([beta] + m))) mn Death or (m, n) [right [[mu].sub.2] (m, recovery of an arrow] (m, n - n) = ([d.sub.1] infected host 1) + [delta])n TABLE 2: Possible changes in the two-population system (4) with the probabilities. Change Probability [DELTA][x.sub.l] = [p.sub.1] = (([b.sub.1] [[1,0].sup.T] - [b.sub.2] (S + I))S + [delta]I) [DELTA]t [DELTA][x.sub.2] = [p.sub.2] = ([d.sub.l] + [[-1,0].sup.T] [d.sub.2] (S + I)) S[DELTA]t [DELTA][x.sub.3] = [p.sub.3] = ([lambda] + [[-1,l].sup.T] [alpha]/ ([beta] + S)) SI [DELTA]t [DELTA][x.sub.4] = [p.sub.4] = ([d.sub.I] + [[0,-l].sup.T] [delta]) I[DELTA]t TABLE 3: Values of parameters used in model dynamics. Parameter Default value Definition ([day.sup.-1]) r Maximal growth rate of 0.2 [18] susceptible host K Carrying capacity 20 [18] [b.sub.1] Intrinsic rate 0.3 [21] [d.sub.1] Crowding coefficient 0.02 [21] [h.sub.[gamma] Half-maximal at a target cell 9 [18] density [lambda] Force of infection through 0.2 [18] contact with infected host [gamma] Force of infection through 0.04 [18] contact with pathogens [d.sub.I] Lysis death rate of infected 2.5 liter [18] host [eta] Pathogens replication factor 115 [18] [mu] Mortalityrate of pathogen 2.2 [18] [delta] Recovery rate of infected host 4 [18]

Printer friendly Cite/link Email Feedback | |

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

Author: | Roy, Priti Kumar; Mondal, Jayanta; Bhattacharyya, Rupa; Bhattacharya, Sabyasachi; Szabados, Tamas |

Publication: | Journal of Applied Mathematics |

Article Type: | Report |

Date: | Jan 1, 2013 |

Words: | 5115 |

Previous Article: | Generalised adaptive harmony search: a comparative analysis of modern harmony search. |

Next Article: | An equilibrium chance-constrained multiobjective programming model with birandom parameters and its application to inventory problem. |

Topics: |