# Accounting for symptomatic and asymptomatic in a SEIR-type model of COVID-19

Jayrold P. Arcede, Randy L. Caga-Anan, Cheryl Q. Mentuda, Youcef Mammeri

► To cite this version:

Jayrold P. Arcede, Randy L. Caga-Anan, Cheryl Q. Mentuda, Youcef Mammeri. Accounting for symptomatic and asymptomatic in a SEIR-type model of COVID-19. Mathematical Modelling of Natural Phenomena, 2020, 15, pp.34. 10.1051/mmnp/2020021. hal-02904005

HAL Id: hal-02904005

https://hal.science/hal-02904005

Submitted on 21 Jul 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ACCOUNTING FOR SYMPTOMATIC AND ASYMPTOMATIC
IN A SEIR-TYPE MODEL OF COVID-19

JAYROLD P. ARCEDE<fnref n="1" />, RANDY L. CAGA-ANAN<fnref n="2" />, CHERYL Q. MENTUDA<fnref n="1" />,<fnref n="3" />
AND YOUCEF MAMMERI<fnref n="3" />,*

**Abstract.** A mathematical model was developed describing the dynamic of the COVID-19 virus over a population considering that the infected can either be symptomatic or not. The model was calibrated using data on the confirmed cases and death from several countries like France, Philippines, Italy, Spain, United Kingdom, China, and the USA. First, we derived the basic reproduction number, $R_0$, and estimated the effective reproduction $R_{eff}$ for each country. Second, we were interested in the merits of interventions, either by distancing or by treatment. Results revealed that total and partial containment is effective in reducing the transmission. However, its duration may be long to eradicate the disease (104 days for France). By setting the end of containment as the day when hospital capacity is reached, numerical simulations showed that the duration can be reduced (up to only 39 days for France if the capacity is 1000 patients). Further, results pointed out that the effective reproduction number remains large after containment. Therefore, testing and isolation are necessary to stop the disease.

**Mathematics Subject Classification.** 92D30, 37N25, 34D20.

Received April 2, 2020. Accepted June 1, 2020.

# 1. INTRODUCTION

In late 2019, a disease outbreak emerged in a city of Wuhan, China. The culprit was a certain strain called Coronavirus Disease 2019 or COVID-19 in brief [13]. This virus has been identified to cause fever, cough, shortness of breath, muscle ache, confusion, headache, sore throat, rhinorrhoea, chest pain, diarrhea, and nausea and vomiting [2, 8]. COVID-19 belongs to the *Coronaviridae* family. A family of coronaviruses that cause diseases in humans and animals, ranging from the common cold to more severe diseases. Although only seven coronaviruses are known to cause disease in humans, three of these, COVID-19 included, can cause a much severe infection, and sometimes fatal to humans. The other two to complete the list were the severe acute respiratory syndrome (SARS) identified in 2002 in China, and the Middle East respiratory syndrome (MERS) originated decade after in Saudi Arabia.

<footnote>Keywords and phrases: COVID-19, SEI_s I_a UR Model, reproduction number, interventions.</footnote>

<footnote n="1">¹ Department of Mathematics, Caraga State University, Butuan City, Philippines.</footnote>

<footnote n="2">² Department of Mathematics and Statistics, MSU-Iligan Institute of Technology, Iligan City, Philippines.</footnote>

<footnote n="3">³ Laboratoire Amiénois de Mathématique Fondamentale et Appliquée, CNRS UMR 7352, Université de Picardie Jules Verne, 80069 Amiens, France.</footnote>

<footnote>* Corresponding author: youcef.mammeri@u-picardie.fr</footnote>
Like MERS and SARS, COVID-19 is a zoonotic virus and believed to be originated from bats transmitted to humans [19]. In comparison with SARS, MERS, the COVID-19 appears to be less deadly. However, the World Health Organization (WHO) reported that it has already infected and killed more people than its predecessors combined. Also, COVID-19 spreads much faster than SARS and MERS. It only took over a month before it surpassed the number of cases recorded by the SARS outbreak in 2012. According to WHO, it only took 67 days from the beginning of the outbreak in China last December 2019 for the virus to infect the first 100 000 people worldwide [12]. As of the 25th of March 2020, a cumulative total of 372 757 confirmed cases, while 16 231 deaths have been recorded for COVID-19 by World Health Organization [14].

Last 30th of January, WHO characterized COVID-19 as Public Health Emergency of International Concern (PHEI) and urge countries to put in place strong measures to detect disease early, isolate and treat cases, trace contacts, and promote social distancing measures commensurate with the risk [15]. In response, the world implemented its actions to reduce the spread of the virus. Limitations on mobility, social distancing, and self-quarantine have been applied. Moreover, health institutions advise people to practice good hygiene to keep from being infected. All these efforts have been made to reduce the transmission rate of the virus.

For the time being, COVID-19 infection is still on the rise. Government and research institutions scramble to seek antiviral treatment and vaccines to combat the disease. Several reports list possible drugs combination to apply, yet it is still unclear which drugs could combat the viral disease and which won't.

Several mathematical models have been proposed from various epidemiological groups. These models help governments as an early warning device about the size of the outbreak, how quickly it will spread, and how effective control measures may be. However, due to the limited emerging understanding of the new virus and its transmission mechanisms, the results are largely inconsistent across studies.

In this paper, we will mention a few models and, in the end, to propose one. Gardner and his team [7] at Center for Systems Science and Engineering, Johns Hopkins University, implemented on a previously published model applied for COVID-19. It is a metapopulation network model represented by a graph with each nodes follows a discrete-time Susceptible-Exposed-Infected-Recovered (SEIR) compartmental model. The model gives an estimate of the expected number of cases in mainland China at the end of January 2020, as well as the global distribution of infected travelers. They believe that the outbreak began in November 2019 with hundreds of infected already present in Wuhan last early December 2019. Wu et al. [16] from WHO Collaborating Centre for Infectious Disease Epidemiology and Control at the University of Hong Kong presented a modeling study on the nowcast and forecast of the 2019-nCoV outbreak at Wuhan. The group used an SEIR metapopulation model to simulate epidemic and found reproductive number $R_0$ around 2.68 (with 95% credible interval 2.47–2.86). Imai et al. [9] estimates $R_0$ around 2.6 with uncertainty range of 1.5–3.5. Zhao et al. [17, 18] found $R_0$ to range from 2.24 to 5.71 based on the reporting rate of cases. If the reporting rate increase 2-fold, $R_0$ = 3.58, if it increase 8-fold, $R_0$ = 2.24. If there is no change in the reporting rate, the estimated $R_0$ is 5.71. Similar to the above authors, Wang et al. [11] employed an infectious disease dynamics model (SEIR model) for modeling and predicting the number of COVID-19 cases in Wuhan. They opined that to reduce $R_0$ significantly, the government should continue implementing strict measures on containment and public health issues. In the same tune as the latter, the model of Danchin et al. [4] also suggests continuing to implement effective quarantine measures to avoid a resurgence of infection. The model consists of five (5) compartments: susceptible, infected, alternative infection, detected, and removed.

Here, we proposed an extension from the classical SEIR model by adding a compartment of asymptomatic infected. We address the challenge of predicting the spread of COVID-19 by giving our estimates for the basic reproductive numbers $R_0$ and its effective reproductive number $R_{eff}$. Afterward, we also assess risks and interventions via containment strategy or treatments of exposed and symptomatic infected.

The rest of the paper is organized as follows. Section 2, outlines our methodology. Here the model was explained, where the data was taken, and its parameter estimates. Section 3 contained the qualitative analysis for the model. Here, we give the closed-form equation of reproductive number $R_0$, then tackling the best strategy to reduce transmission rates. Finally, Section 4 outlines our brief discussion on some measures to limit the outbreak.
## 2. MATERIALS AND METHODS

### 2.1. Confirmed and death data

In this study, we used the publicly available dataset of COVID-19 provided by the Johns Hopkins University [6]. This dataset includes many countries' daily count of confirmed cases, recovered cases, and deaths. Data can be downloaded from https://github.com/CSSEGISandData/COVID-19/tree/master/csse_covid_19.data. These data are collected through public health authorities' announcements and are directly reported public and unidentified patient data, so ethical approval is not required.

### 2.2. Mathematical model

We focus our study on six components of the epidemic flow (Fig. 1), i.e. Susceptible individual (S), Exposed individual (E), Symptomatic Infected individual (I_s), Asymptomatic Infected individual (I_a), Under treatment individual (U), and Removed individual (R). To build the mathematical model, we followed the standard strategy developed in the literature concerning SIR model [1, 5]. We assumed that susceptible can be infected by exposed, symptomatic infected as well as asymptomatic infected individuals. While exposed individuals are generally quarantined, symptomatics are hospitalized. Due to the lack of serological tests to evaluate asymptomatic patients, no treatment is applied to them. To lighten the model, quarantine or hospitalization are not distinguished.

The dynamics is governed by a system of six ordinary differential equations (ODE) as follows

$$S'(t) = -\omega (\beta_e E + \beta_s I_s + \beta_a I_a) \frac{S}{N}$$

$$E'(t) = \omega (\beta_e E + \beta_s I_s + \beta_a I_a) \frac{S}{N} - (\delta_e + \nu_e)E$$

$$I'_s(t) = f \delta_e E - (\gamma_s + \mu_s + \nu_s) I_s$$

$$I'_a(t) = (1 - f) \delta_e E - \gamma_a I_a$$

$$U'(t) = \nu_e E + \nu_s I_s - (\gamma_u + \mu_u)U$$

$$R'(t) = \gamma_s I_s + \gamma_a I_a + \gamma_u U$$

Note that the total living population follows $N'(t) = -\mu_s I_s - \mu_u U$, while death is computed by $D'(t) = \mu_s I_s + \mu_u U$. We assume that there is no new recruit. The parameters are described in Table A of Figure 2.

### 2.3. Parameters estimation

Calibration is made before intervention. Thus it is set $\nu_e = \nu_s = 0$. The model is made up of eleven parameters $\theta = (\omega, \beta_e, \beta_s, \beta_a, \delta_e, f, \gamma_s, \mu_s, \gamma_a, \gamma_u, \mu_u)$ that need to be determined. Given, for n days, the observations $I_{s,obs}(t_i)$ and $D_{obs}(t_i)$, the cost function consists of a nonlinear least square function

$$J(\theta) = \sum_{i=1}^{n} (I_{s,obs}(t_i) - I_s(t_i, \theta))^2 + (D_{obs}(t_i) - D(t_i, \theta))^2,$$

with constraints $\theta \ge 0$, and $0 \le f \le 1$. Here $I_s(t_i, \theta)$ and $D(t_i, \theta)$ denote output of the mathematical model at time $t_i$ computed with the parameters $\theta$. The optimization problem is solved using Approximate Bayesian Computation combined with a quasi-Newton method [3].
<figure><img src="image_2.png" /><figcaption>FIGURE 1. Compartmental representation of the SEI<sub>a</sub>I<sub>s</sub>UR-model. Blue arrows represent the infection flow. Green arrows denote for the treatments. Purple arrow is the death.</figcaption></figure>

3. RESULTS

3.1. Basic and effective reproduction numbers

It is standard to check that the domain

$$
\Omega = \{(S, E, I_s, I_a, U, R) \in \mathbb{R}_+^6; 0 \le S + E + I_s + I_a + U + R \le N(0)\}
$$

is positively invariant. In particular, there exists a unique global in time solution $(S, E, I_s, I_a, U, R)$ in $C(\mathbb{R}_+; \Omega)$
as soon as the initial condition lives in $\Omega$.

Since the infected individuals are in E, Iₐ and Iₛ, the rate of new infections in each compartment (F) and
the rate of other transitions between compartments (V) can be rewritten as

$$
\mathcal{F} = \begin{pmatrix} \omega(\beta_e E + \beta_s I_s + \beta_a I_a) \frac{S}{N} \\ 0 \\ 0 \end{pmatrix}, \quad \mathcal{V} = \begin{pmatrix} (\delta_e + \nu_e)E \\ (\gamma_s + \mu_s + \nu_s)I_s - f\delta_e E \\ \gamma_a I_a - (1-f)\delta_e E \end{pmatrix}.
$$

Thus,

$$
F = \begin{pmatrix} \frac{\omega \beta_e S}{N} & \frac{\omega \beta_s S}{N} & \frac{\omega \beta_a S}{N} \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix},
$$

and

$$
V = \begin{pmatrix} (\delta_e + \nu_e) & 0 & 0 \\ -f\delta_e & \gamma_s + \mu_s + \nu_s & 0 \\ -(1-f)\delta_e & 0 & \gamma_a \end{pmatrix}, \quad V^{-1} = \begin{pmatrix} \frac{1}{\delta_e + \nu_e} & 0 & 0 \\ \frac{f}{\gamma_s + \mu_s + \nu_s} & \frac{1}{\gamma_s + \mu_s + \nu_s} & 0 \\ \frac{1}{\gamma_a} & 0 & \frac{1}{\gamma_a} \end{pmatrix}.
$$

Therefore, the next generation matrix is

$$
FV^{-1} = \begin{pmatrix} \frac{\omega\beta_e S}{(\delta_e + \nu_e)N} + \frac{f\omega\beta_s S}{(\gamma_s + \mu_s + \nu_s)N} + \frac{(1-f)\omega\beta_a S}{\gamma_a N} & \frac{\omega\beta_s S}{(\gamma_s + \mu_s + \nu_s)N} & \frac{\omega\beta_a S}{\gamma_a N} \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}.
$$
A.

<table><thead><tr><th>country/state</th><th>description</th><th>France</th><th>Philippines</th><th>Italy</th><th>Spain</th><th>United Kingdom</th><th>Hubei</th><th>New York</th></tr></thead><tbody><tr><td>J</td><td>relative cost function</td><td>4.70E-03</td><td>2.02E-02</td><td>3.07E-02</td><td>5.94E-03</td><td>4.27E-03</td><td>7.25E-02</td><td>7.96E-04</td></tr><tr><td>ωβ<sub>e</sub></td><td>transmission rate from S to E from contact with E</td><td>0.057899</td><td>0.223627</td><td>0.251883</td><td>0.273061</td><td>0.237637</td><td>0.000000</td><td>0.277300</td></tr><tr><td>ωβ<sub>s</sub></td><td>transmission rate from S to E from contact with I<sub>s</sub></td><td>0.368146</td><td>0.399019</td><td>0.675281</td><td>0.503112</td><td>0.979359</td><td>0.507331</td><td>0.002878</td></tr><tr><td>ωβ<sub>a</sub></td><td>transmission rate from S to E from contact with I<sub>a</sub></td><td>0.471628</td><td>0.186464</td><td>0.000000</td><td>0.000432</td><td>0.598530</td><td>0.000000</td><td>4.263184</td></tr><tr><td>δ<sub>e</sub></td><td>latency rate</td><td>0.269925</td><td>0.143335</td><td>2.872143</td><td>0.572886</td><td>0.183673</td><td>3.445527</td><td>0.135697</td></tr><tr><td>f</td><td>probability of be asymptomatic</td><td>0.387702</td><td>0.084839</td><td>0.426419</td><td>0.490394</td><td>0.051886</td><td>0.598110</td><td>0.159819</td></tr><tr><td><b>1-f</b></td><td><b>probability of be asymptomatic</b></td><td><b>0.612298</b></td><td><b>0.915161</b></td><td><b>0.573581</b></td><td><b>0.509606</b></td><td><b>0.948114</b></td><td><b>0.401890</b></td><td><b>0.840181</b></td></tr><tr><td>γ<sub>s</sub>, γ<sub>a</sub>, γ<sub>u</sub></td><td>recovery rates</td><td>0.103809</td><td>0.122703</td><td>0.063046</td><td>0.044322</td><td>0.383395</td><td>0.033333</td><td>1.787245</td></tr><tr><td>μ<sub>s</sub></td><td>death rate</td><td>0.007093</td><td>0.015090</td><td>0.020413</td><td>0.013520</td><td>0.012764</td><td>0.009753</td><td>0.005178</td></tr><tr><td>R</td><td>effective reproduction number</td><td>4.283329</td><td>3.196568</td><td>3.537937</td><td>4.747048</td><td>2.902201</td><td>7.042639</td><td>4.047891</td></tr><tr><td>median R</td><td>median effective reproduction number</td><td>3.096738</td><td>3.155502</td><td>2.089169</td><td>2.902238</td><td>2.868001</td><td>7.405132</td><td>3.781538</td></tr><tr><td>mean R</td><td>mean effective reproduction number</td><td>3.474858</td><td>3.779976</td><td>3.542378</td><td>3.439437</td><td>3.939954</td><td>7.633783</td><td>4.111556</td></tr></tbody></table>

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>FIGURE 2. A. Parameters calibrated according to data from France, Philippines, Italy, Spain, United Kingdom, Hubei, and New York. B. and C. Calibrated solution (straight line) and data (dots) with respect to day for France, Philippines, Italy, Spain, United Kingdom, Hubei and New York. First is the infected $I_s$ (B.), and the second one is the death $D$ (C.). D. Effective reproduction number with respect to day.</figcaption></figure>

We deduce as in [10] that the basic reproduction number $R_0$ for the Disease Free Equilibrium ($S^*, 0, 0, 0, R^*$), with $N^* = S^* + R^*$, is

$$R_0 := \omega \left( \frac{\beta_e}{\delta_e + \nu_e} + \frac{f \beta_s}{\gamma_s + \mu_s + \nu_s} + \frac{(1 - f) \beta_a}{\gamma_a} \right) \frac{S^*}{N^*} =: R \frac{S^*}{N^*}.$$

This number has an epidemiological meaning. The term $\frac{\omega \beta_e}{\delta_e + \nu_e}$ represents the contact rate with exposed during the average latency period $1/(\delta_e + \nu_e)$. The term $\frac{\omega f \beta_s}{\gamma_s + \mu_s + \nu_s}$ is the contact rate with symptomatic during the average infection period, and the last one is the part of asymptomatic.

In the subsequent, we write DFE when we mean by Disease Free Equilibrium.

**Theorem 3.1.** The DFE ($S^*, 0, 0, 0, R^*$) is the unique positive equilibrium. Moreover it is globally asymptotically stable.

*Proof.* By computing the eigenvalues of the Jacobian matrix, we deduce that if $R_0 < 1$, then DFE is locally asymptotically stable.

Here we will prove that global asymptotic stability is independent that of $R_0$. Indeed, from the last differential equation in our system of ODE, we can deduce that $R$ is an increasing function bounded by $N(0)$. Thus $R(t)$
<figure><img src="image_6.png" /><figcaption>FIGURE 3. A. Boxplot of the posterior distribution computed from France data. B. Effective reproduction number in grey of the posterior distribution, median (=3.096738) in blue straight line, mean (=3.474858) is dotted line. C. Fitted symptomatic infected in grey of the posterior distribution, median in red straight line, mean is dotted line. D. Fitted death in grey of the posterior distribution, median in black straight line, mean is dotted line.</figcaption></figure>

converges to $R^*$ as $t$ goes to $+\infty$. Then integrating over time this equation provides

$$R(t) - R(0) = \int_{0}^{t} \gamma_{s}I_{s}(s) + \gamma_{a}I_{a}(s) + \gamma_{u}U(s) \, ds$$

and

$$R^* - R(0) = \int_{0}^{+\infty} \gamma_{s}I_{s}(s) + \gamma_{a}I_{a}(s) + \gamma_{u}U(s) \, ds,$$

which is finite. Furthermore, $\gamma_s I_s + \gamma_a I_a + \gamma_u U$ goes to 0 as $t \to +\infty$, and each term of this sum does thanks to the positivity of the solution. Adding the two first equations implies that

$$(S + E)' = -(\delta_e + \nu_e)E$$
<table><caption>TABLE 1. Comparison between the maximum number of symptomatic infected and death without control and the solution reducing contact rate to 0, 100% of exposed under treatment, and 100% of symptomatic infected under treatment to reach $T_c = 1$ and $T_c = 1000$. Interventions are assumed to being 53 days after the first confirmed infection.</caption><thead><tr><th>country/state</th><th>France</th><th>Philippines</th><th>Italy</th><th>Spain</th><th>United Kingdom</th><th>Hubei</th><th>New York</th></tr></thead><tbody><tr><td>population</td><td>6.67E+07</td><td>1.03E+08</td><td>6.03E+07</td><td>4.63E+07</td><td>6.49E+07</td><td>5.84E+07</td><td>1.95E+07</td></tr><tr><td>maximum number of symptomatic infected without control</td><td>7.74E+06</td><td>1.75E+06</td><td>9.43E+06</td><td>1.05E+07</td><td>3.49E+05</td><td>2.02E+07</td><td>9.25E+04</td></tr><tr><td>maximum number of dead without control</td><td>1.64E+06</td><td>9.15E+05</td><td>6.15E+06</td><td>5.29E+06</td><td>1.02E+05</td><td>7.92E+06</td><td>8.85E+03</td></tr><tr><td>intervention duration if T'<sub>c</sub> = 1 and full containment starting 53 days after the first confirmed infection</td><td>104</td><td>93</td><td>194</td><td>262</td><td>78</td><td>407</td><td>119</td></tr><tr><td>maximum number of symptomatic infected with full containment</td><td>8.78E+03</td><td>5.75E+02</td><td>1.19E+05</td><td>7.40E+04</td><td>1.16E+04</td><td>2.00E+07</td><td>9.22E+04</td></tr><tr><td>maximum number of dead with full containment</td><td>1.24E+03</td><td>1.94E+02</td><td>4.24E+04</td><td>2.52E+04</td><td>1.96E+03</td><td>6.53E+06</td><td>6.82E+03</td></tr><tr><td>intervention duration if T'<sub>c</sub> = 1000 and full containment starting 53 days after the first confirmed infection</td><td>39</td><td>33</td><td>86</td><td>110</td><td>40</td><td>292</td><td>68</td></tr><tr><td>intervention duration if T'<sub>c</sub> = 1000 and 100% of exposed under treatment starting 53 days after the first confirmed infection</td><td>44</td><td>39</td><td>87</td><td>118</td><td>44</td><td>342</td><td>82</td></tr><tr><td>intervention duration if T'<sub>c</sub> = 1000 and 100% of symptomatic under treatment starting 53 days after the first confirmed infection</td><td>163</td><td>188</td><td>202</td><td>263</td><td>103</td><td>408</td><td>122</td></tr></tbody></table>

