key: cord-0902226-h0o8kwla
authors: Feng, Shanshan; Zhang, Juping; Li, Juan; Luo, Xiao-Feng; Zhu, Huaiping; Li, Michael Y.; Jin, Zhen
title: The Impact of Quarantine and Medical Resources on the Control of COVID-19 in Wuhan based on a Household Model
date: 2022-02-26
journal: Bull Math Biol
DOI: 10.1007/s11538-021-00989-y
sha: 6a0b2d60be8af4abc68225decd410deba7bcccea
doc_id: 902226
cord_uid: h0o8kwla

In order to understand how Wuhan curbed the COVID-19 outbreak in 2020, we build a network transmission model of 123 dimensions incorporating the impact of quarantine and medical resources as well as household transmission. Using our new model, the final infection size of Wuhan is predicted to be 50,662 (95%CI: 46,234, 55,493), and the epidemic would last until April 25 (95%CI: April 23, April 29), which are consistent with the actual situation. It is shown that quarantining close contacts greatly reduces the final size and shorten the epidemic duration. The opening of Fangcang shelter hospitals reduces the final size by about 17,000. Had the number of hospital beds been sufficient when the lockdown started, the number of deaths would have been reduced by at least 54.26%. We also investigate the distribution of infectious individuals in unquarantined households of different sizes. The high-risk households are those with size from two to four before the peak time, while the households with only one member have the highest risk after the peak time. Our findings provide a reference for the prevention, mitigation and control of COVID-19 in other cities of the world.

time for the prevention and control of the epidemic in China and other countries around the world (Chinazzi et al. 2020; Wells et al. 2020 ). However, for Wuhan, because the COVID-19 was an emerging infectious disease, human had little knowledge of COVID-19 and medical resources were not well ready yet. The lack of outpatient medical resources, such as test kits and medical personnel, makes it impossible for all infectious individuals to be timely tested and confirmed, which leads to a longer diagnosis time. The limited number of hospital beds makes it impossible to timely admit and cure all symptomatic infectious individuals, which leads to the prolonged transmission time of the infectious individuals. The lack of centrally quarantined points prevents close contacts from being quarantined in a timely and effective manner, which also leads to the increase in the risk of transmission. In order to overcome the shortage of medical personnel and materials in Wuhan, all other provinces successively sent essential healthcare personnel and medical resources to support Wuhan City in batches.

As an important medical resource, the number of hospital beds is directly related to the isolation and treatment of symptomatic infectious individuals. The more sufficient the number of hospital beds is, the more infectious individuals will be effectively isolated, and the transmission between people will be correspondingly reduced. Therefore, it is of great significance to study the emerging infectious diseases incorporated by the number of hospital beds. For the shortage of hospital beds, Wuhan had been increasing the number of designated hospitals (DHs) and began to build the Huoshenshan Hospital and Leishenshan Hospital on January 23, 2020, and January 25, 2020, respectively, immediately after the lockdown. Moreover, Fangcang Shelter Hospitals (FSHs) have been rapidly established and put into use since February 5, 2020, to treat patients with mild illnesses (Chen et al. 2020) . 6 These measures have provided the most important material basis for the epidemic control in Wuhan. On the other hand, in the absence of an effective vaccine, quarantining close contacts is the most effective method in the prevention and control of emerging infectious diseases. To this end, while building hospitals, Wuhan had been increasing the number of centrally quarantined points. On the premise that there are vacant beds, every confirmed infectious individual is immediately hospitalized for massive isolation. At the same time, their close contacts are checked according to their social networks and are quarantined at home, centrally, or not quarantined. In addition, due to the fact that the COVID-19 occurred during the Spring Festival in China and the strict prevention and control measures implemented, the epidemic showed significant family clustered (see Footnote 5).

