Failure of R0 near the epidemic threshold in the classical SIS model

The primary predictor of a disease outbreak and severity is the basic reproduction number R0, which represents the average number of secondary cases produced by introducing an infected individual into an entirely susceptible population. According to the classical SIS model, a disease with R0 less than one will eventually die out and persist if R0 is greater than one. Using the pair-approximation method, we reconstruct the classical SIS model by explicitly accounting for neighbourhood interactions between susceptible and infected individuals. Specifically, the disease can only be transmitted, with some transmission probability, if a susceptible individual is surrounded by at least one infected individual within its direct neighborhoods. Despite the simplicity of the SIS model present here, results produced by the pair-approximation model deviates significantly from predictions by the mean-field approximation model, particularly near the epidemic threshold. Contrasting the standard SIS model based on the mean-field approach, we find scenarios where the disease dies out even if R_0 is greater than one. We suggest a crucial need for redefining the basic reproduction number on a smaller spatial scale and taking the average R0 over a global scale, rather than applying it globally to an entire population. However, in the realm of more intricate models of infectious diseases, it remains an open question to what extent mean-field approximation predictions diverge from predictions produced by models that consider neighborhood interactions.


Introduction
The basic reproduction number  0 quantifies the average number of secondary infections produced by an infected individual when introduced into a completely susceptible population (Dietz, 1993).The basic reproduction number has long been considered as import predictor of a disease spread and severity, and often plays an important role for supporting decisions of public health officials and policymakers (Guerra et al., 2017;Zhao et al., 2020).Estimation of the basic reproduction number is typically complex as it depends on various biological and environmental factors (Delamater et al., 2019).However, in the simple classical susceptibleinfected-susceptible (SIS) model, the basic reproduction number emerges from just two key parameters and is defined as the average disease transmission rate divided by average recovery rate per an infected individual (Allen, 1994).It is commonly accepted that a disease with  0 less than 1 will eventually die out and persist if  0 is greater than 1.The epidemic threshold, occurring at  0 equal to 1, represents the boundary between the extinction and sustained transmission of the disease.The outcome of a disease transmission away from the epidemic threshold might be intuitively straightforward, depending on the transmission complexity.
However, estimation of  0 for the same disease can vary as the model set up and complexity change.Depending on the method used to calculate  0 , many diseases with an  0 value greater than 1 are predicted to die out and persist when  0 is less than 1, leading to a failure of  0 (Li et al., 2011).The basic reproduction number is primarily defined on a global scale, often disregarding the small local scale on which interactions between susceptible and infected individuals occur.In reality, an infected individual cannot infect a susceptible one if they are spatially distant.Nevertheless, it is commonly assumed by the mean-field approximation method that all individuals within a population do interact with one another, and this simplification has been widely accepted.Network and pair-approximation models of infectious diseases typically account for interactions between neighbouring pairs of individuals (Keeling, 1999;Keeling and Eames, 2005;Payne, 2019).Such models often describe the dynamics of infectious disease more accurately than the mean-field approximation models.The pairapproximation method has been employed to describe systems other than infectious diseases where neighborhood interactions are crucial determinants of the dynamics outcomes (Harada and Iwasa, 1994;Mohammed et al., 2023Mohammed et al., , 2018)).This particular methodology correctly captures systems dynamics predicted lattice-based models, offering mathematical tractability.
Using the pair-approximation method, we reconstruct the classical SIS model by explicitly accounting for neighborhood interactions among susceptible and infected individuals.The pairapproximation approach assumes that a susceptible individual can only get infected if it is surrounded by at least one infected individual present in its direct neighbourhoods.That is a necessary condition for a disease transmission, but the transmission occurs with some probability, following an interaction between susceptible and infected individuals.The core object of this work is to assess whether the  0 concept applies when neighbourhood interactions are taken into account in infectious disease models.Although the SIS model presented here is simple, but results of the pair-approximation model deviates significantly from previous predictions by the mean-field approximation model.Contrasting the SIS model that based mean-field approximation, the pair-approximation model exhibits a scenario where a disease can die out even if  0 is greater than 1.This divergence in the predictions of the two models occurs mainly around the epidemic threshold.More complex models of infection diseases based on the mean-field approximation are thus expected to diverge more significantly from models that explicitly account for neighbourhood interactions.We suggest that estimations of the basic reproduction number should be redefined on a small neighbourhood scale, rather than on a global scale for the entire population.Estimation of the basic reproduction number might not be enough to inform public health officials about the current and future behaviour of the disease outbreak dynamics, and this fundamental concept of  0 can be misleading by omitting some of its important properties such as spatial scales.