and $S + E$ is a nonnegative decreasing function whose derivative tends to zero. Then $E(t) \to_{t \to +\infty} 0$ and $S(t) \to_{t \to +\infty} S^*$. □

This theorem means that the asymptotic behavior does not depend on $\mathcal{R}_0$. For all initial data in $\Omega$, the solution converges to the DFE when time goes to infinity. Nevertheless, to observe initial exponential growth, $\mathcal{R}_0 > 1$ is necessary. Indeed, $S$ is initially close to $N$ such that infected states are given by the linear system of differential equations

$$ \begin{pmatrix} E \\ I_s \\ I_a \end{pmatrix}'(t) = \begin{pmatrix} \omega\beta_e - \delta_e - \nu_s & \omega\beta_s & \omega\beta_a \\ f\delta_e & -(\gamma_s + \mu_s + \nu_s) & 0 \\ (1-f)\delta_e & 0 & -\gamma_a \end{pmatrix} \begin{pmatrix} E \\ I_s \\ I_a \end{pmatrix}. $$

The characteristic polynomial is $P(x) = x^3 + a_2x^2 + a_1x + a_0$, with $a_0 = (\gamma_s + \mu_s + \nu_s)\gamma_a(\delta_e + \nu_s)(1 - \mathcal{R})$. If $\mathcal{R} > 1$, there is at least one positive eigenvalue that coincides with an initial exponential growth rate of solutions.

To better reflect the time dynamic of the disease, the effective reproduction number

$$ \mathcal{R}_{\text{eff}}(t) = \omega \left( \frac{\beta_e}{\delta_e + \nu_s} + \frac{f\beta_s}{\gamma_s + \mu_s + \nu_s} + \frac{(1-f)\beta_a}{\gamma_a} \right) \frac{S(t)}{N(t)}. $$

is represented in Figure 2D and values of $\mathcal{R}$ are computed in Table A of Figure 2.
France

<figure><img src="image_7.png" /><figcaption>FIGURE 4. Comparison of solutions S, E, I_s, I_a, R, D without control in blue, full containment in green, full treatment of symptomatic in red, and full treatment of exposed in cyan for France. Ordinate axis is expressed in log.</figcaption></figure>

**3.2. Model resolution**

To calibrate the model, our simulations start the day of first confirmed infection and finish before interventions to reduce the disease. Therefore $\nu_e$ and $\nu_s$ are assumed to be equal to 0. We assume that the whole population of the country is susceptible to the infection. Seven states with comparable populations are chosen. The objective function $J$ is computed to provide a relative error of order less than $10^{-2}$. In Figure 2, Table A. shows estimated parameters. The rest of the Figure presents the solution and data. Note that the product $\omega\beta_i$ ($i = e, s, a$) is uniquely identifiable but not $\omega$ and $\beta_i$ separately. Figure 3 represents the effective reproduction number, the fitted symptomatic infected and death of the posterior distribution.

**3.3. Strategy to reduce disease to a given threshold**

