Hostname: page-component-7c8c6479df-7qhmt Total loading time: 0 Render date: 2024-03-28T17:21:00.924Z Has data issue: false hasContentIssue false

Real-time estimation of the hospitalization fatality risk of influenza A(H1N1)pdm09 in Hong Kong

Published online by Cambridge University Press:  18 December 2015

J. Y. WONG
Affiliation:
School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong Special Administrative Region, China
P. WU*
Affiliation:
School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong Special Administrative Region, China
E. H. Y. LAU
Affiliation:
School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong Special Administrative Region, China
T. K. TSANG
Affiliation:
School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong Special Administrative Region, China
V. J. FANG
Affiliation:
School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong Special Administrative Region, China
L.-M. HO
Affiliation:
School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong Special Administrative Region, China
B. J. COWLING
Affiliation:
School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong Special Administrative Region, China
*
*Author for correspondence: Dr P. Wu, School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, 21 Sassoon Road, Hong Kong Special Administrative Region, China. (Email: pengwu@hku.hk)
Rights & Permissions [Opens in a new window]

Summary

During the early stage of an epidemic, timely and reliable estimation of the severity of infections are important for predicting the impact that the influenza viruses will have in the population. We obtained age-specific deaths and hospitalizations for patients with laboratory-confirmed H1N1pdm09 infections from June 2009 to December 2009 in Hong Kong. We retrospectively obtained the real-time estimates of the hospitalization fatality risk (HFR), using crude estimation or allowing for right-censoring for final status in some patients. Models accounting for right-censoring performed better than models without adjustments. The risk of deaths in hospitalized patients with confirmed H1N1pdm09 increased with age. Reliable estimates of the HFR could be obtained before the peak of the first wave of H1N1pdm09 in young and middle-aged adults but after the peak in the elderly. In the next influenza pandemic, timely estimation of the HFR will contribute to risk assessment and disease control.

Type
Original Papers
Copyright
Copyright © Cambridge University Press 2015 

INTRODUCTION

During the early stage of an epidemic of a novel influenza strain, timely and reliable estimates of the severity of infections are important for predicting the impact that the influenza viruses will have in the population, and for calibrating the public health response [Reference Wu1]. The hospitalization fatality risk (HFR), a measure of the severity of infection, is defined as the probability of death in cases who required hospitalization for medical reasons [Reference Yu2, Reference Borja-Aburto3]. A well-known complication with the estimation of severity of infection in ‘real time’ during an epidemic is the uncertainty about the final outcomes of some patients that remain in hospital at the time of analysis [Reference Ghani4]. This led to substantial underestimation of the HFR of H7N9 early in the outbreak in China in 2013 [Reference Yu2], and similar issues were also documented in the 2009 influenza pandemic [Reference Garske5] and the 2003 epidemic of SARS [Reference Ghani4, Reference Yip6, Reference Jewell7].

The first objective of this study was to estimate the HFR of pandemic influenza A(H1N1)pdm09 virus infections (H1N1pdm09) during the first wave in Hong Kong in 2009. The second objective was to compare different statistical methods for real-time estimation of the HFR.

METHODS

Sources of data

Data on individual patients with laboratory-confirmed H1N1pdm09 infection from 30 June 2009 to 31 December 2009 were obtained from the Hong Kong Hospital Authority (the e-flu database) [Reference Cowling8]. Epidemiological information was provided in the database, including age, dates of laboratory confirmation, first hospital admission, death and discharge. In this analysis we used data on all patients for whom laboratory confirmation of H1N1pdm09 was obtained from 7 days preceding to 28 days after hospital admission.

Statistical analysis

We used three alternative statistical methods to estimate the HFR in real time, i.e. (1) division of cumulative deaths by cumulative cases (HFR1); (2) division of cumulative deaths by cumulative resolved cases (HFR2); and (3) the method proposed by Garske et al. [Reference Garske5], which is based on the cumulative cases and deaths and a fitted Weibull distribution (HFR3) for the hospitalization-to-death distribution. The three estimators are defined as follows:

$${\rm HF}{{\rm R}_1}\left( t \right)\; = \; \displaystyle{{D\left( t \right)} \over {C\left( t \right)}},$$
$${\rm HF}{{\rm R}_2}\left( t \right)\; = \; \displaystyle{{D\left( t \right)} \over {D\left( t \right) + R\left( t \right)}},$$
$${\rm HF}{{\rm R}_3}\left( t \right){\rm \;} = {\rm \;} \displaystyle{{D\left( t \right)} \over {\mathop \sum \nolimits_{u = 0}^t c\left( u \right)F\left( {t - u} \right)}},$$

where HFR j (t) is the hospitalization fatality risk estimated on day t under method j; D(t) is the cumulative number of deaths on day t; C(t) is the cumulative number of hospitalized cases on day t; R(t) is the cumulative number of recovered cases on day t; c(t) is the number of hospitalized cases on day t; F(.) is the Weibull cumulative distribution function of time from hospitalization to death based on data available at time t.

