Volume 10, Number 6—June 2004
SARS Epidemiology Modeling
To the Editor: To assess the effectiveness of intervention measures during the recent severe acute respiratory syndrome (SARS) pandemic, Zhou and Yan (1) used Richards model, a logistic-type model, to fit the cumulative number of SARS cases reported daily in Singapore, Hong Kong, and Beijing. The key to using mathematical models for SARS epidemiology is understanding the models (2). In the Richards model (1), the function F(S) in the model was described as measuring “the effectiveness of intervention measures.” The parameters in F(S), namely the maximum cases load K and the exponent of deviation a, depict the actual progression of the epidemic as described by the reported data. In other words, the parameter estimates are used to quantify end results of the intervention measures implemented during the outbreaks. Simply put, the all-important question of “what if?” was not answered by their result. To gauge the effectiveness of intervention measures, one should consider a more complicated model with variable maximum case load and growth rate (r) that highlights the time-varying nature of an epidemic and its dependence on the intervention measures implemented during the epidemic.
A more subtle issue deals with the attempt to predict the trend of an epidemic from limited data during early stages of the epidemic which is often futile and sometimes misleading (3). Nevertheless, early prediction of the magnitude of an epidemic outbreak is immeasurably more important than retrospective studies. But how early is too early? Intuitively, the cumulative case curve will always be S-shaped and well-described by a logistic-type model. The essential factor is the time when the inflection of the cumulative case curve occurs, i.e., the moment when a rapid increase in case numbers is replaced by a slower increase. Since the inflection point, approximated by tm (1), dictates the point in time when the rate of increase of cumulative case numbers reaches its maximum, the moment marks the key turning point when the spread of the disease starts to decline. As long as the data include this inflection point and a time interval shortly after, the curve fitting and predicting future case number will be reasonably accurate.
To illustrate this point more precisely, the cumulative SARS case data by onset date in Taiwan were obtained from the SARS databank of Taiwan Center for Disease Control. The data cover the time from February 25, 2003, the onset date of the first confirmed SARS case, to June 15, 2003, the onset date of the last confirmed case; a total of 346 SARS cases were confirmed during the 2003 outbreak in Taiwan (4). The cumulative case data are fitted to the cumulative case function S(t) in Richards model with the initial time t0 = 0 being February 25 and the initial case number S0 = S(0) = 1. Description of the model, as well as the result of the parameter estimation, is shown in the Appendix [PDF - 24 KB - 6 pages]. The estimates for the parameters are r = 0.136 (95% confidence interval [CI] 0.121 to 0.150), K = 343.4 (95% CI 339.7 to 347.1), a = 1.07 (95% CI 0.80 to1.35), and the approximate inflection point at tm = 66.62 (95% CI 63.9 to 69.3) with adjusted r2 >0.998, p <0.0001 for the goodness-of-fit of the model (Figure). The result indicates that the infection occurred on May 3, and the estimate for the maximum case number K = 343.3 is 0.8% off the actual total case numbers.
Moreover, the case number data are sorted by onset date. Given a mean SARS incubation of 5 days (4–6 days) (5), the inflection point for SARS in Taiwan could be traced back to 5 days before May 3, namely April 28. On April 26, the first SARS patient in Taiwan died. Starting April 28, the government implemented a series of strict intervention measures, including household quarantine of all travelers from affected areas (6). In retrospect, April 28 was indeed the turning point of the SARS outbreak in Taiwan.
To address making projections during an ongoing epidemic, we used the same dataset but used various time intervals (all starting February 25) but truncated at various dates around the inflection point of May 3. The resulting parameter estimates are given in the Table of the Appendix [PDF - 24 KB - 6 pages]. For the truncated data ending on April 28 before the inflection, an unreasonable estimate of K = 875.8 was obtained. However, if we use the data ending on May 5, May 10, May 15, and May 20, we obtain estimates of K = 204.9, 253.1, 334.2, and 342.1, respectively. These estimates improve as we move further past the inflection time of May 3 (Figure). Moreover, the last estimate, using data from February 25–May 20 only, produces a 1.1% error from the eventual cumulative case number of 346, with 95% CI of 321.5 to 362.6. This retrospective exercise demonstrates that if the cumulative case data used for predictive purpose during an outbreak contain information on the inflection point and approximately 2 weeks afterwards, the estimate for the total case number can be obtained with accuracy, well before the date of the last reported case. This procedure may be immensely useful for deciding future public health policies although correctly determining the true inflection point during a real ongoing epidemic calls for scrutiny and judicious use of the model, as with all mathematical epidemic models.
- Zhou G, Yan G. Severe acute respiratory syndrome epidemic in Asia. Emerg Infect Dis. 2003;12:1608–10.
- Hsieh YH, Chen CWS. Re: Mathematical modeling of SARS: cautious in all our movements. J Epidem Com Health [serial on the Internet] 2003, [cited 2003 Nov 18]. Available from: http://jech.bmjjournals.com/cgi/eletters/57/6/DC1#66
- Razum O, Becher H, Kapaun A, Junghanss T. SARS, lay epidemiology, and fear. Lancet. 2003;361:1739–40.
- World Health Organization. Summary of probable SARS cases with onset of illness from 1 November 2002 to 31 July 2003. [monograph on the Internet].[cited 2003 Sep 26]. Available from: http://www.who.int/csr/sars/country/table2003_09_23/en/
- World Health Organization. Consensus document on the epidemiology of severe acute respiratory syndrome (SARS). [monograph on the Internet]. [cited 2003 Oct 17]. Available from: http://www.who.int/csr/sars/en/WHOconsensus.pdf
- Lee ML, Chen CJ, Su IJ, Chen KT, Yeh CC, King CC, Use of quarantine to prevent transmission of severe acute respiratory syndrome—Taiwan, 2003. MMWR Morb Mortal Wkly Rep. 2003;52:680–3.
- Figure. SARS cases, Taiwan, 2003, using Richards model; t = real data. A, confirmed cases; B, estimated cases using the truncated data.
Suggested citation for this article: Hsieh YH, Lee JY, Chang HL. SARS epidemiology modeling. Emerg Infect Dis [serial on the Internet] 2004 June [date cited]. Available from: http://wwwnc.cdc.gov/eid/article/10/6/03-1023
In Reply: Our analysis of the dynamics of reported severe acute respiratory syndrome (SARS) clinical cases was conducted in May 2003 during the height of the public panic (1). Our primary goal in that study was to predict “when the epidemic might be brought under control if the current intervention measures were continued” (1). We used the Richards model and successfully predicted the epidemic cessation dates in Beijing, Hong Kong, and Singapore. Our predicted total number of SARS cases was close to the actual number of cases. In addition, we estimated the basic reproductive rate (R0) of SARS infection, and our estimates based on the deterministic model were similar to those based on stochastic models (2,3). Therefore, our analysis provided useful information on the epidemiologic characteristic of SARS infections in three major Asian cities.
Hsieh et al. (4) commented that our article did not address the effect that specific intervention measures might have on the dynamics of SARS infection. Our study was not intended to measure this. As we stated in our article, “the transmission mechanism of the coronavirus that causes SARS and the epidemiologic determinants of spread of the virus are poorly understood.” Any models built on these unknowns are not suitable for assessing the effects of specific intervention measures. A method suggested by Hsieh et al. (4) to merely “consider a more complicated model with variable maximum case load and growth rate” will not answer the question to any extent.
The retrospective analysis of SARS case dynamics in Taiwan by Hsieh et al. (4) found that “as long as the data include this inflection point and time interval shortly after, the curve fitting and predicting future case number will be reasonably accurate.” This notion holds only if the true inflection point is known before an epidemic ends. The main difficulty is how the true inflection point is correctly determined, as noted by Hsieh et al. (4). The time when inflection occurs varies tremendously if truncated data of cumulative SARS case numbers are used. To illustrate this point, we used the cumulative number of reported probable SARS cases in Hong Kong, starting March 17, 2003, but truncated at various dates, and calculated the date when inflection occurred (Table). For example, if the data period from the onset date (March 17, 2003) to the last case reported (June 12, 2003) was used, the date when inflection would occur was estimated as March 19, 2003. If the truncated data ending April 9, April 16, April 30, May 14, and May 28, 2003, were used, the dates when inflection would occur were estimated as April 2, February 7, March 3, March 23, and April 2, 2003, respectively (Table). Clearly, inflection point dates became a moving target as the epidemic progressed. When truncated data ending April 9, April 16, April 30, May 14, and May 28, 2003, were used, the corresponding estimated the maximum numbers of cumulative cases (K) were 1,107, 1,907, 1,819, 1,749, and 1,733, respectively. Estimation of K improved when the data period used for prediction was at least 1 month past the March 19 inflection point obtained from the entire epidemic period. This analysis highlights the difficulty in identifying an optimal inflection point for prediction purposes during an ongoing epidemic when only a partial cumulative case number is available.
We fully agree with Hsieh et al. (4) that the quantitative assessment of the effectiveness of public health intervention measures for SARS is a difficult task for modelers. To make models useful for assessing the effects of specific intervention measures and for predicting the future dynamics during an ongoing epidemic, we need improved knowledge on the transmission mechanisms, pathogenesis, and the epidemiologic determinants of the spread of the virus. Any retrospective analysis of the 2003 SARS epidemic that improves our knowledge of SARS epidemiology is welcome.
- Zhou G, Yan G. Severe acute respiratory syndrome epidemic in Asia. Emerg Infect Dis. 2003;9:1608–10.
- Lipsitch M, Cohen T, Cooper B, Robins JM, Ma S, James L, Transmission dynamics and control of severe acute respiratory syndrome. Science. 2003;300:1966–70.
- Riley S, Fraser C, Donnelly CA, Ghani AC, Abu-Raddad LJ, Hedley AJ, Transmission dynamics of the etiological agent of SARS in Hong Kong: impact of public health interventions. Science. 2003;300:1961–6.
- Hsieh YH, Lee JY, Chang HL. SARS epidemiology and cumulative case curve. Emerg Infect Dis. 2004;10:1165–7.
Please use the form below to submit correspondence to the authors or contact them at the following address:
Guofa Zhou, Department of Biological Sciences, State University of New York, Buffalo NY 14260, USA; fax: 716-645-2975
Comment submitted successfully, thank you for your feedback.
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.
- Page created: June 14, 2011
- Page last updated: June 14, 2011
- Page last reviewed: June 14, 2011
- Centers for Disease Control and Prevention,
National Center for Emerging and Zoonotic Infectious Diseases (NCEZID)
Office of the Director (OD)