 Published:
We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.
A dynamic estimation of the daily cumulative cases during infectious disease surveillance: application to dengue fever
BMC Infectious Diseasesvolume 10, Article number: 136 (2010)
Abstract
Background
In infectious disease surveillance, when the laboratory confirmation of the cases is timeconsuming, there is often a time lag between the number of suspect cases and the number of confirmed cases. This study proposes a dynamic statistical model to estimate the daily number of new cases and the daily cumulative number of infected cases, which was then applied to historic dengue fever data.
Methods
The duration between the date of disease onset and date of laboratory confirmation was assumed to follow a gamma distribution or a nonparametric distribution. A conditional probability of a case being a real case among the unconfirmed cases on a given date was then calculated. This probability along with the observed confirmed cases was integrated to estimate the daily number of new cases and the cumulative number of infected cases.
Results
The distribution of the onsettoconfirmation time for the positive cases was different from that of the negative cases. The daily new cases and cumulative epidemic curves estimated by the proposed method have a lower absolute relative bias than the values estimated solely based on the available dailyconfirmed cases.
Conclusion
The proposed method provides a more accurate realtime estimation of the daily new cases and daily cumulative number of infected cases. The model makes use of the most recent "moving window" of information relative to suspect cases and dynamically updates the parameters. The proposed method will be useful for the realtime evaluation of a disease outbreak when case classification requires a timeconsuming laboratory process to identify a confirmed case.
Background
Timeliness and accuracy of case reporting are two of the most important performance measures when evaluating an infectious disease surveillance system [1–5]. Laboratory confirmation is usually needed for case diagnosis in many infectious diseases. When laboratory confirmation of the diagnosis is timeconsuming, however, there is often a timelag between the onset date of symptoms and the diagnosis date [6]. For example, the median time for confirmation of the meningococcal disease is about 13 days based on the National Notifiable Diseases Surveillance System (NNDSS) dataset for the United States from 1999 to 2001 [7]. Time from disease onset to diagnosis has been also reported to account for most of the delay in case reporting in Korea [8]. Correct estimation of daily cases or dailycumulative infected cases helps the implement of immediate disease control and allows prevention activities for infectious diseases to be instituted [6]. Using a disease surveillance system, one is able to apply statistical methods, such as cumulative sum (CuSum) [9, 10] or autoregressive integrated moving average (ARIMA) [11–15], in order to forecast an epidemic curve or to detect aberrations in disease spread. These estimations are based on having a correct daily number of cases or a dailycumulative number of cases.
An epidemic of dengue fever occurs every year in southern Taiwan. In the year 2002 in particular, there was a largescale epidemic with 5,388 confirmed cases out of totally 15,221 suspect cases [16]. This epidemic continued until March 2003. Surveillance and the control of dengue fever have become one of the most important routine areas of work at the Taiwan Centers for Disease Control (Taiwan CDC) in recent years. By 2008, Taiwan CDC had defined a confirmed case of dengue fever as an acute febrile illness together with one of the following criteria: (1) isolation of dengue virus; (2) demonstration of positive results by realtime reverse transcriptionpolymerase chain reaction (realtime RTPCR); (3) demonstration of positive seroconversion or a fourfold increase in denguespecific IgM or IgG antibody titers in paired serum samples; or (4) demonstration of hightiter denguespecific IgM and IgG antibodies in a single serum specimen [17–19]. When the dengue fever case classification only included confirmed cases by this protocol, the time needed for isolating the agent or measuring the titers for the acute and convalescentphase serum specimens was significant. The result was that there was a gap between the available daily cases or the dailycumulative cases for given a day and the actual final confirmed cases for the same day given that all diagnosis had been completed on that given day.
Assuming that a time cost for laboratory confirmation of diagnosis is sometimes inevitable, daily numbers of infected cases and dailycumulative number of infected cases may be underestimated during an epidemic. The objective of this study was to develop a new method to estimate the number of daily cumulative cases and that this method will be applied to dengue fever in Taiwan, as an example.
Methods
Data sources
Since there are almost no dengue fever cases occurred during the winter in Taiwan, we chose May 1 as the beginning of the dengue epidemic season when estimating the cumulative epidemic curve. The data come from the dengue notification dataset containing suspect cases in Taiwan whose date of onset was from May 1, 2005 to April 30, 2007. All serum samples from suspect cases were sent to the two reference laboratories at the Taiwan CDC in order to further identify if they were positive (dengue fever infected) or negative cases. The reason we retrieved data based on the date of onset rather than the report date was to avoid the influence of lag reporting on the course of the disease. All imported cases of the disease were removed. The variables we used were the date of onset, the date of laboratory confirmation (diagnosis date), and the final confirmed status of each suspect case (a binary variable that is either positive or negative). In this article, we use confirmed dengue cases and positive cases interchangeably. No personal identification information was contained in the dataset.
There were 515 positive cases out of 841 local suspect cases during the 20052006 season and 1,092 positive cases out of 2,360 local suspect cases during the 20062007 season. The median values (interquartilerange; IQR) of the onsettoconfirmation time (OCtime) for positive and negative cases were 9.5 (11) and 20 (14) days respectively for the 20052006 season. The median values (IQR) of the OCtime for positive and negative cases were 7 (9) and 18 (7) days respectively for the 20062007 season. The OCtime for positive cases was, in general, shorter than that for the negative cases. The standard deviations of the OCtime for the positive cases and negative cases were also different. Figure 1 shows that the epidemic started in the late June of 2006, and had a peak around October to November, then had nearly vanished by February 2007, based on the final status of each case.
The proposed method
Confirmation status of the suspect cases
Our proposed method estimates the realtime daily new number of cases and the daily cumulative number of dengue cases; specifically, these numbers of dengue cases are updated daily. Let c be the "current" date when the number of dengue cases is to be estimated. In this study, the date c runs from May 1, 2006 to April 30, 2007. For the i ^{th} reported suspect case counting from the 1^{st} day of the epidemic season, that is May 1 in this study, we define the suspect case's onset date as O _{ i }and the laboratory confirmation date as D _{ i }If D _{ i }> c on date c, the case i does not have a confirmation result as of date c; on the other hand, if D _{ i }≤ c, this case i is either confirmed to be a positive dengue case or has a negative result as of date c. Let the final confirmation status for the i ^{th} suspect case be, where as a positive dengue case, and as a negative case. In the situation where there are unconfirmed suspect cases as of date c, we assigned a probability of being a dengue case, P(i), to those unconfirmed cases (D ≤ c). Then for each suspect case i, the expected final confirmation status on date c, E _{ i }(c), can be written as
The values of P _{ i }(c), and E _{ i }(c)are updated for each case i every day. Without applying the proposed method, one would be only able to observe the case status from the upper part of E _{ i }(c)in equation (1). After E _{ i }(c)is calculated for each suspect case, daily new cases are easily estimated by summing the E _{ i }(c)over all new suspect cases on date c, and cumulative cases can be obtained by summing E _{ i }(c)over all cases from i = 1 to the newest suspect cases on date c.
Estimation of the probability of being dengue case, P _{ i }(c) among unconfirmed suspect cases
P _{ i }(c)is estimated for unconfirmed cases using information from the confirmed cases before date c within one year. Let T _{ i }be the onsettoconfirmation time (OCtime), the time interval between the onset date and the labconfirmation date. The OCtime for the i ^{th} suspect case as of date c, t _{ i }(c), is calculated as follows,
The t _{ i }(c) is the OCtime for confirmed cases and the censored OCtime for unconfirmed cases on date c.
By applying several steps of Bayes' rules, the probability P _{ i }(c) is given by:
To estimate P _{ i }(c) using the information as of date c, we applied the following steps. We first estimated P(Y _{ i }= 1) by calculating the proportion of confirmed positive dengue cases out of the suspect cases using the data with onset date within 1 year before the date c. Based on a parametric approach, we assumed the OCtime for a given case status, P(Ty _{ i }), follow a gamma distribution. Gamma distributions are frequently used to fit timedelay distributions or time event distributions when carrying out disease surveillance analysis [20, 21]. The probability density function of the gamma distribution is , where . The gamma distribution is denoted by with two parameters, the shape parameter α and the scale parameter β, and the mean and variance are αβ and αβ ^{2}, respectively. The values of α and β were estimated and solved by setting up the sample mean and the sample variance of the OCtime equal to αβ and αβ ^{2}, respectively. As mentioned in the previous section, the mean and standard deviation of the OCtime were different between positive and negative cases, we estimated different sets of α and β for the positive dengue cases (Y = 1) and negative cases (Y = 0) separately. A nonparametric approach was also performed in which the probability P(T > t _{ i }(c)Y _{ i }) was simply replaced with the cumulative proportion of confirmed data given their final status. Both the parametric and nonparametric models were based on the data within a 1year "moving window" before date c. The P _{ i } (c) and E _{ i } (c) were also updated everyday.
Evaluation of the proposed methods
To evaluate the performance of the proposed method, we estimated the daily new cases and daily cumulative cases for each calendar date c from May 1, 2006 to April 30, 2007. Four epidemic curves are presented. There are:

(1)
The final status curve, which is the number of dengue cases based on their final confirmation status ("Real data", "gold standard").

(2)
The daily confirmed curve, which is the number of dengue cases based on the confirmed cases as of date c.

(3)
The gammamodel curve, which is the number of dengue cases, estimated using the gamma distribution.

(4)
The nonparametricmodel curve, which is the number of dengue cases, estimated using the nonparametric distribution.
To summarize the magnitude of the bias, we defined the absolute relative bias (ARB) at date c as:
where are the cumulative cases estimated by the proposed methods or by the confirmed cases observed on date c without using the proposed methods, N _{ c }and are the cumulative confirmed cases based on the final status ("real data", "gold standard"). An ARB closer to zero is a more accurate estimate.
All analyses were performed using SAS 9.1.3 software (SAS Institute, Inc., Cary, NC). Special SAS macros for estimating the cumulative cases and daily new cases, based on our proposed model, were developed.
Results
Figure 2 compares the daily new cases estimated by the proposed models, the confirmed curve (confirmed new cases observed on date c without using the proposed models), and the final status curve (confirmed new cases based on final status; the "gold standard"). Since the daily new cases estimated by the proposed methods or the daily new cases observed on date c were different when viewed on different dates, arbitrary view dates of August 1, 2006 (beginning of the epidemic), September 1, 2006 (rising stage, before the peak), October 1, 2006 (rising stage, before but closer the peak) and November 1, 2006 (around the peak), December 1, 2006 (downward stage, after the peak), and January 1, 2007 (end of the epidemic) were chosen to illustrate the results of the estimated daily new cases. Each graph in Figure 2 shows the epidemic curves three weeks before the view dates. When viewed on August 1, most of suspect cases had been labconfirmed before July 15 thus all four curves were close to each others before that date. From July 19 to August 1, the estimated curves by the proposed methods (red dashed lines with triangle symbols by gamma distribution and blue dashed lines with cross symbols by the nonparametric method) were much closer to the final status curve (shaded area) than that by simply observing the dailyconfirmed new cases (purple dashed line). Similar patterns were observed when the results are viewed on September 1, October 1, November 1, and December 1. The observed dailyconfirmed cases usually underestimated the true daily new cases as would be expected, especially within the two weeks before the view date. The curves estimated by gamma distribution or the nonparametric approach were quite similar. However, the daily new cases, as estimated by the proposed method, did not give an accurate estimate towards the end stage of the epidemic, namely when viewed on January 1, 2007.
The epidemic curves in terms of daily cumulative cases are shown on Figure 3. In this figure, the cumulative number of positive cases was updated every day. The two estimated dailycumulative curves by the proposed methods are quite similar to the final status curve before January but again the proposed method does not work well during the end stage of the epidemic. Table 1 compares the ARB of the daily cumulative number of positive cases between the different methods. After the first confirmed positive case appeared on July 6, 2006, the estimates based on the gamma model results in an estimate closer to the real data than the other methods. For other two curves, the nonparametric method performs worst at the end of the epidemic after January 1 and there was about 20 cases higher than the final status curve. The daily confirmed curve was about 50 cases lower than the final status curve during the peak of epidemic.
Figure 4 showed the daily parameter estimates, α and β, of the gamma distributions used to dynamically calculate the daily number of positive cases. The parameter estimates varied from day to day and thus the probability of being a positive case changes. For the negative cases, the parameters had a jump during late September.
Discussion
As noted previously, timeliness and accuracy are the two of most important characteristics when we evaluate an infectious disease surveillance system. Our results show that when an infectious disease required a timeconsuming process for diagnosis, such as the dengue fever using the previously mentioned protocol, the actual daily number of infected cases and cumulative positive cases are potentially underestimated. The proposed method dynamically updates the parameters daily by making use of the most recently available information on suspect cases, and then performed estimates with a lower absolute relative bias than when using observed daily labconfirmed cases only. As shown in Table 1, the proposed method performed a lower median absolute relative bias (ABS range 1.7% ~ 8.5%) than those solely based on daily confirmed cases (ABS range 4.7% ~ 67.7%) between July 6 and December 31. These dates covered the rising stage and around the peak stage which were of public health interest. The proposed method provides a more accurate estimate of the epidemic curves when applied to the dengue fever dataset for Taiwan during the 20062007 season. Based on these results, this approach can be used for the realtime evaluation of the severity of a disease outbreak when case classification requires that a confirmed case involves a timeconsuming process.
In this study, we first established the different distributions for the onsettoconfirmation time of the positive cases and negative cases. Next, either a gamma distribution was assumed in order to estimate the probability of being a confirmed case given cases status in equation (1), or, alternatively, a nonparametric approach was used. We actually experimented with several types of distribution. The estimates using a lognormal distribution were numerically very similar to the results for the gamma distribution. The estimates using a Weibull distribution did not perform as well as the gamma distribution applied in our dengue fever data. From Figure 4, we learn that the shape parameter α changed from 0.5 to 2 and therefore an exponential distribution may not be appropriate. For simplification, we have chosen to present only the results from the gamma distribution as one example of a parametric approach and compare this with a nonparametric approach. As shown in Figure 2 for daily new cases, the differences in the estimates based on parametric approach with Gamma distribution and those with nonparametric approaches were minor. The Figure 3 and Table 1 for cumulative cases showed that a gamma distribution is a more appropriate assumption for the onsettodiagnosis time when estimating the probability of being a positive case using the dengue fever example; nonetheless, the difference between the gamma and the nonparametric method is again only slight except towards the end stage of the epidemic after January 1. The reason that the nonparametric method did not work well after January 1 is because P(Y _{ i }= 1), P(Y _{ i }= 0), and P(Ty _{ i }) had not changed substantially, resulting in a near constant estimate of the daily positive cases.
In practice, any form of the probability of being a positive case can be assumed. It is also not restricted to certain distributions when the models are adapted to different types of infectious disease. When applying this approach to other diseases, researchers should investigate several distributions according to the shape of their data and choose an appropriate one based on some appropriate measures, for instance, those shown in Table 1.
Our method estimated the probability of being a positive case based on the data within a 1year "moving window" before date c and updated P _{ i } (c) and E _{ i } (c) everyday. The epidemic profiles of dengue fever are different from one year to another in Taiwan. Choosing the data from most recent one year was done in order to insure that there was enough information to cover a whole epidemic season. In the early stage of the 20062007 season, the data from the 20052006 season actually contributed more to estimating the daily cases counts. In this study, even the epidemic profiles were not necessarily the same between the 20052006 and 20062007 seasons, the proposed methods performed well.
The study shows that before the first positive case appeared on July 6, the proposed method did not work well and are not that useful (Table 1). Our method worked well after the first positive case appeared during the 20062007 season. Indeed, it needed only four days to be able to consistently estimate the final status curve. In the 20062007 season, Taiwan CDC activated a central command center for intensively dengue epidemic control on October 2. The task of this command center included expanded blood sample collection and it is likely that this resulted in more suspect cases for laboratory confirmation, which might have led to a lower proportion of positive cases. This would influence the estimation of probability of being a positive case over the following few days. As we can see on Figure 4, it also influenced the estimation of the parameters for negative cases. While our manuscript was being prepared, the Taiwan CDC changed their laboratory protocol for dengue fever to one that requires only a single laboratory test for dengue surveillance and control. The result is a substantial reduction in the waiting time for laboratory confirmation. However, confirmation time can never be completely avoided with dengue fever. A situation where a large number of serum specimens are sent for diagnosis at the same time will result in overloading at the laboratory, which might increase the confirmation waiting time. As described previously, the estimation used information based on a "moving window" time period before the estimated date and the parameters of the model are updated everyday. Since the observed confirmed cases counts on date c are always underestimated as long as there is a time lag, our method potentially can be applied while waiting for further investigation of the status of cases.
There are some limitations to our method. Firstly, the approach needs sufficient historical data to be available in order estimate the parameters of the model; therefore our model cannot be applied effectively to an emerging disease, such as SARS or avian flu. Secondly, we used confirmed cases, the dates of onset of which were within 1year before the date estimated and if a case needs more than 1year for diagnosis such a case might never provide any information to the parameter estimation; in such a circumstance a different "moving window" needs to be chosen. Thirdly, when missing diagnosis dates exist, the estimated curve using the nonparametric method cannot converge with the final status curve. There were 134 and 289 cases missing confirmation results for the 20052006 season and the 20062007 season, respectively, at the time that the manuscript was prepared. The nonparametric method estimates by plugging in the cumulative proportion of confirmed data given the final status. As we mention before, at the end stage of epidemic, the probabilities in equation(2) almost remained unchanged. In this study, we could only assume that the proportion of positive cases out of all suspect cases among the missing observations were similar to those having results, which basically assumes that the missing data were missing at random thus ignorable.
Conclusion
When diagnosis of infectious diseases required laboratory confirmation, the time lag between onset and confirmation of a positive cases often exists and case counts are usually underestimated. This study has proposed a statistical method that more accurately estimates the realtime daily new cases and daily cumulative number of infected cases using a dengue fever epidemic as an example. The model makes use of the most recent "moving window" of information on suspect cases and dynamically updated the parameters of the assumed probability distributions. Different parametric or nonparametric distributions of the onsettoconfirmation times can be specified for different infectious diseases. The results show that, after the first confirmed case occurred, the estimated daily new cases or the cumulative case count fit the real data well compared to the daily counts based only on the available confirmed cases; this was done by assuming a gamma distribution for the onset to confirmation times and involved the use of a dynamic oneyear "moving window" of suspected cases when applied to dengue fever outbreaks in Taiwan. This method can be used for the realtime evaluation of a disease outbreak when case diagnosis requires timeconsuming laboratory process.
References
 1.
