Poisson integral type quarantine in a stochastic SIR system

Abstract: We propose a SIR system that includes a Poisson measure term to model the quarantine of infected individuals. An inequality concerning the term representing the transmission rate is given to establish the stochastic stability of the disease free equilibrium. It is further shown that if R0 > 1 then the long-run behavior the system will reside within a neighborhood of the equilibrium in the underlying deterministic version of this system.


Introduction
Many authors have considered continuous time stochastic epidemiological models. For instance, the papers [1] , [2], and [3] analyzed various but similar SIR models perturbed by a Gaussian term in the form of an Itô integral, [4] considered a two-group SIR model with similar Gaussian perturbations, and [5] analyzed a SIS model also perturbed by an Itô integral. However, none of these models have taken into account anomalies that a↵ect the system. Recently, there have been some e↵orts to incorporate anomalies into the dynamic by incorporating a Poisson integral to capture the instantaneous impact. The paper [6] assumed temporary immunity was captured by random jumps in a SIR model, while [7] and [8] both considered general environmental that may instantaneously a↵ect di↵erent SIR models, and [9] assumed a SIS model with a Poisson integral capturing outside e↵ects. Excluding the manuscript in [9], all of the models consider a compensated-Poisson integral. While there is a clear effect on the dynamic, the compensated-Poisson integral disappears in expectation and may not properly capture the phenomena.
Kuske, Gordillo and Greenwood in [3] investigated an SIR model perturbed by multiple and independent Brownian terms, and through a linear approximation, give conditions for a unique stationary measure establishing that an epidemic will occur. Roa, Wang, and Li [1] also considered an SIR, but the authors assume a nonlinear incident rate, and after analyzing the deterministic model, perturb the system twice with a two-dimensional Brownian motion, one for an endemic and the other for the disease free equilibrium. The authors then show that these systems are well defined and give conditions for almost sure convergence to an endemic and the disease free equilibrium. Yu, Jiang, and Shi [4] analyze a Brownian perturbed two-group SIR model with a stationary point at the deterministic endemic and determine su cient conditions for the endemic equilibrium to be stochastically stable. The models examined in [2] and [3] are very close to the models considered in this paper and are explored in detail below.
For a fixed population of size N, we consider the following SIR model of the form where S (t), I(t), and R(t) denote respectively the frequencies of the susceptible, infected, and removed; ergo, S (t) + I(t) + R(t) = 1. The constant µ represents the birth and death rate (newborns are assumed to be susceptible), is the recovery rate for the individuals that are infected, and represents the transmission rate. (Please see [10] for further information.) In this system, when  + µ the disease free equilibrium is globally asymptotically stable. If We investigate a stochastic version of this SIR model under a Gaussian noise of interactions and include quarantine as a random event described by a Poisson process. Take the function j(y) to be continuously di↵erentiable where min j(y) > 0, and max j(y) < 1. Define N(dt, dy) as a Poisson measure with ⌫(·) as its intensity measure, and W(t) as an independent standard Brownian motion. This system is described by the right-continuous stochastic di↵erential equation Since the Itô term is well known and represents the noise in the transmission process, we focus on the Poisson integral and give intuition for this term. Suppose that a specific percentage of the population is quarantined randomly with respect to a Poisson distribution. Denoting this percentage as C and N(·) as a Poisson distribution then the stochastic e↵ect of the quarantine at time t is C · N(t). However, this kind of modeling is not robust enough in capturing the lack of e ciently of placing a subpopulation in quarantine. This is especially true of a population adjusting to an endemic.
To incorporate di↵erent impact sizes into the evolution, we briefly discuss the Poisson process as a temporal and spatial random measure, and interpret an integral under this measure. For further information, please see Applebaum [11]. For a Lévy process X(s) and B 2 B as the intensity measure. If ⌫ R = 1 it is possible for an anomaly to happen an infinite number of times in a finite time interval. Therefore, the assumption that ⌫ R < 1 is quite natural. Taking B (·) to be the identity function over the measurable set B then at time t the total e↵ect is j ( X(s)) B ( X(s)) .
When the deterministic integral Z t 0 Z R j(y)⌫(dy)ds, (ds is the Lebesgue measure) is added to the Poisson integral the expectation of this e↵ect is zero. Specifically, forÑ(ds, dy) := N(ds, dy) ⌫(dy)ds, the expected net e↵ect of the "adjusted" quarantine E [12], Bertoin [13], or Applebaum [11] for further information. This observation will be utilized in the analysis of the dynamic.
For completeness of the model, for the triple (⌦, F , P) we assume the filtration {F t } t2R + is rightcontinuous and F 0 contains all of the null sets of F . Tornatore, Buccellato, and Vetro [2] analyzed System (1.2) when j ⌘ 0 and displayed when < , 2µ the dynamic converged to the disease free equilibria of the noiseless version of the system. The authors displayed through simulations that this inequality is a good threshold. A stochastic SIR with delay was then analyzed and determined that if the same condition holds the process will converge to the disease free equilibria almost surely. Further, Tornatore and Buccellato [14] show that if > + µ + 2 2 then a strong endemic will occur.
Since this seminal paper [2], many authors have further refined their results by considering similar dynamics (see [15], [16], and [17]). When the various parameters were set to zero, the models in [15] and [16] showed that the model given in [2] converges to the disease free equilibria when < +µ Furthermore, for R 0 := + µ the deterministic endemic equilibrium, the authors in [16] showed that the dynamics stay within a neighborhood of R 0 if > + µ, i.e., if R 0 > 1. Interestingly, Chen, Li and Liu [17] showed that the dynamics are stable to the disease free equilibrium when < + µ + 2 2 .
In this paper it will be shown that System (1.2) is well-defined, and we give conditions for stability of the disease free equilibrium for both dynamics, as well as give conditions for System (1.2) to be within a neighborhood of the endemic equilibrium of the underlying deterministic system.
While the methods applied are well-known in the analysis of stochastic di↵erential equations, analysis of the dynamic are not necessarily restricted these and similar methods. For instance, Calatayud, Cortés and Jornet [18] establish an algorithm that will quickly approximate the probability density function. The authors utilize polynomial expansions, the random variable transformation (RVT) technique, and multidimensional integration schemes to derive a solution with an algorithm. The authors display the e cacy of their algorithm with an application to a SIR model. The authors in [19] also focus on the probability density function in a more rigorous manner via a randomized Markov chain SIR-type model.