HFR1 includes all hospitalized patients in the denominator, and will tend to underestimate the true HFR if there are many unresolved cases [Reference Ghani4, Reference Yip6]. HFR2 includes only patients with known outcome (i.e. death/recovery) and may not be able to estimate the HFR in the early stages of an epidemic. HFR3 allows for unresolved outcomes by accounting for the interval between hospital admission and death [Reference Garske5], assuming that the distribution of the interval from hospital admission to death for unresolved patients is consistent with the distribution observed from patients with known outcomes.

Each estimator was applied to the time series of daily hospitalization data from June to December 2009. The data were stratified into three age groups: 20–44 years (young adults), 45–64 years (middle-aged adults) and ⩾65 years (elderly). For each age group, estimates were compared based on the three different statistical methods. In each method, a real-time dataset was reconstructed at different time points incorporating reporting delays in laboratory confirmation, but assuming no reporting delay on outcomes. Estimates using HFR1 and HFR2 were made starting from the first death, while estimates using HFR3 were started after five deaths had occurred because the estimated hospitalization-to-death distribution was unstable with fewer deaths. The HFRs and interval distribution were estimated in a Bayesian framework in order to retain uncertainty in the hospitalization-to-death distribution for HFR3. Uniform priors over the entire range of possible values for all parameters were used for all parameters. Convergence of Markov Chain Monte Carlo algorithms was judged using the R-hat criteria [Reference Brooks9]. Posterior medians and 95% credible intervals (CrI) are reported. All analyses were conducted with R version 3.0.2 (R Foundation for Statistical Computing, Austria).

RESULTS

H1N1pdm09 virus emerged in Hong Kong in May 2009 and caused one major epidemic wave peaking in September–October in that year (Supplementary Fig. S1). A total of 6458 hospitalized patients were included in the e-flu database with laboratory-confirmation from 7 days before to 28 days after admission. Mortality in hospitalized patients was 0·02% (1/4309), 1·29% (13/1005), 3·38% (24/710), and 4·61% (20/434) in the <20, 20–44, 45–64 and ⩾65 years age groups, respectively. HFRs were only estimated for age groups aged ⩾20 years because there was only one death in the <20 years age group.

We examined the hospitalization-to-death intervals using age-specific data from June to December 2009 (Supplementary Fig. S2). The mean hospitalization-to-death intervals were 7·68 days (95% CrI 3·95–23·10), 7·40 days (95% CrI 4·85–12·49) and 12·59 days (95% CrI 8·16–23·54) in the 20–44, 45–64 and ⩾65 years age groups, respectively. However, the Weibull distribution might underestimate the probability of long hospitalization-to-death intervals earlier in the epidemic with fewer deaths (Supplementary Fig. S3).

Point estimates of age-specific HFR1 at the end of 2009 gradually increased with age from 1·36% in the 20–44 years age group to 4·76% in the ⩾65 years age group (Fig. 1). The final HFR3 estimates were very similar, increasing from 1·37% (95% CrI 0·78–2·25) to 4·94% (95% CrI 3·12–7·37) for individuals aged 20–44 to ⩾65 years, respectively (Fig. 1) using a Weibull distribution fitted to the hospitalization-to-death intervals.

Fig. 1. Age-specific real-time estimates of the hospitalization fatality risk from three different methods in Hong Kong, 2009. Black solid lines represent the posterior median, dotted lines show the 95% credible interval, and the grey solid line represents the posterior median of HFR3 on 31 December 2009. HFR, Hospitalization fatality risk.

We applied the three different estimators to the H1N1pdm09 data in real time from August 2009 onwards and found that the HFR estimate was relatively stable (Fig. 1), and HFR estimates were similar for the three methods at different time points (Supplementary Table S1). Estimators that accounted for censoring (HFR2–HFR3) tended to perform better than the HFR1 estimator without adjustment (Fig. 1), and quite consistent estimates of HFR could be obtained from the beginning of September for adults, and early October for the elderly with a much smaller number of hospitalizations.

DISCUSSION

During an ongoing influenza epidemic, final outcomes of some hospitalized cases will not be known at the time of data analysis. This was ignored in some early studies that determined the HFR of H1N1pdm09 and could lead to an inaccurate estimation of the HFR [Reference Ghani4Reference Yip6]. We used two methods, with and without adjustment for censoring, to estimate the HFR of H1N1pdm09 virus infection in 2009 in Hong Kong. Estimators allowing for censoring performed better than models without although we did not find substantial difference in estimates of HFR in real time between these two methods. However, in the SARS epidemic, substantial differences were identified between crude and adjusted estimates [Reference Ghani4].

Previous studies showed that the elderly were at a particular high risk of developing severe illness if infected with H1N1pdm09, indicated by a higher estimate of infection hospitalization risk compared to younger individuals [Reference Riley10, Reference Steens11]. Together with the age-stratified estimates of HFR from our study showing a higher risk of fatality in hospitalized elderly, this suggested that H1N1pdm09 infection was more likely to lead to fatal events in the elderly. Early treatment and better medical support to reduce severe complications after infection for this age group are needed. In addition, we were able to obtain a timely estimate of HFR for adults in early September, before the peak of the first pandemic wave, while such an estimate could only be achieved for the elderly after the peak. This could be due to the smaller number of confirmed deaths in the elderly.