Breen E, Ghebrehewet S, Regan M, Thomson APJ: How complete and accurate is meningococcal disease notification?. Communicable Disease and Public Health. 2004, 7: 334338.
 2.
Chin J: Public health surveillance of AIDS and HIV infections. Bulletin of the World Health Organization. 1990, 68: 529536.
 3.
German RR, Lee LM, Horan JM, Milstein RL, Pertowski CA, Waller MN, GuidelinesWorking Group Centers for Disease Control and Prevention (CDC): Updated guidelines for evaluating public health surveillance systems: Recommendations from the Guidelines Working Group. MMWR. 2001, 50: i31.
 4.
Paulozzi , Mercy J, Frazier L, Annest JL: CDC's National Violent Death Reporting System: background and methodology. Injury Prevention. 2004, 10: 4752. 10.1136/ip.2003.003434.
 5.
Sejvar J, Tangkanakul W, Ratanasang P, Dowell SF, Sangjun N, Bragg S, Ashford D, Tappero J: An outbreak of leptospirosis, Thailandthe importance of the laboratory. Southeast Asian Journal of Tropical Medicine and Public Health. 2005, 36: 289295.
 6.
Jajosky RA, Groseclose SL: Evaluation of reporting timeliness of public health surveillance systems for infectious diseases. BMC Public Health. 2004, 4: 2910.1186/14712458429.
 7.
Yan P, Chen H, Zeng D: Syndromic Surveillance Systems: Public Health and Biodefense. Annual Review of Information Sciences and Technology (ARIST). 2008, 42:
 8.
Yoo HS, Park O, Park HK, Lee EG, Jeong EK, Lee JK, Cho SI: Timeliness of national notifiable diseases surveillance system in Korea: a crosssectional study. BMC Public Health. 2009, 9: 9310.1186/14712458993.
 9.
Cowling BJ, Wong IOL, Ho LM, Riley S, Leung GM: Methods for monitoring influenza surveillance data. Int J Epidemiol. 2006, 35: 13141321. 10.1093/ije/dyl162.
 10.
Gallus G, Chen R: The relative efficiency of the sets and the cusum techniques in monitoring the occurrence of a rare event. Statistics in Medicine. 1988, 7: 997999. 10.1002/sim.4780070911.
 11.
Box G, Jenkins G: Time series analysis forecasting and control. 1976, San Francisco: Holden Day
 12.
Mitsi C, Roch R, Jean L: On the application of integervalued time series models for the analysis of disease incidence. Statistics in Medicine. 1999, 18: 20252039. 10.1002/(SICI)10970258(19990815)18:15<2025::AIDSIM163>3.0.CO;2D.
 13.
Allard R: Use of timeseries analysis in infectious disease surveillance. Bull World Health Organ. 1998, 76: 327333.
 14.
Nobre FF, Monteiro AB, Telles PR, Williamson GD: Dynamic linear model and SARIMA: a comparison of their forecasting performance in epidemiology. Statistics in Medicine. 2001, 20: 30513069. 10.1002/sim.963.
 15.
Williamson GD, Ginner Weatherby H: A monitoring system for detecting aberrations in public health surveillance reports. Statistics in Medicine. 1999, 18: 32833298. 10.1002/(SICI)10970258(19991215)18:23<3283::AIDSIM316>3.0.CO;2Z.
 16.
Number of Reported and Confirmed Cases  By Month. [http://www.cdc.gov.tw/public/Data/7121917451971.pdf]
 17.
World Health Organization: Dengue haemorrhagic fever: diagnosis, treatment and control. 1997, Geneva, WHO, 2
 18.
Huang JH, Liao TL, Chang SF, Su CL, Chien LJ, Kuo YC, Yang CF, Lin CC, Shu PY: Laboratorybased Dengue Surveillance in Taiwan, 2005: A Molecular Epidemiologic Study. Am J Trop Med Hyg. 2007, 77: 903909.
 19.
Shu PY, Chang SF, Yueh YY, Chow L, Chien LJ, Kuo YC, Su CL, Liao TL, Lin TH, Huang JH: Current status of dengue diagnosis at the Center for Disease Control. Dengue Bull World Health Organ/SEARO. 2004, 28: 107117.
 20.
Donnelly CA, Ghani AC, Leung GM, Hedley AJ, Fraser C, AbuRiley S, Raddad LJ, Ho LM, Thach TQ, Chau P, et al: Epidemiological determinants of spread of causal agent of severe acute respiratory syndrome in Hong Kong. The Lancet. 2003, 361: 17611766. 10.1016/S01406736(03)134101.
 21.
Farewell VT, Herzberg AM, James KW, Ho LM, Leung GM: SARS incubation and quarantine times: when is an exposed individual known to be disease free?. Statistics in Medicine. 2005, 24: 34313445. 10.1002/sim.2206.
Prepublication history
The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712334/10/136/prepub
Acknowledgements
This study was supported in part by the Centers for Disease Control, Department of Health, Taiwan, Republic of China (DOH94DC2036 and DOH96DC2038) and was also supported by the Aim for the Top University Plan of National YangMing University. The authors would like to thank the reviewers' excellent comments for improving this manuscript.
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
IFL designed and conducted the study, and finalized the manuscript. PHC participated in the design of the study, performed the statistical analyses, and drafted the manuscript. JHC helped conceive the study, participated in the data collection, gave input to the manuscript, and provided medical advice from the public health perspective. All authors have read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Received
Accepted
Published
DOI
Keywords
 Gamma Distribution
 Positive Case
 Dengue Fever
 Infected Case
 Laboratory Confirmation