Hostname: page-component-76fb5796d-r6qrq Total loading time: 0 Render date: 2024-04-26T16:32:32.882Z Has data issue: false hasContentIssue false
  • English
  • Français

Forecasting the monthly incidence of scarlet fever in Chongqing, China using the SARIMA model

Published online by Cambridge University Press:  21 April 2022

W. W. Wu Open the ORCID record for W. W. Wu [Opens in a new window] ,
Q. Li ,
D. C. Tian ,
H. Zhao ,
Y. Xia ,
Y. Xiong ,
K. Su ,
W. G. Tang ,
X. Chen  and
J. Wang
...Show all authors
W. W. Wu
Affiliation:
Chongqing Center of Disease Control and Prevention, Chongqing 400042, China
Q. Li
Affiliation:
Chongqing Center of Disease Control and Prevention, Chongqing 400042, China
D. C. Tian
Affiliation:
School of Public Health (Shenzhen), Sun Yat-sen University, Shenzhen, Guangdong 510275, China
H. Zhao
Affiliation:
Chongqing Center of Disease Control and Prevention, Chongqing 400042, China
Y. Xia
Affiliation:
Chongqing Center of Disease Control and Prevention, Chongqing 400042, China
Y. Xiong
Affiliation:
Chongqing Center of Disease Control and Prevention, Chongqing 400042, China
K. Su
Affiliation:
Chongqing Center of Disease Control and Prevention, Chongqing 400042, China
W. G. Tang
Affiliation:
Chongqing Center of Disease Control and Prevention, Chongqing 400042, China
X. Chen
Affiliation:
Chongqing Center of Disease Control and Prevention, Chongqing 400042, China
J. Wang
Affiliation:
Chongqing Center of Disease Control and Prevention, Chongqing 400042, China
L. Qi*
Affiliation:
Chongqing Center of Disease Control and Prevention, Chongqing 400042, China
*
Author for correspondence: L. Qi, E-mail: qili19812012@126.com
Rights & Permissions [Opens in a new window]

Abstract

The incidence of scarlet fever has increased dramatically in recent years in Chongqing, China, but there has no effective method to forecast it. This study aimed to develop a forecasting model of the incidence of scarlet fever using a seasonal autoregressive integrated moving average (SARIMA) model. Monthly scarlet fever data between 2011 and 2019 in Chongqing, China were retrieved from the Notifiable Infectious Disease Surveillance System. From 2011 to 2019, a total of 5073 scarlet fever cases were reported in Chongqing, the male-to-female ratio was 1.44:1, children aged 3–9 years old accounted for 81.86% of the cases, while 42.70 and 42.58% of the reported cases were students and kindergarten children, respectively. The data from 2011 to 2018 were used to fit a SARIMA model and data in 2019 were used to validate the model. The normalised Bayesian information criterion (BIC), the coefficient of determination (R2) and the root mean squared error (RMSE) were used to evaluate the goodness-of-fit of the fitted model. The optimal SARIMA model was identified as (3, 1, 3) (3, 1, 0)12. The RMSE and mean absolute per cent error (MAPE) were used to assess the accuracy of the model. The RMSE and MAPE of the predicted values were 19.40 and 0.25 respectively, indicating that the predicted values matched the observed values reasonably well. Taken together, the SARIMA model could be employed to forecast scarlet fever incidence trend, providing support for scarlet fever control and prevention.

Keywords

Forecastinginfectious diseaseSARIMA modelscarlet fevertime series analysis
Type
Original Paper
Information
Epidemiology & Infection , Volume 150 , 2022 , e90
Check for updates
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © Chongqing Center of Disease Control and Prevention, 2022. Published by Cambridge University Press

Introduction