Up to now, a large number of mathematical models have been contributed to study the transmission of COVID-19 and the evaluation of prevention and control measures. These models can be divided into two categories: well-mixed compartment models and heterogeneous mixing network models. Among well-mixed compartment models, the simplest model is to divide individuals into the susceptible, the infectious and the recovered (Roda et al. 2020) . Taking into account the quarantine of close contacts, Tang et al. added compartments for quarantined susceptible and latent individuals to the well-mixed model (Tang et al. 2020) . Li et al. used a simulation model to investigate the role of two types of hospital beds in controlling the epidemic in Wuhan and concluded that FSHs played a key role in bringing the epidemic to an end (Li et al. 2020a) .

Network models consider the contact pattern of human, in which people are regarded as nodes and contacts between people are regarded as links. Xue et al. established a mean-field network model based on node degree to study the spread of the COVID-19 in Wuhan, Toronto and Italy (Xue et al. 2020) . Luo et al. used a network pairwise approximation model to study the transmission of COVID-19 in the Chinese mainland (except Hubei) and Hubei provinces, respectively (Luo et al. 2021) . This model takes into account both the family cluster and the measure of quarantining close contacts, so as to better reproduce the epidemic prevention and control scene during the COVID-19 in China. None of these models takes into account the effect of medical resources. Although Li et al. studied the effect of FSHs on the COVID-19 epidemic in Wuhan and found that FSHs helped slow and eventually stop the COVID-19 in Wuhan (Li et al. 2020a) , they ignored the family cluster and household transmission. In view of the complexity of the COVID-19 prevention and control in Wuhan, none of these models considered the household structure of families of the city due to the complexity.

Considering the significant family cluster transmission of COVID-19, the medical resources in Wuhan and the massive quarantine of close contacts, we build a household model with 123 dimensional to mimic the spread of COVID-19 in Wuhan and to evaluate the impact of medical resources and group quarantine on the spread of COVID-19. The final size of Wuhan is predicted to be 50, 662 (95%CI: 46, 234, 55, 493) , and the epidemic will last until April 25 (95%CI: April 23, April 29). It is worth noticing that the final size in this manuscript means that of reported cases. All simulation results show that the key to a successful control of COVID-19 is firstly to track, test and quarantine close contacts and secondly to increase the number of beds to hospitalize and isolate all the confirmed cases. The two complement to each other and are indispensable.

We collect all the daily reported data on COVID-19 infections, including the confirmed cases (see Fig. 1a ), the recovered cases and death caused by COVID-19 in Wuhan (see Footnote 3). The data are from January 23, 2020, to March 8, 2020. There is a sudden increase in confirmed cases on February 12. This is due to a change in the definition (a) Figure 2 shows the distribution of household size in Hubei Province, China. The household size ranges from one to ten, mostly two and three, and the number of households with seven or more is relatively very small. We will then use the household size distribution of Hubei Province to approximate the distribution of Wuhan. 

After the lockdown of Wuhan City, family members became the main close contacts of the confirmed cases, and the outbreak is characterized by family cluster transmission.

In this scenario, we take households with different sizes as the basic variables and then subdivided these basic variables according to the number of susceptible and infectious within households. There are two ways of transmission: One is within the household, and the other is between the household which is caused by a member of the family going to public places. In view of complex network, Wuhan formed a metapopulation network (as shown in Fig. 3 ). Nodes in this metapopulation network fall into two categories, one representing a household with size from one to ten and the other representing a public place. All public places are regards as a node here. Members in the unquarantined households can move between their homes and the public place. Individuals in the same household are well mixed and are also well mixed in the public place. In addition, households are divided into quarantined households and unquarantined households based on the presence or absence of a confirmed case in the household. The difference between these two households lies in the following: Individuals in a quarantined household cannot move to the public place and can only contact individuals within the household, while for a unquarantined household, one member each day moves to the public place for a period of time every day and then goes home. Individuals' movement between different unquarantined households is ignored in view of strict prevention and control measures. For individuals in Wuhan, this paper divides them into six compartments: susceptible individuals (S) in households, centrally quarantined susceptible (S 0 ), infectious individuals (I ) in households, centrally quarantined infectious individuals (I 0 ), confirmed cases (C) and the removed (recovery or death) (R). Here, an infectious individual is an individual infected with COVID-19 virus and capable of transmitting the virus. I include those infected during the incubation period and those symptomatic infected beyond the incubation period without being confirmed. Although a confirmed individual is also infected with COVID-19 virus, he cannot transmit the virus because The metapopulation network of Wuhan. Nodes in this metapopulation network fall into two categories, one representing a household and the other representing a public place. Households are further divided into two types: quarantined households and unquarantined households based on whether the presence or absence of a confirmed case of being isolated. S in a unquarantined household may be infected by I within the same household at rate β h and may randomly move to the public place and be infected by I from other unquarantined households at rate β g . I in a unquarantined household may become confirmed, then hospitalized and isolated after 1/μ(t) days. Meanwhile, their family members are quarantined at home with probability h, centrally quarantined with probability c, or not quarantined with probability (1 − h − c). For quarantined households, S may only be infected by I within the same households at rate β h ; I may become confirmed, then hospitalized and isolated after 1/μ(t) days, while their family members are still quarantined at home. In centrally quarantined points, there is no contact between individuals, so S 0 will not be infected; I 0 may be confirmed, hospitalized and isolated after 1/μ(t) days. C will be removed after treatment at rate γ . In view of strict prevention and control measures in Wuhan City, all confirmed cases and their close contacts are strictly tracked and recorded even after they are released from isolation or quarantine. So we consider that individuals who are released from isolation or quarantine are no longer involved in transmission. Figure 4 provides a detailed flowchart, and the variables and parameters in the flowchart are explained in Table 1 . According to Fig. 4 , the following household model is established: 

The number of unquarantined households in which the number of S and I is p and q at time t, p + q ∈ [1, 10], p, q ∈ [0, 10]

The number of quarantined households in which the number of susceptible individuals and infectious individuals is p and q at time t, p + q ∈ [1, 9], p, q ∈ [0, 9]

S 0 (t) The number of centrally quarantined susceptible individuals at time t

The number of centrally quarantined infectious individuals at time t

The number of confirmed cases at time t h Household quarantined rate (probability of the family members in a unquarantined household being quarantined at home when

Centrally quarantined rate (probability of the family members in a unquarantined household being centrally quarantined when I become

Transmission rate within household β g Transmission rate in the public place 1/μ(t) The average time from infection to hospitalization (referred to as diagnosis time) of an infectious individual, which depends on medical resources (the number of beds, etc.) γ

Removed rate

In model (1)

(1)], the first and the third terms represent that S are infected by I , including within-household transmission and transmission by I in the public place; the second term represents that I are confirmed and hospitalized; the last term represents that family members are not quarantined when an I of household [S p I q+1 ] is confirmed. For quarantined households [S p I q ] Q [see Eq.

(2)], the first and the third terms represent that S are infected by I in the same household; the second term represents that I are confirmed and hospitalized; the fourth term represents that family members are quarantined at home when an infectious individual of household [S p I q+1 ] is confirmed; the last term represents that family members are quarantined at home when an infectious individual of household [S p I q+1 ] Q is confirmed. The only term in Eq. (3) represents that S in unquarantined households are centrally quarantined when an infectious individual is confirmed. In Eq. (4), the first term represents that I in unquarantined households are centrally quarantined when an infectious individual is confirmed; the second term represents that I are confirmed. In Eq. (5), the first three terms represent that I are confirmed, and the fourth term represents that C are removed.

The parameter 1/μ(t) represents the average time from infection to hospitalization (referred to as diagnosis time) of an infectious individual, which depends on the number of beds, newly symptomatic infectious individuals and confirmed cases, as shown in the following equation:

where B(t) is the number of beds at time t; C n (t) represents the number of newly symptomatic infectious individuals and

1/μ max is the shortest diagnosis time, and it is approximate to the incubation period; 1/μ min is the longest diagnosis time.

, that is, the number of beds is sufficient, I will be confirmed and hospitalized as soon as they become symptomatic, so the diagnosis time is approximately equal to the average incubation period (recorded as the shortest diagnosis time). When B(t) ≤ C(t), that is, the number of beds is extremely short, newly symptomatic infectious individuals have to wait a long time before being confirmed and hospitalized (recorded as the longest diagnosis time). When

, only a fraction of newly symptomatic infectious individuals can be confirmed and hospitalized. In this case, the diagnosis time is between the shortest diagnosis time and the longest diagnosis time.

The basic reproduction number R 0 , which is defined as the average number of secondary cases that one case produces during the course of its infectious period in a totally susceptible population, determines the growing rate of an emerging infectious disease. This definition is for individuals, but the main variables of model (1)-(5) are households. So it is necessary to define the household basic reproduction number. Similar to the definition of R 0 , the household basic reproduction number R * is defined as the average number of secondary households infected by a household with initially one infectious member during the course of its infectious period in a totally susceptible population. The basic reproduction number is defined for the transmission of infection from individual to individual, while the household basic reproduction number for the transmission of infection from household to household. When all households are of size 1, the household basic reproduction number is the same with the basic reproduction number, that is, R * = R 0 . R 0 is usually calculated by the next-generation matrix method (Van den Driessche and Watmough 2002), so we use the same method to calculate R * .

In the early stages of the spread of COVID-19, medical resources are adequate, so μ(t) = μ max . Moreover, there are only unquarantined households [S p I q ] ( p + q ∈ [1, 10], p, q ∈ [0, 10]), so the variables are Adopting the next-generation matrix method (Van den Driessche and Watmough 2002), we have

The first row, diagonal and next diagonal of the block matrix F 55×55 are nonzero matrices; the rest are zero matrices. The diagonal and last diagonal of the block matrix V 55×55 are nonzero matrices, and the rest are zero matrices. In detail,

.

and ρ(A) is the spectral radius of the matrix A. It is difficult to find the explicit expression of R * , but it can be calculated by numerical methods.

Based on the collected data and the MCMC (Markov chain Monte Carlo) method (Gamerman and Lopes 2006) , this section will estimate the parameters of household model, predict the prevalence trend of COVID-19 in Wuhan City and evaluate the effectiveness of prevention and control measures, including quarantine and medical resources.

A total of 152 parameters and initial values of variables are involved in the household model, among which five parameters are derived from the literature, four parameters need to be estimated, and the remaining 143 parameters and initial values need to be pre-calculated.

The five parameters derived from the literature are the shortest diagnosis time (see Footnote 5), the longest diagnosis time (see Footnote 5), the removed rate (Luo et al. 2021) , and the population and the total number of households in Wuhan (see Footnote 13) (see Table 2 ).

The four parameters need to be estimated are the within-household transmission rate β h , the transmission rate in the public place β g , household quarantined rate h and centrally quarantined rate c. Based on the cumulative number of confirmed cases and beds collected from January 23 to March 8 in Wuhan, we adopt the MCMC method (Gamerman and Lopes 2006) for 20,000 iterations with a burn-in of 10,000 iterations to fit the model (1)-(5) to estimate these four parameters. A change in the definition of confirmed cases in Wuhan City on February 12 resulted in a sudden increase in 12,000 confirmed cases (as shown in Fig. 1a) . Under the new definition of confirmed cases, the cumulative number of confirmed cases reported from January 23 to February 11 underestimated the actual number of confirmed cases. Therefore, in the parameter estimation, different weights are applied to the log-likelihood function around February 12 (the weight after February 12 is greater than before February 12) to obtain more realistic results. Now, the parameter values and initial values that need to be calculated are given in the following. Firstly, the number of households H k with a household size of k in Wuhan is calculated according to the probability distribution P(k)(k ∈ [1, 10]) with a household size of k. Considering that the distribution of household sizes in Wuhan is consistent with that in Hubei Province, according to China Statistical Yearbook, the probability distribution of P(k) with a household size of k in Wuhan is shown in Fig. 2 . According to Wuhan Statistical Yearbook, the total number of households in Wuhan is H = 3,108,800, so H k can be calculated by the following equation:

Secondly, the initial values of confirmed cases C(0), I in unquarantined household I (0) and the number of individuals quarantined Q (including quarantined infectious individuals and susceptible individuals) are obtained by simple calculations. C(0) equals the difference between the cumulative number of confirmed cases and the cumulative number of removed individuals (the sum of the number of recovered individuals and deaths) at the initial time. The number of newly confirmed cases in one incubation period can be used to approximate I (0). Although the Wuhan Municipal Health Commission has not announced the number of individuals quarantined, the Hubei Provincial Health Commission announced daily data that individuals are quarantined. 15 According to the ratio of cumulative confirmed cases in Wuhan and Hubei Province, Q can be estimated. Next, we calculate the number of unquarantined households with size k > 0 and I being q > 0. Considering that the larger the population of households with size k is, the greater the number of I is, and letting H I k be the total number of I in all households with size k at the initial time, we have

According to the principle of randomness of infection and clustering of families, we believe that in an unquarantined household with a size of k with infectious individuals, the number of households with different numbers of infectious individuals q(1 ≤ q ≤ k) is equal, which is represented by U (k). For example, when k = 2, the number of unquarantined households with infectious individuals q = 1 is equal to that with q = 2, and both equal to U (2). Then, the total number of infectious individuals in the unquarantined household with size 2 is H I 2 = U (2) + 2U (2). In a similar way, the total number of infectious individuals in the unquarantined household with size k is H E k = U (k) k q=1 q. Thus, when q > 0, the initial values of the unquarantined households behave as [S p I q 

The initial values of quarantined households are calculated in a similar way to that of unquarantined households, which will not be repeated here. All the values of parameters and the initial values are summarized in Tables 2 and 3 .

The formulas involved in Tables 2 and 3 are as follows:

Based on the cumulative number of confirmed cases and beds collected from January 23 to March 8 in Wuhan, we adopt the MCMC method (Gamerman and Lopes 2006) for 20,000 iterations with a burn-in of 10,000 iterations to fit the model (1)-(5) to estimate four parameters (see Table 2 ). A good fitting between the model solution and real data after February 12 is shown in Fig. 5 . On February 12, the definition of confirmed cases in Wuhan changed. Under the new definition of confirmed cases, the reported data before February 12 were lower than the real data, so it is reasonable for the model to overestimate the confirmed cases before February 12. According to Table 2 , withinhousehold transmission rate is larger than the transmission rate in the public place. This is due to more frequent within-household contact and individuals staying within Based on the estimated parameters, under the current prevention and control measures, the household model predicts that the peak size of Wuhan is 17, 356 (95%CI: 16, 222, 18, 571) , that the peak time is on February 15 (95%CI: February 14, February 16), that the final size is 50,662 (95%CI: 46,234, 55,493) and that the time duration is 94 days (95%CI: 92, 98 days) (counting from January 23, 2020). The household basic reproduction number is 11.6540. Compared with the actual situation in Wuhan, the household model does not predict the peak size well. However, the peak time and duration time predicted by the model are only two days away from the actual, and the relative error between the final size and the actual (50,340) is only 0.64%. This indicates that the household model established in this paper is reasonable, and the prediction is more accurate. It is noted that the actual final size of Wuhan did not take into account the imported cases. Figure 6 illustrates the total number and the distribution of I in the unquarantined household with different sizes. Figure 7 illustrates the number of quarantined infectious individuals (including household quarantined infectious individuals and centrally quarantined infectious individuals) and confirmed cases. It is shown that infectious individuals in the unquarantined households are mainly distributed in household with size from 1 to 6. In detail, the high-risk households are mainly with size from 2 to 4 before the peak time, while the households with size 1 become the most at risk after the peak time. After February 3, the number of infectious individuals in the household with size from 2 to 6 drops sharply. At the same time, the number of quarantined infectious individuals and confirmed cases is still increasing. After February 7, the number of quarantined infectious individuals begins to decrease. It is noticed that Huoshenshan Hospital was put into use from February 4 and that FSHs were put into use from February 5. The increase in the number of beds allowed infectious individuals to be confirmed and hospitalized timely. Meanwhile, their close contacts were also (a) Fig. 6 The total number and the distribution of I in the unquarantined households [S p I q ] ( p + q ∈ [1, 10]) The number of individuals Quarantined infectious individuals Confirmed cases timely quarantined. Prompt confirmation of infectious individuals and prompt quarantine of close contacts led to a rapid reduction in the number of infectious individuals in the unquarantined households. The number of infectious individuals in the unquarantined households with size 1 was slow to reduce after February 3. This is mainly because the number of infectious individuals in such household can only be reduced by confirmation. This also reflects the importance of quarantining close contacts.

For emerging infectious diseases, isolating the infectious and quarantining the susceptible can effectively control the spread of the disease from the source. Figure 8 shows the evolution of the cumulative confirmed cases under different quarantine methods for close contacts. Obviously, centralized quarantine of all close contacts is most conducive to disease prevention and control. If close contacts are not quarantined, the infectious disease will continue for 3 years, and eventually about 87% of Wuhan population will be infected. When quarantining close contacts, the centrally quarantined rate and household quarantined rate can characterize the control of the spread of the COVID-19. The greater the sum of these two rates, the more conducive to the control (a) . 8 The effect of quarantining close contacts on the transmission of COVID-19. Subfigure a shows three extreme cases: purple means no close contacts are quarantined; black means all close contacts are quarantined at home; red means all close contacts are centrally quarantined. Subfigure b compares the case where the centralized quarantine rate and the home quarantine rate are both greater than 0, and the sum of the two is 1 and 0.8, respectively (Color figure online) of the spread of the infectious disease is. If the sum of these two rates is a fixed value, the higher the centrally quarantined rate, the more conducive to the control of the spread of disease is. Therefore, in the prevention and control of infectious diseases, the centralized quarantine of close contacts should be achieved to the greatest extent. Even if centralized quarantine of all close contacts cannot be achieved, strict household quarantine for them is also conducive to the prevention and control of infectious diseases.

Massive centralized quarantine of all close contacts is most conducive to disease prevention and control. However, it will take some time to set up centralized quarantine points. In the event of an emerging infectious disease, it will be difficult to put all close contacts under centralized quarantine immediately. Let t 0 indicate the start time of the centralized quarantine for all close contacts (h = 0, c = 1). Figure 9 compares the cumulative number of confirmed cases under five different t 0 . Before reaching peak time, the sooner all close contacts are quarantined, the better the control of the spread of the disease is. However, when t 0 exceeds the peak time, achieving complete centralized quarantine has only a slight impact on the final size and duration time.

Since February 5, 2020, Wuhan has put FSHs with more than 13,000 beds into use. In order to study the necessity of establishing FSHs, this section compares the cumulative confirmed cases and the number of confirmed cases under three cases: whether FSHs are put into use or not, and when the FSHs are established at the same time as city lockdown (as shown in Fig. 10 ). It is shown that the final size with FSHs is reduced by about 17,000. If the establishment time of FSHs was advanced to January 23, then they could be put into use on January 25 under the assumption that it takes 2 days from establishment to the opening of FSHs. This measure reduces the final size by more than 20,000. With the opening of the FSHs and the increase in designated hospitals (including the Huoshenshan and Leishenshan hospitals), the number of beds is increasing. On the one hand, the infectious individuals can be timely confirmed and hospitalized, thus reducing the number of deaths due to disease. On the other hand, the prompt confirmation of infectious individuals contributes to the prompt quarantine of close contacts, thus reducing the number of confirmed cases. Based on the daily reported data on infections, it is estimated that the mortality rate due to disease in Wuhan is about 5% by the ratio of the number of cumulative deaths to the cumulative confirmed cases. In order to study the impact of medical resources (number of hospital beds) on the number of deaths, the final size is calculated with a fixed number of hospital beds. The resulting final size is then multiplied by the mortality rate to give a rough estimate of the number of deaths from disease. Figure 11 illustrates the relationship between the number of beds and deaths. The number of deaths is 2512 according to the final size predicted by the household model. As can be seen from Fig. 11 , when the number of beds is 8150, the number of deaths is approximate to the estimated value. An increase of just 50 beds could prevent about 360 deaths. The more the beds, the fewer the deaths. But when the number of beds exceeds 11,000, the number of 

Considering of the prevention and control measures taken after the city lockdown on January 23, 2020, we built a household model with medical sources and quarantining of close contacts. It is found that about 9% of households with confirmed cases are not effectively quarantined. The household model predicted that the final size of Wuhan is 50, 662 (95%CI: 46, 234, 55, 493) and that the outbreak will not end until April 25 (95%CI: April 23, April 29). The predicted peak time and duration time are only 2 days away from the actual ones, and the relative error of the final scale is only 0.64%. The household basic reproduction number was estimated to be 11.6540. It looks overestimated, but it is actually reasonable. The value of household basic reproduction number depends not only on household sizes but also on the value of the basic reproduction number (Ball et al. 1997; Becker and Dietz 1995; Pellis et al. 2011 ). According to Ref. Li et al. (2020b, c) , Nishiura et al. (2020) and Wu et al. (2020b) , the basic reproduction number of COVID-19 in Wuhan was estimated to be in the interval [2.23, 6.6037]. According to China Statistical Yearbook and Wuhan Statistical Yearbook, the average household size in Wuhan is 2.9373. For the initially infectious household with an infectious member, which is with the average size 2.9373, the other members soon become infectious since R 0 ∈ [2.23, 6.6037]. This process takes a very short time. It is equivalent to all members of the household being infectious at the beginning. Then, each member of the infectious household will infect other households. According to the definition of the household basic reproduction number, a household with the average size 2.9373 may infect 6.5502-19.3970 households. Hence, the estimated household reproduction number was reasonable.

By studying the effects of different times and degrees of quarantine of close contacts on the cumulative confirmed cases, we found that not quarantining close contacts will eventually result in about 87% of infections. Moreover, it is found that putting FSHs into use reduced the final size by about 17,000. If the establishment time of FSHs started on January 23, 2020, immediately after the lockdown, the final size may be reduced by 20,000. If the medical resources are adequate at the time of Wuhan lockdown, the number of deaths may be reduced by nearly 54.26%. It is also shown that the high-risk households are mainly with size from 2 to 4 before the peak time, while the households with size 1 become the most at risk after the peak time. Based on the above analysis results, the keys to the success of COVID-19 prevention and control in Wuhan are to adopt tracking and quarantine measures for close contacts and to open the FSHs and increase DHs for isolating and curing confirmed cases.

Compared with other mathematical models for COVID-19, the dimension of this household model (123) is very high, but the advantages of the model are significant. On the one hand, this model reflects the characteristics of family cluster of COVID-19 which is an extension of the modeling in Li et al. (2020a) . On the other hand, this model can depict the evolution of infectious individuals in household with different sizes and indicate high-risk household type. Moreover, fine division of compartment allows this model to well reproduce the process of COVID-19 transmission and control in Wuhan. The conclusions obtained are more credible. Our findings of the successful experience in Wuhan can help the prevention and control of similar emerging infectious diseases and also provide a good reference for COVID-19 prevention and control in other places of the world.

This paper also has some limitations. Taking into account the strict implementation of Wuhan's epidemic prevention and control measures, this paper believes that people who have been quarantined or isolated (whether susceptible people are released from quarantine or confirmed cases are removed after hospitalization) will no longer participate in transmission. This is a simplification of the household model. When individuals in Wuhan are divided into susceptible S, infectious I , confirmed C, recovered R and dead D, the flowchart is shown in Fig. 12 . The household is divided into quarantined household and unquarantined household based on the presence or absence of quarantined individuals. [S 0 s I 0 i C c R r D d ] is defined as the number of unquarantined households in which the number of S, I , C, R, D is s, i, c, r , d, respectively. Here, (s, i, c, r , d) ∈ {(s, i, c, r , d) ∈ N 5 | s + i + c + r + d ∈ [1, 10]}. [S s 0 I i 0 C c R r D d ] denotes the number of quarantined households in which the number of S, I , C, R, D is s, i, c, r , d, respectively. Here, (s, i, c, r , d) ∈ {(s, i, c, r , d) ∈ N 5 | s +i +c+r +d ∈ [2, 10], s + i ≥ 1, c + r + d ≥ 1}. For this household model with detailed division, its dimension is more than 5600. It can also be used to study the transmission of COVID-19. However, such a high dimension makes it more difficult to select the initial values of the variables and makes the program run slower. Therefore, it is necessary to simplify the model reasonably. The model (1)-(5) is one of the reasonable simplifications of this model. According to China Statistical Yearbook, household sizes in China are divided into 10 categories: 1, 2, 3, 4, 5, 6, 7, 8, 9, and 10 and above. The number of households with size larger than 10 is counted together with those with size 10, and it is small compared to that of other sizes. In addition, considering that gatherings of more than 10 people were prohibited during COVID-19, we set the size of households

Epidemics with two levels of mixing

The effect of household distribution on transmission and control of highly infectious diseases

A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster

Fangcang shelter hospitals: a novel concept for responding to public health emergencies

The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak

Markov chain Monte Carlo: stochastic simulation for Bayesian inference

Clinical features of patients infected with 2019 novel coronavirus in Wuhan

Fangcang shelter hospitals during the COVID-19 epidemic

Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2)

Basic reproduction number and predicted trends of coronavirus disease 2019 epidemic in the mainland of China

Nonpharmaceutical interventions contribute to the control of COVID-19 in China based on a pairwise model

Epidemic growth rate and household reproduction number in communities of households, schools and workplaces

Why is it difficult to accurately predict the COVID-19 epidemic?

Estimation of the transmission risk of 2019-nCov and its implication for public health inverventions

Initial cluster of novel coronavirus (2019-nCoV) infections in Wuhan, China is consistent with substantial human-to-human transmission

Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission

Real-time tentative assessment of the epidemiological characteristics of novel coronavirus infections in Wuhan, China, as at 22

Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study

Impact of international travel and border control measures on the global spread of the novel 2019 coronavirus outbreak

A data-driven network model for the emerging COVID-19 epidemics in Wuhan, Toronto and Italy

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

The authors declare that they have no conflict of interest.

 Fig. 12 The flowchart of household model with more detailed division. The black solid lines represent infection; the black dashed lines are caused by the transformation of C; the yellow lines mean that I is confirmed and that his family members are quarantined; the gray line indicates the release of quarantine (Color figure online) with 10 and above to be 10. We will collect more detailed data and expand household size in the future work. Although the household model established in this paper is for the megacity of Wuhan, it could be extended to study of COVID-19 in other countries and regions. When applied to other provinces and cities in China, medical resources are adequate, so μ(t) should be constant; there is no household quarantine, so h = 0, c = 1 and [S p I q ] Q = 0. When analyzing other countries, the status of family members and the household sizes may vary. The detailed division will result in higher model dimensions which is more complicated; we will keep it as future work.