key: cord-0963406-bn7t1qwj
authors: Inamo, Jun
title: How should we overcome the threat by the pandemic of 2019-nCoV? Epidemic simulation using the SIRS model
date: 2020-05-09
journal: Clin Rheumatol
DOI: 10.1007/s10067-020-05083-7
sha: 2888d4384154787fb8ce110a48e1d9625bc04ad0
doc_id: 963406
cord_uid: bn7t1qwj

nan

Since the first identification in December 2019 in Wuhan, China, a number of cases of novel coronavirus (2019-nCoV)infected patients is increasing around the world. It has been reported that there are high-risk populations with high mortality, such as elderly people and those with comorbidities [1] . However, no specific treatment has been developed yet. Historically, humans have overcome the pandemics of viral infection by preventing secondary infection. Several studies have simulated the spread of 2019-nCoV, and the basic reproduction number (R 0 ) ranges from 2.2 to 3.6 [2, 3] . On the other hand, patients treated with disease-modifying antirheumatic drugs may have insufficient antibody response, implying the possibility of reinfection [4] .

Here, to simulate the efficacy of preventing secondary infection in population with consideration of reinfection, I simulated a scenario of one person with 2019-CoV infection entering into a population of 1000 people with the hypothesis that a part of them lose immunity. I used the SIRS model which is composed of different categories within a population; the susceptible (S), infected (I), recovered (R), and returned-susceptible (S) compartments, respectively. The model can be written as the following ordinary differential equation:

where S, I, and R are proportions of susceptible, infected, and recovered individuals, respectively, β is the transmission rate of infection, γ is recovery rate (inverse of the infectious period), and ξ is returning rate of recovered individuals to the susceptible statue due to loss of immunity. N (= S + I + R) is the total population. The simulation was started with S = 999, I = 1, and R = 0 for the following 100 days after the entry of one person with 2019-nCoV. Considering the uncertainty of parameters in the real-world outside China, I combined parameters based on considerable circumstances. Analyses were performed by R software. First, parameters were set as (β, γ, ξ) = (0.5, 1/20, 0.01) based on assumption that infectious duration is 20 days [5] and 1% of recovered individuals return to the susceptible statue due to loss of immunity. As β decreased 0.5, 0.4, 0.3, 0.1, and 0.05 (Fig. 1 a, b, c, d, e, and f, respectively) , the maximum number of infected people decreased and the peak delayed, suggesting the effectivity of decreasing transmission rate by standard precautions. When β was 0.05, the secondary infection did not occur. R 0 in each situation was 10, 8, 6, 4, 2, and 1, respectively. Next, to evaluate the efficacy of prompt isolation of infected patients, γ was stratified with 1/20, 1/16, 1/12, 1/8, 1/4, and 1/1 with keeping βN and ξ the same. As a result, the maximum number of infected people decreased as γ increased (Fig. 2 a, b, c, d, e, and f, respectively) . The secondary infection did not occur in the case of γ = 1/1. R 0 in each situation was 10, 8, 6, 4, 2, and 1, respectively. These trends were the same as the case ξ = 0 (Supplementary Electronic supplementary material The online version of this article (https://doi.org/10.1007/s10067-020-05083-7) contains supplementary material, which is available to authorized users. This study suffers from several limitations. First, since the clinical characteristics of 2019-nCoVare largely unknown, the value of parameters used in this simulation may not be proper. Second, the SIRS model is a simple model based on some assumptions; the population is closed (no vital dynamics such as birth and death) and the course of the infection is short (emergent outbreak). Also, the possibility that asymptomatic patients could be infectious to others is not considered in this model [6] . Third, because the SIR model assumes a homogeneous infection network between individuals, the heterogeneity of individuals such as disease status was not considered. Forth, other than standard precautions, β is also affected by virologic characteristics and factors of hosts. Therefore, the result of this study cannot apply to data in the real-world directly. However, this study demonstrated that outbreaks can be prevented by performing appropriate precautionary measures even in populations with consideration of the possibility of reinfection.

Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72 314 cases from the Chinese Center for Disease Control and Prevention

Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: a data-driven analysis in the early phase of the outbreak

Feng Z (2020) Early transmission dynamics in

Immunogenicity and safety of pneumococcal vaccination in patients with systemic sclerosis

SARS-CoV-2 viral load in upper respiratory specimens of infected patients

Transmission of 2019-nCoV infection from an asymptomatic contact in Germany

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations