key: cord-0953363-haszvpis
authors: Hsieh, Ying-Hen; Lee, Jen-Yu; Chang, Hsiao-Ling
title: SARS Epidemiology Modeling
date: 2004-06-03
journal: Emerg Infect Dis
DOI: 10.3201/eid1006.031023
sha: 03e352c35635f746ff112e5a766eefe10efb59b4
doc_id: 953363
cord_uid: haszvpis

nan

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 mathemati-cal 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.

Predicting the trend of an epidemic from limited data during early stages of the epidemic 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 t m (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 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. Emerging Infectious Diseases • www.cdc.gov/eid • Vol. 10, No. 6, June 2004 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 (R 0 ) 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 epidemiological 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) (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 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 one 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.

*State University of New York, Buffalo, New York, USA

To the Editor: A novel coronavirus caused more than 8,000 proba-ble cases of severe acute respiratory syndrome (SARS) worldwide (1, 2) during the 2003 outbreak. Before the etiologic agent was identified, the diagnosis of SARS was made according to a set of clinical-epidemiologic criteria as suggested by the Centers for Disease Control and Prevention (CDC) (1-3) . These criteria remained important in the initial diagnosis and prompt isolation of patients because the overall sensitivity of initial reverse transcriptase-polymerase chain reaction (RT-PCR) testing for SARSassociated coronavirus (SARS CoV) RNA on upper respiratory specimens ranged from approximately 60% to 70% (though sensitivity improved with a second test) (4, 5) . In a SARS screening clinic at the Prince of Wales emergency department, the positive predictive value (PPV) of these criteria was estimated to be 54% (95% CI 39% to 69%) (6) . The relative importance of the clinical versus epidemiologic criteria had not been evaluated. By using paired serologic testing to determine SARS-CoV infection (3), we evaluated the relative importance of the clinical-epidemiologic diagnostic criteria during an outbreak.

Patients with a diagnosis of SARS, and who were admitted to one of five regional hospitals in Hong Kong for isolation and treatment from March 4 to June 6, 2003, were included in this retrospective analysis. Probable SARS case-patients were those who met the CDC clinical criteria for severe respiratory illness of unknown etiology (3), and met the epidemiologic criterion for exposure in either a close or a possible contact. Close contact was defined as caring for, living with, or having direct contact with body fluids of a probable SARS patient (e.g., working in the same medical ward or staying in the same household) within 10 days of initial symptoms. Because Hong Kong was the documented SARS transmission site from February 1 to July 11, 2003, a modified epidemiologic criterion of possible contact was adopted. Possible contact was defined as staying or working in the same hospital compound, or residing in the same building where case clusters of SARS had been reported, within 10 days of symptoms onset.

Laboratory testing of paired immunoglobulin (Ig) G antibody to SARS-CoV was used to determine infection (7). Positive serologic evidence of infection was defined as a four-fold rise in antibody titer or detection of antibody in convalescentphase serum. Seronegativity was defined as absence of antibody in convalescent-phase serum obtained >21 days after symptom onset (3) . Seronegativity in this defined time frame (>21 days -serum collected before July 11, 2003, and beyond 28 days) excluded the diagnosis of SARS (3). Samples from patients showing nonspecific fluorescent signals were considered negative for SARS-CoV infection. RT-PCR was performed on clinical specimens (respiratory, fecal) from all patients (1, (3) (4) (5) .

Demographic and laboratory parameters and history of close contact were compared between the seropositive and seronegative groups. Student t test was used to analyze continuous variables. A p value of <0.05 was considered statistically significant. Odds ratio (OR) and 95% confidence interval (CI) were calculated for categorical variables.

During the study period, 475 patients were hospitalized with probable SARS. One hundred patients were excluded because their serologic results were either missing (n = 37) or they died before day 21 of illness (no convalescent-phase serum, n = 63). Three hundred seventy-five patients were included in the analyses; 353 (94.1%) patients were serology-positive for SARS-CoV. Two hundred sixty-three of the 353 patients (74.5%) had a 4-fold increase in antibody titers, and 90 of the 353 patients

Age and variant Cretzfeldt-Jakob disease

Comparative epidemiology of scrapie outbreaks in individual sheep flocks

Severe acute respiratory syndrome epidemic in Asia

Re: Mathematical modeling of SARS: cautious in all our movements

SARS, lay epidemiology, and fear

World Health Organization. Summary of probable SARS cases with onset of illness from 1

World Health Organization. Consensus document on the epidemiology of severe acute respiratory syndrome (SARS)

Use of quarantine to prevent transmission of severe acute respiratory syndrome-Taiwan

Severe acute respiratory syndrome epidemic in Asia

Transmission dynamics and control of severe acute respiratory syndrome

Transmission dynamics of the etiological agent of SARS in Hong Kong: impact of public health interventions

SARS epidemiology and cumulative case curve