Scarlet fever is a contagious disease caused by group A streptococcus pyogenes (GAS), primarily affecting school children aged between 2 and 10 years [Reference Mahara1]. Patients and scarlet fever streptococcus carriers are the main infectious sources of scarlet fever. GAS is transmitted from human to human via direct contact, droplet spread from people with pharyngeal colonisation or carriage and contaminated fomites in addition to being reportedly food borne [Reference Ralph and Carapetis2]. The incubation period of scarlet fever is 2 to 10 days, and transmission period usually lasts 10 to 21 days without appropriate treatment [Reference Lamden3]. The symptoms of the disease include cute fever, sore throat, and diffuse bright red rash and obvious desquamation after rash regression [Reference Zhong and Wang4]. The incidence rate of scarlet fever is high in kindergartens and primary schools, but extremely low in children under the age of 2 because they are protected by maternal antibodies [Reference Mahara5]. Scarlet fever has been a feared killer during the eighteenth and nineteenth centuries in Europe, but it disappeared in twentieth century because of widely use of antibiotics as well as advancements in health care and nutrition [Reference Andrey and Posfay-Barbe6]. However, in past decade, scarlet fever had resurged in several countries and districts, including Hong Kong and South Korea in Asia, and UK in Europe [Reference Lamagni7–Reference Kim and Cheong9]. The causes of resurgence of scarlet fever may include microbial, host and meteorological factors [Reference Liu10]. Due to lack of vaccines for streptococcus pyogenes infection, the resurgence of scarlet fever has become a public health problem of global concern [Reference Ryu and Chun11].

Accurately modelling and forecasting the scarlet fever epidemic holds great importance as it helps the authorities make the necessary arrangements and allows a timely response. Many statistical models have been used for forecasting and early warning of disease outbreaks. One of the most popular model is the autoregressive integrated moving average model (ARIMA), which has been used to describe the trends of many diseases in epidemiological studies, including Haemorrhagic fever, H5N1 avian influenza and COVID-19 [Reference Li12–Reference Singh14]. Seasonal autoregressive integrated moving average (SARIMA) model considers the overall trend as well as seasonal variation, which was widely used in time series modelling [Reference Wang15]. Since scarlet fever presents typical seasonal characteristics, a SARIMA model may be more effective than ARIMA model in predicting the scarlet fever epidemic. However, no research has focused on the prediction of scarlet fever epidemics that considered specific seasonal trends observed in Chongqing, Southwest China. Therefore, in this study, we aimed to develop a prediction model using SARIMA model for scarlet fever in Chongqing, China.

Material and methods

Study region

With the latitude of 28°10′~32°13′ and longitude of 89°35′~103°04′, Chongqing is the largest municipality with over 32 million registered inhabitants in Southwest of mainland China (Fig. 1). The city of Chongqing has a subtropical monsoon humid climate characterised by a hot, humid summer and a cold, dry winter. The annual average temperature and precipitation in the city are 17°C and 1100 mm, respectively.

Fig. 1. Geographical location of Chongqing, China.

Data collection

All scarlet fever cases diagnosed in hospitals were reported to ‘National Disease Reporting Information System’ (NDRIS) within 24 h after diagnoses since 2004, according to the requirement of the law of the People's Republic of China on Prevention and Treatment of Infectious Diseases. The scarlet fever cases included clinical diagnosed and laboratory-confirmed cases. The diagnosis criteria for scarlet fever were provided by National Health Commission of the People’ Republic of China in 2008. A clinical diagnosed scarlet fever case was defined as a patient with fever, angina, rash, which rapid detection of group A Streptococcus test positive. A laboratory-confirmed case must be based on an aetiological examination. In Chongqing, scarlet fever can only be diagnosed in hospitals authorised by the health commission. We retrieved monthly reported scarlet fever cases information during January 2011 to December 2019 from NDRIS, including basic characteristics and the date of symptoms onset of cases. Additionally, the map of Chongqing was downloaded from National Earth System Science Data Center (http://www.geodata.cn/).

SARIMA model construction

