Data-Based Analysis, Modelling and Forecasting of the novel Coronavirus (2019-nCoV) outbreak

Abstract

Since the first suspected case of coronavirus diesease-2019 (COVID-19) on December 1st, 2019, in Wuhan, Hubei Province, China, a total of 40,235 confirmed cases and 909 deaths have been reported in China up to February 10, 2020, evoking fear locally and internationally. Here, based on the publicly available epidemiological data (WHO, CDC, ECDC, NHC and DXY) for Hubei from January 11 to February 10, 2020, we provide estimates of the main epidemiological parameters, i.e. the basic reproduction number (R₀) and the infection, recovery and mortality rates, along with their 90% confidence intervals. As the number of infected individuals, especially of those with asymptomatic or mild courses, is suspected to be much higher than the official numbers, which can be considered only as a sample of the actual numbers of infected and recovered cases in the population, we have repeated the calculations under a second scenario that considers twenty times the number of confirmed infected cases and forty times the number of recovered, leaving the number of deaths unchanged. Our computations and analysis were based on a mean field Susceptible-Infected-Recovered-Dead (SIRD) model. Thus, based on the SIRD simulations, the estimated average value of R₀ in both scenarios was found to be 2.5, the one computed considering the period from the 11th of January until the 18th of January, using the official counts of the confirmed cases was found to be 4.6, while the one computed under the second scneario was found to be 3.2. Furthermore, on the estimated parameters from both scenarios, we provide tentative three-week forecasts of the evolution of the outbreak at the epicenter. Our forecasting flashes a note of caution for the presently unfolding outbreak in China. Based on the official counts for confirmed cases, the simulations suggest that the cumulative number of infected will surpass 80,000 (as a lower bound) and could reach 160,000 (with an upper bound of 330,000) by February 29. Regarding the number of deaths, simulations forecast that on the basis of the up to the 10th of February data and estimations of the actual number of infected and recovered in the population, the death toll might exceed 7,000 by February 29. However, our analysis further reveals a significant decline of the mortality rate to which various factors may have contributed, such as the severe control measures taken in Hubei, China (e.g. quarantine and hospitalization of infected individuals), but mainly it is due to the fact that the actual number of infected and recovered cases in the population is estimated to be more than twenty times higher than reported, thus resulting in a much lower mortality rate, which according to our computations is of the order of 0.15%.