Next Article in Journal
Analysis of the Possibility of Transport Mode Switch: A Case Study for Joinville Students
Previous Article in Journal
Framework for Risk Identification of Renewable Energy Projects Using Fuzzy Case-Based Reasoning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

COVID-19 Healthcare Planning: Predicting Mortality and the Role of the Herd Immunity Barrier in the General Population

by
Julio Emilio Marco-Franco
1,*,
Natividad Guadalajara-Olmeda
1,
Silvia González-de Julián
1 and
David Vivas-Consuelo
2
1
Research Centre for Economics Engineering, Universitat Politècnica de València, 46022 Valencia, Spain
2
Research Unit for Health Economics and Management, Universitat Politècnica de València, 46022 Valencia, Spain
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(13), 5228; https://doi.org/10.3390/su12135228
Submission received: 11 May 2020 / Revised: 22 June 2020 / Accepted: 24 June 2020 / Published: 27 June 2020

Abstract

:
Using a mathematical model for COVID-19 incorporating data on excess of mortality compared to the corresponding period of the previous year obtained from the daily monitoring of mortality in Spain (MoMo), the prediction of total number of casualties in Spain for the first outbreak has been computed. From this figure, and following a stepwise meta-analysis of available reports, the case fatality rate (CFR) and the infectious case fatality rate (IFR) for the outbreak have been estimated. As the impact of age on these rates is notable, it is proposed to include an age-related adjusted fatality ratio in future comparative analyses between studies, calculated by adjusting the results by risk ratio to a reference age band (e.g., 60–69). From the casualty figures, and the corresponding CFR and IFR ratios, the forecast of serologically positive cases in the general Spanish population has been estimated at approximately 1% (0.87–1.3%) of the samples. If the data are confirmed by the ongoing study of the Carlos III Institute, until a vaccine is found, the immunity acquired in the general population after the infectious outbreak is far from the 65–70% herd immunity required as a barrier for COVID-19.

1. Introduction

In December 2019, China reported 27 cases of pneumonia of unknown origin. The agent causing this pneumonia was identified as a new virus named SARS-CoV-2. The disease associated with this virus has been called COVID-19. On March 11, the World Health Organization (WHO) declared the pandemic as global. Since the beginning of the epidemic to the date of this paper (April 20), more than two and a half million cases have been reported worldwide, and more than 200,000 in Spain.
Unlike other pandemics that originated in less economically developed nations, COVID-19 has emerged in several of the world′s largest economies, spreading across the globe causing enormous disruption to economic, working, and social life. To analyze the impact on healthcare economics and policies, accurate and reliable data are essential.
As has been recognized by the health authorities themselves, the official data on COVID-19 are incomplete and inconsistent. Reported values are far from what may be obtained with mathematical modeling [1,2]. For that reason, we have elaborated our data based on more reconcilable information, the mortality. According to the daily monitoring of mortality in Spain report from the Ministry of Justice (MoMo), there has been an unexpected increment of mortality figures when compared with the same period in 2019 [3]. Given that there are currently no other epidemic diseases, we assume that this increment corresponds to undetected cases of COVID-19, directly or indirectly. Based on this hypothesis, we estimated the expected casualties, CFR (case fatality rate), IFR (infectious case fatality rate), and the percentage of population which may result in positive serological testing for COVID-19 in Spain, using a mathematical model, and a meta-analysis with easy obtainable data from previous reports. This mathematical model is applied to the infectious outbreak occurring in Spain starting on March 9. This proceeding may be especially useful for making estimates in the event of a disease outbreak in the coming months. This estimation is used as a guide to consider the convenience of health policies based on herd immunity in Spain.

2. Materials and Methods

