Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

When will the battle against novel coronavirus end in Wuhan: a SEIR modeling analysis

Kangkang Wan, Jing Chen, View ORCID ProfileChangming Lu, Lanlan Dong, Zhicheng Wu, View ORCID ProfileLianglu Zhang
doi: https://doi.org/10.1101/2020.02.16.20023804
Kangkang Wan
1Wuhan Ammunition Life-Tech Co., Ltd., Wuhan, Hubei, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jing Chen
2Clinics of Oilcrops Research Institute, CAAS, Wuhan, Hubei, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Changming Lu
1Wuhan Ammunition Life-Tech Co., Ltd., Wuhan, Hubei, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Changming Lu
Lanlan Dong
1Wuhan Ammunition Life-Tech Co., Ltd., Wuhan, Hubei, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhicheng Wu
1Wuhan Ammunition Life-Tech Co., Ltd., Wuhan, Hubei, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lianglu Zhang
1Wuhan Ammunition Life-Tech Co., Ltd., Wuhan, Hubei, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Lianglu Zhang
  • For correspondence: zhangll{at}whammunition.com
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Data/Code
  • Preview PDF
Loading

Abstract

Recent outbreak of 2019-nCoV in Wuhan raised serious public health concerns. By February 15, 2020 in Wuhan, the total number of confirmed infection cases has reached 37,914, and the number of deaths has reached 1123, accounting for 56.9% of the total confirmed cases and 73.7% of the total deaths in China. People are eager to know when the epidemic will be completely controlled and when people’s work and life will be on the right track. In this study we analyzed the epidemic dynamics and trend of 2019-nCoV in Wuhan by using the data after the closure of Wuhan city till February 12, 2020 based on the SEIR modeling method. The optimal parameters were estimated as R0=1.44 (interquartile range: 1.40-1.47),TI=14 (interquartile range: 14-14) and TE=3.0 (interquartile range: 2.8-3.1). Based on these parameters, the number of infected individuals in Wuhan city may reach the peak around February 19 at about 45,000 people. Once entering March, the epidemic would gradually decline, and end around the late March. It is worth noting that the above prediction is based on the assumption that the number of susceptible population N = 200,000 will not increase. If the epidemic situation is not properly controlled, the peak of infected number can be further increased and the peak time will be a little postponed. It was expected that the epidemic would subside in early March, and disappear gradually towards the late March.

Introduction

Wuhan is the largest city in central China with a total population of more than 11 million[1]. The epidemic of 2019-nCoV pneumonia has been raging in the whole country, especially in Hubei province for nearly a month. In late December 2019, 67 cases of 2019-nCoV pneumonia were reported in Wuhan[2]. In order to prevent the further spread of 2019-nCoV, Wuhan began to close the city from 10:00 on January 23, banning all vehicles from entering and leaving the city. Tens of thousands of medical staff, soldiers and people from all walks of life have been involved in the campaign.

The spread of the epidemic has caused a huge threat to people’s health and life safety, at the same time, it has a serious impact on China’s social life and national economy. By February 15, 2020, the total number of confirmed cases has reached 37914, and the number of deaths has reached 1123 in Wuhan, accounting for 56.9% of the total confirmed cases and 73.7% of the total deaths in China[3]. With the increase of medical staff from all over the country, the opening of several large novel hospitals, and the adoption of anti epidemic measures, more patients can get efficient and timely treatment. The number of confirmed cases increased sharply on February 12 and 13, while the total number of suspected cases decreased gradually[3].

People are eager to know when the epidemic will be completely controlled and when people’s work and life will be on the right track. In order to help the public to understand the future trend of the epidemic, we analyzed the epidemic dynamic and trend of 2019-nCoV in Wuhan city by using the SEIR modeling method based on the actual data and published references.

Methods

1. Epidemic transmission model

The SEIR model is a classical epidemic model for the flows of people between four states: susceptible (S), exposed (E), infected (I), and recovery (R). Each of those variables represents the number of people in those groups. The relationship among the four groups is elucidated in Figure 1, where β1 is the probability of S to E after I contacts S, γ1 is the probability of E to I, and γ2 is the probability of I to R. Since 2019-nCoV is also infectious in the incubation period, we introduced parameter β2 here to represent the probability of S to E after E contact S. We used the “susceptible – exposed – infected – recovered” model[4] to describe the prevalent characteristics of 2019-nCoV in Wuhan.

Fig. 1.
  • Download figure
  • Open in new tab
Fig. 1. Relationship among four groups according to SEIR model