Pair-approximation model of disease transmission
Pair approximation is a method of constructing ordinary differential equation models, describing the dynamics of neighbouring pairs, the global and local dynamics of a given population (Harada and Iwasa, 1994).It provides an explicit description for neighborhood interactions among individuals and accurately capture the spatiotemporal dynamics predicted by lattice-based models (Mohammed et al., 2023).
Here, we use the pair-approximation method to reconstruct the classical SIS model by explicitly considering neighbourhood interactions among susceptible and infected individuals (Keeling, 1999;Payne, 2019).We denote by   and   the probabilities that a randomly chosen individual is susceptible and infected, respectively.These probabilities represent the global frequencies .
(/) and (/) of susceptible and infected individuals, where  and  are the number of susceptible and infected individuals in a population of size  =  + , with   +   = 1.The disease is transmitted locally when a susceptible individual comes into a direct contact with an infected individual present in its own neighborhood.The pair-approximation method describes the average local frequency of infected individuals present in the neighborhood of a susceptible individual.The number of nearest-neighbouring sites of a susceptible individual is quantified by the positive integer .The average local frequency of infected individuals in the nearest neighborhoods of a susceptible individual is quantified by  | , representing the conditional probability that a randomly chosen individual in the neighborhood of a susceptible individual is infected (Harada and Iwasa, 1994).The local frequency of infected individuals is the number of infected individuals in the neighborhood of a susceptible individual, averaged over the number of nearest-neighbouring sites .A susceptible individual can be infected if at least one of neighbours in the  sites is infected.
The disease transmission rate is denoted as .The infected individual recovers at rate  in a period of time 1/.In the SIS model, a susceptible individual is assumed to become susceptible again after recovery.
The differential equations governing the global disease dynamics are described by ̇=  |   −   . (2) The first term in Eqns.( 1) and ( 2) describes the neighborhood interaction between susceptible and infected individuals.The last term describes the overall recovery rate.The conditional probability  | needs to be specified.It is mathematically defined as  | =   /  (Harada and Iwasa, 1994), where the pair   =   represents the joint probability that a randomly chosen neighbouring pair of individuals are susceptible and infected.Notably, there are four possibilities of choosing a neighbouring pair of individuals at random:   ,   ,   , and   , with   + 2  +   = 1.
Following (Keeling, 1999;Payne, 2019), we construct differential equations governing the dynamics of neighbouring pairs.The following transitions contribute to the dynamics of the   pair: 1.A pair of infected individuals   contributes positively to the pair   if one of the infected individuals in the   pair has been recovered 2. A pair of susceptible individuals   contributes positively to the pair   if one of the susceptible individuals in the   pair has been infected by one of its direct neighbour from outside the pair, i.e., from the remaining nearest-neighbouring sites ( − 1). 3. A pair of susceptible and infected individuals   contributes negatively to the pair   if the susceptible individual in the   pair has been infected by its direct neighbour in the pair 4. A pair of susceptible and infected individuals   contributes negatively to the pair   if the susceptible individual in the   pair has been infected by one of its direct neighbour from outside the pair, i.e., from the remaining nearest-neighbouring sites ( − 1). 5. A pair of susceptible and infected individuals   contributes negatively to the pair   if the susceptible individual in the   pair has been recovered Using the pair-approximation approach (Harada and Iwasa, 1994;Keeling, 1999;Mohammed et al., 2023), the dynamics of the   pair are thus described by (3) To reduce the model complexity, we replace the neighbouring pair   by (1 −   − 2  ), and we only need to describe the   pair dynamics to construct an equation for the average local frequency  | .In Eqn.
(3), the first term indicates the transition from   to either   or   , thus the factor 2. In the subsequent term, the quantity / represents disease transmission rate per infected individual present in the neighborhood of a susceptible individual, and the transition from   to either   or   will require an infected individual from the remaining neighbouring sites ( − 1), thus the multiplication by the conditional probability  | .This conditional probability would require the computation of higher-order frequencies, and would not allow to obtain a closed system of differential equations.Therefore, pair-approximation method ignores the effect of the neighbour-of-the neighbour (it assumes correlation is weaker with non-direct neighbours, i.e., correlation between sites decreases exponentially with the distance between them, (see, e.g., Harada and Iwasa, 1994)  The following transitions contribute to the dynamics of the   pair: 1.A pair of susceptible and infected individuals   contributes positively to the pair   if the susceptible individual in the pair   has been recovered 2. A pair of susceptible individuals   contributes negatively to the pair   if any of the susceptible individuals in the pair   has been infected by its direct neighbour from outside the pair, i.e., from the remaining nearest-neighbouring sites ( − 1).
The dynamics of the neighbouring pair   are thus described by Using the definition of the local frequency of infected individuals  | =   /  , and taking the time derivative of both sides, the dynamics of the average local frequency of infected individuals in the neighborhoods of a susceptible individual are described by Equations ( 1) -( 5) represent the reconstructed SIS model under neighborhood interactions among susceptible and infected individuals present in the population.Equations ( 1) and ( 2) corresponds to the mean-field approximation model when  | is replaced by   .

The basic reproduction number
As defined in the traditional SIS model (Daley and Gani, 2001), the basis reproduction number is given as  0 = /.According to the mean-field SIS model, a disease with  0 < 1 will die out and persist if  0 > 1.The bifurcation point occurs at  0 = 1, separating the extinction and persistence regions of the disease. .

Numerical analysis
The numerical simulation of the pair-approximation and the mean-field approximation SIS models has been implemented in MATLAB.Parameter values used to simulate the model are provided in Table 1.

Results
We showcased the SIS model simulation for the global frequencies of susceptible and infected individuals,   and   , respectively.The dynamics of neighbouring pairs (  and   ) and the average local dynamics of infected individuals present in the neighborhood of a susceptible individual,  | .We compared results produced by the pair-approximation model with previous predictions by the classical mean-field SIS model.
In the pair-approximation model, the number of nearest-neighbouring sites of a susceptible individual  reflects the potential number of susceptible and infected individuals present within the direct neighbourhood of a susceptible individual.In the mean-field approximation model, a susceptible individual has the same likelihood of interacting with all infected individuals present within a population.That is a very unlikely situation and can affect predictions of disease outbreaks.However, it is intuitive to assume that a susceptible individual can typically get infected when interacts with an infected individual present within its direct four neighbouring sites (i.e.,  = 4).
We first restricted our model to four nearest-neighbouring sites.We found that results produced by the pair-approximation model significantly deviated from predictions by the mean-field approximation.This deviation occurred, particularly, close to the epidemic threshold,  0 = 2.5 and  0 = 1.17 (Fig. 1b and c).For large and small values of  0 (10 and 0.8), corresponding to the sustained transmission and extinction of the disease, the two models produced very similar results (Fig. 1a and d).According to the classical mean-field SIS model, the disease can spread out and persist when  0 = / > 1, and this principle has been widely accepted for this specific model.However, results of the proposed pair-approximation model demonstrated scenarios where the disease dies out even if  0 > 1 (Fig. 1c; and Fig. 3a ).This divergence is likely driven by the consideration of neighborhood interactions between susceptible and infected individuals in the pair-approximation model.In fact, this is the primary difference between the mean-field and the pair-approximation models presented in this work.
Furthermore, we simulated the dynamics of the reconstructed SIS model for a large number of nearest-neighboring sites ( = 100) of a susceptible individual.This indicates that a susceptible individual can potentially interact with 100 individuals in its surroundings, an unlikely situation.Some of these 100 individuals can be infected and can transmit the disease to this susceptible individual at transmission rate .By considering a very large number of individuals in the neighborhoods of a susceptible individual, the pair-approximation model corresponds to the mean-field approximation model.Therefore, our results demonstrated almost an identical behaviour of the disease outbreak dynamics, produced by the pair-approximation and the meanfield approximation models (Fig. 2; and Fig. 3b).The two models approached the same longterm attractor of the disease, with a small difference in the initial transient phase.
. Fig. 2. The temporal dynamics of the global frequencies of susceptible (  ) and infected (  ) individuals under four scenarios of the  0 value (10, 2.5, 1.17, and 0.8).The description, initial conditions, and parameters used are the same as in the caption of figure 1, and  = 100.
We systematically assessed the influence of the basic reproduction number on the epidemic size and compared the corresponding results of the two models.Away from the epidemic threshold ( 0 = 1), the outcome of the disease transmission is intuitively straightforward.Of a particular interest is to predict the disease transmission near the epidemic threshold.Our findings underscored significant differences between the results of the two models (Fig. 3a), especially around the epidemic threshold.At small to moderate values of  0 , the epidemic size predicted by the pair-approximation model is lower than that predicted by the mean-field approximation model (Fig. 3a).More precisely, the epidemic size curve produced by the pair-approximation model is shifted forward, extending the disease extinction region.Notably, this forward shift emerged only under the .consideration of four nearest-neighbouring individuals ( = 4).The results of the two models produced the exact dependence of the epidemic size on the basic reproduction number  0 when the number of nearest-neighbouring sites is chosen to be very large,  = 100 (Fig. 3b).
The temporal dynamics for the average local frequency of infected individuals within the neighbourhood a susceptible individual, and the dynamics of neighbouring pairs are shown for a number of nearest-neighbouring sites  = 4 (Fig. 4) and  = 100 (Fig. 5).
Results demonstrated that the average local frequency of infected individuals decreases as the basic reproduction number  0 increases (Fig. 4 and 5).A small number of infected individuals are present in the neighbourhoods of each susceptible individual when the basic reproduction number is small.Thus, the probability of choosing an infected individual from the neighbourhood of a susceptible individual is low.The average local frequency of infected individuals will also reflect the probability of choosing a neighbouring pair of susceptible and infected individuals.The probability of randomly choosing a neighbouring pair of susceptible and infected   depends on the absolute difference between the global probabilities   and   .
The probability of choosing the pair   is small when the absolute difference is high, and vice versa.The probability of randomly choosing the neighbouring pair   typically increases as the global frequency   increases.Thus, a small values of  0 will lead to higher probabilities of a neighbouring pair of susceptible individuals (Fig. 4).However, under the consideration of very large number of nearest-neighbouring sites, the probability of choosing an infected individual from the neighbourhood of a susceptible individual is relatively high when  0 is slightly greater than one (Fig. 5c).

Conclusion
The basic reproduction number  0 has been used for several decades as an indicator of whether the disease will die out or persist.It is widely acknowledged that the disease is likely to die out if  0 < 1, while persist if  0 > 1.In this work, we reconstructed the classical SIS model using the pair-approximation method (Harada and Iwasa, 1994;Keeling, 1999;Payne, 2019), and presented a theoretical evidence that a disease with  0 > 1 dies out in a simple SIS model.Despite this simplicity of the model, results produced by the pair-approximation model deviated pronouncedly from predictions by the mean-field approximation model.The advantage of the pair-approximation approach over the commonly-used mean-field approximation is that it explicitly accounts for neighborhood interactions among susceptible and infected individuals.However, a recent pair-approximation SIS model produced results in alignment with predictions by the classical mean-field SIS model, yet diverging from our findings in terms of the dependence of epidemic size on  0 (Payne, 2019).That model, while informative, solely focused on dynamics of pairs and did not consider the average local frequency dynamics of infected individuals; thus, overlooked the effect of the number of nearest-neighboring sites of a susceptible individual on the disease transmission dynamics.In the realm of more intricate models of infectious diseases, it remains an open question to what extent mean-field approximation predictions deviate from predictions produced by models that consider neighborhood interactions (Keeling, 1999).We suggest a crucial need for redefining the basic reproduction number on a smaller spatial scale and taking the average  0 over a global scale, rather than applying it globally to an entire population.This redefinition may depend on the average number of infected individuals within the neighbourhoods of a susceptible individual (Keeling, 1999).This fundamental epidemiological concept of  0 might otherwise be misleading (Li et al., 2011). .

Fig. 3 .
Fig. 3.The equilibrium dynamics of the epidemic size (y-axis) as a function the basic reproduction number  0 (x-axis), predicted by the pair-approximation (PA) model (bright curve) and the mean-field approximation (MFA) model (dark curve).Initial conditions and parameter values are the same as in the caption of figure 1.

Fig. 5 .
Fig.5.The temporal dynamics of the average local frequency ( | , pink curve), neighbouring pairs of susceptible and infected individuals (  , grey curve), and neighbouring pairs of susceptible individuals (  , cyan curve) predicted by the pair-approximation model under four scenarios of the  0 value (10, 2.5, 1.17, and 0.8).Initial conditions and parameter values are the same as in the caption of figure 4.
and approximates  | ≈  | and  | ≈  | .That is, correlations with non-nearest neighbours are approximately reconstructed from nearest neighbour correlations.Similar explanations apply to the other terms in Eqn.(3).

Table 1 :
Variable and parameter descriptions, and numerical values used in the model.