The classical ARIMA (p, d, q) model is a frequently used technique for forecasting time series data. It has been proved as an effective and useful forecasting tool [Reference Wang16]. The SARIMA model is an extension of ARIMA model, which contains seasonal characteristics of time series [Reference Lin17]. The basic form of a SARIMA model is represented as SARIMA (p, d, q) (P, D, Q) s, where AR is the autoregressive model and MA is the moving average model; d and D represent the non-seasonal and seasonal differencing, respectively; p, q, P and Q denotes the orders of non-seasonal and seasonal autoregressive and moving average order, respectively; s denotes the length of seasonal period of the observed series, defined as 12 [Reference Xu18]. The formula of the SARIMA (p, d, q) (P, D, Q) s model is:

$$\phi ( B ) \Fcy ( {B^s} ) ( {\Delta^d\Delta_s^D y_t} ) = \theta ( B ) \Theta ( {B^s} ) \varepsilon _t$$
$$E( {\varepsilon_t} ) = 0, \;\;Var( {\varepsilon_t} ) = \sigma _\varepsilon ^2 , \;\;E( {\varepsilon_t\varepsilon_s} ) = 0, \;\;s\ne t$$
$$E( {y_t\varepsilon_t} ) = 0, \;\;\forall _s < t$$
$$\phi ( B ) = 1-\alpha _1B-\alpha _2B^2- \ldots \alpha _pB^p$$
$$\theta ( B ) = 1-\beta _1B-\beta _2B^2- \ldots \beta _qB^q$$
$$ \Phi \Fcy ( {B^s} ) = 1-\alpha _1B^s-\alpha _2B^{2s}- \ldots \alpha _PB^{Ps}$$
$$ \Theta ( {B^s} ) = 1-\beta _1B^s-\beta _2B^{2s}- \ldots \beta _QB^{Qs}$$

Where B denotes the backward shift operator, y t represents the actual number of scarlet fever cases at time t, and ɛ t denotes the residual errors from the SARIMA model. ϕ(B) is the non-seasonal p order autoregressive coefficient polynomial, and θ(B) is the non-seasonal q order autoregressive coefficient polynomial. $\Phi \Fcy ( {B^s} )$ denotes the seasonal P order autoregressive coefficient polynomial, Θ (B s) denotes the seasonal Q order autoregressive coefficient polynomial.

The entire time series data was divided into two parts: a training set consisting of data collected from 1 January 2011 to 31 December 2018 to fit the models, and a validation set collected from 1 January to 31 December 2019 to validate the fitted models. The process of the SARIMA model involves four steps [Reference Xu18]: Firstly, a time series graph of monthly scarlet fever cases from 2011 to 2019 was drawn to test its stationary intuitively. Then, non-seasonal difference or seasonal difference should be used to transform the time series to achieve stationarity if the original data is not stationary. The stationarity was checked by the autocorrelation function (ACF) and partial autocorrelation function (PACF) of the differenced series. Secondly, the orders of p, q, P and Q were identified by using the ACF and PACF. The least square method was used to estimate the model parameters. Thirdly, the Bayesian information criterion (BIC), the coefficient of determination (R 2) and the root mean squared error (RMSE) were used to evaluate the goodness-of-fit of the SARIMA models. The model with the lowest BIC and RMSE values, the highest R 2 value was regarded as the optimal model. The Ljung-Box Q test, ACF and PACF of the residuals were used to diagnose whether the residuals sequence was a white noise sequence. Finally, the optimal model was used for a prospective prediction, and the RMSE and mean absolute per cent error (MAPE) were used to assess the accuracy of the model.

$${\rm RMSE} = \sqrt {\;\displaystyle{1 \over n}\mathop \sum \nolimits ( {Y_{{\rm obs}}-Y_{{\rm pred}}} ) \;^2\;} $$
$$\! {\rm MAPE} = \displaystyle{1 \over n}\mathop \sum \nolimits{\displaystyle{{\vert Y_{{\rm obs}}-Y_{{\rm pred}}\vert } \over {Y_{{\rm obs}}}}}$$

Where Y obs and Y pred represent the observed and predicted numbers of scarlet fever cases, respectively, and n represents the number of months for forecasting.