To temporarily reduce the value of $R_{\text{eff}}$, three strategies regarding different interventions are being compared. The first one consists in reducing the number of contacts $\omega$. Full containment is translated by $\omega = 0$, and no containment by $\omega = 1$. The second strategy is expressed by treatment of symptomatic infected. It is translated by modifying the value of $\nu_s$. The treatment of exposed is the third one. Parameters $\omega, \nu_e$, and $\nu_s$ vary such that intervention is carried out from 0 to 100% of susceptible, exposed, and symptomatic respectively. Given a
Philippines

<figure><img src="image_8.png" /><figcaption>FIGURE 5. Comparison of solutions S, E, I_s, I_a, R, D without control in blue, full containment in green, full treatment of symptomatic in red, and full treatment of exposed in cyan for the Philippines. Ordinate axis is expressed in log.</figcaption></figure>

critical infection threshold $T_c$, the model runs until the time $t_c$ such that

$$E(t_c) + I_s(t_c) + I_a(t_c) \leq T_c.$$

To juxtapose the benefit of the intervention, we assume that interventions start 53 days after the first confirmed infection. We remind that the first infection in France was confirmed on January 24th, 2020, and containment begins on March 17th, 2020. Comparison between three strategies can be found in Table 1 and Figures 4–5. Without intervention to control the disease, the maximum number of symptomatic infecteds varies from $3.49 \times 10^5$ to $2.02 \times 10^7$. The maximum number of deaths totals from $8.85 \times 10^3$ to $7.92 \times 10^6$. We observe that any intervention strongly reduces the number of dead. Concerning France, Philippines, Italy, Spain, and the United Kingdom, when containment is fully respected and when the sum of infecteds is reduced to 1, the maximum number of symptomatic infecteds and deaths has been cut sharply, of order $10^3$. It varies now from $5.75 \times 10^2$ to $7.04 \times 10^4$ and $1.94 \times 10^2$ to $2.52 \times 10^4$ respectively. To wait from 104 to 407 days is the price to pay. On the contrary, for the states of Hubei and New York, 53 days to intervene seems to be already too late. We can also see in Table 1 that treating only the symptomatic does not reduce the duration.

Note that when the intervention ends at time $t_c$, the number of susceptible $S(t_c)$ is large so that the effective reproduction number $R_{eff}$ is larger than 1.
## France

A.

<figure><img src="image_9.png" /></figure>

B.

<figure><img src="image_10.png" /></figure>

C.

<figure><img src="image_11.png" /></figure>

## Philippines

A.

<figure><img src="image_12.png" /></figure>

B.

<figure><img src="image_13.png" /></figure>

C.

<figure><img src="image_14.png" /></figure>

FIGURE 6. A. Comparison of the maximum number of dead. B. Comparison of the maximum number of symptomatic infected. C. Comparison of the intervention duration to reach $T_c$ = 1000 with respect to percentage of containment (green), treatment of symptomatic (red), and treatment of exposed (cyan) for France and Philippines.

Figure 6 compares the maximum number of dead and symptomatic but infected individuals, as well as the intervention duration, to reach $T_c$ = 1000 varying from 0 to 100%. Containment is the most efficient when it is respected by more than 76% in France, 63% in the Philippines. Beyond that, treating the exposed is the best choice. We also observe that the intervention duration becomes long below 89% in France, 82% in the Philippines. This can be understood by too little susceptibility to achieve recovery but enough for the disease to persist.

## 4. DISCUSSION

Without intervention, we observe in Figures 4-5 that the number of susceptible $S$ is decreasing; most of the individuals are recovering, which generates population immunity. It translates that the disease free equilibrium is globally asymptotically stable. Nevertheless, the price to pay is high, the number of deaths being excessive. As presented in Figure 2, the effective reproduction is decreasing and points out that control has to be done as fast as possible.

The other important information is that, as discovered by Danchin et al. [4], an alternative transmission way may occur. Here, it is due to the proportion of asymptomatic infected individuals that is not negligible, as shown in Table 1.

Finally, with the little knowledge about COVID-19 nowadays, decreasing transmission, i.e. $\beta_e, \beta_s, \beta_a$, is the preferred option. The simplest choice consists in reducing contact between individuals. Table 1 and Figures 4-6 show that total and partial containment do indeed drastically reduce the disease. However, the
duration of containment may be too long and then impracticable especially if we aim at totally eradicating the infection ($T_c = 1$). Instead, to stop the containment as soon as the capacity of the hospitals has been reached could be privileged. When this criterion is set to 1000 patients ($T_c = 1000$), the duration goes from 104 to 39 days for France. A similar reduction in duration is also obtained for other countries. Again, we see that the earlier the intervention, the more effective it is. Due to the high number of susceptible, it is worth noting that the effective reproduction number remains large after containment. Screening tests, especially to carry out exposed individuals, are then necessary to be carried out, and the positive individuals are quarantined.

***Acknowledgements.*** The authors would like to thank the anonymous reviewers for their insightful comments and their prompt reply.

## REFERENCES

[1] F. Brauer and C. Castillo-Chavez, Mathematical Models in Population Biology and Epidemiology. Springer, New York, NY (2012).

[2] N. Chen, M. Zhou, X. Dong, J. Qu, F. Gong, Y. Han, Y. Qiu, J. Wang, Y. Liu, Y. Wei, J. Xia, T. Yu, X. Zhang and L. Zhang, Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in wuhan, china: a descriptive study. *The Lancet* **395** (2020) 507–513.

[3] K. Csilléry, M.G. Blum, O.E. Gaggiotti and O. François, Approximate bayesian computation (abc) in practice. *Trends Ecol. Evol.* **25** (2010) 410–418.

[4] A. Danchin, T.W.P. Ng and G. Turinici, A new transmission route for the propagation of the sars-cov-2 coronavirus. *medRxiv*, 2020.

[5] O. Diekmann and J.A.P. Heesterbeek, Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation. John Wiley, Chichester (2000).

[6] E. Dong, H. Du and L. Gardner, An interactive web-based dashboard to track covid-19 in real time. *The Lancet Infectious Diseases* **20** (2020) 533–534.

[7] L. Gardner, Modeling the spreading risk of 2019-ncov. Technical report, Center for Systems Science and Engineering, Johns Hopkins University (2020).

[8] D.S. Hui, E.I. Azhar, T.A. Madani, F. Ntoumi, R. Kock, O. Dar, G. Ippolito, T.D. Mchugh, Z.A. Memish, C. Drosten, A. Zumla and E. Petersen. The continuing 2019-ncov epidemic threat of novel coronaviruses to global health - the latest 2019 novel coronavirus outbreak in wuhan, china. *Int. J. Infet. Dis.* **91** (2020) 264–266.

[9] N. Imai, I. Dorigatti, A. Cori, C. Donnelly, S. Riley and N.M. Ferguson, Report 2: Estimating the potential total number of novel coronavirus cases in wuhan city, china. Technical report, Imperial College London (2020).

[10] P. van den Driessche and J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. *Math. Biosci.* **180** (2000) 29–48.

[11] H. Wang, Z. Wang, Y. Dong et al., Phase-adjusted estimation of the number of coronavirus disease 2019 cases in wuhan, china. *Cell Discovery* **6** (2020) 10.

[12] World Health Organization, COVID-19 (2020).

[13] World Health Organization, Naming the coronavirus disease (COVID-19) and the virus that causes it (2020).

[14] World Health Organization, Novel Coronavirus (2019-nCoV) Situation. Report 64 (2020).

[15] World Health Organization, Statement on the Second Meeting of the International Health Regulations. (2005). Emergency Committee Regarding the Outbreak of Novel Coronavirus (2019-nCoV) (2020).

[16] J.T. Wu, K. Leung and G.M. Leung, Nowcasting and forecasting the potential domestic and international spread of the 2019-ncov outbreak originating in wuhan, china: A modelling study. *The Lancet* **395** (2020) 689–697.

[17] S. Zhao, Q. Lin, J. Musa Ran, S.S. Yang, G Wang., W. Lou, D.Y. Gao, L. Yang, D. He and M.H. Wang, 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. *Int. Journal Infect. Dis.* **92** (2020) 214–217.

[18] S. Zhao, S.S. Musa, Q. Lin, J. Ran, G. Yang, W. Wang, Y. Lou, et al., Estimating the unreported number of novel coronavirus (2019-ncov) cases in china in the first half of january 2020: A data-driven modelling analysis of the early outbreak. *J. Clin. Med.* **9** (2020) 388.

[19] P. Zhou, X. Yang, X. Wang, et al., A pneumonia outbreak associated with a new coronavirus of probable bat origin. *Nature* **579** (2020) 270–273.