 Research article
 Open Access
 Published:
Agespecific rate of severe and critical SARSCoV2 infections estimated with multicountry seroprevalence studies
BMC Infectious Diseases volumeÂ 22, ArticleÂ number:Â 311 (2022)
Abstract
Background
Knowing the agespecific rates at which individuals infected with SARSCoV2 develop severe and critical disease is essential for designing public policy, for infectious disease modeling, and for individual risk evaluation.
Methods
In this study, we present the first estimates of these rates using multicountry serology studies, and public data on hospital admissions and mortality from early to mid2020. We combine these under a Bayesian framework that accounts for the high heterogeneity between data sources and their respective uncertainties. We also validate our results using an indirect method based on infection fatality rates and hospital mortality data.
Results
Our results show that the risk of severe and critical disease increases exponentially with age, but much less steeply than the risk of fatal illness. We also show that our results are consistent across several robustness checks.
Conclusion
A complete evaluation of the risks of SARSCoV2 for health must take nonfatal disease outcomes into account, particularly in young populations where they can be 2 orders of magnitude more frequent than deaths.
Background
The SARSCoV2 pandemic had impacts of historic proportion in both public health and society. Remarkably, there is considerable uncertainty regarding the full spectrum of health effects of SARSCoV2 infection. On the one hand, some of the effects of SARSCoV2 are relatively well understood, such as the infection fatality rate (IFR) and its dependence on age. Three different metaanalyses now exist that estimate the agestratified IFR of SARSCoV2 using multicountry seroprevalence studies [1,2,3]. These studies document an exponential increase of the IFR with age and show considerable agreement on their estimated IFRs by agestratum. On the other hand, the rate of less extreme infection outcomes, and their dependence on age, remains uncertain despite being similarly important for public health. Examples of this are the rate of severe infections (Infectionsevere rate, ISR), which we define as infections resulting in hospitalization or outofhospital death, and the rate of critical infections (Infectioncritical rate, ICR), which we define as infections resulting in admission to intensive care unit (ICU) or outofICU death.
But despite their relevance to analyzing the development of the pandemic and for future planning, estimates of ISR and ICR for these ages using multicountry data are still missing from the literature (see estimates of the infectionhospitalization rate, for France [4, 5], for Denmark [6], for Indiana, USA [7], for Connecticut, USA [8], for Qatar [9], and a modelbased analysis with early nonserological pandemic data [10]). To fill this gap, we present a metaanalysis of the agestratified rates of severe and critical disease of SARSCoV2 across several locations, combining seroprevalence studies from early to mid 2020 with public data on the numbers of agestratified hospitalizations, ICU admissions, and deaths.
Results
We analyzed locations with seroprevalence studies that were either listed in the metaanalysis of Levin et al. [1], to which we refer for further details, or the studies providing their own agestratified rates of infectionhospitalization rate. We included in the analysis 15 locations with serosurveys (11 using representative samples and 4 using convenience samples), and 2 locations with comprehensive testing and contact tracing (see methods section M1 for a description of this classification of locations). Together, these locations represent 5% of the worldâ€™s population. Furthermore, to account for outof hospital and outof ICU deaths, which are common among the elderly, we computed the number of severe cases as the number of hospitalizations plus outofhospital deaths, and the number of critical cases as the number of ICU admissions plus outofICU deaths.
The estimated probability of severe, critical, and fatal disease outcomes (ISR, ICR, and IFR, respectively) are shown for each age and location as colored points in Fig.Â 1, using a logtransformed vertical axis. As expected from the reports in previous analyses of IFR [1,2,3] and of the infectionhospitalization ratio [4], the three outcome ratios show an approximately exponential increase in risk with respect to age, which becomes a homogeneous linear effect on the logscale. Thus, for each outcome rate (ISR, ICR, and IFR) we fitted Bayesian logistic regression models with a linear age effect on the logit scale (this effect becomes nonlinear in the risk scale). We used logistic regression because it is a commonly used model for disease outcomes that approximate the loglinear rateage patterns observed in our study. The logistic regression model had an intercept and ageslope shared across locations, plus locationspecific random effects on both the slope and the intercept of the regression to account for the heterogeneity between locations. We also accounted for the uncertainty of the seroprevalence estimates (through the specification of the prior distribution), and the sampling (binomial) variability of the observed outcomes.
The logistic regression models fitted the data well for the three outcome rates (Fig.Â 1, although see the discussion of section S1 in Additional file 1 about a possible deviation from the trend for the youngest ages), and the patterns observed were similar across locations. Importantly, the slope of IFR with respect to age (0.133, 95% credibility interval: [0.123â€“0.143]) was higher than the slope of ICR (0.099, [0.089â€“0.108]) and of ISR (0.076, [0.067â€“0.083]), indicating that the risks of severe and critical disease are more evenly distributed across ages than the risk of death. The models also had different intercepts, of âˆ’Â 12.9 [âˆ’Â 13.6, âˆ’Â 12.3] for IFR, âˆ’Â 9.9 [âˆ’Â 10.3, âˆ’Â 9.4] for ICR and âˆ’Â 7.3 [âˆ’Â 7.7, âˆ’Â 6.8] for ISR, reflecting the difference of 1 order of magnitude between ICR and IFR for the youngest ages, and of 2 orders of magnitude between ISR and IFR. As an example, according to our estimates, people in the 20â€“25Â years old range are on average 779 [467â€“1223] times less likely to die from COVID19 than 70â€“75Â year old people, 133 [83â€“202] times less likely to develop critical disease, and 38 [27â€“51] times less likely to develop severe disease.
Predicted risk levels by age (and the corresponding 95% credibility intervals) are shown in Table 1 (see Table S1 for finer age stratification, and Table S2 for the estimated model parameters in Additional file 1). Also, we verified that our estimates are robust to the correction for outofhospital and outofICU deaths (Additional file 1: Figs. S1, S2), the ages used to fit the model (Additional file 1: Fig. S3), the method of estimating SARSCoV2 infections (Additional file 1: Fig. S4), the date of outcome data collection (Additional file 1: Figs. S5, S6), and the delay between the epidemic wave and the seroprevalence study (Additional file 1: Fig. S7, to control for seroreversion). More details about these comprehensive robustness analyses are provided in the Additional file 1.
Next, we validated our estimates by estimating the ISR and ICR of SARSCoV2 indirectly using a novel ratioofratios approach. We start from the agespecific IFR reported in the three different metaanalyses [1,2,3], which were not used in the analysis of Fig.Â 1. Because the IFR is the expected ratio between deaths and infections, we can estimate the ISR as the ratio IFR/SFR, where SFR is the ratio between deaths and severe infections (severe fatality rate). We approximated the agespecific SFR by fitting a Bayesian logistic regression model to published data of COVID19 hospital mortality (Fig.Â 2A, data sources listed in Table 4), which is the ratio between inhospital deaths and hospitalizations. The approximation of SFR by hospital mortality assumes that all deaths occur in hospitals, which is expected to hold well for all but the oldest age bins (see Additional file 1: Section S2, Fig. S2). Then, we estimated the agespecific ISR by taking the ratio between the IFRs and the SFRs. We applied the same procedure to estimate ICR, using the ICU mortality of COVID19 patients. The values of the parameters obtained by fitting the model to hospital and ICU mortality are shown in Additional file 1: Table S3, and the agespecific estimates are shown in Additional file 1: Table S4.
We note that we used the same hospital and ICU mortality data for this analysis and for correcting for outofhospital and outofICU deaths in Fig.Â 1. Although this means that the two analyses share some data in common, the results of the regression of Fig.Â 1 are similar when performed on the uncorrected data (Additional file 1: Figs. S1, S2), supporting the use of this validation method.
The estimates from the indirect method are shown by the colored points in Fig.Â 2B. To aid comparison, we show in black the lines obtained from the fit to serology data of Fig.Â 1. Firstly, we see that our IFR estimates and the IFR estimates from the three metaanalyses are very similar (Fig.Â 2B, right). Second, we note that the estimates of ISR and ICR obtained from seroprevalence and disease outcome data with the direct method (Fig.Â 1) are in close agreement with the estimates obtained with the indirect method (with the largest differences being with Brazeau et al. [3] for the younger ages; however, this study also reports estimates different than those reported in the other two studies).
Conclusions
In conclusion, we present estimates of the rates of severe and critical SARSCoV2 infections during the first half of 2020. These are the first estimates based on multicountry seroprevalence data, which we combine in a rigorous way using Bayesian methods, to account for the uncertainty of each study as well as temporal and geographical heterogeneity. We find that while young and middleaged individuals had low rates of fatal infection, they had much higher rates of severe and critical infection, emphasizing the need to consider these disease outcomes in these populations. The estimates presented here are an important reference of the health impacts of COVID19 during 2020, as well as an important baseline over which to build more updated estimates, by combining them with estimates of the relative change in risk across locations and time.
Methods
The data sets included in this study come from locations where agestratified seroprevalence studies have been performed (see M1), plus locations with agestratified prevalence coming from exhaustive contract tracing (see M2). For each of these locations, we searched for agestratified data on Hospitalizations and ICU admissions (see M3). We used these two sources of data to estimate the agestratified rates of severe (ISR) and critical (ICR) SARSCoV2 infections for each of the locations. We used this data to fit Bayesian randomeffects logistic regression models for each of the outcomes (see M4). We also searched for studies reporting agestratified mortality for COVID19 patients admitted to the hospital or to the ICU (see M5). We used this data to fit a Bayesian randomeffects logistic regression model to obtain the agespecific hospital and ICU mortality for COVID19. This regression was then combined with estimates of agespecific IFR extracted from the literature, to estimate the ISR and ICR through a ratioofratios method (see M6). The regression of hospital and ICU mortality was also used to correct the hospital and ICU data described in M3 for outofhospital and outofICU deaths (see M7). All data and code are available online (see M8).
(M1) Data from seroprevalence studies
We used a curated list of seroprevalence studies released prior to 18 September 2020 that is presented in Levin et al. [1]â€”a systematic review and a metaanalysis. The list is restricted to developed countries; we refer the reader to Levin et al. [1] for an exhaustive list of other existing studies, and the criteria used for excluding seroprevalence studies from their final analysis. The locations used by Levin et al. [1] can be divided into three groups: those with representative seroprevalence studies, those with convenience seroprevalence studies, and those with comprehensive testing and tracing. Representative seroprevalence studies refer to those in which the population included in the serosurvey aims to be a representative sample of the population. Convenience seroprevalence studies are those that perform the serosurvey over samples that are conveniently available but not necessarily representative of the population, like blood donor samples. Locations with comprehensive testing and tracing are defined by Levin et al. [1] as locations that up to the date of interest, had over 300 tests performed for each detected case (in these locations, we corrected for underascertainment following Levin et al. [1]).
We then searched for agestratified hospitalization and ICU data to match the representative seroprevalence studies listed in Additional file 1: Appendix Tables I1 and I3, and the convenience seroprevalence studies listed in Appendix Table I2 from Levin et al. [1]. From the 11 representative seroprevalence studies included in those lists, we were able to find agestratified hospitalization or ICU data for 7 locations (England; France; Ireland; Netherlands; Spain; Atlanta, USA; Geneva, Switzerland) and we failed to find such data for the 4 remaining locations (Italy; Portugal; Indiana, USA; Salt Lake City, USA). From the 4 convenience seroprevalence studies listed, we were able to find hospitalization or ICU data for all three locations (Ontario, Canada; Sweden; Belgium; New York, USA).
In addition to the seroprevalence studies used in Levin et al. [1], we included the seroprevalence study carried out in Iceland up to April 4 2020, which reports the results of a representative sample of the population [46], and the representative seroprevalence studies from Indiana, Connecticut and Denmark, which were used to estimate the ISRs in the existing literature [6,7,8].
Furthermore, in three cases we also changed the use of some seroprevalence studies with respect to Levin et al. to match them to the available hospitalization and ICU data. The first case is the seroprevalence study from France, which offers data by region. We were only able to find the agestratified hospitalization data for the region of ÃŽledeFrance; therefore, we only used seroprevalence data from this region. The second case is the New York seroprevalence study, where we could only find hospitalization data for New York City but not for New York State; thus, we only used the seroprevalence for New York City. For Ontario, Canada, we could only find agestratified hospitalization and ICU data up to July 31 2020, and so we used the seroprevalence report for this date, which is different from the report date used by Levin et al. [1]. Table 2 summarizes the final list of seroprevalence studies included in our analysis.
Matching agebins
The agebins reported by each of the studies did not always match the agebins in the corresponding hospitalization and ICU reports. Therefore, in some cases we extrapolated or interpolated the seroprevalence estimates obtained for a given agebin into a different agebin. For example, for New York City, seroprevalence was reported for the 18â€“34Â years old age range, but hospitalization data was reported for the 18â€“44Â year old age range. Therefore, to make use of this hospitalization data, we assumed that the proportion of seropositive individuals in the 18â€“44Â years old range is the same as the proportion for the 18â€“34Â year old age range. All such decisions were contrasted with other available data, and agreed upon by the two authors. Furthermore, these assumptions are all documented in the publicly available analysis code.
Correcting for test characteristics
The positive rate of a test depends on disease prevalence and on the test characteristics. Most of the seroprevalence estimates used were already corrected for test characteristics. For the results that were not corrected for test characteristics, we used the GladenRogan formula (Rogan and Gladen 1978) to adjust the estimates as follows:
(M2) Countries with comprehensive tracing included in the analysis
Following [1], we also included in our analysis two countries (Republic of Korea and New Zealand) with comprehensive tracing programs where the number of infections detected through testing are thought to approximate the total number of infections accurately. As in the original Levin et al. study, we corrected the prevalence estimates for these countries using the agespecific ratio between the number of infections estimated through seroprevalence and the number of positive tests in Iceland (Gudbjartsson et al. 2020).
(M3) Hospitalizations, ICU admissions, and deaths data
We obtained the agestratified hospitalizations, ICU, and death data in relevant government websites of the locations, using google search, and looking for relevant regionwide studies. We selected the data reports that were closest to the end of the serosurvey date. The list of data sources and the end dates for their cumulative outcome numbers are shown in Table 3.
Also, as described above for the serology data, in some cases we interpolated or extrapolated some data for these disease outcomes, or we combined different data sources with incomplete data (e.g., agedistribution of an outcome from one source, with the total count of the outcome from another source) to obtain the data for these outcomes with the appropriate age bins. For example, for Belgium, we were only able to find the unstratified number of cumulative ICU admissions at the desired date of May 8th, 2020, but we were able to find the age distribution for cumulative ICU admissions up to June 14th, 2020. Therefore, we distributed the cumulative ICU admissions of May 8th across age strata, following the distribution from June 14th. As mentioned above, all these decisions were agreed upon by the authors, and they are thoroughly documented in the publicly available analysis code.
(M4) Estimation of outcome probabilities with serology and outcome data
We fitted Bayesian logistic regression models to the serology and outcome data. We describe the model for severe SARSCoV2 outcome; the same model was fitted to severe, critical, and fatal disease outcomes.
Let \(\{{y}_{la},{x}_{la}\}\) represent the number of severe SARSCov2 infections \({(y}_{la})\) experienced among \({x}_{la}\) individuals infected with SARSCoV2, at the location l (lâ€‰=â€‰1, â€¦, L), for the age stratum a (aâ€‰=â€‰1, â€¦, A_{l}) of the l^{th} location. The ISR for this locationstratum is defined as \({\theta }_{R,la}=E[\frac{{y}_{la}}{{x}_{la}}]\).
Bayesian likelihood
The probability of \({y}_{la}\) given \({\theta }_{R,la}\) and the number of infections\({(x}_{la})\), is given by the Binomial likelihood\(p({y}_{la}) \propto {\theta }_{R,la}^{{y}_{la}}(1{\theta }_{R,la}{)}^{{x}_{la}{y}_{la}}\). Assuming conditional independence across locations and strata, and taking\(y=({y}_{11},{y}_{12}, ...,{y}_{L{A}_{L}})\), and \({\theta }_{R}=({\theta }_{R,11},{\theta }_{R,12}, ...,{\theta }_{R,L{A}_{L}})\) we have \(p({y\theta }_{R}) \propto {\prod }_{l=1}^{l=L} {\prod }_{a=1}^{a={A}_{l }}{\theta }_{R,la}^{{y}_{la}}(1{\theta }_{R,la}{)}^{{x}_{la}{y}_{la}}\)
Modeling the number of SARSCov2 cases from seroprevalence data
The selected seroprevalence studies provide agestratified estimates (and SE) of disease prevalence. Rather than assuming that prevalence was known with complete certainty, we used the reported point estimates and SE to specify a Beta prior for prevalence for each location. Specifically, we used the reported prevalence and its SE to estimate (through first and second moments matching) the shape parameters of the Beta distribution used for each location. Then, prevalence was modeled as \({\theta }_{P,la}\sim Beta({\alpha }_{1,la},{\alpha }_{2,la})\), where \({\alpha }_{1,la}\) and \({\alpha }_{2,la}\) are locationagestratum specific shape parameters. Then, the number of cases was defined as \({x}_{la}={N}_{la}\times {\theta }_{P,la}\), where \({N}_{la}\) is the size of the population at location l and age stratum a.
Modeling severity rates using randomeffects logistic regression
Infectionseverity rates were modeled using a logit of the form.
Therefore, \({\theta }_{R,la}=\frac{{e}^{{\eta }_{la}}}{1+{e}^{{\eta }_{la}}}\) where \({\eta }_{la}=[\mu +{u}_{l}] + [\beta +{b}_{l}]\times ag{e}_{la}\).
Above, \(\mu\) and \(\beta\) are the average intercept and slopes across locations, and \({u}_{l}\) and \({b}_{l}\) are locationspecific random effects on the intercept and the slope, respectively, with prior distribution \(p(u,b{\sigma }_{\mu }^{2},{\sigma }_{b}^{2})={\Pi }_{l=1}^{L}N({u}_{l}0,{\sigma }_{u}^{2})N({b}_{l}0,{\sigma }_{b}^{2})\). The shared intercept (\(\mu\)) and regression coefficient (\(\beta\)) were assigned flat priors, and the standard deviations for the random effects, \({\sigma }_{u}\) and \({\sigma }_{b}\) were assigned gamma priors with shape and rate parameters equal to 4.
For each age stratum a at location l, the value of \(ag{e}_{la}\) used for fitting corresponded to the median age of the stratum. For age strata with an open upper bound (e.g. 70â€‰+â€‰age), we used 90Â years as the upper bound of the stratum.
The posterior distribution of the model described above does not have a closed form; therefore, we used Monte Carlo Markov Chain (MCMC) methods to generate samples from the posterior distribution for all the model unknowns \(\{\mu ,\beta ,x,u,b,{\sigma }_{\mu }^{2},{\sigma }_{b}^{2}\}\). We used 4 chains with 2500 iterations each. A script that implements the above model in Stan [64] is available in the online code.
Prediction of outcome rates
We used the samples of the posterior distribution to generate posterior samples for the infection severity rates for specific ages using the inverselogit function: \({\theta }_{R}(age{)}_{s}=\frac{{e}^{{\mu }_{s}+age\times {\beta }_{s}}}{1+{e}^{{\mu }_{s}+age\times {\beta }_{s}}}\), where \(s\) is an index for the sample from the posterior distribution. We then used these samples to estimate the posterior means and posterior credibility regions reported in Figs.Â 1 and 2. We report the severity rates for age intervals by estimating the rate of the mean age of the interval.
The predicted outcome rates obtained from the model fit are shown in Additional file 1: Table S1, and the mean and credible intervals for the main model parameters are shown in Additional file 1: Table S2.
(M5) Hospital and ICU mortality data
Our robustness analysis was based on an indirect estimator (a ratioofratios) of ISR and ICR. To derive this estimator, we used mortality data from hospitalized and ICU SARSCov2 patients. We searched in the literature for reports on agestratified mortality of patients admitted to the hospital or the ICU with a COVID19 diagnosis. We also used the data sources from Table 3 that provided mortality numbers for hospitalized or ICU patients. We identified 8 ICU mortality reports and 8 mortality hospital reports with agestratified data, which were either published studies in the literature or public reports from official organisations. The reports used are listed in Table 4.
(M6) Indirect estimation of ISR and ICR using IFR and hospital mortality data
To validate the estimates obtained with the data and methods described above, we used an alternative source of data and a different estimation method to obtain agespecific ISR and ICR. Specifically, we combined agespecific reports of IFRs from the literature with the hospital and ICU mortality data listed in Table 4 to obtain the ISR and ICR using a ratioofratios method, as explained below.
Let \(IF{R}_{sa}\) be the expected ratio between deaths and infections estimated in a study s (sâ€‰=â€‰1, â€¦, S) for age stratum a (aâ€‰=â€‰1, â€¦, A_{s}) and let \(SF{R}_{a}\) be the expected ratio between deaths and severe COVID19 cases for age stratum a (aâ€‰=â€‰1, â€¦, A_{s}). Then, we have that the estimated ISR for age stratum a estimated from study s is \(IS{R}_{sa}=\frac{IF{R}_{sa}}{SF{R}_{a}}\). Thus, by estimating the values of \(SFR\) for different ages, we can use agespecific IFR values reported in the literature to obtain estimates of agespecific ISR.
To approximate the agespecific SFR, we fitted a Bayesian logistic regression to agestratified hospital death for COVID19 patients. Let \(\{{d}_{la},{h}_{la}\}\) represent the number of deaths \({(d}_{la})\) among \({h}_{la}\) individuals hospitalized with COVID19 for the age stratum a (aâ€‰=â€‰1, â€¦, A_{l}) in location l. The hospital mortality for this locationstratum is defined as \({\theta }_{HM,la}=E\left[\frac{{d}_{la}}{{h}_{la}}\right]\).
To estimate \({\theta }_{HM}(age)\), we used Bayesian randomeffects logistic regressions, like the one described in section M4, to the hospital death data. The only difference with the procedure in M4 is that, in this case, the denominators \({h}_{la}\) were known, and thus we directly used these fixed \({h}_{la}\) values (unlike the \({x}_{la}\) from M4, for which a distribution over possible values was obtained using seroprevalence estimates).
We use \({\theta }_{HM}(a)\) as our estimate of \(SF{R}_{a}\). These two quantities are equal if we assume that all deaths occur in the hospital (note that our definition of severe case, stated in the main text, is a case that results in either hospital admission or outofhospital death). As discussed in section S3 and shown in Additional file 1: Figs. S1, S2, outofhospital deaths make only a very small fraction of severe cases for all but the oldest agestrata. Also, we find that outofhospital deaths make up a minority of the deaths for all but the oldest ages (analysis not shown).
Then, to account for the uncertainty of the \(IF{R}_{sa}\) estimates in our estimations, we fitted a Beta distribution to the mean and credible interval of each \(IF{R}_{sa}\) through moment matching, to obtain \(IF{R}_{sa}\sim Beta({\alpha }_{1,sa},{\alpha }_{2,sa})\) (for Brazeau et. al. (2020) we only used the point estimates since credible intervals on the mean estimates are not reported).
Finally, we estimated \(IS{R}_{sa}=\frac{IF{R}_{sa}}{SF{R}_{a}}\) by generating samples from the posterior distribution of each \(SF{R}_{a}\) (obtained from the Bayesian logistic regression model) and from the Beta distribution fitted for each \(IF{R}_{sa}\). In total, we generated 50.000 samples of this ratio for each \(IS{R}_{sa}\).
The same procedure was applied to estimate the \(IC{R}_{sa}\), by fitting the model to ICU death data. The estimated hospital and ICU mortality rates obtained from these models are shown in Additional file 1: Table S3, and the parameters obtained from fitting the model are shown in Additional file 1: Table S4.
(M7) Correction for outofhospital and outofICU deaths
Some COVID19 deaths occur outside of the ICU, or outside of the hospital. This happens when the patient prognosis is poor, such as in elderly and frail patients, and it may be accentuated when health systems are operating at high occupancy. This phenomenon is particularly notable in our data for some locations and ages, where the number of reported deaths is larger than the number of reported ICU admissions (in some cases by more than one order of magnitude).
Our definitions of severe and critical COVID19 outcomes include these outofhospital and outofICU deaths, besides hospitalizations and ICU admissions. Therefore, we obtained the number of severe cases by adding to hospitalizations the number of outofhospital deaths. Likewise, we obtained the number of critical cases by adding to the number of ICU patients the number of outofICU deaths. For some locations, we could obtain data on the outofhospital and outofICU deaths, but for other locations this data was absent, and so we estimated it using the death data.
Let \({y}_{la}\) be the cumulative number of hospitalizations for a location l and age stratum a, for which no outofhospital death data is available. Also, let \({{m}_{la}^{tot}}\) be the total number of deaths reported for this location and age stratum. First, we obtained the expected number of inhospitaldeaths, \({m}_{la}^{h}\), by combining the number of hospitalizations with the expected hospital mortality for this age, \({\theta }_{HM}(a)\) (fitted as described in section M6), \({m}_{la}^{h}={{y}_{la}\times {\theta }_{HM}(a)}\). Then, we obtain the expected number of outofhospital \({m}_{la}^{ooh}\) deaths by subtracting from the total number of deaths, \({m}_{la}^{ooh}={m}_{la}^{tot}{m}_{la}^{h}\) (setting \({m}_{la}^{ooh}\) to 0 if the result is negative).
The same procedure is performed for the ICU data to obtain the number of critical cases.
See sections Additional file 1: Supplementary S3 and Figs. S1, S2 for an analysis showing the effect of this correction method on the data, and the robustness of the results to the removal of this correction.
Discussion
In this work, we present the first estimates of ISR and ICR of SARSCoV2 obtained through a metaanalysis of serology studies from early to mid 2020. Our estimates show that, like the IFR, the ISR and ICR increase exponentially with age; however, the rate of increase in the risk of severe and critical disease outcomes with age is smaller than the rate of increase in lethality, which is in agreement with previous studies [5,6,7,8,9,10,11]. However, previous studies show considerable variability, probably due to the uncertainty in serology estimates, differences in local reporting protocols, and geographical variability in the impacts of COVID19. Thus, this analysis presents the most up to date estimation and comparison of these rates, summarizing the best available evidence of several locations with a Bayesian approach. Our simple Bayesian regression analysis that takes into account several sources of uncertainty and this novel ratioofratio methods constitute two complementary methods that may aid future work on estimating these parameters under the changing nature of the pandemic. For example, these methods can be extended so as to estimate the outcome rates for the new variants of SARSCoV2.
Furthermore, we provide extensive validation of our estimates (see Additional file 1). First, we performed several robustness analyses, controlling for various potential sources of bias in estimates. In Additional file 1: Fig. S3 we show that despite adverse outcomes concentrating in the older ages (giving them more statistical weight), the estimates for younger ages are robust to excluding older ages from the regression. In Additional file 1: Fig. S4 we show that our estimates are also robust to excluding the locations where prevalence was estimated from nonrepresentative samples, which have increased risk of bias. In Additional file 1: Fig. S6, we show that our estimates are robust to excluding the locations with the fastest changing epidemics at the time of data collection, and are thus robust to the choice of dates for outcome data collection. Then, in Additional file 1: Fig. S7 we show that our results are robust to excluding the locations with the longest delays between epidemic wave and seroprevalence study, which are the most susceptible to seroreversion, and thus our estimates are not strongly affected by seroreversion. Finally, besides these comprehensive robustness analyses, a highlight of our study is the validation of our results with an independent estimation method, based on the ratioofratios approach (Fig.Â 2).
Our results are highly relevant for aspects of COVID19 modeling, such as estimating the number of unreported infections from hospital and ICU data, allowing to better estimate the present levels of natural immunity [4, 12, 13]; prediction of the effects of publichealth policies implemented along the pandemic [14]; evaluating policy decisions such as vaccine allocation [15, 16]; the prediction of health outcomes in countries with high or low vaccination rates that account for the agedistribution of each country [17, 18]. Particularly, our estimates are important for analyzing the risk of COVID19 for younger populations. These populations have very low risks of death, but as seen in our estimates, the risk of severe disease can be 2 orders of magnitude larger than the risk of death, and thus severe and critical outcomes are essential to properly characterize the risk of these populations. One illustrative example is the discussion around vaccination of young individuals against COVID19. The FDA estimates that the rate of mRNAvaccineinduced myocarditis, a side effect which is mild in some cases but which can result in severe outcomes, is of 1/5000 for males in ages between 16â€“17, the population at highest risk [19] (in line with reports from other locations such as Israel [20]). Although doing a riskâ€“benefit analysis for adolescent vaccination is very complicated, and outside of the scope of this work, it is notable that while our estimated IFR is 0.12 [0.07â€“0.20] times the rate of vaccine induced myocarditis (i.e. 8 times smaller), our estimated ISR is 12 [7â€“19] times larger this rate. Thus, although such a direct comparison has many limitations, and should not be taken as a riskâ€“benefit analysis, it shows how radically the conclusions of risk evaluation for young individuals depends on the disease outcome being considered.
Importantly, we note that the dynamic nature of the COVID19 pandemic makes any estimates of outcome rates transient, since those rates are expected to change in space and time as new variants emerge and social behavior and medical practices change. As such, generalization of our estimates across time and space requires caution. For example, substantial drops in hospital and ICU mortality between the first and second waves have been reported in developed countries [21,22,23,24,25,26], partly due to improvements in care, although this is accompanied by considerable geographical heterogeneity [27]. On the other hand, the emergence of variants of concern was associated with increased rates of hospital mortality [28,29,30] and severe, critical and fatal cases [31,32,33,34,35,36,37] for early variants (Alpha, Gamma, Delta), and with a decrease in disease severity for the posterior variant Omicron [38]. Recent changes in disease severity are also due to the introduction of effective vaccines, with reductions in the rates of hospitalization and death of over 90% for the BNT162b2 vaccine [39]. More recently, oral antiviral medications have become available that can reduce the risk of hospitalization from COVID19 close to 90% in patients at high risk of developing severe COVID19 [40]. However, the effects of all of these changes on disease severity have been estimated in relative terms, rather than in absolute changes in risk. The estimates of severe and critical disease that we provide in this work can thus serve as the baseline to estimate the absolute risks of COVID19 after such changes (e.g. the risk of severe disease for vaccinated individuals, or for unvaccinated individuals in the presence of Omicron). One particularly important problem where we may wish to calculate such absolute risks of COVID19 severe and critical disease is to anticipate the effects of variants with potential for immune escape [41], or the effects of waning of vaccine effectiveness [42].
Finally, we note that the current work has some limitations that should be considered. One limitation is that our estimates rely on publicly available data on the number of hospitalizations, ICU admissions, and deaths, which may not match exactly the real number of severe and critical cases. Hospitalizations can underestimate the number of severe cases if a health system is overwhelmed, and some severe cases are not admitted or are not reported properly. On the other hand, the number of COVID19related hospitalizations may overestimate the number of severe cases if some people are admitted without severe COVID19 infection, for example due to an abundance of caution in lowoccupancy situations, or due to an incidental positive test at the time of admission in highprevalence situations. Another limitation stems from the fact that the protocols for testing, admitting, and treating patients can vary between locations, and may depend on the strain of the healthcare system. For example, a large COVID19 wave may induce a higher rate of severe disease, due to limited treatment capacity, and conversely to a smaller rate of hospital admissions due to limited resources. The underlying health of each population, and quality of care will also determine the outcome rates as specific locations. For example, race and socioeconomic status have been reported as risk factors for critical COVID19 [43, 44], and striking differences in COVID19 IFR between developed and developing countries have been reported [45]. In line with this, our statistical analysis shows that there is significant variability in the adverse outcome rates in between locations (in our models this is captured by the standard deviation of the intercept and slope of the logISR and the logICRage regression, see Additional file 1: Table S2). Thus, although we rigorously take variability between countries into account in our estimates, caution is required when extrapolating our estimates of the average ISR and ICR to specific locations.
Availability of data and materials
All data and code used in this project are available at https://github.com/dherrera1911/estimate_covid_severity.
Abbreviations
 IFR:

Infection fatality rate
 ICR:

Infectioncritical rate
 ISR:

Infectionsevere rate
 FDA:

Food and Drug Administration
 ICU:

Intensive Care Unit
 SFR:

Severe fatality rate
References
Levin AT, Hanage WP, OwusuBoaitey N, Cochran KB, Walsh SP, MeyerowitzKatz G. Assessing the age specificity of infection fatality rates for COVID19: systematic review, metaanalysis, and public policy implications. Eur J Epidemiol. 2020;35(12):1123â€“38.
Oâ€™Driscoll M, Dos Santos GR, Wang L, Cummings DAT, Azman AS, Paireau J, et al. Agespecific mortality and immunity patterns of SARSCoV2. Nature [Internet]. 2020 Nov 2 [cited 2020 Nov 5]; Available from: http://www.nature.com/articles/s4158602029180.
Brazeau N, Verity R, Jenks S, Fu H, Whittaker C, Winskill P, et al. Report 34: COVID19 infection fatality ratio: estimates from seroprevalence [Internet]. Imperial College London; 2020 Oct [cited 2020 Nov 2]. Available from: http://spiral.imperial.ac.uk/handle/10044/1/83545.
HozÃ© N, Paireau J, Lapidus N, Tran Kiem C, Salje H, Severi G, et al. Monitoring the proportion of the population infected by SARSCoV2 using agestratified hospitalisation and serological data: a modelling study. Lancet Public Health. 2021;6(6):e408â€“15.
Lapidus N, Paireau J, LevyBruhl D, de Lamballerie X, Severi G, Touvier M, et al. Do not neglect SARSCoV2 hospitalization and fatality risks in the middleaged adult population. Infect Dis Now. 2021;51(4):380â€“2.
Espenhain L, Tribler S, SvÃ¦rke JÃ¸rgensen C, Holm Hansen C, Wolff SÃ¶nksen U, Ethelberg S. Prevalence of SARSCoV2 antibodies in Denmark: nationwide, populationbased seroepidemiological study. Eur J Epidemiol. 2021. https://doi.org/10.1007/s10654021007968.
Menachemi N, Dixon BE, WoolsKaloustian KK, Yiannoutsos CT, Halverson PK. How many SARSCoV2infected people require hospitalization? Using random sample testing to better inform preparedness efforts. J Public Health Manag Pract. 2021;27(3):246â€“50.
Mahajan S, Caraballo C, Li SX, Dong Y, Chen L, Huston SK, et al. SARSCoV2 infection hospitalization rate and infection fatality rate among the noncongregate population in Connecticut. Am J Med. 2021;134(6):812816.e2.
Seedat S, Chemaitelly H, Ayoub HH, Makhoul M, Mumtaz GR, Al Kanaani Z, et al. SARSCoV2 infection hospitalization, severity, criticality, and fatality rates in Qatar. Sci Rep. 2021;11(1):18182.
Verity R, Okell LC, Dorigatti I, Winskill P, Whittaker C, Imai N, et al. Estimates of the severity of coronavirus disease 2019: a modelbased analysis. Lancet Infect Dis [Internet]. 2020 Mar 30 [cited 2020 Apr 3]; Available from: http://www.sciencedirect.com/science/article/pii/S1473309920302437.
Blackburn J, Yiannoutsos CT, Carroll AE, Halverson PK, Menachemi N. Infection fatality ratios for COVID19 among noninstitutionalized persons 12 and older: results of a randomsample prevalence study. Ann Intern Med. 2021;174(1):135â€“6.
Irons NJ, Raftery AE. Estimating SARSCoV2 infections from deaths, confirmed cases, tests, and random surveys. Proc Natl Acad Sci [Internet]. 2021 Aug 3 [cited 2021 Jul 28];118(31). Available from: https://www.pnas.org/content/118/31/e2103272118.
Russell TW, Golding N, Hellewell J, Abbott S, Wright L, Pearson CAB, et al. Reconstructing the early global dynamics of underascertained COVID19 cases and infections. BMC Med. 2020;18(1):332.
Davies NG, Barnard RC, Jarvis CI, Russell TW, Semple MG, Jit M, et al. Association of tiered restrictions and a second lockdown with COVID19 deaths and hospital admissions in England: a modelling study. Lancet Infect Dis. 2021;21(4):482â€“92.
Hjorleifsson KE, Rognvaldsson S, Jonsson H, Agustsdottir AB, Andresdottir M, Birgisdottir K, et al. Reconstruction of a largescale outbreak of SARSCoV2 infection in Iceland informs vaccination strategies. Clin Microbiol Infect [Internet]. 2022 Feb 16 [cited 2022 Feb 25];0(0). Available from: https://www.clinicalmicrobiologyandinfection.com/article/S1198743X(22)000854/fulltext.
Matrajt L, Eaton J, Leung T, Brown ER. Vaccine optimization for COVID19: who to vaccinate first? Sci Adv. 2020;7(6):eabf1374.
Davies NG, Kucharski AJ, Eggo RM, Gimma A, Edmunds WJ, Jombart T, et al. Effects of nonpharmaceutical interventions on COVID19 cases, deaths, and demand for hospital services in the UK: a modelling study. Lancet Public Health [Internet]. 2020 Jun 2 [cited 2020 Jun 25]; Available from: http://www.sciencedirect.com/science/article/pii/S246826672030133X.
Sandmann FG, Davies NG, Vassall A, Edmunds WJ, Jit M, Sun FY, et al. The potential health and economic value of SARSCoV2 vaccination alongside physical distancing in the UK: a transmission modelbased future scenario analysis and economic evaluation. Lancet Infect Dis. 2021;21(7):962â€“74.
FDA. August 23, 2021 Summary basis for regulatory actionâ€”comirnaty [Internet]. [cited 2021 Sep 30]. Available from: https://www.fda.gov/media/151733/.
Mevorach D, Anis E, Cedar N, Bromberg M, Haas EJ, Nadir E, et al. Myocarditis after BNT162b2 mRNA vaccine against Covid19 in Israel. N Engl J Med. 2021;385:2140â€“9.
Anesi GL, Jablonski J, Harhay MO, Atkins JH, Bajaj J, Baston C, et al. Characteristics, outcomes, and trends of patients with COVID19related critical illness at a learning health system in the United States. Ann Intern Med. 2021;174(5):613â€“21.
Asch DA, Sheils NE, Islam MN, Chen Y, Werner RM, Buresh J, et al. Variation in US hospital mortality rates for patients admitted with COVID19 during the first 6 months of the pandemic. JAMA Intern Med. 2021;181(4):471â€“8.
AznarGimeno R, PaÃ±oPardo JR, Esteban LM, LabataLezaun G, EsquillorRodrigo MJ, Lanas A, et al. Changes in severity, mortality, and virus genome among a Spanish cohort of patients hospitalized with SARSCoV2. Sci Rep. 2021;11(1):18844.
Navaratnam AV, Gray WK, Day J, Wendon J, Briggs TWR. Patient factors and temporal trends associated with COVID19 inhospital mortality in England: an observational study using administrative data. Lancet Respir Med. 2021;9(4):397â€“406.
Doidge JC, Gould DW, FerrandoVivas P, Mouncey PR, Thomas K, ShankarHari M, et al. Trends in intensive care for patients with COVID19 in England, Wales, and Northern Ireland. Am J Respir Crit Care Med. 2021;203(5):565â€“74.
Armstrong RA, Kane AD, Kursumovic E, Oglesby FC, Cook TM. Mortality in patients admitted to intensive care with COVID19: an updated systematic review and metaanalysis of observational studies. Anaesthesia. 2021;76(4):537â€“48.
Kadri SS, Sun J, Lawandi A, Strich JR, Busch LM, Keller M, et al. Association between caseload surge and COVID19 survival in 558 U.S. Hospitals, March to August 2020. Ann Intern Med. 2021;174(9):1240â€“51.
Freitas ARR, Beckedorff OA, Cavalcanti de LPG, Siqueira AM, Castro de DB, Costa da CF, et al. The emergence of novel SARSCoV2 variant P.1 in Amazonas (Brazil) was temporally associated with a change in the age and sex profile of COVID19 mortality: a population based ecological study. Lancet Reg Health Am. 2021;1:100021.
Jassat W, Mudara C, Ozougwu L, Tempia S, Blumberg L, Davies MA, et al. Increased mortality among individuals hospitalised with COVID19 during the second wave in South Africa 2021. https://doi.org/10.1101/2021.03.09.21253184v1.
Khedar RS, Mittal K, Ambaliya HC, Mathur A, Gupta JB, Sharma KK, et al. Greater COVID19 severity and mortality in hospitalized patients in second (delta variant) wave compared to the first: single centre prospective study in India. 2021. https://doi.org/10.1101/2021.09.03.21263091v1.
Challen R, BrooksPollock E, Read JM, Dyson L, TsanevaAtanasova K, Danon L. Risk of mortality in patients infected with SARSCoV2 variant of concern 202012/1: matched cohort study. BMJ. 2021;10(372):n579.
Davies NG, Jarvis CI, Edmunds WJ, Jewell NP, DiazOrdaz K, Keogh RH. Increased mortality in communitytested cases of SARSCoV2 lineage B.1.1.7. Nature. 2021;593(7858):270â€“4.
Fisman DN, Tuite AR. Progressive increase in virulence of novel SARSCoV2 variants in Ontario, Canada 2021. https://doi.org/10.1101/2021.07.05.21260050v3.
Fisman DN, Tuite AR. Agespecific changes in virulence associated with SARSCoV2 variants of concern. 2021; https://doi.org/10.1101/2021.09.25.21264097v1.
Grint DJ, Wing K, Williamson E, McDonald HI, Bhaskaran K, Evans D, et al. Case fatality risk of the SARSCoV2 variant of concern B.1.1.7 in England, 16 November to 5 February. Eurosurveillance. 2021;26(11):2100256.
Ong SWX, Chiew CJ, Ang LW, Mak TM, Cui L, Toh MPH, et al. Clinical and virological features of SARSCoV2 variants of concern: a retrospective cohort study comparing B.1.1.7 (Alpha), B.1.315 (Beta), and B.1.617.2 (Delta) [Internet]. Rochester, NY: Social Science Research Network; 2021 Jun [cited 2021 Sep 30]. Report No.: ID 3861566. Available from: https://papers.ssrn.com/abstract=3861566.
Tuite AR, Fisman DN, Odutayo A, Bobos P, Allen V, Bogoch II, et al. COVID19 Hospitalizations, ICU admissions and deaths associated with the new variants of concern [Internet]. Ontario COVID19 Science Advisory Table; 2021 Mar [cited 2021 Sep 30]. Available from: https://covid19sciencetable.ca/sciencebrief/covid19hospitalizationsicuadmissionsanddeathsassociatedwiththenewvariantsofconcern.
Iuliano AD. Trends in disease severity and health care utilization during the early omicron variant period compared with previous SARSCoV2 High Transmission Periodsâ€”United States, December 2020â€“January 2022. MMWR Morb Mortal Wkly Rep [Internet]. 2022 [cited 2022 Feb 23];71. Available from: https://www.cdc.gov/mmwr/volumes/71/wr/mm7104e4.htm.
Andrews N, Tessier E, Stowe J, Gower C, Kirsebom F, Simmons R, et al. Duration of protection against mild and severe disease by COVID19 vaccines. N Engl J Med. 2022;386(4):340â€“50.
Hammond J, LeisterTebbe H, Gardner A, Abreu P, Bao W, Wisemandle W, et al. Oral nirmatrelvir for highrisk, nonhospitalized adults with COVID19. N Engl J Med. 2022.
Collie S, Champion J, Moultrie H, Bekker LG, Gray G. Effectiveness of BNT162b2 vaccine against omicron variant in South Africa. N Engl J Med. 2022;386(5):494â€“6.
Feikin DR, Higdon MM, AbuRaddad LJ, Andrews N, Araos R, Goldberg Y, et al. Duration of effectiveness of vaccines against SARSCoV2 infection and COVID19 disease: results of a systematic review and metaregression. The Lancet [Internet]. 2022 Feb 21 [cited 2022 Feb 23];0(0). Available from: https://www.thelancet.com/journals/lancet/article/PIIS01406736(22)001520/fulltext.
Quan D, Luna Wong L, Shallal A, Madan R, Hamdan A, Ahdi H, et al. Impact of race and socioeconomic status on outcomes in patients hospitalized with COVID19. J Gen Intern Med. 2021;36(5):1302â€“9.
Patel AP, Paranjpe MD, Kathiresan NP, Rivas MA, Khera AV. Race, socioeconomic deprivation, and hospitalization for COVID19 in English participants of a national biobank. Int J Equity Health. 2020;19(1):114.
Levin AT, OwusuBoaitey N, Pugh S, Fosdick BK, Zwi AB, Malani A, et al. Assessing the burden of COVID19 in developing countries: systematic review, metaanalysis, and public policy implications. medRxiv. 2021;396:10258.
Gudbjartsson DF, Norddahl GL, Melsted P, Gunnarsdottir K, Holm H, Eythorsson E, et al. Humoral immune response to SARSCoV2 in Iceland. N Engl J Med. 2020;383(18):1724â€“34.
Waterfield T, Watson C, Moore R, Ferris K, Tonry C, Watt A, et al. Seroprevalence of SARSCoV2 antibodies in children: a prospective multicentre cohort study. Arch Dis Child. 2021;106(7):680â€“6.
Ward H, Atchison C, Whitaker M, Ainslie KEC, Elliott J, Okell L, et al. SARSCoV2 antibody prevalence in England following the first peak of the pandemic. Nat Commun. 2021;12(1):905.
Carrat F, de Lamballerie X, Rahib D, BlanchÃ© H, Lapidus N, Artaud F, et al. Seroprevalence of SARSCoV2 Among Adults in Three Regions of France Following the Lockdown and Associated Risk Factors: A Multicohort Study [Internet]. Rochester, NY: Social Science Research Network; 2020 Oct [cited 2021 Jan 6]. Report No.: ID 3696820. Available from: https://papers.ssrn.com/abstract=3696820.
Ireland Health Service Executive. Preliminary report of the results of the Study to Investigate COVID19 Infection in People Living in Ireland (SCOPI): a national seroprevalence study, JuneJuly 2020 [Internet]. 2020 [cited 2021 Jul 4]. Available from: https://www.hpsc.ie/az/respiratory/coronavirus/novelcoronavirus/scopi/.
Vos ERA, den Hartog G, Schepp RM, Kaaijk P, van Vliet J, Helm K, et al. Nationwide seroprevalence of SARSCoV2 and identification of risk factors in the general population of the Netherlands during the first epidemic wave. J Epidemiol Commun Health. 2020;75:489â€“95.
PollÃ¡n M, PÃ©rezGÃ³mez B, PastorBarriuso R, Oteo J, HernÃ¡n MA, PÃ©rezOlmeda M, et al. Prevalence of SARSCoV2 in Spain (ENECOVID): a nationwide, populationbased seroepidemiological study. The Lancet. 2020;396(10250):535â€“44.
Biggs HM. Estimated community seroprevalence of SARSCoV2 antibodiesâ€”two Georgia Counties, April 28â€“May 3, 2020. MMWR Morb Mortal Wkly Rep [Internet]. 2020 [cited 2021 Jan 20];69. Available from: https://www.cdc.gov/mmwr/volumes/69/wr/mm6929e2.htm.
Rosenberg ES, Tesoriero JM, Rosenthal EM, Chung R, Barranco MA, Styer LM, et al. Cumulative incidence and diagnosis of SARSCoV2 infection in New York. Ann Epidemiol. 2020;1(48):2329.e4.
Ontario Public Health. COVID19 Seroprevalence in Ontario: July 4 to July 31, 2020 [Internet]. 2020 p. 9. Available from: https://www.publichealthontario.ca//media/documents/ncov/epi/2020/10/covid19episeroprevalenceinontariojuly31.pdf?la=en.
Sweden Public Health Authority. Veckorapport om covid19, vecka 24 [Internet]. 2020 [cited 2021 Jul 4]. Available from: https://www.folkhalsomyndigheten.se/globalassets/statistikuppfoljning/smittsammasjukdomar/veckorapportercovid19/2020/covid19veckorapportvecka24_final.pdf.
PerezSaez J, Lauer SA, Kaiser L, Regard S, Delaporte E, Guessous I, et al. Serologyinformed estimates of SARSCoV2 infection fatality risk in Geneva, Switzerland. Lancet Infect Dis [Internet]. 2020 Jul 14 [cited 2021 Jan 17];0(0). Available from: https://www.thelancet.com/journals/laninf/article/PIIS14733099(20)305843/abstract.
Herzog S, Bie JD, Abrams S, Wouters I, Ekinci E, Patteet L, et al. Seroprevalence of IgG antibodies against SARS coronavirus 2 in Belgiumâ€”a serial prospective crosssectional nationwide study of residual samples. medRxiv. 2020;142:251.
Mahajan S, Srinivasan R, Redlich CA, Huston SK, Anastasio KM, Cashman L, et al. Seroprevalence of SARSCoV2Specific IgG antibodies among adults living in Connecticut: postinfection prevalence (PIP) study. Am J Med. 2021;134(4):526534.e11.
Chishinga N, Gandhi NR, Onwubiko UN, Telford C, Prieto J, Smith S, et al. Characteristics and risk factors for hospitalization and mortality among persons with COVID19 in Atlanta Metropolitan Area. Infect Diseases (except HIV/AIDS). 2020. https://doi.org/10.1101/2020.12.15.20248214.
Molenberghs G, Faes C, Aerts J, Theeten H, Devleesschauwer B, Sierra NB, et al. Belgian COVID19 mortality, excess deaths, number of deaths per million, and infection fatality rates (8 Marchâ€“9 May 2020). medRxiv. 2020. https://doi.org/10.1101/2020.06.20.20136234.
Eythorsson E, Helgason D, Ingvarsson RF, Bjornsson HK, Olafsdottir LB, Bjarnadottir V, et al. Clinical spectrum of coronavirus disease 2019 in Iceland: population based cohort study. BMJ. 2020;2(371):4 4529.
Park HY, Lee JH, Lim NK, Lim DS, Hong SO, Park MJ, et al. Presenting characteristics and clinical outcome of patients with COVID19 in South Korea: A nationwide retrospective observational study. Lancet Reg Healthâ€”West Pac [Internet]. 2020 Dec 1 [cited 2021 Jan 5];5. Available from: https://www.thelancet.com/journals/lanwpc/article/PIIS26666065(20)300614/abstract.
Carpenter B, Gelman A, Hoffman MD, Lee D, Goodrich B, Betancourt M, et al. Stan: a probabilistic programming language. J Stat Softw. 2017;76(1):1â€“32.
Cummings MJ, Baldwin MR, Abrams D, Jacobson SD, Meyer BJ, Balough EM, et al. Epidemiology, clinical course, and outcomes of critically ill adults with COVID19 in New York City: a prospective cohort study. The Lancet. 2020;395(10239):1763â€“70.
COVIDICU Group on behalf of the REVA Network and the COVIDICU Investigators. Clinical characteristics and day90 outcomes of 4244 critically ill adults with COVID19: a prospective cohort study. Intensive Care Med. 2021;47(1):60â€“73.
Ranzani OT, Bastos LSL, Gelli JGM, Marchesi JF, BaiÃ£o F, Hamacher S, et al. Characterisation of the first 250 000 hospital admissions for COVID19 in Brazil: a retrospective analysis of nationwide data. Lancet Respir Med. 2021;9(4):4Â 7418.
Oliveira E, Parikh A, LopezRuiz A, Carrilo M, Goldberg J, Cearras M, et al. ICU outcomes and survival in patients with severe COVID19 in the largest health care system in central Florida. PLoS ONE. 2021;16(3):e0249038.
GonzÃ¡lezDambrauskas S, VÃ¡squezHoyos P, Camporesi A, DÃazRubio F, PiÃ±eresOlave BE, FernÃ¡ndezSarmiento J, et al. Pediatric critical care and COVID19. Pediatrics. 2020;146(3):e20201766.
PrataBarbosa A, LimaSetta F, dos Santos GR, Lanziotti VS, de Castro REV, de Souza DC, et al. Pediatric patients with COVID19 admitted to intensive care units in Brazil: a prospective multicenter study. J Pediatr (Rio J). 2020;96(5):582â€“92.
Richardson S, Hirsch JS, Narasimhan M, Crawford JM, McGinn T, Davidson KW, et al. Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized with COVID19 in the New York City area. JAMA. 2020;323(20):2052.
Karagiannidis C, Mostert C, Hentschker C, Voshaar T, Malzahn J, Schillinger G, et al. Case characteristics, resource use, and outcomes of 10 021 patients with COVID19 admitted to 920 German hospitals: an observational study. Lancet Respir Med. 2020;8(9):853â€“62.
Salje H, Kiem CT, Lefrancq N, Courtejoie N, Bosetti P, Paireau J, et al. Estimating the burden of SARSCoV2 in France. Science [Internet]. 2020 May 13 [cited 2020 May 13]; Available from: https://science.sciencemag.org/content/early/2020/05/12/science.abc3517.
Docherty AB, Harrison EM, Green CA, Hardwick HE, Pius R, Norman L, et al. Features of 20 133 UK patients in hospital with COVID19 using the ISARIC WHO clinical characterisation protocol: prospective observational cohort study. BMJ. 2020;22(369):m1985.
Berenguer J, Ryan P, RodrÃguezBaÃ±o J, JarrÃn I, CarratalÃ J, PachÃ³n J, et al. Characteristics and predictors of death among 4035 consecutively hospitalized patients with COVID19 in Spain. Clin Microbiol Infect. 2020;26(11):1525â€“36.
Maquilon C, Gongora J, Antolini M, Alvarado B, Valdes N, Benavente A, et al. Risk factors on admission and condition at discharge of 529 consecutive COVID19 patients at a tertiary care center in Santiago, Chile [Internet]. In Review; 2020 Sep [cited 2021 Jun 15]. Available from: https://www.researchsquare.com/article/rs71187/v1.
Acknowledgements
We thank ElÃas EyÃ¾Ã³rsson, Gideon MeyerowitzKatz, Nicholas Davies and Nana OwusuBoaitey for useful discussions and comments on this work.
Funding
The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.
Author information
Authors and Affiliations
Contributions
DHE conceived the study, collected the data, wrote the original draft and wrote the analysis code. DHE and GC designed the data analysis, interpreted the results, and edited the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable. Our research only utilized publiclyavailable information from published articles, preprints, and official government sources.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1
: Description of data: Robustness checks of the main analysis.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
HerreraEsposito, D., de los Campos, G. Agespecific rate of severe and critical SARSCoV2 infections estimated with multicountry seroprevalence studies. BMC Infect Dis 22, 311 (2022). https://doi.org/10.1186/s12879022072620
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12879022072620
Keywords
 SARSCoV2
 COVID19
 Severity
 Critical disease
 Metaanalysis
 Serology