The statistical analysis was performed with SPSS 22.0 (SPSS Inc., USA), and the map was created with ArcGIS software (version 10.2, ESRI Inc., Redlands, CA, USA). P < 0.05 was considered as statistical significance.

Results

Descriptive analysis

A total of 5073 scarlet fever cases were reported in Chongqing from 2011 to 2019, with an average annual incidence of 1.87 per 100 000 (ranged from 1.16 per 100 000 to 2.72 per 100 000). The male to female ratio was 1.44, with 2997 male cases and 2076 female cases. The median age of the cases was 6 years old and 81.86% were children aged 3 to 9 years old. Regarding to the population classification, students have the highest proportion (42.70%, 2166), followed by children in kindergarten (42.58%, 2160) and scattered children (13.50%, 685) (Table 1).

Table 1. The demographic characteristics of scarlet fever cases in Chongqing, 2011–2019

a Scattered children refer to children who do not reach the age of 3 years old to go kindergarten, or are taken care of by their family members.

Temporal pattern

The variation of yearly and monthly distributions of scarlet fever cases was shown in Figures 2 and 3. Semiannual peaks of scarlet fever were observed during the study period, with peaks occurred from April to June and October to December each year (Fig. 2). The first peaks were much higher than the second peaks. The incidence of scarlet fever was highest in 2015 (2.72 per 100 000) (Fig. 3). Figure 4 showed that the annual incidence rate among children aged 3–9 years was increasing gradually, with the incidence rate peaked in 2015. The incidence rates of other age groups remain stable.

Fig. 2. Monthly scarlet fever cases from 2011 to 2019 in Chongqing, China.

Fig. 3. Annual incidence and number of scarlet fever cases reported in Chongqing, 2011–2019.

Fig. 4. The incidence of scarlet fever by age group in Chongqing.

Spatial pattern

In general, the variation of geographical distribution of scarlet fever indicated that main urban districts had higher incidence than other regions. The top five districts reported most cases were Yubei (828 cases), Shapinba (599 cases), Jiulongpo (449 cases), Nanan (398 cases), Hechuan (350 cases). A total of 2624 scarlet fever cases were reported in these five districts, accounting for 51.72% of all cases (Fig. 5).

Fig. 5. Spatial distribution of scarlet fever incidence in Chongqing, 2011–2019.

SARIMA model

As shown in Figure 6, the number of scarlet fever cases in Chongqing fluctuated with time in 2011–2019 years, which showing an obvious seasonal periodicity. We performed a first-order non-seasonal difference and a first-order seasonal difference with a period of 12 to eliminate the trends and seasonal effects, respectively. Figures 7 and 8 indicated that the differenced time series tend to be stationary.

Fig. 6. Sequence diagram of monthly scarlet fever incidence in Chongqing, 2011–2019.

Fig. 7. Sequence diagram of scarlet fever incidence after a first-order non-seasonal difference and first – order seasonal difference.

Fig. 8. Autocorrelation function (ACF) and partial autocorrelation function (PACF) of monthly scarlet fever cases after a first-order non-seasonal difference and first – order seasonal difference.

According to the characteristics of ACF and PACF distribution (Fig. 8), we selected four models provisionally, and excluded the models in which the residuals were not likely to be white noise. Based on the results of the goodness-of-fit test statistics, we confirmed the optimal SARIMA (3, 1, 3) (3, 1, 0) 12 model, which had the lowest BIC (BIC = 6.054), lowest RMSE (RMSE = 15.81) and relatively high R 2 (R 2 = 0.789). Estimation of the SARIMA model parameters and the testing results were presented in Table 2 and S1 in Supplementary material. Figure 9 showed that the observed values and fitted values of the model matched well. Furthermore, the ACF and PACF of the residuals of SARIMA (3, 1, 3) (3, 1, 0)12 fell within the random confidence interval, indicating that the residuals did not deviate from a zero-mean white noise process (Fig. 10). Finally, the results of the Ljung-Box Q test (Q = 14.663, P = 0.101) for the model suggest that the residual series was a white noise sequence. The results of Goodness-of-fit test demonstrated that the performance of the fitted SARIMA (3, 1, 3) (3, 1, 0)12 model was reasonably well.

Fig. 9. Observed and forecasted values of the SARIMA (3, 1, 3) (3, 1, 0)12 model (UCL, upper confidence limit; LCL, low confidence limit).

Fig. 10. Autocorrelation function (ACF) and partial autocorrelation function (PACF) of the residuals series of the SARIMA (3, 1, 3) (3, 1, 0)12 model.

Table 2. Goodness of fits for the SARIMA models corresponding to different choices of p, q and P, Q which had passed the white noise test

Finally, the SARIMA (3, 1, 3) (3, 1, 0) 12 model was applied to forecast the monthly incidences of scarlet fever cases for 12 months from January to December in 2019. Figure 11 showed the comparison of observed and forecasted scarlet fever cases by the model, and the observed cases were within 95% confidence interval (CI) of the forecasted values. Table 3 showed the values of observed and predicted in 2019. The RMSE and MAPE between the predicted and actual values were 19.40, and 0.25 respectively, which verified the feasibility and effectiveness of the SARIMA (3, 1, 3) (3, 1, 0)12 model.

Fig. 11. Sequence diagram of observed and predicted values of scarlet fever in 2019 (UCL, upper confidence limit; LCL, low confidence limit).

Table 3. Comparison of observed and forecasted scarlet fever from January to December in 2019 by the SARIMA (3, 1, 3) (3, 1, 0)12 model

Discussion

Over the past decade, the incidence of scarlet fever has an anomalous increase in some Asian and European countries and regions, and this trend had exerted a marked influence on Chinese population since 2011 [Reference Chen19]. It is necessary to fully understand the epidemiological characteristics and predict the incidence of scarlet fever to provide information and evidence for its prevention and control.

Even though strategies and measurements, such as symptoms surveillance in the gate of kindergartens and primary schools, timely isolation of cases and daily disinfecting of environment, were brought up for intervening the transmission of scarlet fever, the incidence of scarlet fever increased sharply from 2004 to 2019 in China and other countries [Reference Andrey and Posfay-Barbe6, Reference Lu20, Reference Kim, Kim and Bae21]. Previous studies indicated that meteorological factors and air pollution factors may lead to an increase in scarlet fever incidence, such as ambient particulate matter, nitrogen dioxide, temperature and relative humidity were significantly associated with scarlet fever incidence [Reference Cheng22–Reference Liu24]. Moreover, some studies suggested that the rapid increase trend may be associated by the overuse of antibiotics, which would bring mutation and drug resistance to Streptococcal species [Reference You25]. However, the specific mechanism was remained unclear, and more research studies are needed in the future.

The time series of scarlet fever exhibited semiannual peaks in Chongqing, with a large peak during April to June and a small peak during October to December. The two peaks were coincided with local school semesters, which was similar to other provinces in China [Reference Duan23, Reference Lu26]. Therefore, it is crucial to enhance scarlet fever control and prevention measures in schools.

The incidence of scarlet fever was higher in boys, with a male-to-female ratio of 1.44. This observation was similar to the findings in Korea [Reference Kim, Kim and Bae27], mainland China [Reference Mahara5, Reference Zhang28]. Consistent with previous studies, males were more frequently infected by scarlet fever than females, this phenomenon may be related to some factors including male susceptibility at the host genetic level, host immune status, behaviour pattern, and it may also be due to reporting bias. The incidence of scarlet fever was highest among 3–9 years of age during 2011–2019. A study in Zhejiang, China, also reported similar finding from 2004 to 2018 [Reference Lu20]. In addition, several studies in Hong Kong [Reference Andrey and Posfay-Barbe6] and England [Reference Mahara5] indicated that the most prevalent ages of scarlet fever were 3–5, and 3–7, respectively. Therefore, kindergartens and primary schools are suggested to be the focus of monitoring and control of scarlet fever. It may be helpful to strengthen the monitoring and reporting of scarlet fever cases, so as to achieve early detection, early diagnosis and early isolation.

In this study, we found that the SARIMA (3, 1, 3) (3, 1, 0)12 model was the best-fitted mathematical model for forecasting monthly scarlet fever incidence based on the surveillance data from January 2011 to December 2018 in Chongqing, China. The model was then used to predict the incidence of scarlet fever in Chongqing in 2019. The results showed that the actual monthly incidence was within 95% CI of the predicted values, and the RMSE and MAPE between the predicted values and the actual values was 19.40 and 0.25 respectively, indicating that the model has reasonable accuracy in predicting the incidence of scarlet fever in Chongqing. Compared with the actual values, there had slight difference in predicted value in 2019, but the change trends of the two were basically the same. Due to the incidence of infectious disease was affected by complex factors, the time series model cannot predict it absolutely accurately, but the predicted value of the model can still accurately estimate the possible increase in the future.

The SARIMA model requires at least 7 or 8 seasonal cycles to estimate seasonal parameters [Reference Wu29]. We used nine seasonal cycles of incidence of scarlet fever, which meet the requirements of applying the SARIMA model. However, the epidemic of scarlet fever was affected by climates, environments and other factors [Reference Cheng22–Reference Liu24], and the change of any factor may lead to variety of the epidemic law of the disease. Therefore, the single constructed prediction model cannot be used as a permanent prediction tool. The prediction model should be adjusted according to the changes of the actual situation. Besides, further studies should explore variables that mainly affect the incidence of scarlet fever, including environmental factors, host factors and pathogen factors, and then incorporate them into the prediction model [Reference Dechuan30].

To the best of our knowledge, this is the most comprehensive study of scarlet fever up to now in Chongqing, China. The findings can be helpful for the control and prevention of scarlet fever epidemics in this area.

Several limitations of this study should be mentioned when interpreting the predicting results. First, the data were obtained from NDRIS, and the subclinical or asymptomatic cases of scarlet fever without seeking medical consultation were not reported in the system, therefore the reported data may underestimate the incidence. Second, the SARIMA model was only applicable to short-term prediction, and the model needs to be dynamically updated with new data to ensure the accuracy and stability of the prediction. Third, meteorological and environments factors such as temperature, rainfall, and PM10 may affect scarlet fever epidemics, but we were unable to incorporate these data into our study. Finally, we did not have data for emm gene typing of group A streptococcus in Chongqing.

In conclusion, the study provides key epidemic characteristics of scarlet fever in Chongqing, China. The SARIMA model can effectively predict the short-term scarlet fever incidence. The prediction model could be helpful in developing timely disease control and prevention strategies for scarlet fever in Chongqing, such as properly allocating health resources.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0950268822000693.

Acknowledgements

This work was supported by Chongqing Health Commission Program (No.2019GDRC014), Chongqing Science and Technology Bureau Program (No. CSTC2021 jscx-gksb-N0005).

Conflict of interest

We declare that all authors have no conflict of interest.

Ethical standards

This study was approved by Ethical Committee of Chongqing Center for Disease Control and Prevention.

Data availability statements

The datasets used and analysed during the current study is available from the corresponding author Li Qi (E-mail: ) on reasonable request.

Footnotes

*

These authors contributed equally to this work

References

