With many countries experiencing a resurgence in COVID-19 cases, it is important to forecast disease trends to enable effective planning and implementation of control measures. This study aims to develop Seasonal Autoregressive Integrated Moving Average (SARIMA) models using 593 data points and smoothened case and covariate time-series data to generate a 28-day forecast of COVID-19 case trends during the third wave in Malaysia. SARIMA models were developed using COVID-19 case data sourced from the Ministry of Health Malaysia's official website. Model training and validation was conducted from 22 January 2020 to 5 September 2021 using daily COVID-19 case data. The SARIMA model with the lowest root mean square error (RMSE), mean absolute percentage error (MAE) and Bayesian information criterion (BIC) was selected to generate forecasts from 6 September to 3 October 2021. The best SARIMA model with a RMSE = 73.374, MAE = 39.716 and BIC = 8.656 showed a downward trend of COVID-19 cases during the forecast period, wherein the observed daily cases were within the forecast range. The majority (89%) of the difference between the forecasted and observed values was well within a deviation range of 25%. Based on this work, we conclude that SARIMA models developed in this paper using 593 data points and smoothened data and sensitive covariates can generate accurate forecast of COVID-19 case trends.
This paper aims to develop an automated web application to generate validated daily effective reproduction numbers (Rt) which can be used to examine the effects of super-spreading events due to mass gatherings and the effectiveness of the various Movement Control Order (MCO) stringency levels on the outbreak progression of COVID-19 in Malaysia. The effective reproduction number, Rt, was estimated by adopting and modifying an Rt estimation algorithm using a validated distribution mean of 3.96 and standard deviation of 4.75 with a seven-day sliding window. The Rt values generated were validated using thea moving window SEIR model with a negative binomial likelihood fitted using methods from the Bayesian inferential framework. A Pearson's correlation between the Rt values estimated by the algorithm and the SEIR model was r = 0.70, p