Data source
We obtained data from the official website of Health Commission of Sichuan Province (http://wsjkw.sc.gov.cn/). The number of mumps cases per month was reported from January 2017 to July 2020. Data used in this study are available in the Supplementary Materials.
Non-pharmaceutical interventions against COVID-19
The first confirmed COVID-19 case from Wuhan in Sichuan was reported on January 21, 2020, and emergency response to major public health events in Sichuan upgraded to level 1 (the highest level) on January 24, 2020. A series of non-pharmaceutical interventions against COVID-19 were implemented, including social distancing, contact tracking, personal protection, travel restriction, etc. First, school holidays were extended, and schools were not allowed to open without permission of the local government, especially the kindergarten and primary school. Second, mass gatherings and major public events were canceled, and people were encouraged to stay at home as much as possible. Third, personal protective behaviors became normalized, such as hand washing and face masks. Additionally, the local government organized environmental disinfection in public places. Continuous rigorous social and public health measures could effectively decrease the spread of mumps.
Model construction
We modelled mumps incidence per month in Sichuan using the seasonal autoregressive integrated moving average (SARIMA) model, ARIMA (p, d, q) × (P, D, Q) S given by:
where p is the AR order, d is the degree of differencing, q is the MA order; P is the seasonal AR order, D is the degree of seasonal differencing, Q is the seasonal MA order, and S is the length of the seasonal period; B is the backshift operator.
Figure 1 showed the modelling process. We divided all data into three intervals: training interval (01/2017-06/2019), testing interval (07/2019-12/2019), and projecting interval (01/2020-12/2020). The stationarity of time series was judged by the Augmented Dickey-Fuller (ADF) test, and white noise was judged by the Ljung-Box test. We estimated d and D based on the time series diagram and ADF test. Next, we identified parameters p, d, Q, D by the autocorrelation function (ACF) and partial autocorrelation function (PACF), and optimized parameters by Akaike information criterion (AIC). Finally, residuals were diagnosed as white noise by the Ljung-Box test. All analyses were done using R software version 4.0.2. Code used in this study is available in the Supplementary Materials.
Model prediction and assessment
We projected the monthly incidence of mumps from July 2019 to December 2020. The projecting precision of the model was assessed by calculating the mean absolute percentage error (MAPE) between the reported cases and projected cases from July 2019 to December 2019. We compared the monthly number of mumps cases reported in 2020 against the average from corresponding periods in the preceding years (2017-2019) and the projected values by the SARIMA model. We performed a monthly one-sided (less) paired difference t-test using R version 4.0.2[6]. We calculated the percent change with 95% confidence interval (95%CI) in mumps cases reported in 2020 compared with the projected values using the Wilson method[11].