Mahara, G et al. (2016) The association between environmental factors and scarlet fever incidence in Beijing region: using GIS and spatial regression models. International Journal of Environmental Research and Public Health 13, 1083.CrossRefGoogle ScholarPubMed
Ralph, AP and Carapetis, JR (2013) Group a streptococcal diseases and their global burden. Current Topics in Microbiology and Immunology 368, 1–27.Google Scholar
Lamden, KH (2011) An outbreak of scarlet fever in a primary school. Archives of Disease in Childhood 96, 394–397.CrossRefGoogle Scholar
Zhong, H and Wang, W (2020) Mathematical modelling for scarlet fever with direct and indirect infections. Journal of Biological Dynamics 14, 767–787.CrossRefGoogle ScholarPubMed
Mahara, G et al. (2016) Spatiotemporal pattern analysis of scarlet fever incidence in Beijing, China, 2005–2014. International Journal of Environmental Research and Public Health 13, 131.CrossRefGoogle Scholar
Andrey, DO and Posfay-Barbe, KM (2016) Re-emergence of scarlet fever: old players return? Expert Review of Anti-infective Therapy 14, 687–689.CrossRefGoogle ScholarPubMed
Lamagni, T et al. (2018) Resurgence of scarlet fever in England, 2014–16: a population-based surveillance study. The Lancet Infectious Diseases 18, 180–187.CrossRefGoogle ScholarPubMed
Lee, CF, Cowling, BJ and Lau, EHY (2017) Epidemiology of reemerging scarlet fever, Hong Kong, 2005–2015. Emerging Infectious Diseases 23, 1707–1710.CrossRefGoogle ScholarPubMed
Kim, JH and Cheong, HK (2018) Increasing number of scarlet fever cases, South Korea, 2011–2016. Emerging Infectious Diseases 24, 172–173.CrossRefGoogle ScholarPubMed
Liu, Y et al. (2018) Resurgence of scarlet fever in China: a 13-year population-based surveillance study. The Lancet Infectious Diseases 18, 903–912.CrossRefGoogle ScholarPubMed
Ryu, S and Chun, BC (2018) Investigation of scarlet fever outbreak in a kindergarten. Infection & Chemotherapy 50, 38–42.CrossRefGoogle Scholar
Li, S et al. (2016) Time series analysis of hemorrhagic fever with renal syndrome: a case study in Jiaonan County, China. PLoS One 11, e0163771.CrossRefGoogle ScholarPubMed
Kane, MJ et al. (2014) Comparison of ARIMA and random forest time series models for prediction of avian influenza H5N1 outbreaks. BMC Bioinformatics 15, 276.CrossRefGoogle ScholarPubMed
Singh, S et al. (2020) Forecasting daily confirmed COVID-19 cases in Malaysia using ARIMA models. Journal of Infection in Developing Countries 14, 971–976.CrossRefGoogle ScholarPubMed
Wang, H et al. (2018) Time-series analysis of tuberculosis from 2005 to 2017 in China. Epidemiology and Infection 146, 935–939.CrossRefGoogle ScholarPubMed
Wang, Y et al. (2019) Development and evaluation of a deep learning approach for modeling seasonality and trends in hand-foot-mouth disease incidence in mainland China. Scientific Reports 9, 8046.CrossRefGoogle ScholarPubMed
Lin, Y et al. (2015) Application of an autoregressive integrated moving average model for predicting injury mortality in Xiamen, China. BMJ Open 5, e008491.CrossRefGoogle ScholarPubMed
Xu, Q et al. (2017) Forecasting the incidence of mumps in Zibo City based on a SARIMA model. International Journal of Environmental Research and Public Health 14, 925.CrossRefGoogle ScholarPubMed
Chen, H et al. (2019) Epidemiological study of scarlet fever in Shenyang, China. BMC Infectious Diseases 19, 1074.CrossRefGoogle ScholarPubMed
Lu, Q et al. (2019) Analysis of epidemiological characteristics of scarlet fever in Zhejiang Province, China, 2004–2018. International Journal of Environmental Research and Public Health 16, 3454.CrossRefGoogle Scholar
Kim, J, Kim, JE and Bae, JM (2019) Incidence of scarlet fever in children in Jeju Province, Korea, 2002–2016: an age-period-cohort analysis. Journal of Preventive Medicine and Public Health = Yebang Uihakhoe Chi 52, 188–194.CrossRefGoogle ScholarPubMed
Cheng, W et al. (2020) The association between ambient particulate matters, nitrogen dioxide, and childhood scarlet fever in Hangzhou, Eastern China, 2014–2018. Chemosphere 246, 125–826.CrossRefGoogle Scholar
Duan, Y et al. (2016) Impact of meteorological changes on the incidence of scarlet fever in Hefei City, China. International Journal of Biometeorology 60, 1543–1550.CrossRefGoogle ScholarPubMed
Liu, Y et al. (2020) Exposure to air pollution and scarlet fever resurgence in China: a six-year surveillance study. Nature Communications 11, 4229.CrossRefGoogle ScholarPubMed
You, Y et al. (2018) Scarlet fever epidemic in China caused by Streptococcus pyogenes serotype M12: epidemiologic and molecular analysis. EBioMedicine 28, 128–135.CrossRefGoogle ScholarPubMed
Lu, JY et al. (2019) Effect of meteorological factors on scarlet fever incidence in Guangzhou City, Southern China, 2006–2017. Science of the Total Environment 663, 227–235.CrossRefGoogle Scholar
Kim, J, Kim, JE and Bae, JM (2019) Incidence of scarlet fever in children in Jeju Province, Korea, 2002–2016: an age-period-cohort analysis. Journal of Preventive Medicine and Public Health = Yebang Uihakhoe Chi 52, 188–194.CrossRefGoogle ScholarPubMed
Zhang, Q et al. (2017) Spatiotemporal epidemiology of scarlet fever in Jiangsu Province, China, 2005–2015. BMC Infectious Diseases 17, 596.CrossRefGoogle Scholar
Wu, J et al. (2007) Prediction of incidence of notifiable contagious diseases by application of time series model. Journal of Mathematical Medicine 20, 90–92.Google Scholar
Dechuan, K et al. (2019) Application of exponential smoothing model in predicting incidence of scarlet fever in Shanghai. Disease Surveillance 34, 932–936.Google Scholar
Figure 0

