A Time-dependent SIR model for COVID-19 with Undetectable Infected Persons

2020·Arxiv

Abstract

Abstract

In this paper, we conduct mathematical and numerical analyses to address the following important questions for COVID-19: (Q1) Is it possible to contain COVID-19? (Q2) If COVID-19 can be contained, when will be the peak of the epidemic, and when will it end? (Q3) How do the asymptomatic infections affect the spread of disease? (Q4) If COVID-19 cannot be contained, what is the ratio of the population that needs to be infected in order to achieve herd immunity? (Q5) How effective are the social distancing approaches? (Q6) If COVID-19 cannot be contained, what is the ratio of the population infected in the long run? For (Q1) and (Q2), we propose a time-dependent susceptible-infected-recovered (SIR) model that tracks two time series: (i) the transmission rate at time t and (ii) the recovering rate at time t. Such an approach is not only more adaptive than traditional static SIR models, but also more robust than direct estimation methods. Using the data provided by the National Health Commission of the People’s Republic of China (NHC) [1], we show that the one-day prediction errors for the numbers of confirmed cases are almost less than 3%. Also, the turning point, defined as the day that the transmission rate is less than the recovering rate, is predicted to be Feb. 17, 2020. After that day, the basic reproduction number, known as the value at time t, is less than 1. In that case, the total number of confirmed cases is predicted to be around 80, 000 cases in China under our model. For (Q3), we extend our SIR model by considering two types of infected persons: detectable infected persons and undetectable infected persons. Whether there is an outbreak in such a model is characterized by the spectral radius of a matrix that is closely related to the basic reproduction number . We plot the phase transition diagram of an outbreak and show that there are several countries, including South Korea, Italy, and Iran, that are on the verge of COVID-19 outbreaks on Mar. 2, 2020. For (Q4), we show that herd immunity can be achieved after at least fraction of individuals being infected and recovered from COVID-19. For (Q5) and (Q6), we analyze the independent cascade (IC) model for disease propagation in a random network specified by a degree distribution. By relating the propagation probabilities in the IC model to the transmission rates and recovering rates in the SIR model, we show two approaches of social distancing that can lead to a reduction of

Index Terms—COVID-19, SARS-CoV-2, 2019-nCoV, Coronavirus, Time-dependent SIR model, asymptomatic infection, herd immunity, superspreader, independent cascade, social distancing.

I. INTRODUCTION

At the beginning of December 2019, the first COVID-19 victim was diagnosed with the coronavirus in Wuhan, China. In the following weeks, the disease spread widely in China mainland and other countries, which causes global panic. The virus has been named “SARS-CoV-2,” and the disease it causes has been named “coronavirus disease 2019 (abbreviated “COVID-19”). There have been 80, 151 people infected by the disease and 2, 943 deaths until Mar. 2, 2020 according to the official statement by the Chinese government. To block the spread of the virus, there are some strategies such as city-wide lockdown, traffic halt, community management, social distancing, and propaganda of health education knowledge that have been adopted by the governments of China and other countries in the world.

Unlike the Severe Acute Respiratory Syndrome (SARS) and other infectious diseases, one problematic characteristic of COVID-19 is that there are asymptomatic infections (who have very mild symptoms). Those asymptomatic infections are unaware of their contagious ability, and thus get more people infected. The transmission rate can increase dramatically in this circumstance. According to the recent report from WHO [2], only 87.9% of COVID-19 patients have a fever, and 67.7% of them have a dry cough. If we use body temperature as a means to detect COVID-19 infected cases, then more than 10% of infected persons cannot be detected.

Due to the recent development of the epidemic, we are interested in addressing the following important questions for COVID-19:

(Q1) Is it possible to contain COVID-19? Are the commonly used measures, such as city-wide lockdown, traffic halt, community management, and propaganda of health education knowledge, effective in containing COVID-19?

(Q2) If COVID-19 can be contained, when will be the peak of the epidemic, and when will it end?

(Q3) How do the asymptomatic infections affect the spread of disease?

(Q4) If COVID-19 cannot be contained, what is the ratio of the population that needs to be infected in order to achieve herd immunity?

(Q5) How effective are the social distancing approaches, such as reduction of interpersonal contacts and canceling mass gatherings in controlling COVID-19?

(Q6) If COVID-19 cannot be contained, what is the ratio of the population infected in the long run?

For (Q1), we analyze the cases in China and aim to predict how the virus spreads in this paper. Specifically, we propose using a time-dependent susceptible-infected-recovered (SIR) model to analyze and predict the number of infected persons and the number of recovered persons (including deaths). In the traditional SIR model, it has two time-invariant variables: the transmission rate and the recovering rate . The transmission rate means that each individual has on average contacts with randomly chosen others per unit time. On the other hand, the recovering rate indicates that individuals in the infected state get recovered or die at a fixed average rate . The traditional SIR model neglects the time-varying property of and , and it is too simple to precisely and effectively predict the trend of the disease. Therefore, we propose using a time-dependent SIR model, where both the transmission rate and the recovering rate are functions of time t. Our idea is to use machine learning methods to track the transmission rate and the recovering rate , and then use them to predict the number of the infected persons and the number of recovered persons at a certain time t in the future. Our time-dependent SIR model can dynamically adjust the crucial parameters, such as and , to adapt accordingly to the change of control policies, which differs from the existing SIR and SEIR models in the literature, e.g., [3], [4], [5], [6], and [7]. For example, we observe that city-wide lockdown can lower the transmission rate substantially from our model. Most data-driven and curve-fitting methods for the prediction of COVID-19, e.g., [8], [9], and [10] seem to track data perfectly; however, they are lack of physical insights of the spread of the disease. Moreover, they are very sensitive to the sudden change in the definition of confirmed cases on Feb. 12, 2020 in the Hubei province. On the other hand, our time-dependent SIR model can examine the epidemic control policy of the Chinese government and provide reasonable explanations. Using the data provided by the National Health Commission of the People’s Republic of China (NHC) [1], we show that the one-day prediction errors for the numbers of confirmed cases are almost less than 3% except for Feb. 12, 2020, which is unpredictable due to the change of the definition of confirmed cases.

For (Q2), the basic reproduction number , defined as the number of additional infections by an infected person before it recovers, is one of the commonly used metrics to check whether the disease will become an outbreak. In the classical SIR model, is simply as an infected person takes (on average) days to recover, and during that period time, it will be in contact with (on average) persons. In our time-dependent SIR model, the basic reproduction number is a function of time, and it is defined as . If , the disease will spread exponentially and infects a certain fraction of the total population. On the contrary, the disease will eventually be contained. Therefore, by observing the change of with respect to time or even predict in the future, we can check whether certain epidemic control policies are effective or not. Using the data provided by the National Health Commission of the People’s Republic of China (NHC) [1], we show that the turning point (peak), defined as the day that the basic reproduction number is less than 1, is predicted to be Feb. 17, 2020. Moreover, the disease in China will end in about 6 weeks after its peak in our (deterministic) model if the current contagious disease control policies are maintained in China. In that case, the total number of confirmed cases is predicted to be around 80, 000 cases in China under our (deterministic) model.

For (Q3), we extend our SIR model to include two types of infected persons: detectable infected persons (type I) and undetectable infected persons (type II). With probability (resp. ), an infected person is of type I (resp. II), where . Type I (resp. II) infected persons have the transmission rate (resp. ) and the recovering rate (resp. ). The basic reproduction number in this model is

In practice, type I infected persons have a lower transmission rate than that of type II infected persons (as type I infected persons can be isolated). For such a model, whether the disease is controllable is characterized by the spectral radius of a matrix. If the spectral radius of that matrix is larger than 1, then there is an outbreak. On the other hand, if it is smaller than 1, then there is no outbreak. One interesting result is that the spectral radius of that matrix is larger (resp. smaller) than 1 if the basic reproduction number in (1) is larger (resp. smaller) than 1. The curve that has the spectral radius equal to 1 is known as the percolation threshold curve in a phase transition diagram [11]. Using the historical data from Jan. 22, 2020 to Mar. 2, 2020 from the GitHub of Johns Hopkins University [12], we extend our study to some other countries, including Japan, Singapore, South Korea, Italy, and Iran. Our numerical results show that there are several countries, including South Korea, Italy, and Iran, that are above the percolation threshold curve, and they are on the verge of COVID-19 outbreaks on Mar. 2, 2020.

The British prime minister, Boris Johnson, once suggested having a sufficiently high fraction of individuals infected by COVID-19 and recovered from the disease to achieve herd immunity. To address the question in (Q4), we argue that herd immunity corresponds to the reduction of the number of susceptible persons in the SIR model, and herd immunity can be achieved after at least fraction of individuals being infected and recovered from the COVID-19.

For (Q5), we consider two commonly used approaches for social distancing: (i) allowing every person to keep its interpersonal contacts up to a fraction of its normal contacts, and (ii) canceling mass gatherings. For the analysis of social distancing, we have to take the social network (and its network structure) into account. For this, we consider the independent cascade (IC) model for disease propagation in a random network specified by a degree distribution . The IC model has been widely used for the study of the influence maximization problem in viral marketing (see, e.g., [13]). In the IC model, an infected node can transmit the disease to a neighboring susceptible node (through an edge) with a certain propagation probability. Repeatedly continuing the propagation, we have a subgraph that contains the set of infected nodes in the long run. By relating the propagation probabilities in the IC model to the transmission rates and recovering rates in the SIR model, we show two results for social distancing: (i) for the social distancing approach that allows every person to keep its interpersonal contacts up to (on average) a fraction a of its normal contacts, the basic reproduction number is reduced by a factor of , and (ii) for the social distancing approach that cancels mass gatherings by removing nodes with the number of edges larger than or equal to , the basic reproduction number is reduced by a factor of

For (Q6), there is a piece of solid evidence for an outbreak in the State of New York when Andrew Cuomo, the governor of New York State, said on Apr. 23, 2020, that 13.9% of a group of 3, 000 people tested positive for COVID-19 antibodies. In the SIR model with a stationary transmission rate and a stationary recovering rate, it is well-known (see, e.g., [11]) that the ratio of the population infected in the long run, denoted by r, can be computed from the fixed point equation:

However, the SIR model does not take the network structure into account. To see the effect of the degree distribution to the ratio of the population infected in the long run, we consider the IC model for disease propagation in a random network generated by the configuration model with the degree distribution . We show that if , then a certain proportion of the population will be infected. Moreover, the ratio of the population infected in the long run can also be computed from a fixed point equation. When the degree distribution is a Poisson degree distribution, it reduces to (2). Our numerical results show that (2) is a conservative estimate of the ratio of the population infected in the long run in comparison with real networks that have power-law degree distributions.

In Table I, we provide a list of notations that are used in this paper.

The rest of the paper is organized as follows: In Section II, we propose the time-dependent SIR model. We then extend the model to the SIR model with undetectable infected persons in Section III. In Section IV, we consider the independent cascade model for disease propagation in a random network specified by a degree distribution. In Section V, we conduct several numerical experiments to illustrate the effectiveness of our models. In Section VI, we put forward some discussions and suggestions to control COVID-19. The paper is concluded in Section VII.

Table I: List of notations

II. THE TIME-DEPENDENT SIR MODEL

A. Susceptible-infected-recovered (SIR) Model

In the typical mathematical model of infectious disease, one often simplify the virus-host interaction and the evolution of an epidemic into a few basic disease states. One of the simplest epidemic model, known as the susceptible-infected-recovered (SIR) model [11], includes three states: the susceptible state, the infected state, and the recovered state. An individual in the susceptible state is one who does not have the disease at time t yet, but may be infected if one is in contact with a person infected with the disease. The infected state refers to an individual who has a disease at time t and may infect a susceptible individual potentially (if they come into contact with each other). The recovered state refers to an individual who is either recovered or dead from the disease and is no longer contagious at time t. Also, a recovered individual will not be back to the susceptible state anymore. The reason for the number of deaths is counted in the recovered state is that, from an epidemiological point of view, this is basically the same thing, regardless of whether recovery or death does not have much impact on the spread of the disease. As such, they can be effectively eliminated from the potential host of the disease [14]. Denote by S(t), X(t) and R(t) the numbers of susceptible persons, infected persons, and recovered persons at time t. Summing up the above SIR model, we believe it is very similar to the COVID-19 outbreak, and we will adopt the SIR model as our basic model in this paper.

In the traditional SIR model, it has two time-invariant variables: the transmission rate and the recovering rate . The transmission rate means that each individual has on average contacts with randomly chosen others per unit time. On the other hand, the recovering rate indicates that individuals in the infected state get recovered or die at a fixed average rate . The traditional SIR model neglects the time-varying property of and . This assumption is too simple to precisely and effectively predict the trend of the disease. Therefore, we propose the time-dependent SIR model, where both the transmission rate and the recovering rate are functions of time t. Such a time-dependent SIR model is much better to track the disease spread, control, and predict the future trend.

B. Differential Equations for the Time-dependent SIR Model

For the traditional SIR model, the three variables S(t), X(t) and R(t) are governed by the following differential equations (see, e.g., the book [11]):

We note that

where n is the total population. Let and be transmission rate and recovering rate at time t. Replacing and by and in the differential equations above yields

The three variables S(t), X(t) and R(t) still satisfy (3). Now we briefly explain the intuition of these three equations. Equation (4) describes the difference of the number of

susceptible persons S(t) at time t. If we assume the total population is n, then the probability that a randomly chosen person is in the susceptible state is S(t)/n. Hence, an individual in the infected state will contact (on average) people in the susceptible state per unit time, which implies the number of newly infected persons is (as there are X(t) people in the infected state at time t). On the contrary, the number of people in the susceptible state will decrease by . Additionally, as every individual in the infected state will recover with rate , there are (on average) people recovered at time t. This is shown in (6) that illustrates the difference of R(t) at time t. Since three variables S(t), X(t) and R(t) still satisfy (3), we have

which is the number of people changing from the susceptible state to the infected state minus the number of people changing from the infected state to the recovered state (see (5)).

C. Discrete Time Time-dependent SIR Model

Due to the COVID-19 data is updated in days [1], we revise the differential equations in (4), (5), and (6) into discrete time difference equations:

Again, the three variables S(t), X(t) and R(t) still satisfy (3).

In the beginning of the disease spread, the number of confirmed cases is very low, and most of the population are in the susceptible state. Hence, for our analysis of the initial stage of COVID-19, we assume , and further simplify (8) as follows:

From the difference equations above, one can easily derive and of each day. From (9), we have

Given the historical data from a certain period , we can measure the corresponding by using (11) and (12). With the above information, we can use machine learning methods to predict the time varying transmission rates and recovering rates.

D. Tracking Transmission Rate β(t) and Recovering Rate γ(t) by Ridge Regression

In this subsection, we track and predict and by the commonly used Finite Impulse Response (FIR) filters in linear systems. Denote by and the predicted transmission rate and recovering rate. From the FIR filters, they are predicted as follows:

where J and K are the orders of the two FIR filters (0 < , and are the coefficients of the impulse responses of these two FIR filters.

There are several widely used machine learning methods for the estimation of the coefficients of the impulse response of an FIR filter, e.g., ordinary least squares (OLS), regularized least squares (i.e., ridge regression), and partial least squares (PLS) [15]. In this paper, we choose the ridge regression as our estimation method that solves the following optimization problem:

where and are the regularization parameters.

E. Tracking the Number of Infected Persons ˆX(t) and the Number of Recovered Persons ˆR(t) of the Time-dependent SIR Model

In this subsection, we show how we use the two FIR filters to track and predict the number of infected persons and the number of recovered persons in the time-dependent SIR model. Given a period of historical data , we first measure by (11) and (12). Then we solve the ridge regression (with the objective functions in (15) and (16) and the constraints in (13) and (14)) to learn the coefficients of the FIR filters, i.e., and . Once we learn these coefficients, we can predict and at time by the trained ridge regression in (13) and (14).

Denote by (resp. ) the predicted number of infected (resp. recovered) persons at time t. To predict and at time t = T, we simply replace and by and in (9) and (10). This leads to

To predict and for t > T, we estimate and by using (13) and (14). Similar to those in (17) and (18), we predict and as follows:

The detailed steps of our tracking/predicting method are outlined in Algorithm 1.

We note that this deterministic epidemic model is based on the mean-field approximation for X(t) and R(t). Such an approximation is a result of the law of large numbers. Therefore, when X(t) and R(t) are relatively small, the mean-field approximation may not be as accurate as expected. In those cases, one might have to resort to stochastic epidemic models, such as Markov chains.

III. THE SIR MODEL WITH UNDETECTABLE INFECTED PERSONS

According to the recent report from WHO [2], only 87.9% of COVID-19 patients have a fever, and 67.7% of them have a dry cough. This means there exist asymptomatic infections. Recent studies in [7] and [16] also pointed out the existence of the asymptomatic carriers of COVID-19. Those people are unaware of their contagious ability, and thus get more people infected. The transmission rate can increase dramatically in this circumstance.

To take the undetectable infected persons into account, we propose the SIR model with undetectable infected persons in this section. We assume that there are two types of infected persons. The individuals who are detectable (with obvious symptoms) are categorized as type I infected persons, and the asymptomatic individuals who are undetectable are categorized as type II infected persons. For an infected individual, it has probability to be type I and probability to be type II, where . Besides, those two types of infected persons have different transmission rates and recovering rates, depending on whether they are under treatment or isolation or not. We denote and as the transmission rate and the recovering rate of type I at time t. Similarly, and are the transmission rate and the recovering rate for type II at time t.

A. The Governing Equations for the SIR Model with Undetectable Infected Persons

Now we derive the governing equations for the SIR model with two types of infected persons. Let (resp. ) be the number of type I (resp. type II) infected persons at time t. Similar to the derivation of (8), (9) in Subsection II-C, we assume that in the initial stage of the epidemic and split X(t) into two types of infected persons. We have the following difference equations:

where , and are constants. It is noteworthy that those constants can also be time-dependent as we have in Section II. However, in this section, we set them as constants to focus on the effect of undetectable infected persons. Rewriting (21) and (22) in the matrix form yields the following matrix equation:

where . Let A be the transition matrix of the above system equations, i.e.,

It is well-known (from linear algebra) such a system is stable if the spectral radius (the largest absolute value of the eigenvalue) of A is less than 1. In other words, and will converge gradually to finite constants when t goes to infinity. In that case, there will not be an outbreak. On the contrary, if the spectral radius is greater than 1, there will be an outbreak, and the number of infected persons will grow exponentially with respect to time t (at the rate of the spectral radius).

B. The Basic Reproduction Number

To further examine the stability condition of such a system, we let

Note that is simply the basic reproduction number of a newly infected person as an infected person can further infect on average (resp. ) persons if it is of type I (resp. type II) and that happens with probability (resp. ). In the following theorem, we show that there is no outbreak if and there is an outbreak if . Thus, in (25) is known as the percolation threshold for an outbreak in such a model [11].

Theorem 1. If , then the spectral radius of A in (24) is less than 1 and there is no outbreak of the epidemic. On the other hand, if , then the spectral radius of A in (24) is larger than 1 and there is an outbreak of the epidemic.

Proof. (Theorem 1)

First, we note that and are recovering rates and they cannot be larger than 1 in the discrete-time setting, i.e., it takes at least one day for an infected person to recover. Thus, the matrix A is a positive matrix (with all its elements being positive). It then follows from the Perron-Frobenius theorem that the spectral radius of the matrix is the larger eigenvalue of the matrix.

Now we find the larger eigenvalue of the matrix A. Let I be the identify matrix and

Let

and

It is straightforward to show that the two eigenvalues of are

and

Note that . In view of (26), the larger eigenvalue of the transition matrix A is .

If , we know that , and . Thus, we have from (28) that . In view of (30), we conclude that

This shows that and the spectral radius of A is less than 1.

On the other hand, if , then and we have from (30) that

This shows that and the spectral radius of A is larger than 1.

The relation between the system parameters and the phase transition will be shown in Subsection V-G.

C. Herd Immunity

Herd immunity is one way to resist the spread of a contagious disease if a sufficiently high fraction of individuals are immune to the disease, especially through vaccination. One interesting strategy, once suggested by Boris Johnson, the British prime minister, is to have a sufficiently high fraction of individuals infected by COVID-19 and recovered from the disease to achieve herd immunity. The question is, what will be the fraction of individuals that need to be infected to achieve herd immunity for COVID-19.

To address such a question, we note that herd immunity corresponds to the reduction of the number of susceptible persons in the SIR model. In our previous analysis, we all assume that every person is susceptible to COVID-19 at the early stage and thus . For the analysis of herd immunity, we assume that there is a probability h that a randomly chosen person is susceptible at time t. This is equivalent to that fraction of individuals are immune to the disease. Under such an assumption, we then have

In view of the difference equation for X(t) in (8), we can rewrite (21)-(23) to derive the governing equations for herd immunity as follows:

In comparison with the original governing equations in (21)-(23), the only difference is the change of the transmission rate of type I (resp. type II) from to (resp. from to ). Thus, herd immunity effectively reduces the transmission rates by a factor of h. As a direct consequence of Theorem 1, we have the following corollary.

Corollary 2. For a contagious disease modeled by our SIR model with two types of infected persons that has in (25) greater than 1, herd immunity can be achieved after at least fraction of individuals being infected and recovered from the contagious disease, where

IV. THE INDEPENDENT CASCADE (IC) MODEL FOR DISEASE PROPAGATION IN NETWORKS

Our analysis in the previous section does not consider how the structure of a social network affects the propagation of a disease. There are other widely used policies, such as social distancing, that could not be modeled by our SIR model with undetectable infected persons in Section III. To take the network structure into account, in this section, we consider the independent cascade (IC) model for disease propagation. The IC model was previously studied by Kempe, Kleinberg, and Tardos in [13] for the influence maximization problem in viral marketing. In the IC model, there is a social network modeled by a graph G = (V, E), where V is the set of nodes, and E is the set of edges. An infected node can transmit the disease to a neighboring susceptible node (through an edge) with a certain propagation probability. As there are two types of infected persons in our model, we denote by (resp. ) the propagation probability that a type I (resp. type II) infected node transmits the disease to an (immediate) neighbor of the infected node. Once a neighboring node is infected, it becomes a type I (type II) infected node with probability (resp. and it can continue the propagation of the disease to its neighbors. Continuing the propagation, we thus form a subgraph of G that contains the set of infected nodes in the long run. Call such a subgraph the infected subgraph. One interesting question is how one controls the spread of the disease so that the total number of nodes in the infected subgraph remains small even when the total number of nodes is very large.

A. The Infected Tree in the Configuration Model

The exact network structure, i.e., the adjacency matrix of the network G, is in general very difficult to obtain for a large population. However, it might be possible to learn some characteristics of the network, in particular, the degree distribution of the nodes. The configuration model (see, e.g., the book [11]) is one family of random networks that are specified by degree distributions of nodes. A randomly selected node in such a random network has degree k with probability . The edges of a node are randomly connected to the edges of the other nodes. As the edge connections are random, the infected subgraph appears to be a tree (with high probability) if one follows an edge of an infected node to propagate the disease to the other nodes in such a network. The tree assumption is one of the most important properties of the configuration model. Another crucial property of the configuration model is the excess degree distribution. The probability that one finds a node with degree k + 1 along an edge connected to that node is

Thus, excluding the edge coming to the node, there are still k edges that can propagate the disease. This is also the reason why is called the excess degree distribution. Note that the excess distribution is in general different from the degree distribution . They are the same when is the Poisson degree distribution. In that case, the configuration model reduces to the famous Erd¨os-R´enyi random graph.

As the infected subgraph is a tree in the configuration model, we are interested to know whether the size of the infected tree is finite. We say that there is no outbreak if the size of the infected tree of an infected node is finite with probability 1. Let (resp. ) be the probability that the size of the infected tree of a type I (resp. type II) node is finite via a specific one of its neighbors. Then,

To see the intuition of (37), we note that either the neighbor is infected or not infected. It is not infected with probability . On the other hand, it is infected with probability . Then with probability (resp. ), the infected neighbor is of type I (reps. type II). Also, with probability , the neighbor has additional k edges to transmit the disease. From the tree assumption, the probability that these k edges all have finite infected trees is (resp. ) if the infected neighbor is of type I (reps. type II). The equation in (38) follows from a similar argument.

Let

be the moment generating function of the excess degree distribution. Then we can simplify (37) and (38) as follows:

From (40) and (41), we can solve and by starting from and

It is easy to show (by induction) that and . Thus, they converge to some fixed point solution and of (40) and (41).

B. Connections to the Previous SIR Model

Now we show the connections to the SIR model in Section III by specifying the propagation probabilities and .

Suppose that one end of a randomly selected edge is a type I node. Then this type I node will infect persons on average from the SIR model in Section III. Since the average excess degree is , the average number of neighbors infected by this type I node is . In order for a type I node to infect the same average number of nodes in the SIR model in Section III, we have

Similarly for , we have

With the propagation probabilities and specified in (44) and (45), we have the following stability result.

Theorem 3. For the IC model (for disease propagation) in a random network constructed by the configuration model, suppose that the propagation probabilities and are specified in (44) and (45). Then the size of the infected tree is finite with probability 1 if

Under such a condition, there is no outbreak.

Proof. (Theorem 3)

Let and . It suffices to show that u = e is the unique solution for the system of equations in (40) and (41) if . We prove this by contradiction. Suppose that there is a solution of (40) and (41) that either or when .

Since the moment generating function in (39) is a convex function, we have from the first order Taylor’s expansion for and that

Note that . Replacing (46) and (47) into (40) and (41), we have

This can be further simplified by using (44) and (45). Thus, we have

Let

and

We now rewrite (50) in the following matrix form:

It is straightforward to see that the two eigenvalues of B are

Moreover, the eigenvector corresponding to the eigenvalue is

Recursively expanding (53) for m times yields

Since , both and converge to the zero vectors as . Letting in (54) yields . This contradicts to the assumption that either or .

C. Social Distancing

Social distancing is an effective way to slow down the spread of a contagious disease. One common approach of social distancing is to allow every person to keep its interpersonal contacts up to (on average) a fraction a of its normal contacts (see, e.g., [17], [18]). In our IC model, this corresponds to that every node randomly disconnects one of its edges with probability .

As in the previous subsection, we let (resp. ) be the probability that the size of the infected tree of a type I (resp. type II) node is finite via a specific one of its neighbors. Then,

To see (55), note that a neighboring node of an infected node can be infected only if (i) the edge connecting these two nodes is not removed (with probability ), and (ii) the disease propagates through the edge (with the propagation probability ). This happens with probability . Then with probability (resp. ), the infected neighbor is of type I (reps. type II). Also, with probability , the neighbor has additional k edges to transmit the disease. From the tree assumption, the probability that these k edges all have finite infected trees is (resp. ) if the infected neighbor is of type I (reps. type II). The equation in (56) follows from a similar argument.

In comparison with the two equations in (37) and (38), we conclude that this approach of social distancing reduces the propagation probabilities and to and , respectively. As a direct consequence of Theorem 3, we have the following corollary.

Corollary 4. Suppose that a social distancing approach allows every person to keep its interpersonal contacts up to (on average) a fraction a of its normal contacts. For the IC model (for disease propagation) in a random network constructed by the configuration model, the size of the infected tree is finite with probability 1 if

Under such a condition, there is no outbreak.

Another commonly used approach of social distancing is canceling mass gatherings. Such an approach aims to eliminate the effect of “superspreaders” who have lots of interpersonal contacts. For this, we consider a disease control parameter and remove nodes with the number of edges larger than or equal to in our IC model. Analogous to the derivation of (37) and (38), we have

To see (58), we note that a type I infected person only infects a finite number of persons along an edge if (i) the disease does not propagate through the edge (with probability ), (ii) the disease propagates through the edge and the neighboring node is removed (with probability ), and (iii) the disease propagates through the edge and the neighboring node only infects a finite number of persons (with probability ). The argument for (59) is similar.

Analogous to the stability result of Theorem 3, we have the following stability result for a social distancing approach that cancels mass gatherings.

Theorem 5. Consider a social distancing approach that cancels mass gatherings by removing nodes with the number of edges larger than or equal to . For the IC model (for disease propagation) in a random network constructed by the configuration model, suppose that the propagation probabilities and are specified in (44) and (45). Then the size of the infected tree is finite with probability 1 if

Under such a condition, there is no outbreak.

Proof. (Theorem 5)

As in the proof of Theorem 3, we let and . It suffices to show that u = e is the unique solution for the system of equations in (58) and (59) if the inequality in (60) is satisfied. We prove this by contradiction. Suppose that there is a solution of (58) and (59) that either or when the inequality in (60) is satisfied.

Since is a convex function for and , we have . It then follows from (58) and (59) that

Writing these two equations in the matrix form and using (44) and (45) yields

where B is the matrix in (51). Note that is simply the expected excess degree. Following the same argument as that in Theorem 3, one can easily show that and when the inequality in (60) is satisfied. This contradicts to the assumption that either or .

Unfortunately, it is difficult to obtain an explicit expression for to prevent an outbreak in (60). For this, we will resort to numerical computations in the next section.

D. The Ratio of the Population Infected in the Long Run

Andrew Cuomo, the governor of New York State, said on Apr. 23, 2020, that 13.9% of a group of 3, 000 people tested positive for COVID-19 antibodies. It is even higher in the City of New York. Such a test implies that a certain proportion of the population in the State of New York were infected (and recovered), and that is a piece of solid evidence for an outbreak in the State of New York. In all our previous studies, we have been focusing on how to prevent an outbreak. If the disease cannot be contained, we ask the question of what is the ratio of the population infected in the long run.

For the SIR model, the ratio of the population infected in the long run, denoted by r, is defined as

If the SIR model has the stationary transmission rate and the stationary recovering rate , then the basic reproduction number is simply . It is well-known (see, e.g., [11]) that r in the stationary SIR model satisfies the fixed point equation:

As discussed in the previous section, the SIR model does not take the network structure into account. To see the effect of the degree distribution to the ratio of the population infected in the long run, we consider the IC model for disease propagation in a random network generated by the configuration model with the degree distribution . In Theorem 3, we already showed that there is no outbreak if . In this subsection, we will show that if , then a certain proportion of the population will be infected. This is stated in the following theorem.

Theorem 6. For the IC model (for disease propagation) in a random network constructed by the configuration model, suppose that the propagation probabilities and are specified in (44) and (45). Let (resp. ) be the moment generating function of the degree distribution (resp. excess degree distribution). If

then there is a nonzero probability r that a randomly selected node is infected in the long run (after the propagation of the disease in the IC model). The probability r can be computed by the following equation:

where

and v is the probability that one end node of a randomly selected edge is infected by one of its neighbors that is not on the other end of the selected edge. Moreover, the probability v is the unique solution in (0, 1] of the following fixed point equation:

Proof. (Theorem 6) We first show (67). Using (66), we can rewrite (67) as follows:

To explain (68), let us consider a node, say node x, that is on one end of a randomly selected edge (see Figure 1 for an illustration). Excluding the neighbor on the other end of the randomly selected edge, we call the remaining neighbors of x its excess neighbors (like the excess degree distribution). The probability on the left-hand side of (68) is simply the probability that node x is not infected by one of its excess neighbors. Suppose that there are k excess neighbors of x. An excess neighbor of x, say node y, is on the other end of an edge and it is not infected with probability v. If node y is not infected, then y cannot infect x. On the other hand, if y is infected, then with probability (resp. is of type I (resp. II) and it does not infect x with the probability (resp. ). Thus, the probability that node x is not infected from node y is . As there are k (excess) neighbors, node x is not infected is

(in the configuration model). Since the probability that there are k excess neighbors of node x is , averaging over the excess degree distribution yields (68).

Similarly, we have

This then leads to (65).

With and being specified in (44) and (45) and in (1), we have

Figure 1: An illustration of the derivation of (68). Node x is on one end of a randomly selected edge, and node y is an (excess) neighbor of x. Here the number of excess neighbors k is 4.

To find r, we need to solve v from the fixed point equation in (67). For this, we let and rewrite (67) as follows:

Since and is a convex function,

Analogous to the argument in the proofs of Theorem 3 and Theorem 5, it is easy to see from (70) and (71) that (and thus v = 0) if . Moreover, there is a unique (and thus v > 0) for the fixed point equation in (71) if . Such a solution can be solved iteratively by using the initial value 0 for as described in (42).

In particular, for the ER model, the degree distribution is the Poisson degree distribution with mean c, i.e.,

The excess degree distribution of the Poisson degree distribution is the same as the degree distribution, i.e., . In this case, we have r = v and r is the unique solution in (0, 1] of the fixed point equation

when . This is exactly the same as (64) for the SIR model with the stationary transmission rate and the stationary recovering rate . The fixed point equation in (72) has a very intuitive explanation. The left-hand side is the probability that a randomly selected node is not infected, while the right-hand side is the probability that an infected node does not infect any of its neighbors as the degree distribution is Poisson.

V. NUMERICAL RESULTS

A. Dataset

In this section, we analyze and predict the trend of COVID-19 by using our time-dependent SIR model in Section II and the SIR model with undetectable infected persons in Section III. For our analysis and prediction of COVID-19, we collect our dataset from the National Health Commission of the People’s Republic of China (NHC) daily Outbreak Notification [1]. NHC announces the data as of 24:00 the day before. We collect the number of confirmed cases, the number of recovered persons, and the number of deaths from Jan. 15, 2020 to Mar. 2, 2020 as our dataset. The confirmed case is defined as the individual with positive real-time reverse transcription polymerase chain reaction (rRT-PCR) result. It is worth noting that in the Hubei province, the definition of the confirmed case has been relaxed to the clinical features since Feb. 12, 2020, while the other provinces use the same definition as before.

B. Parameter Setup

For our time-dependent SIR model, we set the orders of the FIR filters for predicting and as 3, i.e., J = K = 3. The stopping criteria of the model is set to . Since the numbers of infected persons before Jan. 27, 2020 are too small to exhibit a clear trend (which may contain noises), we only use the data after Jan. 27, 2020 as our training data for predicting and .

We use the scikit-learn library [19] (a third-party library of Python 3) to compute the ridge regression. The regularization parameters of predicting and are set to 0.03 and respectively. Since the transmission rate is nonnegative, we set it to 0 if it is less than 0. Then, we use Algorithm 1 to predict the trend of COVID-19.

C. Time Evolution of the Time-dependent SIR Model

In Figure 2, we show the time evolution of the number of infected persons and the number of recovered persons. The circle-marked solid curves are the real historical data by Mar. 2, 2020, and the star-marked dashed curves are our prediction results for the future. The prediction results imply that the disease will end in 6 weeks, and the number of the total confirmed cases would be roughly 80, 000 if the Chinese government remains their control policy, such as city-wide lockdown and suspension of works and classes.

In Figure 3, we show the measured and from the real historical data. We can see that decreases dramatically, and increases slightly. This is a direct result of the Chinese government that tries to suppress the transmission rate by city-wide lockdown and traffic halt. On the other hand, due to the lack of effective drugs and vaccines for COVID-19, the recovering rate grows relatively slowly. Additionally, there is a definition change of the confirmed case on Feb. 12, 2020 that makes the data related to Feb. 11, 2020 have no reference value. We mark these data points for and with the gray dashed curve.

In an epidemic model, one crucial question is whether the disease can be contained and the epidemic will end, or whether there will be a pandemic that infects a certain fraction of the total population n. To answer this, one commonly used metric is the basic reproduction number that is defined as the average number of additional infections by an infected person before it recovers. In the classical SIR model, is simply as an infected person takes (on average) days to recover, and during that period time, it will be in contact with

Figure 2: Time evolution of the time-dependent SIR model of the COVID-19. The circle-marked solid curve with dark orange (resp. green) color is the real number of infected persons X(t) (resp. recovered persons R(t)), the star-marked dashed curve with light orange (resp. green) color is the predicted number of infected persons (resp. recovered persons).

Figure 3: Measured transmission rate and recovering rate of the COVID-19 from Jan. 15, 2020 to Feb. 19, 2020. The two curves are measured according to (12) and (11) respectively.

(on average) persons. In our time-dependent SIR model, the basic reproduction number is a function of time, and it is defined as . If , the disease will spread exponentially and infects a certain fraction of the total population n. On the contrary, the disease will eventually be contained. Therefore, by observing the change of with respect to time or even predicting in the future, we can check whether certain epidemic control policies are effective or not.

In Figure 4, we show the measured basic reproduction number , and the predicted basic reproduction number . The blue circle-marked solid curve is the measured and the purple star-marked dashed curve is the predicted (from Feb. 15, 2020). It is clear that has decreased dramatically since Jan. 28, 2020, and it implies that the control policies work in China. More importantly, it shows that the turning point is Feb. 17, 2020 when . In the following days after Feb. 17, 2020, X(t) will decrease exponentially, and that will lead to the end of the epidemic in China. Our model predicts precisely that will go less than 1 on Feb. 17, 2020 by 3 days in advance (Feb. 14, 2020). The results show that our model is very effective in tracking the characteristics of and .

Figure 4: Basic reproduction number of the time-dependent SIR model of the COVID-19 in China. The circle-marked solid curve with blue color is the based on the given data from Jan. 27, 2020 to Feb. 20, 2020, the star-marked dashed curve with purple color is the predicted based on the data from Jan. 27, 2020 to Feb. 15, 2020, and the dashed line with red color is the percolation threshold 1 for the basic reproduction number.

D. One-day Prediction

To show the precision of our model, we demonstrate the prediction results for the next day (one-day prediction) in Figure 5. It contains the predicted number of infected persons (orange star-marked dashed curve), the predicted number of recovered persons (green star-marked dashed curve), and the real number of infected and recovered persons (dark orange and dark green circle-marked solid curves) every day. The unpredictable days due to the change of the definition of the confirmed case on Feb. 12, 2020. are marked as gray. The predicted curves are extremely close to the measured curves (obtained from the real historical data). In this figure, we also annotate the predicted number of infected persons and the predicted number of recovered persons on Mar. 3, 2020.

We further examine our prediction accuracy in Figure 6. The error rates are all within except for the predicted number of recovered persons on Feb. 1, Feb. 3, and Feb. 5, 2020. The gray dashed curve stands for the unpredictable points due to the change of definition of the confirmed case. However, from the prediction results after Feb. 16, 2020, we find that our model can still keep tracking and accurately and overcome the impact of the change of the definition.

Figure 5: One-day prediction for the number of infected and recovered persons. The unpredictable points due to the change of definition of the confirmed case are marked as gray. The circle-marked solid curve with dark orange (resp. green) color is the real number of infected persons X(t) (resp. recovered persons R(t)), the star-marked dashed curve with light orange (resp. green) color is the predicted number of infected persons (resp. recovered persons ).

Figure 6: Errors of the one-day prediction of the number of infected and recovered persons. The unpredictable points due to the change of definition of the confirmed case on Feb. 12, 2020 are marked as the gray dash curve.

E. Connections to the Wuhan City Lockdown

In this subsection, we show the connections between the epidemic prevention policies issued by the Chinese government and the historical data of the time-varying transmission rates.

As shown in Figure 3, the disease has been gradually controlled in China as time goes on. Excluding the small number of cases (X(t) < 500) before Jan. 21, 2020 (that causes the curve to fluctuate a lot), it is notable that increases gradually then drops dramatically during Jan. 23, 2020 to Jan. 28, 2020, and it reaches the peak point on Jan. 26, 2020, which coincides with the trends of the moving out in Wuhan during the Chunyun (Spring Festival travel season) [20] in Figure 7. Especially, the emigration trend is almost the same as the during Jan. 20, 2020 to Jan. 25, 2020. We speculate that people rushed into the public transportation system when the announcement of the Wuhan city lockdown was out, which significantly increases the contact among people and speeds up the spread of the virus. As a result of that, the transmission rate increases substantially. As pointed out in Zhong’s study [21], the median of the incubation period of COVID-19 is 3 days among 1099 valid confirmed cases, which makes the emigration trend aligns with if the extra 3 days are taken into account. Finally, the disease is gradually getting under control after the lockdown. The basic reproduction number is less than 1, i.e., since Feb. 17, 2020.

Figure 7: Trends of moving out during the Chunyun (Spring Festival travel season) in Wuhan city. The vertical axis represents the ratio between the number of people leaving the city and the resident population in Wuhan. The orange curve shows the ratio in 2019, while the blue curve shows the ratio in 2020. We redraw this figure by our-self from https://qianxi.baidu.com/ [20].

F. Basic Reproduction Numbers of Several Other Countries

In addition to the dataset for China, we also measure the basic reproduction number on Mar. 31, 2020 for several countries from the datasets in [12]. This is shown in the last column of Table II. As the data for the cumulative numbers of recovered persons for these countries are noisy, we also show the estimated under various assumptions of the average time to recover . The values for the five countries, including United States of America, the United Kingdom, France, Iran, and Spain are very high. On the other hand, it seems that Italy is gaining control of the spread of the disease after the Italian government announces the lockdown and forbids the gatherings of people on Mar. 10, 2020. Also, both Germany and Republic of Korea are capable of controlling the spread of the disease.

G. The Effects of Type II Infected Persons

In this subsection, we show how undetectable (type II) infected persons affect the epidemic. In particular, we are in-

Table II: The estimated under various assumptions of the average time to recover () from COVID-19, and the measured on Mar. 31, 2020.

terested in addressing the question of whether the existence of undetectable infected persons (type II) can cause an outbreak.

To carry out our numerical study, we need to fix some variables in the system of difference equations in (21)-(23). For the transmission rate (resp. recovering rate) of type I infected persons, i.e., (resp. ), we set it to be the measured (resp. ) on Mar. 1, 2020 in China. The rationale behind this is that type I infected persons were detected and they were under treatment and isolation after Mar. 1, 2020 in China. Also, as there is no medicine for COVID-19, we may assume that these two types of infected persons have the same recovering rate, i.e, . In view of the system equation in (21)-(23), there are still two free variables and , where is the probability that an infected person is of type II and is the transmission rate of type II.

In Figure 8, we illustrate how and affect the outbreak of the COVID-19. Such a figure is known as the phase transition diagram in [11]. The black curve in Figure 8 is the curve when the spectral radius of the transition matrix A in (24) equals to 1. This curve represents the percolation threshold of COVID-19. If and fall above the black curve (in the orange zone), then there will be an outbreak. On the contrary, if and fall below the black curve (in the yellow zone), then there will not be an outbreak. As shown in Figure 8, we would like to point out the importance of detecting an infected person. As long as more than 90% of those infected persons can be actually detected and properly isolated and treated, it is possible to contain the spread of the disease even if the transmission rate of type II infected persons, i.e., , is as high as 0.7. On the other hand, suppressing the transmission rate of type II infected persons can also be effective in controlling the disease while the detection rate is not that high. For example, wearing masks and washing hands can be an effective epidemic prevention mechanism to reduce .

In the following experiments, we extend our study to other countries, including Japan, Singapore, South Korea, Italy, and Iran. We collect the historical data from Jan. 22, 2020 to Mar. 2, 2020 from the GitHub of Johns Hopkins University [12]. For these countries, the transmission rates and the recovering rates (measured from the time-dependent SIR

Figure 8: Phase transition diagram of an outbreak with respect to and . The black curve is the percolation threshold. The orange area means the disease will be an outbreak while the yellow area means the disease is under control.

model in Section II) during the initial period with a rapid increase of the number of confirmed cases can be viewed as and . This is because during that period of time, there is no epidemic prevention intervention, and all the infected persons are basically not detected. It is interesting to note that different countries might have different and . On the other hand, only 87.9% of COVID-19 cases have a fever from the report of WHO [2]. If we use body temperature as a means to detect an infected person, then only 87.9% of COVID-19 cases can be detected. For this, we set .

With and specified in the previous study for China, we plot the phase transition diagram in Figure 9 in terms of the two variables and . Again, the black curve is the curve when the spectral radius of the transition matrix A in (24) equals to 1. Such a curve represents the percolation threshold of a COVID-19 outbreak. If and fall above the black curve (in the orange zone), then there will be an outbreak. On the contrary, if and fall below the black curve (in the yellow zone), then there will not be an outbreak. The countries with large confirmed cases, including Japan, Singapore, South Korea, Italy, and Iran, are marked in Figure 9. From Figure 9, we observe that both Singapore and Japan are below the percolation threshold. But Japan is much closer to the percolation threshold. On the other hand, both South Korea and Italy are above the percolation threshold, and they are on the verge of a potential outbreak on Mar. 2, 2020. These two countries must implement epidemic prevention policies urgently. Not surprisingly, on Mar. 10, 2020, the Italian government announces the lockdown and forbids the gatherings of people. It is worth mentioning that there are two marks for Iran in the Figure. The one above the percolation threshold is the one without adding the number of deaths into the number of recovered persons. The other one below the percolation threshold is the one that adds the number of deaths into the number of recovered persons. For some unknown reason, the death rate in Iran is higher than the other countries. The high death rate seems to prevent an outbreak in Iran.

Figure 9: Phase transition diagram of an outbreak with respect to and . The black curve is the percolation threshold. The orange area means the disease will be an outbreak, while the yellow area means the disease is under control.

H. The Effects of Social Distancing

In this subsection, we show the numerical results of the social distancing approach that cancels mass gatherings by removing nodes with the number of edges larger than or equal to . As shown in Theorem 5, the basic reproduction number is reduced by a factor of

the excess degree distribution of . For this experiment, we use the dataset collected by [22] from Facebook. This dataset represents the verified Facebook page (with blue checkmark) networks of the artist category. The blue checkmark means Facebook has confirmed that an account is the authentic presence of the public figure, celebrity, or global brand it represents. Each node in the network represents the page, and edges between two nodes are mutual likes among them. This dataset is composed of 50, 515 nodes and 819, 306 edges. Some other properties are listed as follow: mean degree 32.4, max degree 1, 469, diameter 11, and clustering coefficient 0.053. In Figure 10, we show the log-log plots of the degree distribution and the excess degree distribution of this dataset. The degree distribution appears to be a (truncated) Pareto distribution with the exponent 1.69 (the slope in the figure).

Figure 10: The degree distribution and the excess degree distribution of the Facebook dataset.

In Figure 11, we plot the reduction ratio

a function of . The ratio is between 0 and 1, and it is monotonically increasing in . Using the values (on Mar. 31, 2020) in the last column of Table II, we also show that the minimum s to prevent an outbreak in Italy, U.S., and South Korea are 63, 195, and 435, respectively; moreover, the affected fraction of tail distributions are 13.1%, 2.2%, and 0.4%, respectively. In particular, if canceling mass gathering is the only measure used for controlling COVID-19 in the U.S. with the value of 12.59 on Mar. 31, 2020, then one can prevent an outbreak by “removing” all the nodes with degrees larger than or equal to 63 (in the Facebook dataset), and the removal might affect 13.1% of the nodes in the Facebook dataset.

� ��� ��� ���� ���� ���� ������ ��� ��� k0

Figure 11: The reduction ratio

The minimum s to prevent an outbreak in Italy, U.S., and South Korea are 63, 195, and 435, respectively.

I. The Ratio of the Population Infected in the Long Run In this subsection, we show the numerical results of the ratio of the population infected in the long run in Subsection

IV-D. As shown in Theorem 6, if , then there is a nonzero probability r that a randomly selected node is infected in the long run. In Figure 12, we plot the infected probability r as a function of for various degree distributions. For this numerical experiment, we use the same Facebook dataset in the previous subsection. Note that the mean degree and the mean excess degree of this Facebook dataset are 32.4 and 155.6, respectively. We also generate three random networks: two from the Erd¨os-R´enyi (ER) model [23] with the mean degrees c = 32 and c = 155, and one from the Barab´asiAlbert (BA) model [24] with the mean degree c = 32. The numbers of nodes are all set to 50, 000. The degree distribution of the BA model (generated by using the linear preferential attachment rule) is known to follow the asymptotic power-law distribution with for large k’s. The mean excess degree of the BA model is 84.05.

From Figure 12, one can see that there is an outbreak when . The probability r is increasing in . Also, as shown in this figure, the two curves of the two ER models (even with different mean degrees) overlap with each other. That means that r is independent of the mean degree when the degree distribution is Poisson (as shown in (72) at the end of Subsection IV-D).

On interesting observation is that the infected probability r of the ER model is larger than the infected probabilities of the other networks with the power-law degree distributions, i.e., the Facebook dataset and the BA model. This is because we relate the propagation probability by in (70) (to ensure that the average number of additional infections by an infected person in the IC model is the same as that in the SIR model). Thus, a network with a large mean excess degree has a low propagation probability (when is fixed). That leads to a smaller ratio of the population infected in the long run. To explain this further, we note that the probability of having a super spreader (a node with a very large degree) in a network with a power-law degree distribution is higher than that of the ER model. A super spreader is capable of infecting a large number of its first neighbors. However, the neighbors of a super spreader tend to have relatively low degrees (contacts) than that of a super spreader. As such, it makes the disease more difficult to spread further in the IC model. On the other hand, the probability of having a super spreader in the ER model is exponentially small. Since the infected subgraph appears to be a tree (if one follows an edge of an infected node to propagate the disease to the other nodes in the IC model), the disease may spread slowly in the early stage of the epidemic in the ER model than that in the BA model. However, if , the disease may continue to propagate in the ER model to a larger fraction of the total population in the long run. For example, when , the percentage of infected persons is more than 90% in the ER model.

VI. DISCUSSIONS AND SUGGESTIONS

As of Mar. 2, 2020, it seems that COVID-19 has been gradually controlled in China since the prevention policies, such as city-wide lockdown was issued in China in Jan. 2020. Although the policies such as traffic halt, small community

Figure 12: The probability r that a randomly selected node is infected in the long run as a function of for various degree distributions.

management, and city-wide lockdown can effectively reduce the transmission rates and ; however, these relatively extreme policies not only restrict the right of personal freedom but also affect the normal operation of society. These extreme policies forced several companies and factories to halt production, which impacts all sectors of the economy. Therefore, to strike a balance between the prevention of disease and ensuring the normal operation of society is crucial, and it is important to suggest the so-called “optimal” control policies.

For this, we would like to put forward some discussions and suggestions for controlling the spread of COVID-19 based on the observation made from the results of our system models. From our results in Theorem 1, Corollary 2, Theorem 3, Corollary