Two different mathematical procedures are used. The first one starts on the date with more than 30 new casualties per day and extends to the calculated point with less than 30 new deaths per day, assuming neither significant changes in the evolution process nor rebounds and a lockdown of 45 days. The number of deaths at a given time t (today) is the sum of the past infections weighted by their probability of death, where the probability of death depends on the number of days since infection.
Using a previous reported methodology [4], the expected number of deaths, 𝑑𝑡, on a given day t, is given by the following discrete sum:
d t = τ = 0 t 1 c τ π t τ
where 𝑐τ is the number of new infections on day τ and where π t τ is the probability of death for day t for those getting infected on day τ; it may be discretized via
π s = s 0.5 s + 0.5 π m ( τ ) d τ
for s = t–τ for s = 2, 3, …ending when the daily casualties are <30, and
π 1 = 0.5 1.5 π m ( τ ) d τ
In forecasting new outbreak data at a time t, i(t), it could be interesting to consider the use of time-related exponential growth rate (r) models [5].
i ( t ) = i 0   e r t
where i0 is the expected number of infected cases at time t = 0.
The cumulative incidence I(t) is the integral of i(t) over the period 0 – t
I ( t ) = 0 t i ( s ) d s = 0 t i 0   e r ( t τ )   d ( t τ ) = i 0 ( e r t 1 ) / r
The cumulative incidence may be adjusted to the date of report by a factor u dependent on the parameters of the delay distribution. For estimating the distribution of time delay from onset of disease to death, the authors have used correct truncation and modeling of a log-normal distribution [6].
In case of a log-normal distribution f ( t ; θ ) , with parameters θ d { a d ; b d } , the factor u ( r ; θ d ) is the multiplying parameter for adjustment of I(t) by date of report t, to the time from onset to death. The factor u results from [5]:
u ( r ; θ d ) = 0 e x p ( r s ) f ( s ; θ d ) d s
It is also possible to evaluate the effects of lockdown with transmission models using a Bayesian framework and jointly infer parameters, as have done the French Pasteur group [7] (p.12), in our case for 45 days of lockdown. Other models may analyze the serial interval (the time between symptom onset of a primary and secondary cases). However, due to the uncertainties about real infected cases, this mathematical approach has been used only as a test to confirm if the approximate data obtained with the meta-analysis were consistent.

3. Results