Fig. 1. Geographical location of Chongqing, China.

Figure 1

Table 1. The demographic characteristics of scarlet fever cases in Chongqing, 2011–2019

Figure 2

Fig. 2. Monthly scarlet fever cases from 2011 to 2019 in Chongqing, China.

Figure 3

Fig. 3. Annual incidence and number of scarlet fever cases reported in Chongqing, 2011–2019.

Figure 4

Fig. 4. The incidence of scarlet fever by age group in Chongqing.

Figure 5

Fig. 5. Spatial distribution of scarlet fever incidence in Chongqing, 2011–2019.

Figure 6

Fig. 6. Sequence diagram of monthly scarlet fever incidence in Chongqing, 2011–2019.

Figure 7

Fig. 7. Sequence diagram of scarlet fever incidence after a first-order non-seasonal difference and first – order seasonal difference.

Figure 8

Fig. 8. Autocorrelation function (ACF) and partial autocorrelation function (PACF) of monthly scarlet fever cases after a first-order non-seasonal difference and first – order seasonal difference.

Figure 9

Fig. 9. Observed and forecasted values of the SARIMA (3, 1, 3) (3, 1, 0)12 model (UCL, upper confidence limit; LCL, low confidence limit).

Figure 10

Fig. 10. Autocorrelation function (ACF) and partial autocorrelation function (PACF) of the residuals series of the SARIMA (3, 1, 3) (3, 1, 0)12 model.

Figure 11

Table 2. Goodness of fits for the SARIMA models corresponding to different choices of p, q and P, Q which had passed the white noise test

Figure 12

Fig. 11. Sequence diagram of observed and predicted values of scarlet fever in 2019 (UCL, upper confidence limit; LCL, low confidence limit).

Figure 13

Table 3. Comparison of observed and forecasted scarlet fever from January to December in 2019 by the SARIMA (3, 1, 3) (3, 1, 0)12 model

Supplementary material: File

Wu et al. supplementary material

Wu et al. supplementary material

Download Wu et al. supplementary material(File)
File 39 KB