This is an ordinary differential equation model, described by the following equations: Embedded Image Among which, S(t), E(t), I(t) and R(t) represent the number of people in the group of the susceptible, the exposed, the infected and the recovered on the day t, respectively. N is the total number of possible contact people, which is assumed to be fixed and N = S +E + I + R.

The simulation used the fourth-order Runge-Kutta algorithm[5] to solve it numerically, with a step size fixed at 0.01 for R0, and at 0.1 for TE and TI.

2 Estimation of parameters for the model

Parameters β1, β2, γ1 and γ2 were estimated according to the reference[6] using the formula below: Embedded Image

Among which, R0 is the basic reproduction number, TI is the time of infectious period, TE is the time of incubation period, and β1, β2, γ1 and γ2 share the same meanings as in Figure 1.We set the range of TE to 1-7 days, TI to 1-14 days, R0 to 1-5 days according to the reference[7].

The optimized parameters of TE, TI and R0 were determined on the condition of RMSE being the smallest. Embedded Image

Where I and Î indicate the real and simulated number of the infected, respectively while the R and Embedded Image indicate the real and simulated number of the recovery in a specific time, respectively.

Calculate the real number of people (I and R) and the smallest RMSE of Î and Embedded Image obtained by the model, and choose the TE, TI and R0 with the smallest RMSE as the optimal model parameters.

In order to avoid over-fitting in parameter optimization, we took 80% data each time, and repeated this process 100 times, and finally took the average value of R0, TI and TE of 100 times as the optimal parameter values.

3 Data source

The data were collected from the official website of Hubei Provincial Health Committee (http://wjw.hubei.gov.cn/)[3], and shown in Table 1. We used the data of 22 days from January 22 to February 12 when Wuhan city was shut down and all the public transportation was suspended.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 1.

Number of the infected and the recovery at specific date in Wuhan

For construction of the model, data of 22 days were divided into two stages. The first stage is from January 23 to February 7, and the second stage is from February 8 to February 12. During the second stage, Wuhan city took a number of measures, including timely diagnosis, timely treatment and effective isolation of the infected population, which will have an important impact on the parameters of the model.

4 Initial parameter settings

To establish the model, we firstly estimated the parameters of the susceptible (S), the exposed (E), the infected (I) and the recovery(R) based on the latest data issued on February 12:

N = 200000, which is the total number of potential close contacts in Wuhan on February 12.

S=N-I, in which S is the number of the susceptible and I is the number of the infected.

I0 = 425, which is the number of susceptible individuals at the beginning of the model run.

E0 = 426, which is the number of exposed individuals at the beginning of the model run

R0 = 28, which is the number of recovered individuals at the beginning of the model run

Results and Analysis

Epidemic prediction based on SEIR model

The epidemic of the novel coronavirus pneumonia in Wuhan was studied by SEIR modeling. The results showed that, at the time when Wuhan was closed, the number of initially infected individuals was I0 = 425, the number of initially exposed individuals was E0 = 426, and the number of initially recovered patients was R0 = 28.

Next, we separated the data into two stages: January 22-February 7 and February 8-February 12. In the first stage, TI = 14(interquartile range:14-14), TE = 3.0 (interquartile range : 2.8-3.1), R0 = 1.44 (interquartile range : 1.40-1.47) (Figure 2). The data showed that the infectious time of the infected person (I) is 14 days, and the incubation period is about 3 days, which is close to the data (Interquartile range: 3.0-7.2) estimated in the reference [8]. The propagation base R0 of this study is 1.44, which is significantly lower than the R0 estimated by other papers before the closure of Wuhan[9-11].

Fig 2.
  • Download figure
  • Open in new tab
Fig 2. Distribution of the R0 (A), TI (B) and TE (C) for 100 random samplings

In the second stage (after February 8), we set the number of susceptible population N to be fixed at 200000, and the infection cycle of infected population decreased from 14 days to 4 days, i.e. TI = 4, which was estimated according to the data of 5 days from February 8 to February 12, so as to get the epidemic development trend of 90 days since January 22, including the number of infected people, the number of latent people and the number of recovered people (Figure 3, Table 3). The results showed that the number of infected people increased slowly in the early stage (January 22-january 31), but during February 1-February 12, the number of infected people increased rapidly, which is expected to peak around February 19, reaching about 45000 people. Subsequently, the number of infections will decrease. Once entering March, the epidemic would gradually decline, and the epidemic would end around the end of March. It is worth noting that the above prediction is based on the assumption that the number of susceptible population N = 200000 will not increase.

View this table:
  • View inline
  • View popup
Table 2.

Distribution of the R0 (A), TI (B) and TE (C) for 100 random samplings

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 3.

Expected number of cases for infected, recovery and exposed

Fig 3.
  • Download figure
  • Open in new tab
Fig 3. Epidemic trend of 2019-nCoV within 90 days after the closure of Wuhan

In Figure 3 and Figure 4, red line indicates the trend of cumulative infection number over time, blue line is the trend of cumulative rehabilitation number over time, and green line is the trend of cumulative latent number over time. Vertical dash line indicates the peak time of cumulative infection number.

Fig 4.
  • Download figure
  • Open in new tab
Fig 4. Epidemic trend of 2019-nCoV within 90 days after the closure of Wuhan city assuming the number of susceptible people n = 300000 (A) and 400000 (B)

If the epidemic situation is not properly controlled, the number of susceptible population will continue to increase on the basis of current N = 200000. If the number of susceptible population increases to N = 300000, and other parameters remain unchanged, the peak number can increase to 75000, and the epidemic peak time will also be postponed at around February 21 (Figure 4). If it is increased to N = 400000 and other parameters remain unchanged, the peak number can be increased to 100000, and the epidemic peak will be postponed to around February 22 (Figure 4). Even in both cases, the epidemic would subside in early March, and disappear gradually towards the late March.

Discussion

Although some modeling studies on the epidemiological characteristics of 2019-nCoV epidemic have been reported so far, they had some limitations, such as the data come from the early stage of the epidemic. Due to the rapid change of the epidemic situation and the closure of Wuhan on January 23, many parameters related to the model have also changed, which affect the applicability and reliability of the model. This study used the latest 2019 nCoV data in Wuhan area, analyzed the epidemiological characteristics of 2019 nCoV epidemic after in Wuhan city was shut down. Compared with other studies, the R0 value produced in this study is smaller, indicating that the closure and subsequent measures have played an important role in the spread of the epidemic.

The infection time index (TI) obtained in this study was higher than that of SARS[12] and MERS[13], but lower than that of 2019-nCoV in literatures[14] reported earlier. This result may be related to the sudden outbreak of the epidemic, the lack of medical resources for early response, and the failure of timely diagnosis and treatment of infected patients. A large number of mild patients and asymptomatic virus carriers were not isolated in time. The incubation period (TE) is about 3 days, which is close to the data in the reference[14].

According to the latest reported data, the cumulative number of people infected on February 13 and 14 was 35991 and 37914 respectively, which is close to the number predicted by our estimation (Table 2). According to this study, the number of infected people will reach the peak in February 19 at about 45000 infected individuals.

Conclusions

With the implementation of more follow-up measures, including strict restrictions on people going out, accelerating the treatment of infected individuals, and clinical trials of new drugs, the development of 2019-nCoV epidemic in Wuhan will be effectively controlled, and the number of infected individuals will gradually decrease. It was expected that the epidemic would subside in early March, and disappear gradually towards the late March. If the epidemic situation is not properly controlled, the peak of infected number can be further increased and the peak time will be a little postponed.

Data Availability

The data were collected from the official website of Hubei Provincial Health Committee (http://wjw.hubei.gov.cn/)

http://wjw.hubei.gov.cn/

Declarations

Ethics approval and consent to participate

The ethical approval or individual consent was not applicable.

Availability of data and materials

All data and materials used in this work were publicly available.

Consent for publication

Not applicable.

Funding

This work was not funded.

Disclaimer

The funding agencies had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; or decision to submit the manuscript for publication.

Conflict of Interests

The author declared no competing interests.

Reference

  1. 1.↵
    Wee, Rolando Y. “The Largest Cities In China.” World Atlas, Jul. 9, 2019, worldatlas.com/articles/20-biggest-cities-in-china.html.
  2. 2.↵
    W. H. Organization, Coronavirus. World Health Organization, Available: https://www.who.int/health-topics/coronavirus (cited January 19, 2020).
  3. 3.↵
    The official website of Hubei Provincial Health Committee (http://wjw.hubei.gov.cn/) (cited February 19, 2020).
  4. 4.↵
    Pastor-Satorras R, Castellano C, Mieghem PV, et al. Epidemic processes in complex networks. Review of Modern Physics 2015, 87(3): 925–979.
    OpenUrl
  5. 5.↵
    Wallinga J, Lipsitch M, How generation intervals shape the relationship between growth rates and reproductive numbers. Proceedings of the Royal Society B: Biological Sciences, 2007, 274(1609):599–604.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    Y Yang, Q Lu, M Liu, et al. Epidemiological and clinical features of the 2019 novel coronavirus outbreak in China. medRxiv 2020.02.10.20021675; doi :https://doi.org/10.1101/2020.02.10.200216.
  7. 7.↵
    Qun Li, Xuhua Guan, Peng Wu et al., Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia, New Engl J Med, January 29, 2020, DOI: 10.1056/NEJMoa2001316
    OpenUrlCrossRefPubMed
  8. 8.↵
    Y Yang, Q Lu, M Liu, et al. Epidemiological and clinical features of the 2019 novel coronavirus outbreak in China, medRxiv, 2020, January 29, 2020 DOI: 10.1056/NEJMoa2001316
    OpenUrlCrossRef
  9. 9.↵
    T Zhou, Q Liu, Z Yang, et al. Preliminary prediction of the basic reproduction number of the Wuhan novel coronavirus 2019-nCoV,arXiv, 2020 - arxiv.org.
  10. 10.
    JM Read, JRE Bridgen, DAT Cummings,Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions, medRxiv 2020.01.23.20018549; doi : https://doi.org/10.1101/2020.01.23.20018549
  11. 11.↵
    Z Cao, Q Zhang, X Lu, et al. Incorporating Human Movement Data to Improve Epidemiological Estimates for 2019-nCoV. medRxiv 2020.02.07.20021071; doi : https://doi.org/10.1101/2020.02.07.20021071
  12. 12.↵
    Donnelly CA, Ghani AC, Leung GM, et al. Epidemiological determinants of spread of causal agent of severe acute respiratory syndrome in Hong Kong. Lancet 361, 2003: 1761–1766.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    Assiri A, McGeer A, Perl TM, et al., Hospital outbreak of Middle East respiratory syndrome coronavirus. New Engl J Med 2013;369(5):407–16
    OpenUrlCrossRefPubMedWeb of Science
  14. 14.↵
    JT Wu, K Leung, GM Leung, Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study, The Lancet, 2020, https://doi.org/10.1016/S0140-6736(20)30260-9
Back to top
PreviousNext
Posted February 18, 2020.
Download PDF
Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
When will the battle against novel coronavirus end in Wuhan: a SEIR modeling analysis
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
When will the battle against novel coronavirus end in Wuhan: a SEIR modeling analysis
Kangkang Wan, Jing Chen, Changming Lu, Lanlan Dong, Zhicheng Wu, Lianglu Zhang
medRxiv 2020.02.16.20023804; doi: https://doi.org/10.1101/2020.02.16.20023804
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
When will the battle against novel coronavirus end in Wuhan: a SEIR modeling analysis
Kangkang Wan, Jing Chen, Changming Lu, Lanlan Dong, Zhicheng Wu, Lianglu Zhang
medRxiv 2020.02.16.20023804; doi: https://doi.org/10.1101/2020.02.16.20023804

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Epidemiology
Subject Areas
All Articles
  • Addiction Medicine (349)
  • Allergy and Immunology (668)
  • Allergy and Immunology (668)
  • Anesthesia (181)
  • Cardiovascular Medicine (2648)
  • Dentistry and Oral Medicine (316)
  • Dermatology (223)
  • Emergency Medicine (399)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (942)
  • Epidemiology (12228)
  • Forensic Medicine (10)
  • Gastroenterology (759)
  • Genetic and Genomic Medicine (4103)
  • Geriatric Medicine (387)
  • Health Economics (680)
  • Health Informatics (2657)
  • Health Policy (1005)
  • Health Systems and Quality Improvement (985)
  • Hematology (363)
  • HIV/AIDS (851)
  • Infectious Diseases (except HIV/AIDS) (13695)
  • Intensive Care and Critical Care Medicine (797)
  • Medical Education (399)
  • Medical Ethics (109)
  • Nephrology (436)
  • Neurology (3882)
  • Nursing (209)
  • Nutrition (577)
  • Obstetrics and Gynecology (739)
  • Occupational and Environmental Health (695)
  • Oncology (2030)
  • Ophthalmology (585)
  • Orthopedics (240)
  • Otolaryngology (306)
  • Pain Medicine (250)
  • Palliative Medicine (75)
  • Pathology (473)
  • Pediatrics (1115)
  • Pharmacology and Therapeutics (466)
  • Primary Care Research (452)
  • Psychiatry and Clinical Psychology (3432)
  • Public and Global Health (6527)
  • Radiology and Imaging (1403)
  • Rehabilitation Medicine and Physical Therapy (814)
  • Respiratory Medicine (871)
  • Rheumatology (409)
  • Sexual and Reproductive Health (410)
  • Sports Medicine (342)
  • Surgery (448)
  • Toxicology (53)
  • Transplantation (185)
  • Urology (165)