Stochastic extinction of the disease
o as the simplex. To ensure the process is well-defined, we show that System (1.2) is invariant to the simplex. In other words, given an initial condition that lies in simplex, System (1.2) will stay in the simplex for all finite time.
Denoting L as the infinitesimal generator for System (1.2), L has the form where g x is a twice di↵erentiable function.
Proposition 2.1. For all finite t, given that Proof. We will invoke Theorem 2.1 in Meyn and Tweedie [20] and apply a norm-like function (V : 3 ! R + that is measurable and V x ! 1 when x tends to the boundary) such that, for some c > 0, LV  cV.
Define V x = log(x 1 ) log(x 2 ) log(x 3 ). Then for x 2 3 LV x = 1 Clearly, there exists a positive constant c such that LV x  cV x . ⇤ To determine conditions for convergence to the disease free equilibrium, we use the stochastic Lyapunov method. Examining the process, one can see that the only stationary position is the point e 1 = (1, 0, 0). Noticing Therefore there exists a positive constant k such that L f  k f . Theorem 4 and Remark 2 in [21] (page 325) yields that for an ✏ > 0 there exists a neighborhood of e 1 , say U, such that for Fix an arbitrary ✏ > 0 and define ⌧ ✏ = inf which will assert our theorem. Define ✏ 0 = sup n S ⌧ ✏ o . Since the process does not hit the boundary in finite time, we have that ✏ 0 < 1 a.s. For the function f above, L f x  ↵ 1 x 1 2 , for some ↵ > 0.

Conditions for an epidemic
Following the method applied in [16], we will show that if the deterministic condition R 0 > 1 holds, i.e., + µ > 1, then the dynamic remains within a neighborhood of the deterministic epidemic equilibrium. Although the result in [17] might seem to contradict this result, since < + µ + 2 2 is the condition needed for almost sure convergence to the disease free equilibria, this does not mean this equilibrium is not in the neighborhood given by [16]. What we can say for certain is that when the disease free equilibria is not in this neighborhood, an endemic will occur. Define as the deterministic endemic equilibrium. Since our process is in the simplex, we will show that process is "near" , which would imply that the process is in a neighborhood The result is very similar to [16], where the jump perturbation adds to the radius of the ball around the deterministic endemic equilibria. In particular, the term k J displays how when a random percentage of the infected population is consistently quarantined this decreases the impact to the general population.
Clearly, when both perturbations are small, an endemic will occur. Proof.
, (which is similar to the function defined in [16]). Then by Itô's extended lemma we have First we simplify and manipulate the drift, then we will simplify and manipulate the Poisson integrand. Notice that k For the Poisson integrands, we have Z  The last equality comes from the integrand being bounded and ⌫ R < 1. Taking (3.1) and the deterministic measure in (3.2) we see that

Thus the Poisson integral is
By the Strong Law of Large Numbers for local martingales, we have almost surely that we are able to conclude that

Conclusion
In this paper quarantine of infected individuals was modeled with a Poisson measure. Due to the imperfections of human interactions, we have given a model that better represents the e↵ect of quarantine. After the derivation of the model, we determined that this dynamic was well defined. In Theorem 2.1 we showed that if the inequality < min , 2µ holds then the dynamic will converge almost surely to the disease free equilibrium. The condition for the disease free equilibrium maybe counterintuitive since the term in the right side of the inequality, shows that quarantine does not fully a↵ect the dynamic. However, what this tells us is that the quarantine does not have a complete e↵ect on the dynamic of the disease and, in some instances, may only slightly dampen the impact to the population.
In Theorem 3.1, computational conditions for an endemic where given. The condition for an endemic is rather encompassing, as it includes all possible strengths of recurrence. The result shows the complexity the quarantine adds to the dynamic encompassed in the complex term k J . The term k J displays how the quarantine has an e↵ect on the entire populace, and in particular how a large quarantine does not completely stop the evolution of the disease and how a small quarantine has a larger e↵ect than expected.    Figure 1(c), shows intuition behind the random quarantine on how and this term severely a↵ects the dynamic and increases the size of the neighborhood of the deterministic endemic equilibrium, as shown in Theorem 3.1.

Conflict of interest
The authors declare no conflicts of interest in this paper.