According to information from MoMo [3], unexpected versus real death values reported for the period from March 17 to April 18 were respectively 25,907 and 63,676 (all-death causes), representing 40.7% of the total (Table 1). The expected value of 37,769 is consistent with the National Institute of Statistics of Spain (INE) 2018 [8] report of a daily average of 1172 deaths. COVID-19 deaths reported by the Ministry of Health (corresponding to April 18) were 20,043 cases (Ministry of Health, Spain, daily release information: https://www.mscbs.gob.es/profesionales/saludPublica/ccayes/alertasActual/nCov-China/home.htm (accessed April 19, 2020)).
With the adjusted values of casualties of 45 days, starting on March 9, a prediction curve for the estimated period of pandemic has been elaborated (Figure 1).
The evolution follows a right skewed curve (skewness = 0.53, Kurtosis = −1.10) with mean = 9.4058 and standard error of the mean (SEM) = 1.003517915. This study is point-estimate focused on cumulative data (casualties) at the end of the outbreak (defined in this study as less than 30 new daily deaths in Spain).
If the total predicted period of this (first) outbreak of 69 days is correct (Figure 1), the adjusted overall total mortality in Spain would be about 30,568 cases (27,307–33,830). This would correspond to overall casualties for this outbreak of 649 per million population (PMP) (0.0649%). [9]
The estimation for regions (CCAA), based on the relative ratio of reported casualties PMP, is presented in Table 2, including total number of casualties and case fatality rate adjusted PMP. This computation assumes a CFR equally distributed among the regions, and only dependent on population, something unproved at present time.
For the meta-analysis procedure to estimate seroprevalence, it is necessary to determine the proportion of symptomatic versus asymptomatic patients. The two key figures, CFR and IFR, are forecasted based on previously available reports.
The report from the Diamond Princess cruise ship, where an outbreak occurred and which was quarantined from January 20 to February 29, 2020 [10] is a very revealing paper. From a total of 3711 people on board (passengers and crew), 705 became sick (19.0%), and seven died (a casualty in the sample of 0.003%). On February 20, 2020, 3063 PCR tests were performed with a positive result of a total of 634 people (20.7%), with 476 of them over 60 years. Of the 634 confirmed cases, approximately half of them where asymptomatic (306). Of these 634 cases, 313 cases were female. The age pattern was: 6 aged 0–19 years (0.94%), 152 aged 20–59 years (23.98%) and 476 aged 60 years and older (75.08%), with a genetic diversity from 28 countries [11]. The data were statistically modeled, and observations treated as survival data with right censoring. The probability of being asymptomatic once infected and the infection time for each case were estimated using a Hamiltonian Monte Carlo algorithm. The estimated total number of true asymptomatic cases was 113.3 (95% credible interval: 98.2–128.3%) and the estimated asymptomatic proportion among all infected cases 17.9% (95% credible interval: 15.5–20.2%). The results are summarized in Table 3.
This reported probability of death following the development of symptoms (1.35%) is very close to that published in China, analyzing 79,394 cases, with 1.4% (0.9–2.1%) [12] and to the results (1.6%) in another report [13]. Starting from the death figure, this allows us to estimate the number of symptomatic cases. However, as seen, the crude mortality percentage is about 4 times higher than that previously found for the general population of 0.06%. This underlines the importance of age band. As mentioned, more than 75% of the PCR-detected cases were aged 60 or older. All 7 deaths reported from the Diamond Princess cruise were patients of 70 years or more [14].
In the second study used in the meta-analysis, one from China with 72,314 case records [15], CFR increased to 8% in patients aged 70–79 years, and 14.8% in patients aged 80 or more years, with an overall CFR of 2.3%. Other papers report CFR values around 5% [7,16,17]. Once more, the age pyramid is of paramount importance in mortality (Table 4), as is access to ICU, particularly in aged people.
In the stepwise process, it must be considered that age, gender, and comorbidity, particularly cardiovascular, play important roles in the final CFR result. Oke and co-workers reported data from the Italian Health Ministry scientific adviser (Professor Walter Ricciardi) indicating that 88% of Italian death certificates related to COVID-19 included at least one pre-morbidity and frequently two or three [18]. Consequently, the evaluation of age, gender, and comorbidity profile becomes crucial when comparing different data series. The Italian series dated March 26, including 73,780 cases [19], provide CFR information closer to our Mediterranean society (Table 5):
As to better include age band in future estimations, the report from The Centre for Evidence-Based Medicine [18] reproduces the whole series of the Italian report [19] describing the statistics as «a grouped-binomial logistic regression with log-link function with main effects for age-band and sex (no two-way interaction terms). Deviance statistic is 30.9 on 6 degrees of freedom» (para.16) providing a table of risk ratio (Table 6) taking age band 60–69 as a reference.

4. Discussion

During the COVID-19 outbreak, data reported by the authorities has proved to be inconsistent. We have estimated an under-reporting of the number of deaths of 29%, close to what has been found in the UK with inconsistencies of about 24% [20], and in other countries [21,22]. A median time delay of 13 days from illness onset to death (17 days with right truncation) [6], and the median basic reproduction number (Ro) 4–6, not far from (2–5), found in other SARS viruses such as the Singapore outbreak [5,23,24], have been reported. The casualties we found for the period studied (649 PMP) are also consistent with other reports [9], and the characteristics of the curve, including a period of about 10 weeks, is consistent with an RNA virus pattern [25].
In cases of incomplete information, such as in the COVID-19 outbreak, death rate may provide the more reliable information to begin with, but one of the important points when comparing different fatality ratios is to analyze data adjusted to the corresponding age band. In this regard, we propose to include an additional reference index, the age-adjusted case fatality ratio (aaCFR), based on risk ratio, setting an age band (e.g., 60–69) as reference. Taking the age bands into consideration, the estimation of the evolution of casualties may be more precise. Once the number of casualties has been determined with the mathematical model, the meta-analysis using available data from the literature allows the estimation of CFR, IFR, and seroprevalence.
Based on the above-mentioned reports, the CFR in Spain could be 4–7%, which is half to one third of the 15% CFR reported by WHO for SARS [26] (p.10), but according to the report, «global case-fatality ratio of 11% was recorded at the end of the outbreak». Consequently, it may be that both SARS outbreaks are not that different in fatality rate.
In the study of the Diamond Princess cruise ship, the IFR result of about half the CFR [10,11], (CFR 2.3% (CI 95% 0.75–5.3%), and IFR 1.2% (CI 95% 0.38–2.7%)), data from Wuhan, and other reports, including WHO [27,28], allow us to forecast that IFR is about half of CFR. This ratio is supported by a computation of data, using Bayesian Markov-chain Monte Carlo methodology, in an age-stratified CFR and IFR model, which resulted in an (adjusted) IFR/CFR ratio close to 0.5 (0.478) [29]. The IFR could be analyzed by predicting attack-rate for age groups [4], but as the main interest here was to move from mortality data to overall population affected in order to evaluate the gross number of possible infected patients, this age-band analysis is not essential.
It is difficult to make a comparison with influenza A (H1N1), as a review of 77 CFR estimations from 50 studies showed a substantial heterogeneity in ranking, from less than 1 to more than 10,000 deaths per 100,000 cases or infections [30]. The official report of the Spanish Surveillance flu group computed a CFR of 0.43 deaths per 1000 cases for the 2009 (H1N1) pandemic [31].
The number of infected patients, assuming most of them will develop herd immunity (natural immunity)—something far from being proved—could be a gross indirect index of the extension and severity of future outbreaks.
Approximating mostly to the Italian report, a crude CFR of 10% over the 30,568 casualties estimated by mathematical modeling at the end of the outbreak will represent a crude IFR value of 5% at most, or around 0.6 million infected patients including both symptomatic and asymptomatic cases (1.3% of a Spanish population of 47.1 million), a figure about half the lowest range of the prediction by the Imperial College Report [4] for Spain. If the WHO lower estimation of CFR/IFR is considered (1/3), then the value is 0.4 million infected patients (0.87% tests will result positive).
Another indirect and approximate estimate of the highest percentage of seropositive cases can be obtained based on the number of hospitalized patients in Spain, whose percentage in relation to the total number of declared cases is 55% [32]. Let us assume that this percentage is not the result of a health policy and protocols in Spain (or Italy) different from those of China and other countries, but rather the result of incomplete information on cases. Let us assume that the number of hospitalized patients corresponds in reality to only 15–20% of the infected population, in line with the WHO report [28]. Over an estimation of about 225,000 cases estimated at the end of this outbreak (in the terms of less than 30 new deaths per day as mentioned), this will represent approximately 123,000 hospitalized patients; if this value is only 15–20% of total symptomatic cases, it will mean 819–615 thousand cases. Taking the highest value and the highest ratio of symptomatic versus asymptomatic cases (0.5), that extreme limit would give an estimation of about 1.2 million patients who have come into contact with the virus (either with or without symptoms), and assuming all their tests will be positive, this upper value estimation represents only 2.6% of the population.
Consequently, serologic analysis is expected to show immunity about 0.87–1.3% of the population, a value close to other preliminary studies, such as the one from Stanford University in Santa Clara with 3324 cases and a result of 1.5% (exact binomial 95% CI 1.1–2.0%). Their test performance specificity was 99.5% (95% CI 99.2–99.7%) and sensitivity s 82.8% (95% CI 76.0–88.4%). The unweighted prevalence adjusted for test performance characteristics was 1.2% (95% CI 0.7–1.8%). After weighting for the population demographics of Santa Clara County, the prevalence was 2.8% (95% CI 1.3–4.7%), using bootstrap to estimate confidence bounds [33]. In our case, less than 3%.
A value about 2–3% positive results in testing the general population has been suggested as a realistic result by the General Director of WHO, Dr Tedros Adhanom Ghebreyesus [27], and this seems to be congruent with the first results obtained in the Netherlands study. Higher percentages (14%), such as reported in the German study with a limited sample of 500 subjects in Heinsberg [9], have been criticized as possible false caveats. The low rate of casualties (0.37%) reported in that study is also to be noted, far from the overall reported fatality rate in Germany of 2% [34]. There is also a surprising ratio of infected people (2%) versus those with antibodies (14%), discrepant with other reports (as commented above) suggesting asymptomatic patients to be about 1/3 or 1/2 of the number of cases with clinical symptomatology.

5. Conclusions

The COVID-19 pandemic outpaced capabilities to control and record cases and deaths, and reporting of data has been inconsistent in many official records [1,2,35]. Clinical assumptions have been made without serological confirmation and many casualties in nursing or home care institutions have not been recorded as related to COVID-19. For these reasons, our study is based on mortality data compared to the period of the previous year [3] and combined with an epidemiological mathematical model which allows estimating the number of cases and deaths at the end of this (first) outbreak. Then, using this result and the relative ratios from other studies, a meta-analysis has been developed to estimate seroprevalence. The heterogeneity of the data in different reports may be, in part, related to the different incidence of the disease according to age bands [13], and an adjustment for relative risk for future reports is proposed in the future. As an important point in health planning is to estimate whether herd immunity may represent a barrier to COVID-19 in Spain, our results based on a meta-analysis contrasted with a mathematical model prediction of approximately 1% of positive tests in the general population (data that in the future may be contrasted with the results of the ongoing study carried out by the Instituto Carlos III). Even considering an extreme value of 2.6%, obtainable by another indirect method, our weighted range (1% to 3%), in line with that reported by the WHO, is very far from the 65–70% necessary for herd immunity to be an effective barrier for COVID-19. As stated by Australian epidemiologist Gideon Meyerowitz–Katz, without vaccination, there is no place for herd immunity as a barrier for COVID-19 [36].

Author Contributions

J.E.M.-F. and D.V.-C.; methodology, N.G.-O.; software, J.E.M.-F.; validation, J.E.M.-F. and D.V.-C.; formal analysis, N.G.-O.; investigation, J.E.M.-F.; resources, J.E.M.-F.; data curation, N.G.-O.; writing—original draft preparation, S.G.-d.J.; writing—review and editing, S.G.-d.J.; visualization, J.E.M.-F.; supervision, D.V.-C.; project administration, D.V.-C; funding acquisition, D.V.-C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We would like to thank John Wright for help with English editing. To Javier Marco-Franco MD, deceased by Covid-19 on May 14, 2020, IN MEMORIAM and to the rest of the health professionals who also fell by the wayside in the fight against the pandemic.

Conflicts of Interest

The authors declare no conflict of interest. The authors declare that they have no affiliations with or involvement in any organization or entity with any financial interest in the subject matter or materials discussed in this manuscript.

References

  1. Baker, T. Mathematical modeling draws more accurate picture of coronavirus cases 2020. Available online: https://www.eurekalert.org/pub_releases/2020-04/mcog-mmd041020.php (accessed on 20 June 2020).
  2. Red Nacional de Vigilancia Epidemiológica-[Spanish epidemiological surveillance network]. Informe sobre la situación de COVID-19 en España-[Report on the situation of COVID-19 in Spain]. 2020. Available online: https://www.isciii.es/QueHacemos/Servicios/VigilanciaSaludPublicaRENAVE/EnfermedadesTransmisibles/Documents/INFORMES/Informes%20COVID-19/Informe%20nº%2022.%20Situación%20de%20COVID-19%20en%20España%20a%2013%20de%20abril%20de%202020.pdf (accessed on 20 June 2020).
  3. Ministerio de Justicia [Spanish Ministry of Justice]. MoMo-Vigilancia de los excesos de mortalidad por todas las causas [Surveillance of excess mortality from all causes (MoMo)]; Ministerio de Justicia [Spanish Ministry of Justice]: Madrid, Spain, 2020; Available online: https://www.isciii.es/QueHacemos/Servicios/VigilanciaSaludPublicaRENAVE/EnfermedadesTransmisibles/MoMo/Documents/informesMoMo2020/MoMo_Situacion%20a%2011%20de%20mayo_CNE.pdf (accessed on 20 June 2020).
  4. Flaxman, S.; Mishra, S.; Gandy, A.; Unwin, J.T.; Coupland, H.; Mellan, T.A.; Zhu, H.; Berah, T.; Eaton, J.W.; Guzman, P.N.P.; et al. Estimating the number of infections and the impact of non-pharmaceutical interventions on COVID-19 in 11 European countries. Imp. Coll. Lond. 2020, 1–35. [Google Scholar] [CrossRef]
  5. Jung, S.; Akhmetzhanov, A.R.; Hayashi, K.; Linton, N.M.; Yang, Y.; Yuan, B.; Kobayashi, T.; Kinoshita, R.; Nishiura, H. Real-time estimation of the risk of death from novel coronavirus (covid-19) infection: Inference using exported cases. J. Clin. Med. 2020, 9, 523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Linton, N.M.; Kobayashi, T.; Yang, Y.; Hayashi, K.; Akhmetzhanov, A.R.; Jung, S.; Yuan, B.; Kinoshita, R.; Nishiura, H. Incubation period and other epidemiological characteristics of 2019 Novel Coronavirus infections with right truncation: A Statistical analysis of publicly available case data. J. Clin. Med. 2020, 9, 538. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Salje, H.; Kiem, C.T.; Lefrancq, N.; Courtejoie, N.; Paireau, J.; Andronico, A.; Hoze, N.; Richet, J.; Dubost, C.; Salje, H.; et al. Estimating the Burden of SARS-CoV-2 in France to Cite this Version: HAL Id: Pasteur-02548181; Institute Pasteur: Paris, France, 2020; Available online: https://hal-pasteur.archives-ouvertes.fr/pasteur-02548181/document (accessed on 20 June 2020).
  8. Instituto Nacional de Estadística (INE) [National Institute of Statistics of Spain]. Defunciones según la Causa de Muerte año 2018. Notas de prensa. Defunciones por causa muerte 2018 [Death by cause death, 2018] 2019, 2018, 1–8. Available online: https://www.ine.es/prensa/edcm_2018.pdf (accessed on 20 June 2020).
  9. Streeck, H.; Hartmann, G.; Exner, M.; Schmid, M. Vorläufiges Ergebnis und Schlussfolgerungen der COVID-19 Case-Cluster-Study (Gemeinde Gangelt) [Preliminary results and conclusions of the COVID-19 case group study (Gangelt municipality); University Hospital Bonn: Bonn, Germany, 2020; Available online: https://www.land.nrw/sites/default/files/asset/document/zwischenergebnis_covid19_case_study_gangelt.pdf (accessed on 20 June 2020).
  10. Rajgor, D.D.; Lee, M.H.; Archuleta, S.; Bagdasarian, N.; Quek, S.C. The many estimates of the COVID-19 case fatality rate. Lancet Infect. Dis. 2020, 3099, 30244. [Google Scholar] [CrossRef] [Green Version]
  11. Mizumoto, K.; Kagaya, K.; Zarebski, A.; Chowell, G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Eurosurveillance 2020, 25, 2000180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wu, J.T.; Leung, K.; Bushman, M.; Kishore, N.; Niehus, R.; de Salazar, P.M.; Cowling, B.J.; Lipsitch, M.; Leung, G.M. Estimating clinical severity of COVID-19 from the transmission dynamics in Wuhan, China. Nat. Med. 2020. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Riou, J.; Hauser, A.; Counotte, M.J.; Althaus, C.L. Adjusted age-specific case fatality ratio during the COVID-19 epidemic in Hubei, China, January and February 2020. medRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  14. Russell, T.W.; Hellewell, J.; Jarvis, C.I.; van Zandvoort, K.; Abbott, S.; Ratnayake, R.; Flasche, S.; Eggo, R.M.; Edmunds, W.J.; Kucharski, A.J. Estimating the infection and case fatality ratio for coronavirus disease (COVID-19) using age-adjusted data from the outbreak on the Diamond Princess cruise ship, February 2020. Eurosurveillance 2020, 25, 6–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Wu, Z.; McGoogan, J.M. Characteristics of and important lessons from the Coronavirus Disease 2019 (COVID-19) Outbreak in China: Summary of a report of 72314 Cases from the Chinese Center for Disease Control and Prevention. JAMA J. Am. Med. Assoc. 2020, 323. [Google Scholar] [CrossRef] [PubMed]
  16. Chen, N.; Zhou, M.; Dong, X.; Qu, J.; Gong, F.; Han, Y.; Qiu, Y.; Wang, J.; Liu, Y.; Wei, Y.; et al. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study. Lancet 2020, 395, 507–513. [Google Scholar] [CrossRef] [Green Version]
  17. Wang, X.; Ma, Z.; Ning, Y.; Chen, C.; Chen, R.; Chen, Q.; Zhang, H.; Li, C.; He, Y.; Wang, T.; et al. Estimating the case fatality ratio of the COVID-19 epidemic in China. medRxiv 2020. [Google Scholar] [CrossRef]
  18. Oke, J.; Heneghan, C. Global Covid-19 Case Fatality Rates; Centre for Evidence-Based Medicine, Oxford University: Oxford, UK, 2020; Available online: https://www.cebm.net/covid-19/global-covid-19-case-fatality-rates/ (accessed on 20 June 2020).
  19. Italian Task Force COVID-19. Epidemia COVID-19 [Covid-19 epidemic]. 2020. Available online: https://www.epicentro.iss.it/coronavirus/bollettino/Bollettino-sorveglianza-integrata-COVID-19_26-marzo%202020.pdf (accessed on 20 June 2020).
  20. Weaver, M. UK coronavirus death toll reaches 1,789 amid data reporting concerns |World news| The Guardian. Available online: https://www.theguardian.com/world/2020/mar/31/uk-coronavirus-death-toll-reaches-1789-amid-data-reporting-concerns (accessed on 7 May 2020).
  21. London Business School. How Many People Really Die From Covid-19? Lessons from Italy. Forbes, 6 April 2020. [Google Scholar]
  22. Wu, J.; McCann, A.; Peltier, E. Coronavirus-missing-deaths. Available online: https://www.nytimes.com/interactive/2020/04/21/world/coronavirus-missing-deaths.html (accessed on 10 May 2020).
  23. Lipsitch, M.; Cohen, T.; Cooper, B.; Robins, J.M.; Ma, S.; James, L.; Gopalakrishna, G.; Chew, S.K.; Tan, C.C.; Samore, M.H.; et al. Transmission dynamics and control of severe acute respiratory syndrome. Science 2003, 300, 1966–1970. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Sanche, S.; Lin, Y.T.; Xu, C.; Romero-Severson, E.; Hengartner, N.; Ke, R. High Contagiousness and rapid spread of severe acute respiratory syndrome coronavirus 2. Emerg. Infect. Dis. 2020, 26. [Google Scholar] [CrossRef] [PubMed]
  25. Holmes, E.C. The Evolution and Emergence of RNA Viruses; Oxford University Press: Oxford, UK, 2010; Volume 16, ISBN 9780199211135. [Google Scholar]
  26. World Health Organization. Consensus Document on The Epidemiology of Severe Acute Respiratory Syndrome (SARS); WHO Publishing: Geneva, Switzerland, 2003; pp. 1–47. Available online: https://www.who.int/csr/sars/WHOconsensus.pdf?ua=1 (accessed on 7 May 2020).
  27. Boseley, S. WHO warns that few have developed antibodies to Covid-19. The Guardian. 2020. Available online: https://www.theguardian.com/society/2020/apr/20/studies-suggest-very-few-have-had-covid-19-without-symptoms (accessed on 20 June 2020).
  28. World Health Organization (WHO) Tech Guid #2. Health Systems Respond to COVID-19 Technical Guidance #2 Creating surge capacity for acute and intensive care Recommendations for the WHO European Region (6 April 2020). JAMA 2020, 323, 1239–1242. [Google Scholar] [CrossRef]
  29. Verity, R.; Okell, L.C.; Dorigatti, I.; Winskill, P.; Whittaker, C.; Imai, N.; Cuomo-Dannenburg, G.; Thompson, H.; Walker, P.; Fu, H.; et al. Estimates of the severity of COVID-19 disease. medRxiv 2020. [Google Scholar] [CrossRef]
  30. Wong, J.Y.; Kelly, H.; Ip, D.K.M.; Wu, J.T.; Leung, G.M.; Cowling, B.J. Case fatality risk of influenza a (H1N1pdm09): A systematic review. Epidemiology 2013, 24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Cámara-Larrauri, A.; Jiménez-Jorge, S.; Simón-Méndez, L.; de Mateo-Ontañón, S. Vigilancia de la pandemia de la gripe 2009 (H1N1) en España [Surveillance of 2009 flu (H1N1) pandemic in Spain]. Rev. Esp. Salud Publica 2010, 24, 569–588. [Google Scholar]
  32. Centro de coordinación de alertas y emergencias sanitarias; R98-[Spanish centre of health alerts and emergencies. R98]. Actualización no 98. Enfermedad por el coronavirus (COVID-19)- [Update No 98. Coronavirus (Covid-19)]; Centro de coordinación de alertas y emergencias sanitarias: Madrid, Spain, 2020; Available online: https://www.mscbs.gob.es/profesionales/saludPublica/ccayes/alertasActual/nCov-China/documentos/Actualizacion_98_COVID-19.pdf (accessed on 7 May 2020).
  33. Bendavid, E.; Mulaney, B.; Sood, N.; Shah, S.; Ling, E.; Bromley-Dulfano, R.; Lai, C.; Weissberg, Z.; Saavedra, R.; Tedrow, J.; et al. COVID-19 Antibody Seroprevalence in Santa Clara County, California. medRxiv 2020. [Google Scholar] [CrossRef]
  34. Regalado, A. Blood tests show 14% of people are now immune to covid-19 in one town in Germany. MIT Technology Review, 9 April 2020. [Google Scholar]
  35. Wu, Z.; McGoogan, J.M. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China. JAMA 2020, 323, 1239–1242. [Google Scholar] [CrossRef] [PubMed]
  36. Meyerowitz-Katz, G. Here’s Why Herd Immunity Won’t Save Us From The COVID-19 Pandemic. Science Alert, 30 March 2020. [Google Scholar]
Figure 1. COVID-19: evolution of adjusted (standardized) casualties per million population (PMP). Authors’ forecast based on first 45 days of outbreak starting March 9.
Figure 1. COVID-19: evolution of adjusted (standardized) casualties per million population (PMP). Authors’ forecast based on first 45 days of outbreak starting March 9.
Sustainability 12 05228 g001
Table 1. Spain: Under-estimation of causalities (March 17 to April 18) in official report.
Table 1. Spain: Under-estimation of causalities (March 17 to April 18) in official report.
Total CasesExpected CasesUnexpectedReported COVID-19Report Under-Estimation
63,67637,76925,90720,043586429%
Source: Authors’ computation with information from MoMo, Ministry of Justice [3,8].
Table 2. Spain: estimation of regional overall outbreak casualties per million population (PMP) over a period of 69 days.
Table 2. Spain: estimation of regional overall outbreak casualties per million population (PMP) over a period of 69 days.
Regions (CCAA)Total CasualtiesCFR PMPRegions (CCAA)Total CasualtiesCFR PMP
Andalucía1501178.7Valencia1563312.7
Aragón939708.6Extremadura556523.4
Asturias328321.8Galicia519192.4
Baleares231192.5Madrid10,4421562.0
Canarias17377.8Melilla331.3
Cantabria241414.5Murcia167112.0
Castilla La Mancha30711506.7Navarra568870.7
Castilla y León2196913.9Pais Vasco1621743.1
Cataluña6019790.0La Rioja4151.320.3
Ceuta562.4
Spain (Overall)30,568649.0
Source: Authors’ computation with information from INE [8] and MoMo, Ministry of Justice [3].
Table 3. Epidemiological data from Diamond Princess cruise.
Table 3. Epidemiological data from Diamond Princess cruise.
Sample SizePCR PositiveOver InfectedOver PCR +Over Sample
True AsympSymp.IFRCFRCasualties
21%18%82%1.11%1.35%0.23%
(n)3063630113.3516.7 7
Source: Authors’ computation with data from [10,11]. IFR: infectious fatality rate. CFR: case fatality rate.
Table 4. Epidemiological data from Chinese Centre for Disease Control and Prevention.
Table 4. Epidemiological data from Chinese Centre for Disease Control and Prevention.
AgeCasualtiesCases% CFR (95% CI)
≤9 years04160
10 to 19 years15490.18 [0.03–1.02]
20 to 49 years6319,7900.32 [0.25–0.41]
50 to 59 years13010,0081.3 [1.1–1.5]
60 to 69 years30985833.6 [3.2–4.0]
70 to 79 years31239188.0 [7.2–8.9]
≥80 years208140814.8 [13.0–16.7]
Overall102344,4152.30
Source: Authors’ computation with data from Wu & McGoogan [15].
Table 5. Epidemiological data from Italy on 26 March 2020.
Table 5. Epidemiological data from Italy on 26 March 2020.
Total Cases
Age Band# Cases% Age Band# Death% Death/ Age Band% CFR
0–94240.6-0.00.00
10–195100.7-0.00.00
20–2927133.7-0.00.00
30–3049596.8170.20.34
40-49916712.6671.00.73
50–5914,33519.72433.61.70
60–6913,14918.176111.25.79
70–7914,09019.4240335.317.05
80–8910,92915.0270239.724.72
≥9025173.56088.924.16
Total72,793100.06801100.09.34
Males
Age Band# Cases% Gender# Death% Death/ Gender% CFR
0–924457.5-0.00.00
10–1926151.2-0.00.00
20–29120344.3-0.00.00
30–30246549.71482.40.57
40–49459750.14973.11.07
50–59799855.819078.22.38
60–69875566.660679.66.92
70–79930966.1184676.81.83
80–89619556.7180866.929.18
≥9087734.827344.931.13
Total41,90457.6478670.411.42
Females
Age Band# Cases% Gender# Death% Death/ Gender% CFR
0–918042.5-0.00.00
10–1924948.8-0.00.00
20–29151055.7-0.00.00
30–30249450.3317.60.12
40-49457049.91826.90.39
50–59633744.25221.40.82
60-69439433.415420.23.50
70–79478133.955523.111.61
80-89473443.389433.118.88
≥90164065.233454.920.37
Total3088942.4201029.66.51
Source: Authors’ recalculation of data from records with complete information (987 cases excluded) published by the Higher Institute of Health (Istituto Superiore di Sanità, ISS). CFR case fatality ratio [19].
Table 6. COVID-19 risk ratio for gender and age band.
Table 6. COVID-19 risk ratio for gender and age band.
CategoryRisk Ratio95% CI
Age 30–390.06[0.038–0.10]
Age 40–490.14[0.11–0.17]
Age 50–590.31[0.27–0.35]
Age 60–69 (reference)1-
Age 70–792.95[2.7–3.2]
Age 80–894.47[4.1–4.8]
Age 90+4.83[4.4–5.3]
Female1-
Male1.66[1.58–1.74]
Source: Centre for Evidence-Based Medicine research with data from Italian Task Force COVID-19 [19].

Share and Cite

MDPI and ACS Style

Marco-Franco, J.E.; Guadalajara-Olmeda, N.; González-de Julián, S.; Vivas-Consuelo, D. COVID-19 Healthcare Planning: Predicting Mortality and the Role of the Herd Immunity Barrier in the General Population. Sustainability 2020, 12, 5228. https://doi.org/10.3390/su12135228

AMA Style

Marco-Franco JE, Guadalajara-Olmeda N, González-de Julián S, Vivas-Consuelo D. COVID-19 Healthcare Planning: Predicting Mortality and the Role of the Herd Immunity Barrier in the General Population. Sustainability. 2020; 12(13):5228. https://doi.org/10.3390/su12135228

Chicago/Turabian Style

Marco-Franco, Julio Emilio, Natividad Guadalajara-Olmeda, Silvia González-de Julián, and David Vivas-Consuelo. 2020. "COVID-19 Healthcare Planning: Predicting Mortality and the Role of the Herd Immunity Barrier in the General Population" Sustainability 12, no. 13: 5228. https://doi.org/10.3390/su12135228

APA Style

Marco-Franco, J. E., Guadalajara-Olmeda, N., González-de Julián, S., & Vivas-Consuelo, D. (2020). COVID-19 Healthcare Planning: Predicting Mortality and the Role of the Herd Immunity Barrier in the General Population. Sustainability, 12(13), 5228. https://doi.org/10.3390/su12135228

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop