 Research article
 Open Access
 Published:
Time series analysis of demographic and temporal trends of tuberculosis in Singapore
BMC Public Health volume 14, Article number: 1121 (2014)
Abstract
Background
Singapore is an intermediate tuberculosis (TB) incidence country, with a recent rise in TB incidence from 2008, after a fall in incidence since 1998. This study identified population characteristics that were associated with the recent increase in TB cases, and built a predictive model of TB risk in Singapore.
Methods
Retrospective time series analysis was used to study TB notification data collected from 1995 to 2011 from the Singapore Tuberculosis Elimination Program (STEP) registry. A predictive model was developed based on the data collected from 1995 to 2010 and validated using the data collected in 2011.
Results
There was a significant difference in demographic characteristics between resident and nonresident TB cases. TB risk was higher in nonresidents than in residents throughout the period. We found no significant association between demographic and macroeconomic factors and annual incidence of TB with or without adjusting for the populationatrisk. Despite growing nonresident population, there was a significant decrease in the nonresident TB risk (p < 0.0001). However, there was no evidence of trend in the resident TB risk over this time period, though differences between different demographic groups were apparent with ethnic minorities experiencing higher incidence rates.
Conclusion
The study found that despite an increasing size of nonresident population, TB risk among nonresidents was decreasing at a rate of about 3% per year. There was an apparent seasonality in the TB reporting.
Background
Tuberculosis (TB) is a leading cause of death globally [1]. Approximately 5% of those infected progress to active disease within the first two years, but the majority experience a state of prolonged latency with a lifetime risk of 1015% eventually progressing to active disease [2]. The risk of progression is affected by immune status and is much higher among immunecompromised individuals [2]. The situation is complicated further by the spread of multidrugresistant (MDR) TB [3], a particular problem in South East Asia.
Since the World Health Organization (WHO) declared TB to be a global health emergency in 1993 [4], there has been a substantial progress on many indicators [3]. The most significant achievement was that global TB mortality reduced by 45% in 2012 [3]. However, the rate of decrease of global TB incidence (2% per annum) remains quite low [3]. In this context, the study of TB prevalence in Singapore is important for the following two reasons. As an economically developed city state (GDP per capita $52,000) [5], Singapore, lies on an important maritime transit route, is a transit hub for travel in Asia and AsiaPacific and a significant global destination of investment and labour. However, it is also a part of the region that accounts for 29% of global TB incidence [3]. Contrasting the global trends of most developed nations, Singapore has witnessed a resurgence of TB incidence since 2008.
TB incidence rates among Singapore residents declined from 300 (per 100,000 population) in the early 1960s to 106 (per 100,000 population) in the mid1980s [6, 7]. In 1997, the Singapore Tuberculosis Elimination Programme (STEP) was launched to address a decadelong stagnation in TB rates at 50–55 (per 100,000 population) between 1987 to 1997. Thereafter, TB incidence rates among Singapore residents declined from 58 (per 100,000 population) in 1998 to an historic low of 35 (per 100,000 population) in 2006 and 2007. However, the TB incidence rates began to rise around 39 per 100,000 population since 2008. The reasons for this surge remain unclear since there has been no association study on the influence of the demographic characteristics on TB in Singapore.
In this research, our aim was to assess temporal dynamics in TB incidence in Singapore at both short (monthly) and large (annual) scales, and to investigate the association of several population level demographic factors with annual TB incidence. For the annual analysis. the primary methodology we used was the standard autoregressive integrated moving average (ARIMA) regression model [8]. Additionally, using monthly data from 1995–2011, we assessed whether the burgeoning nonresident population led to an increase in TB risk in Singapore. We proposed using the seasonal autoregressive moving average (SARIMA), ARIMA models with periodic components, to predict the temporal trends of the more volatile monthly TB risk among residents and nonresidents in Singapore and detect seasonality.
Methods
Study setting
Singapore is a citystate in the middle of the Malay archipelago approximately 100 Km from the equator. Its population of 5.3 million (2013) has grown rapidly over the last decade due to the influx of foreigners who now constitute approximately 28% of the total population, comprising 3.8 million residents (citizens and permanent residents) and 1.5 million nonresidents (i.e. those who hold longterm visit or work passes) [5]. The proportion of the latter group has increased markedly since 2005, with most migrants coming from South and East Asian countries.
Data
This study retrospectively analyzed all newly laboratoryconfirmed TB cases reported to Singapore Tuberculosis Elimination Program (STEP) registry from 1995 to 2011. Laboratory confirmation was based on positive M. tuberculosis culture. Under the Infectious Disease Act in Singapore, all suspected and confirmed TB cases are to be notified to the STEP registry of the TB Control Unit (TBCU) within 72 hours of starting TB treatment and/or laboratory confirmed results. The TBCU is the national unit for treatment of TB patients, contact investigation, preventive therapy, and educational training of health care workers in Singapore.
Population data were obtained from Singapore Department of Statistics and Ministry of Health [9, 10].
Statistics and analysis
Annual TB cases analysis: 1995–2011
We fitted time series (ARIMA) regression models for the yearly time series counts of TB cases to establish whether a relationship exists between population characteristics and yearly TB cases over time. We used the standard Box Jenkins approach [11] to fit the ARIMA models, characterized by 3 main parameters: an autoregressive (AR) term, indicating the strength of relationship between successive years’ incidence, a moving average (MA) term, that checks for the dependence of yearly incidence on current and past model residuals, and a differencing term, that is usually applied to make the data stationary, when the time series has a long term trend [11]. Further as the annual TB time series exhibited change of variance over the time horizon, yearly notifications were log transformed, which had the effect of standardizing the scale of the error terms.
Following stand and procedure, sample autocorrelation (ACF) and partial autocorrelation function (PACF) plots were used to identify the order of the MA and AR terms included in the ARIMA model [12, 13]. Specifically, if we observe statistically significant autocorrelation in the sample ACF plot at q time lags and significant sample PACF at p lags then, an ARIMA model with pAR components and qMA components is proposed. The final model (order) selection is made using a combination of goodness of fit and Akaike’s Information Criterion (AIC). Here we the determined the goodness of fit using the residual autocorrelation function to ensure that no additional autocorrelation was present in each fitted model, and models were compared using Akaike’s Information Criterion (AIC). The model with the lowest value of the AIC was selected to analyze yearly TB cases. We developed separate univariate models, each encompassing a different lag period controlling the effect of potential predictor (from lag 0, that is the current year, to the immediate preceding 4 years). Population variables included numbers of visitors’ arrivals, the sizes of the elderly population, the nonresident population, and the resident population broken down by ethnicity (Chinese, Indian, Malay and Other ethnicity), Gross Domestic Product (GDP) rates, population density and HIV notification rates. Population density was the midyear population estimate divided by the country’s area (square kilometers) which grew during the study period due to land reclamation. We also adjusted for the populationatrisk (resident or nonresident population) as the denominator for demographic variables that did not already account for them, such as GDP, HIV notifications and visitor’s arrivals.
Midyear population estimates from the Department of Statistics were used to generate crude incidence and agestandardized incidence rates [5]. In the direct method of standardization, all agestandardized incidence rates were derived by applying the categoryspecific incidence rates of each population to the standard world population [14].
Monthly TB risk analysis: 1996–2011
The limitation of the above annual analysis is that the available length of annual time series is small and so the analysis is underpowered we therefore supplement it with a shorter, monthly time step. However we have sufficient monthly temporal observations. In this section we describe the method used to investigate monthly temporal risk of TB.
The TB risk among resident or nonresident population was defined as Y_{t} = (c_{t}/n_{t}) × 100 000, where c_{t} denotes the monthly count of TB reported cases (residents or nonresidents) and n_{t} is assumed to be populationatrisk (number of resident or nonresident population) in month t. This definition takes into account the dynamics of the population size in explaining TB incidence. Population figures were derived from annual population data [5].
To investigate the monthly time series of TB risk we use classical time series decomposition [11] that is conceptually,
In our models we have used linear and cosine components to model the trend and cyclic components while seasonality and irregular components are modelled using SARIMA models (see [12]).
To account for the observed periodicity in the time series, Y_{t}, we introduced a combination of cosine components. This is a well known approach for addressing deterministic periodic effects in a time series data (see for example [Ch 5, 30] [15]) , and Wiener and Khintchine [13] have shown, in a seminal result, that a stationary time series can be completely written as a combination of sine and cosine components and error. From an application point of view this technique of modelling the trend and cyclical components can be perceived as a regression of the data Y_{t} on periodic covariates, which are functions of time, while the stochastic dependence of the data, on its past including seasonality, is modelled by fitting SARIMA models to the residuals of this regression.
Plots of incidence over time showed volatility in the yearly and monthly TBreporting indicating a possible heteroskedasticity of variance. We converted the raw data to its natural logarithm to stabilize the variance of the time series and improve the estimation.
Explicit formulae for modelling the trend of resident and nonresident monthly TB risk are given below. We use the subscripts “R” and “NR” to indicate model components for residents and nonresidents respectively.
Resident model
Nonresident model
In each model β _{ 0 } is the intercept term and β _{ 1 } is the long term linear rate of change of TB risk per month. The principal reason for proposing the linear temporal component is to assess if monthly TB risk overall follows an increasing or decreasing trend, β _{ k } are the coefficients of the k ^{th} time period T _{ k }. All periods T_{k} are greater than 1 year. Based on sample autocorrelation and periodogram analysis of log(Y _{ t } ) we proposed periods, T _{ k }. The corresponding coefficients β are then estimated and tested for statistical significance, using the method of least squares. Specific details on k and the coefficients β for each model are given in the Results section and Table 1. The ‘residual’ correlated noise, e _{ t } , were modelled using SARIMA models. Details of SARIMA models are given in Table 2. Initial investigation suggested that the models were both appropriate and parsimonious.
The long latency of TB is frequently observed as a major challenge in constructing epidemiological models for TB [16]. Therefore, we fitted similar linear models for TB risk to adjust for different latent periods of lag 0–4 years by changing the susceptible population n_{t}, denominator of Y _{ t }, to be the population size in the appropriately lagged year. For example, to account for a latent period of 1 year, we divided TB count for a particular time point by the population at risk during previous year.
The prediction models were developed using monthly reported TB data between 1996 and 2010 and validated by predicting the data for 2011.
The goodness of fit was examined using AIC, Mean Squared Prediction Error (MSPE) and residual autocorrelation function. A lower AIC value and MSPE indicate a better fit of the model. The bestfitting model was used to build the subsequent predictive model. The prediction models were fit using the generic function auto.arima( ) in the forecast package of R [17]. Approval was obtained from the institutional review board (DSRB 2009/00232) to conduct this study.
Results
There were 40,046 laboratoryconfirmed TB cases reported to the STEP registry between 1995 and 2011. Table 3 shows the demographic characteristics of TB cases by residential status. Among both residents and nonresidents, most TB cases were aged 15–64 years old. Compared to resident TB cases, nonresident cases were more likely to be younger, belong to other (i.e. not belonging to any of the three main ethnic groups in Singapore, Chinese, Indian and Malay) ethnic groups and to be in prison. As residents and nonresidents have substantially different demographic features and health issues, we did separate time series analyses for the two groups. Figure 1 shows the increasing trend of resident and nonresident TB cases corresponding to increasing population size and density. The proportion of nonresident TB cases contributing to the total case burden increased from 25.5% in 1995 and 28.9% in 2004 to 47.7% in 2011.
Predictors of annual TB cases
On the basis of parameter estimation and goodness of fit test statistics, plausible ARIMA models were used for time series regression analysis. For residents, an ARIMA (2, 0, 1) model was selected whereas for nonresidents, an ARIMA (1, 1, 1) was preferred. We found no significant association between demographic factors and annual incidence of TB with or without adjustment of populationatrisk.
Annual TB risk analysis
TB incidence rates were consistently higher among nonresidents than residents over the period studied (Figure 2a). Among residents, the TB incidence rate was higher among those aged ≥65 years and who belong to Malay and Other ethnic groups (Figure 2b,c). Agestandardized incidence rates (ASIR) among residents were significantly lower in 2011 (ASIR: 36.5 per 100,000 population, 95% CI = 34.838.3 per 100,000 population) by 24% compared with 1995 (ASIR: 48.5 per 100,000 population, 95% CI = 4650.9 per 100,000 population) (Figure 2d). Despite a drop between 2004 and 2007, the ASIR has been increasing since 2008. Compared with the 2007 ASIR (29.9 per 100,000 population, 95% CI = 28.331.6 per 100,000 population), there was a 22% increase by 2011 (36.5 per 100,000 population, 95% CI = 34.838.3 per 100,000 population). When stratified by ethnic groups among residents, Malay and Other ethnicities had significantly higher ASIR over years (Additional file 1: Figure S1).
Monthly TB risk prediction
We noted that the TB incidence rates continued to be higher among nonresidents than residents. The annual time series analysis also indicated that TB incidence among nonresidents continued to grow with an increasing population size. However, from a public health perspective, it is important to assess the dynamics of TB incidence relative to the growth in population size (TB risk).
This study found that TB risk among the nonresident population was significantly linearly decreasing at a rate of 3% per year (monthly linear estimate = −0.3%, 95% CI = −0.3%, 1.9%, p < 0.0001) (Table 1).
However, we found no increasing or decreasing linear trend in the resident TB risk. Both resident and nonresident TB risk time series had significant longterm deterministic periodicities. We modelled these patterns using standard trigonometric components (detailed in Table 1). For residents we observe that the monthly TB risk time series has periodicities of approximately 18 months, 4 years, 8 years. While for the nonresident group we observe that periodicity coincides with the length of the time series. Additional file 2: Figure S2 shows the fitted curves superimposed on the resident and nonresident TB risk. From a purely epidemiology point of view, we would not wish to ascribe too much importance to the particular values of periodicity beyond its contribution to goodnessoffit of the model. However the statistical significance of the periodic components underscore that monthly risk (and hence the monthly reported cases) have multiple long term periodicities and further research is needed to understand if there are socioeconomic factors influencing these periodic effects.
The linear models with suspected latent periods of 0–4 years among residents and nonresidents were found to be similar. In particular, for each lagged model, the rate of change of TB risk, Y_{t}, continued to be statistically significantly negative among all lags of latent period. Due to the similarity of the above lagged linear models, we reported the analyses with lag 0 and 4 (Table 1). Next, we investigated the autocorrelation among monthly residuals of the above linear models. The sample autocorrelation plots for the residuals of linear models for both resident and nonresident TB risk show a slow decaying periodic nature with significant autocorrelation at lags of multiples of 12 months (Additional file 3: Figure S3 and Additional file 4: Figure S4). On average, for both Residents and NonResidents the peak of TB risk is observed in the month of July with relatively more cases are reported in March, July and October, though we would like to point out that this seasonal difference is not substantial. These observations in a time series warrant the use of SARIMA models [13].
Based on AIC and MSPE values, we selected SARIMA models to be the best fit for monthly TB risk of residents and nonresidents (Table 2).
Additional file 5: Figures S5 and Additional file 6: Figure S6 show the sample autocorrelation and partial autocorrelation of residuals of the final SARIMA models. Sample autocorrelation of residuals are within the bounds of 95% confidence interval with no significant autocorrelation between residuals at different lag times, and the models passed the BoxLjung test (p = 0.63) on autocorrelation of the final residuals and KolmogorovSmirnov test for normality (p = 0.58). Thus, we selected an additive combination of a deterministic linear model and a stochastic SARIMA for resident and nonresident TB risk. This analysis detected a 12month stochastic seasonality and a long term deterministic periodicity in the reporting patterns of resident and nonresident TB risk.
Note
Different models and modeling techniques can be used for making inference on the same time series data. The periodic nature of sample autocorrelation function plots for residuals of the trend adjusted time series of TB risk for residents as well as nonresidents depicted in Additional file 3: Figure S3 and Additional file 4: Figure S4 indicate that one could also use Generalized Autoregressive Moving Average models [18].
Discussion
Although the TB incidence rate in Singapore is lower than most other Asian countries, it is higher than most other developed countries [3]. We observed that TB incidence rates were higher among the elderly, Malay, the ‘Other’ ethnic group (a mix of other Asian ethnicities, Europeans and those of mixed racial heritage) and nonresidents. In Singapore, the resident elderly population had almost doubled by 2011 compared to 1995 [9] and the higher TB incidence in the elderly group could partly explain the increase in TB incidence in the resident population since 2008. Among all resident ethnic groups, the Malay group had the highest increase in ASIR. A previous study showed that this ethnic group in Singapore has clinical and socioeconomic vulnerability to TB infection and transmission [19].
Although HIV notification rates have been rising in Singapore, the increasing trend of HIV did not have a significant influence on yearly TB cases. This is possibly because the number of HIV cases is low especially in comparison to the overall TB burden in Singapore [20]. In addition, nonresidents applying for work permits in Singapore are screened for HIV on arrival and are deported if they subsequently are diagnosed to be infected, limiting the impact of HIV on TB in the nonresident population.
Since around 2005, liberalization of immigration policy resulted in a marked increase in population and economic growth with a rapid influx of foreignborn immigrants in Singapore. The nonresident population of Singapore has grown to about 37% of the country’s total labor force in 2012 [21]. Majority of immigrants come from high incidence TB countries such as India and China. An increasing trend of nonresident TB cases contributing to the overall proportion of TB cases over the years could suggest that mass immigration from high TB incidence countries is increasingly contributing to burden of TB in recent years in Singapore, similar to other lowincidence countries such as the US, Canada and the UK [22]. Given plausible interactions between residents and nonresidents, TB infections in the nonresident population might have led to a subsequent transmission to residents, and a recent DNA fingerprinting study in Singapore also suggested a crosstransmission between residents and nonresidents [23]. Consistent with this hypothesis, the increasing trend of TB incidence among residents since 2008 seemed to follow a similar rise among nonresidents after 2005. However, our study did not find a direct influence of increasing number of nonresident population on resident TB cases. Importantly, when adjusted for populationatrisk, the TB risk among nonresident population was found to be decreasing at a rate of about 3% per year despite the burgeoning size of the nonresident population.
We found both short term seasonal and long term periodic variations in TB reporting. Seasonal variation has been reported in the countries with distinct climatic seasons [24–28], but its presence in Singapore is curious as this city state is located very close to equator and has little climatic variability across the year. This study highlights a need for further studies to identify the factors explaining the seasonality of TB reporting patterns in Singapore. This study has shown that SARIMA model could be useful in the shortterm prediction of TB in Singapore.
Under the implementation of National Tuberculosis Program in Singapore, nonresidents can also access free DOTS treatment since 2010 in Singapore [29]. This could be a plausible factor behind a decreasing trend of TB risk among nonresidents. Additionally, it could also be explained due to a much greater influx of employment pass holders (professionals and white collar workers) than work permit holders (manual laborers and domestic workers) in the recent years [30]. Findings from this study may permit better understanding of the underlying timeseries dynamics of TB risk among different demographic groups and review TB control measures at the national scale in Singapore. In ongoing research, we have identified that TB even within a city state of Singapore is a heterogeneous multilevel event influenced by several demographic, socioeconomic, spatiotemporal and environmental factors. Thus, strategies for intervention need to be more specific for any demographic group and not based solely on overall incidence.
A clear strength of the current study is the completeness and consistency of the data: Singapore, as a city that is also a country, has a clearly demarcated borderline and population catchment, while legislation means that any case of TB has to be notified by the attending clinician to the ministry of health. Over the course of the study period, all such notifications have been shared with the national TB registry that has collated these data consistently. Another strength is the rapidness of changes to the population, which potentially would allow large effects of demographics to be seen.
However there are also limitations. These could be broadly classified as either systemic or epidemiological. A major limitation is that the study period is relatively short for annual analysis. The implication of this is that there is limited power to detect smaller demographic effects, and especially to identify independent effects of demographic changes at an annual level, given that many indicators have risen in tandem with each other. The monthly analysis is an attempt to identify effects at other temporal scales. While this limits the ability to relate TB incidence to longer scale changes, such as demographic indices for which only yearly data are available, it does permit estimation of secular changes and withinyear changes, modelled as linear and sinusoidal terms, respectively. Use of monthly data naturally increases the data size and makes a more powerful analysis than the yearly analysis by itself.
First, the ecological nature of the study limited the ability to draw conclusions at the individual level due to ecological bias. A further difficulty is the complicated epidemiology of TB, with a mixture of immediate onset of symptoms in some cases and a long latent period in others. This means that any demographic changes with a genuine effect on TB infection may have their influence dispersed over time, further lowering the power. Another weakness is the absence of socioeconomic and demographic information on nonresidents, beyond their total numbers, as these are not published by the government. Additionally, nonresident TB cases considerd in this study included both longstay (dependent pass, work pass, student pass and long term social visit pass) and shortstay (shortterm social visit pass with a maximum of 89 days) nonresidents. We expected shortstay nonresidents might have contributed to the overall nonresident TB cases to a small degree and exclusion of these cases might not have affected our results to a large extent. Lastly, as we could only obtain reliable data from 1995, we could not reliably assess the impact of STEP (which was formed in 1997) itself. Since the data are from a passive surveillance system, the possible biases in disease reporting and potential underreporting of TB cases might influence the precision of the analysis. However, laboratoryconfirmed tuberculosis is notifiable by law in Singapore, and is moreover notified by both the laboratory and the physicianincharge. There was also no major difference in data collection and extraction procedures used throughout the study period, hence it is very unlikely that there were significant errors. Nonetheless, a significant proportion of TB cases are culturenegative, and discarding these cases may have affected the analyses. Additionally, careful assessment of the accuracy and comprehensiveness is required when interpreting the result of the analysis since TB notifications may not be considered as a direct measure of TB incidence [31]. Because of the long latency of TB, a longer period of time is required to determine the impact of increasing nonresident TB incidence on the resident population. Future molecular epidemiological studies are needed to differentiate recent infection or reactivation of TB between resident and nonresident population over time.
Conclusion
The study found that despite an increasing size of nonresident population, TB risk among nonresidents was decreasing at a rate of about 3% per year. There was an apparent seasonality in the TB reporting patterns and SARIMA model could be useful in predicting the shortterm trend of TB risk.
Abbreviations
 TB:

Tuberculosis
 ARIMA:

Autoregressive integrated moving average
 STEP:

Singapore tuberculosis elimination program
 TBCU:

TB control unit
 GDP:

Gross domestic product
 ACF:

Autocorrelation function
 PACF:

Partial autocorrelation function
 AIC:

Akaike information criterion
 DOTS:

Directly observed treatment short course
 MSPE:

Mean square prediction error
 SARIMA:

Seasonal autoregressive integrated moving average
 ASIR:

Agestandardized incidence rates
 AR:

Autoregressive
 MA:

Moving average
 D:

Differencing.
References
 1.
Tuberculosis; Fact Sheet. [http://www.who.int/mediacentre/factsheets/fs104/en/]
 2.
Maher D: The natural history of Mycobacterium tuberculosis infection in adults. Tuberculosis: A comprehensive clinical reference. Edited by: Schaaf HS, Zumla A. 2009, UK: Elsevier Health Sciences, 129132.
 3.
Global tuberculosis report 2013. [http://apps.who.int/iris/bitstream/10665/91355/1/9789241564656_eng.pdf]
 4.
Tuberculosis Fact Sheet. [http://www.who.int/mediacentre/factsheets/who104/en/print.html]
 5.
Key Annual Indicators. [http://www.singstat.gov.sg]
 6.
Chee C, Wang YT: TB control in Singapore: where do we go from here?. Singapore Med J. 2012, 53: 236
 7.
Communicable Disease Surveillance Report in Singapore 2010. [http://www.moh.gov.sg/content/dam/moh_web/Publications/Reports/2012/Communicable%20Diseases%20Surveillance%20in%20Singapore%202011/Full%20Version%20of%20the%20Report%202011.pdf]
 8.
Permanasari AE, Rambli DRA, Dominic PDD: Performance of univariate forecasting on seasonal diseases: the case of tuberculosis. Software Tools and Algorithms for Biological Systems. 2011, New York: Springer, 171179.
 9.
Statistics; Latest Data. [http://www.singstat.gov.sg/statistics/latest_data.html#12]
 10.
Communicable Diseases Surveillance in Singapore 2011. [http://www.moh.gov.sg/content/moh_web/home/Publications/Reports/2012/_communicable_diseasessurveillanceinsingapore2011.html]
 11.
Box GEP GMJ: Time series analysis: forecasting and control. 1976, San Francisco: Holden Day, 181218.
 12.
Brockwell P, Davis R: Time Series: Theory and Methods. Springer Series in Statistiscs. 1991, New York: Springer Verlag
 13.
Shumway RH, Stoffer DS: Time Series Analysis and its applications. 2011, Springer Texts in Statistics
 14.
Doll R, Payne P, Waterhouse J: Cancer Incidence in Five Continents: A Technical Report. Book Cancer Incidence in Five Continents: A Technical Report. 1966, City: Springer Verlag (for UICC)
 15.
Davis RA, Dunsmuir WTM, Streett SB: Observation driven models for poisson counts. Biometrika. 2003, 90 (4): 777790. 10.1093/biomet/90.4.777.
 16.
Magombedze G, Dowdy D, Mulder N: Latent tuberculosis: models, computational efforts and the pathogen's regulatory mechanisms during dormancy. Front Bioeng Biotechnol. 2013, 1:
 17.
Hyndman R: R package version 5.4 forecast: Forecasting functions for time series and linear models. [http://cran.rproject.org/web/packages/forecast/index.html]
 18.
Benjamin MA, Rigby RA, Stasinopoulos DM: Generalized autoregressive moving average model. J Am Stat Assoc. 2003, 98 (461):
 19.
Lim L, Enarson D, Reid A, Satyanarayana S, Cutter J, Kyi Win K, Chee C, Wang Y: Notified tuberculosis among Singapore residents by ethnicity, 2002–2011. Public Health Action. 2013, 3: 311316. 10.5588/pha.13.0055.
 20.
Update on the HIV/AIDS situation in Singapore. 2013, [http://www.moh.gov.sg/content/moh_web/home/statistics/infectiousDiseasesStatistics/HIV_Stats/updateonthehivaidssituationinsingapore20130.html]
 21.
Labour Force. 2013, [http://stats.mom.gov.sg/Pages/LabourForceinSingapore2013LabourForce.aspx]
 22.
MacPherson D, Gushulak B: Balancing prevention and screening among international migrants with tuberculosis: population mobility as the major epidemiological influence in lowincidence nations. Public Health. 2006, 120: 712723. 10.1016/j.puhe.2006.05.002.
 23.
Lim LK, Sng LH, Win W, Chee CB, Hsu LY, Mak E, Earnest A, Ong ME, Cutter J, Wang YT: Molecular epidemiology of mycobacterium tuberculosis complex in Singapore, 2006–2012. PLoS One. 2013, 8: e8448710.1371/journal.pone.0084487.
 24.
Li XX, Wang LX, Zhang H, Du X, Jiang SW, Shen T, Zhang YP, Zeng G: Seasonal variations in notification of active tuberculosis cases in China, 2005–2012. PLoS One. 2013, 8: e6810210.1371/journal.pone.0068102.
 25.
Douglas A, Strachan D, Maxwell J: Seasonality of tuberculosis: the reverse of other respiratory diseases in the UK. Thorax. 1996, 51: 944946. 10.1136/thx.51.9.944.
 26.
Kumar V, Singh A, Adhikary M, Daral S, Khokhar A, Singh S: Seasonality of tuberculosis in Delhi India: a time series analysis. Tuberc Res Treat. 2014, 2014:
 27.
Luquero F, SanchezPadilla E, SimonSoria F, Eiros J, Golub J: Trend and seasonality of tuberculosis in Spain, 19962004 [Short Communication]. Int J Tuberc Lung Dis. 2008, 12: 221224.
 28.
Rios M, Garcia J, Sanchez J, Perez D: A statistical analysis of the seasonality in pulmonary tuberculosis. Eur J Epidemiol. 2000, 16: 483488. 10.1023/A:1007653329972.
 29.
Health: News; Free TB treatment.
 30.
Foreign Workforce Numbers. [http://mom.gov.sg/statisticspublications/others/statistics/Pages/ForeignWorkforceNumbers.aspx]
 31.
WHO Western Pacific Region: Western Pacific Surveillance and Response; Using tuberculosis surveillance data for informed programmatic decisionmaking. [http://www.wpro.who.int/wpsar/volumes/04/1/2013_ED_Nishikiori/en/index.html]
Prepublication history
The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712458/14/1121/prepub
Acknowledgments
The authors appreciate the support of DukeNUS/SingHealth Academic Medicine Research Institute, the medical editing assistance of Taara Madhavan (Associate, Clinical Sciences, DukeNUS Graduate Medical School) and the advice of Dr Jeffery Cutter (Director, Communicable Diseases Division, Ministry of Health Singapore).
Financial support for this work was provided by a grant from the Centre for Infectious Disease Epidemiology and Research, Saw Swee Hock School of Public Health, National University of Singapore, Singapore.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
All the authors have neither commercial nor personal associations or any sources of support that might pose a competing interest in the subject matter or materials discussed in this manuscript.
Authors’ contributions
Conceived the study objectives and methodology; HLY, AE, SD, WW. Analyzed and interpreted the data; HLY, AE, SD, ARC, WW. Drafted the manuscript; HLY, AE, SD, CCBE, ARC, WW. All authors revised it critically for important intellectual content and approved and read the final manuscript.
Win Wah, Sourav Das contributed equally to this work.
Electronic supplementary material
12889_2014_7217_MOESM3_ESM.tiff
Additional file 3: Figure S3: Autocorrelation (ACF) and partial autocorrelation (PACF) plots of fitted residuals (Residents). (TIFF 859 KB)
12889_2014_7217_MOESM4_ESM.tiff
Additional file 4: Figure S4: Autocorrelation (ACF) and partial autocorrelation (PACF) plots of fitted residuals (Nonresidents). (TIFF 858 KB)
12889_2014_7217_MOESM5_ESM.tiff
Additional file 5: Figure S5: Autocorrelation and partial autocorrelation of final SARIMA model residuals (Residents). (TIFF 888 KB)
12889_2014_7217_MOESM6_ESM.tiff
Additional file 6: Figure S6: Autocorrelation and partial autocorrelation of final SARIMA model residuals (Nonresidents). (TIFF 898 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
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 https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://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
Wah, W., Das, S., Earnest, A. et al. Time series analysis of demographic and temporal trends of tuberculosis in Singapore. BMC Public Health 14, 1121 (2014). https://doi.org/10.1186/14712458141121
Received:
Accepted:
Published:
Keywords
 Timeseries analysis
 Tuberculosis