There are several limitations of this study. First, HFR was assumed to be constant over the course of the epidemic. However, this assumption may be violated if there are changes in case ascertainment or improvements in treatment. Further work could explore methods for estimation of changes in the HFR over time, and the optimal timing and sample size necessary to obtain reliable HFR estimates. When moving from the containment phase to mitigation phase in the 2009 pandemic, there was change in hospital admission guidelines and it is possible that the earliest estimates would bias downwards because mild cases were also admitted to hospitals for isolation during the containment phase [Reference Wu1, Reference Garske5]. Therefore, we only analysed patients admitted to hospitals after the containment period (i.e. 1 May 2009 to 29 June 2009). Improvement in treatment of influenza would also decrease HFR estimates over time [Reference Yip6]. Second, HFR estimates might be biased by the choice of distribution functions for the hospitalization-to-death interval. The choice of distribution is uncertain and it may be censored [Reference Garske5]. Finally, the lack of information on the date of official announcement for each confirmed case prevented us from fully simulating the real-time situation during the epidemic, therefore we used the best available data instead as an alternative, i.e. dates of laboratory confirmation, instead.

SUPPLEMENTARY MATERIAL

For supplementary material accompanying this paper visit http://dx.doi.org/10.1017/S0950268815003179.

ACKNOWLEDGEMENTS

We acknowledge the Hospital Authority Strategy and Planning Division, Quality and Safety Division, and Information Technology Division and the Center for Health Protection for the collation of the e-flu database. We thank Edward Goldstein and Joseph Wu for helpful discussions.

This study was funded by the Harvard Center for Communicable Disease Dynamics from the National Institute of General Medical Sciences (grant no. U54 GM088558), the Area of Excellence Scheme of the University Grants Committee of Hong Kong (AoE/M-12/06), the Health and Medical Research Fund, Food and Health Bureau, Government of the Hong Kong Special Administrative Region (grant no. 14131432), and a commissioned grant from the Health and Medical Research Fund, Food and Health Bureau, Government of the Hong Kong Special Administrative Region. The funding bodies had no role in study design, data collection and analysis, preparation of the manuscript, or the decision to publish.

DECLARATION OF INTEREST

B.J.C. has received research funding from MedImmune Inc. and Sanofi Pasteur, and consults for Crucell NV. The other authors report no conflict of interest.

References

REFERENCES

1. Wu, JT, et al. Estimating infection attack rates and severity in real time during an influenza pandemic: analysis of serial cross-sectional serologic surveillance data. PLoS Medicine 2011; 8: e1001103.CrossRefGoogle ScholarPubMed
2. Yu, H, et al. Human infection with avian influenza A H7N9 virus: an assessment of clinical severity. Lancet 2013; 382: 138145.Google Scholar
3. Borja-Aburto, VH, et al. Epidemiological characterization of a fourth wave of pandemic A/H1N1 influenza in Mexico, winter 2011–2012: age shift and severity. Archives of Medical Research 2012; 43: 563570.CrossRefGoogle ScholarPubMed
4. Ghani, AC, et al. Methods for estimating the case fatality ratio for a novel, emerging infectious disease. American Journal of Epidemiology 2005; 162: 479486.Google Scholar
5. Garske, T, et al. Assessing the severity of the novel influenza A/H1N1 pandemic. British Medical Journal 2009; 339: b2840.Google Scholar
6. Yip, PS, et al. A comparison study of realtime fatality rates: severe acute respiratory syndrome in Hong Kong, Singapore, Taiwan, Toronto and Beijing, China. Journal of the Royal Statistical Society: Series A 2005; 168: 233243.Google Scholar
7. Jewell, NP, et al. Non-parametric estimation of the case fatality ratio with competing risks data: an application to severe acute respiratory syndrome (SARS). Statistics in Medicine 2007; 26: 19821998.Google Scholar
8. Cowling, BJ, et al. The effective reproduction number of pandemic influenza: prospective estimation. Epidemiology 2010; 21: 842846.Google Scholar
9. Brooks, S, et al. Handbook of Markov Chain Monte Carlo. Boca Raton: Chapman & Hall/CRC; 2011.Google Scholar
10. Riley, S, et al. Epidemiological characteristics of 2009 (H1N1) pandemic influenza based on paired sera from a longitudinal community cohort study. PLoS Medicine 2011; 8: e1000442.Google Scholar
11. Steens, A, et al. Age-dependent patterns of infection and severity explaining the low impact of 2009 influenza A (H1N1): evidence from serial serologic surveys in the Netherlands. American Journal of Epidemiology 2011; 174: 13071315.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. Age-specific real-time estimates of the hospitalization fatality risk from three different methods in Hong Kong, 2009. Black solid lines represent the posterior median, dotted lines show the 95% credible interval, and the grey solid line represents the posterior median of HFR3 on 31 December 2009. HFR, Hospitalization fatality risk.

Supplementary material: File

Wong supplementary material S1

Wong supplementary material

Download Wong supplementary material S1(File)
File 153 KB