Next Article in Journal
KCS-FCnet: Kernel Cross-Spectral Functional Connectivity Network for EEG-Based Motor Imagery Classification
Next Article in Special Issue
Editorial on Special Issue “Medical Data Processing and Analysis”
Previous Article in Journal
A Deep Learning Framework for Cardiac MR Under-Sampled Image Reconstruction with a Hybrid Spatial and k-Space Loss Function
Previous Article in Special Issue
Accuracy Analysis of Deep Learning Methods in Breast Cancer Classification: A Structured Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improvement of Time Forecasting Models Using Machine Learning for Future Pandemic Applications Based on COVID-19 Data 2020–2022

by
Abdul Aziz K Abdul Hamid
1,2,
Wan Imanul Aisyah Wan Mohamad Nawi
1,
Muhamad Safiih Lola
1,3,*,
Wan Azani Mustafa
4,5,*,
Siti Madhihah Abdul Malik
1,
Syerrina Zakaria
1,
Elayaraja Aruchunan
6,
Nurul Hila Zainuddin
7,
R.U. Gobithaasan
1,3 and
Mohd Tajuddin Abdullah
8,9
1
Faculty of Ocean Engineering Technology and Informatics, Universiti Malaysia Terengganu, Kuala Nerus 21030, Terengganu, Malaysia
2
Special Interest Group on Applied Informatics and Intelligent Applications (AINIA), Universiti Malaysia Terengganu, Kuala Nerus 21030, Terengganu, Malaysia
3
Special Interest Group on Modeling and Data Analytics (SIGMDA), Universiti Malaysia Terengganu, Kuala Nerus 21030, Terengganu, Malaysia
4
Faculty of Electronic Engineering & Technology, Pauh Putra Campus, Universiti Malaysia Perlis (UniMAP), Arau 02600, Perlis, Malaysia
5
Centre of Excellence for Advanced Computing, Pauh Putra Campus, Universiti Malaysia Perlis (UniMAP), Arau 02600, Perlis, Malaysia
6
Faculty of Science, Institute of Mathematical Sciences, Universiti Malaya, Kuala Lumpur 50603, Malaysia
7
Mathematics Department, Faculty of Science and Mathematics, Universiti Pendidikan Sultan Idris, Tanjong Malim 53900, Perak Darul Ridzuan, Malaysia
8
Faculty of Fisheries and Food Science, Universiti Malaysia Terengganu, Kuala Nerus 21030, Terengganu, Malaysia
9
Fellow Academy of Sciences Malaysia, Level 20, West Wing Tingkat 20, Menara MATRADE, Jalan Sultan Haji Ahmad Shah, Kuala Lumpur 50480, Malaysia
*
Authors to whom correspondence should be addressed.
Diagnostics 2023, 13(6), 1121; https://doi.org/10.3390/diagnostics13061121
Submission received: 22 January 2023 / Revised: 17 February 2023 / Accepted: 20 February 2023 / Published: 15 March 2023
(This article belongs to the Special Issue Medical Data Processing and Analysis)

Abstract

:
Improving forecasts, particularly the accuracy, efficiency, and precision of time-series forecasts, is becoming critical for authorities to predict, monitor, and prevent the spread of the Coronavirus disease. However, the results obtained from the predictive models are imprecise and inefficient because the dataset contains linear and non-linear patterns, respectively. Linear models such as autoregressive integrated moving average cannot be used effectively to predict complex time series, so nonlinear approaches are better suited for such a purpose. Therefore, to achieve a more accurate and efficient predictive value of COVID-19 that is closer to the true value of COVID-19, a hybrid approach was implemented. Therefore, the objectives of this study are twofold. The first objective is to propose intelligence-based prediction methods to achieve better prediction results called autoregressive integrated moving average–least-squares support vector machine. The second objective is to investigate the performance of these proposed models by comparing them with the autoregressive integrated moving average, support vector machine, least-squares support vector machine, and autoregressive integrated moving average–support vector machine. Our investigation is based on three COVID-19 real datasets, i.e., daily new cases data, daily new death cases data, and daily new recovered cases data. Then, statistical measures such as mean square error, root mean square error, mean absolute error, and mean absolute percentage error were performed to verify that the proposed models are better than the autoregressive integrated moving average, support vector machine model, least-squares support vector machine, and autoregressive integrated moving average–support vector machine. Empirical results using three recent datasets of known the Coronavirus Disease-19 cases in Malaysia show that the proposed model generates the smallest mean square error, root mean square error, mean absolute error, and mean absolute percentage error values for training and testing datasets compared to the autoregressive integrated moving average, support vector machine, least-squares support vector machine, and autoregressive integrated moving average–support vector machine models. This means that the predicted value of the proposed model is closer to the true value. These results demonstrate that the proposed model can generate estimates more accurately and efficiently. Compared to the autoregressive integrated moving average, support vector machine, least-squares support vector machine, and autoregressive integrated moving average–support vector machine models, our proposed models perform much better in terms of percent error reduction for both training and testing all datasets. Therefore, the proposed model is possibly the most efficient and effective way to improve prediction for future pandemic performance with a higher level of accuracy and efficiency.

1. Introduction

The city of Wuhan in Hubei Province, China, made history as the first point of the spread of the coronavirus disease (COVID-19) due to severe respiratory syndrome. On January 31, the World Health Organization (WHO) declared for the first time that COVID-19 is a “public health emergency of international concern” [1]. The virus was originally thought to have come from a fish market in Wuhan. On 11 January 2020, the gene sequence that China openly shared through personal contacts fueled its rapid spread, with a total of 9,129,146 confirmed cases, including 473,797 deaths worldwide as of 24 June 2020 [2]. However, as of 1 May 2021, COVID-19 has infected more than 151 million people and caused three million deaths worldwide. Countries such as the USA, Brazil, Russia, Spain, UK, Italy, France, Germany, China, India, Iran, and Pakistan have been hit the hardest by COVID-19. The first cases of COVID-19 reported in Malaysia on 2 January 2020 were detected by Chinese tourists entering the country from Singapore [3]. Only single-digit daily cases were reported in the initial phase, but this increased to 235 by 26 March [4]. The number of daily cases in Malaysia continued to increase exponentially by reaching around 20,000 in August 2021. The Malaysian government declared the implementation the Movement Control Order (MCO) from 18 March to 3 May 2020, the Conditional MCO (CMCO) from 4 May to 9 June 2020, and the Recovery MCO (RMCO) from 10 June 2020 to 31 March 2021. All travel and socio-economic activities (religious and cultural gatherings were not allowed) have been restricted across the country to keep new infections at bay and avoid overloading the country’s healthcare system during this time. All government and private offices and educational institutions, including transportation hubs, have been closed, citizens have been ordered to stay at home, and interstate travel has been banned, with fines of up to RM 10,000 for violators.
Since the WHO declared the COVID-19 outbreak a pandemic, not only governments from around the world, but also dedicated medical institutions have made many efforts to find vaccines and treatments to control the spread of the virus. Statisticians and public health scientists have also performed extensive statistical modelling, especially regarding the forecasting of COVID-19 cases, to help the health system prevent the contagion catastrophe. In this scenario, the ability to most effectively determine the growth rate at which the epidemic is spreading is very critical to counterattack and help governments, social planning, and policy making to accurately address the epidemic. Therefore, the motivation behind this research compared to the existing research is to (i) develop the most accurate and efficient predictive model related to the spread of COVID-19 in Malaysia and (ii) compare the performance of this new model with autoregressive integrated moving average (ARIMA), support vector machine (SVM), least-squares support vector machine (LSSVM), autoregressive integrated moving average-support vector machine ARIMA–SVM, and autoregressive integrated moving average–least-squares support vector machine models (ARIMA–LSSVM).
During the pandemic, many studies have been conducted using various mathematical and statistical models to predict the spread of the COVID-19 pandemic. One of the most popular time-series prognostic models for analyzing and predicting disease spread is the ARIMA (p, d, q) model [5,6,7]. Predicting new daily cases of COVID-19 was a difficult task, as cases increased daily. In the first wave, the pattern of COVID-19 cases has been continuously increasing for a period and then decreasing. However, for the second wave, it appears to be picking up again, and some of the COVID-19 cases are difficult to predict. In this scenario, some researchers predict the pattern of COVID-19 using ARIMA [8,9,10,11,12,13,14]. However, the ARIMA model has a limitation in that it typically can only handle a linear time-series data structure [15]. ARIMA model approximations are insufficient to pose a time-series prediction obstacle for researchers, especially for nonlinear patterns [16]. Despite its superior performance, the classification performance of Support Vector Machines (SVMs) and the generalizability of the classifier are often affected by the dimension or number of feature variables used, as mentioned by Lee [17]. As a result of the development of Vector Machines models, this process will be able to provide the most accurate and efficient result in each prediction case. SVMs, first introduced in 1995 by Vladimir Vapnik [18] in the field of statistical learning theory and structural risk minimization, have proven useful in a variety of prediction problems and classifications. SVMs could also manage or address difficulties such as non-linearity, local minimum, and high dimension where the ARIMA model could not [15,19,20,21]. SVM models have recently been used to handle problems such as nonlinear, local minimum, and high dimension. SVM can even guarantee higher accuracy for long-term predictions compared to other computational approaches in many practical applications. However, the single SVM model as a single ARIMA model also has some limitations, as the SVM model can only handle non-linear data and not linear data. With the limitations of a single ARIMA and SVM model, as well as an in-depth analysis of time-series prediction, hybrid approaches have become the best approach to overcome both limitations, and have a very significant impact in many areas due to their dynamic nature and higher level of predicting accuracy, efficiency, and precision. This approach is crucial because of the problems encountered in time-series forecasting, where almost all real time series contain linear and nonlinear correlation patterns between the data. Recently, the hybridization of prediction methods has been used with great success to achieve higher prediction accuracy [15,16,19,20,22,23,24,25,26].
Regarding the spread of COVID-19, the hybrid time-series model approach is crucial for predicting the impact of the COVID-19 outbreak, and has proven successful in predicting COVID-19 [27,28,29,30,31,32,33]. This study aims to (a) propose the ARIMA–LSSVM hybrid model approach to achieve better forecast results when it is able to produce the best estimator, i.e., produce small error terms; additionally, it aims to (b) examine the performance of the proposed models by comparing them to ARIMA and SVM models using three daily cases of COVID-19 data in Malaysia, that is, daily new positive cases, daily new deaths, and daily new recovered cases. Despite recent advances in time series and on COVID-19, the modelling process does not include COVID-19 cases specifically in Malaysia to help authorities manage the spread of this outbreak by producing more efficient, more accurate data, and more accurate forecasting results.
This study makes a significant contribution to the field of pandemic prediction and prevention by introducing novel approaches to dealing with COVID-19 data. Rather than relying on traditional methods, this research utilizes evidence-based prediction techniques, which have been shown to be more accurate and efficient. The use of these intelligent forecasting models enables local health authorities to create more precise and effective preventive measures, especially in the face of future outbreaks.
This study is particularly innovative in its use of hybrid forecasting models by machine learning for Malaysia’s future pandemics, such as avian flu or novel coronavirus strains. According to Moore [34], the scenario is for the next possible new pandemic of avian influenza virus strain H7N9 or a novel coronavirus. The predictive models developed are more precise, accurate, and efficient in anticipating the dynamic spread of the virus. This approach has been tested on real-world data, including daily new cases, daily new death cases, and daily new recovered cases of COVID-19, making it a valuable tool for public health officials and researchers. This research also has significant implications for future outbreaks, particularly in countries with tropical rainforests such as Malaysia. By predicting the spread of COVID-19 early on, this model can help policymakers build better healthcare facilities, take legislative action, and avoid economic losses. While a vaccine is now available, this model remains useful in accurately forecasting and preventing the impact of future pandemics, including those caused by new virus strains.
This study’s innovative and evidence-based methods make a valuable contribution to pandemic prediction and prevention, providing significant insights that can be used to mitigate the impact of future outbreaks. The implications of this research extend to public health authorities, policymakers, and researchers worldwide, offering powerful tools for mitigating the devastating effects of pandemics. The remainder of this paper is structured as follows. Materials and Methods goes into detail about the method we used to develop our proposed model. The hybrid ARIMA–SVM model used in this study is then briefly described. The results and discussion present the performance of our proposed model based on three known COVID-19 case datasets. Finally, we wrap up the article and make suggestions for future research.

2. Materials and Methods

2.1. ARIMA Modelling

The ARIMA (p, d, q) autoregressive integrated moving average model is one of the families in time-series forecasting that is widely used for time-series forecasting series datasets due to its flexibility with different time categories [16]. It also explicitly considers several standard patterns in time-series analysis, allowing for a powerful and easy-to-use way to produce accurate time-series forecasts. However, limitations may occur due to the existence of assumptions of a linear form that represents a linear relationship between the future value of the time series with the current value, the past value, and random noise in the model [15,16,17,21,26]. In the ARIMA model, p and q are the numbers of the autoregressive and moving average terms, and they are always listed in the order of the model, while d is the integer representing the differential order. The ARIMA model type with mean μ is represented mathematically as follows:
Y t = β + θ 1 Y t 1 + θ 2 Y t 2 + + θ p Y t p + e t 1 e t 1 2 e t 2 q e t q
where Y t   and   e t are the actual value and the random error at time t, respectively. Both are assumed to be independently and identically distributed (iid) with a mean 0 and a constant variance of σ 2 ; θ i   ( i = 1 , 2 , , q ) and j ( j = 0 , 1 , 2 , , q ) are the model parameters that need to be predicted.

2.2. Support Vector Machines Modelling

The Support Vector Machine (SVM) introduced by Vladimir Vapnik [18], which incorporates statistical learning theory, can handle larger dimensional data better, even with a small number of instances generalizability. Because the models select boundary support vectors from the input data, they process the data quickly. The SVM regression function is written as follows.
For linear and regressive dataset { x i , Y i } the function is formulated as follows:
f ( x ) = W T x + b
The coefficient W and b are estimated by minimizing.
1 2 W T W + C 1 n i 1 n L ε ( y i , f ( x i ) )
where ε is called the ε-intensive loss function and is formulated as follows:
ε ( Y , f ( x ) ) = {   0   if   | Y f ( x ) | e | Y f ( x ) |       others
Equation (3) can be transformed to the following constrained formulation by introducing positive slack variables ξ and ξ i * :
m i n 1 2 W T W + C i = 1 n ( ξ i + ξ i * ) W x i + b i Y i ξ + ξ i * W x i b i + Y i ξ + ξ i * ξ i , ξ i * 0 i = 1 , 2 , , N
We always use dual theory to convert the above formula into a convex quadratic programming problem when solving it. Adding the Lagrange Equation (5) results in the following term:
m i n 1 2 i , j = 1 n ( α i * α i ) ( α j * α j ) α i T α j i = 1 n α i * ( Y i e ) α i ( Y i + e )
Subject to
i = 1 Y ( α i α i * ) = 0   ,       α i , α i * [ 0 , C ]
When a dataset cannot be regressed linearly, we map it to a high dimension feature space and regress it linearly. The following is the formulation:
m i n 1 2 i , j = 1 n ( α i * α i ) ( α j * α j ) φ ( x i ) T φ ( x j ) i = 1 n α i * ( Y i e ) α i ( Y i + e )
Subject to
i = 1 n ( α i α i * ) = 0   ,       α i , α i * [ 0 , C ]
Let K ( X i , X j ) = { φ ( X i ) · φ ( X j ) } = φ T ( X j ) φ ( X i ) ; K ( x , x ) is the inner product of feature space and is called kernel function. Any symmetric function that satisfies Mercer condition can be used as Kernel function [19]. The Gaussian kernel function is specified in this study.
K ( x i , x j ) = e x p ( | | x i x i | | 2 / ( 2 σ 2 ) )
SVMs were used to estimate the nonlinear behaviour of the forecasting dataset because Gaussian kernels perform well under general smoothness assumptions [22].

2.3. Least-Square Support Vector Machines Modelling

The Least-Squares Support Vector Machines (LSSVM) proposed by Suykens and Vandewalle [35] is a modification of the standard SVM. LSSVM formulates the training process by solving linear problem quicker than SVM through quadratic programming. Additionally, this model is also more time efficient when analysing huge data. Consider a given training set { ( x i ,   y j ) ,   i = 1 , 2 , n } with x i R n as input data and y i R as output data. LSSVM defines the regression function as:
min J   ( ω , e ) = 1 2 ω T ω + C 2 i = 1 n e i 2
Subject to
Y i = ω T φ ( x i ) + b + e i ;     i = 1 , 2 , 3 n
where ω is the weight vector; y is the regularization parameter where it determines the trade-off between the training error minimization and smoothness of the estimated function; e is the approximation error; φ (.) is the nonlinear function; and b is the bias term. Constrained optimization of Equation (9) can be translated to unconstrained optimization by constructing Lagrange function. This can be obtained by using Karush–Kuhn–Tucker (KKT) condition, where it partially differentiates with respect to ω , b , e ,   and   φ ( . ) :
[ 0 I v t I v Ω + c 1 I ] [ b a ] = [ 0 Y ]
where Y = [ Y 1 , , Y n ] T ; I v = [ 1 , , 1 ] T ; a = [ a 1 , , a n ] T ; Ω = { Ω i j |   i , j = 1 n }   ,   Ω i j = φ ( x i ) T φ ( x i ) = K ( x i ,   x j ) ;   K ( . ) is the Radial Basis Function (RBF) kernel function that obtains a and b by calculating linear operations.

2.4. Proposed Hybrid Model

Despite the various time-series models presented, the accuracy, efficiency, and precision of time-series forecasts are becoming crucial for many decision-making processes today. However, these factors do not appear in ARIMA and SVM models. This is also the main reason why the time-series forecasting model is crucial, more demanding, and dynamic, as well as actively researched in many fields of study. ARIMA and SVM models have also prevailed in their linear or nonlinear domains [15,25,26]. However, none of these are generic principles that can be generalized to all situations. Therefore, a hybrid approach using both linear and non-linear modelling capabilities is recommended. This approach is mainly proposed to improve the overall prediction effectiveness. Therefore, there is no research on how to improve the effectiveness of predictive models created in Malaysia, particularly in the case of COVID-19.
There are two reasons for using hybrid models in this study. First, a single ARIMA and SVM model may not be sufficient to identify all the time series’ characteristics. The second assumption is that one or both cannot recognise the actual data generation process. This study’s hybrid models were built in two stages. Part I discusses linear autocorrelation composition, and Part II discusses nonlinear components. Thus,
Y t = t + N t
where t and N t are denoted as the linear composition and the nonlinear component, respectively. Based on the data, these two parts must be approximated. Part I focuses on linear modelling, which employs the ARIMA model to model the linear composition. The model from the first model included residuals, which are nonlinear interactions that cannot be modelled by a linear model or possibly a linear relationship. Thus,
t = [ i = 1 p θ i z t i i = 1 p j e t j ] + ε t = ^ t + ε t
Let ε i denoted as the residual from the linear model at time t. Then,
ε t = Y t ^ t
where ^ t is the predicted value for time t from the estimated relationship in (1), with ε t as the residual at time t from the linear model. The residual dataset after ARIMA fitting will only contain non-linear relationships that can be represented by a linear model [15]. The first stage results, which include forecast values and residuals from linear modelling, are then used in Part II.
Following Part II, the emphasis is on nonlinear modelling, where LSSVM is used to predict the nonlinear connection that occurs in residuals of linear modelling and original data. Then, the residual can be calculated using LSSVM by modelling various configurations as follows:
Part II focuses on nonlinear modelling, and LSSVM is used to model the nonlinear (possibly linear) relationship that occurs in residuals of linear modelling as well as original data. The residual can then be calculated using LSSVM by modelling different configurations as follows:
ε t = f   ( ε t 1 , ε t 2 , ε t n ) + e t
ε t = f   ( ε t 1 , ε t 12 ) + e t
Y t = f   ( Y t 1 , Y t 12 ,     ^ t ) + e t
Y t = f   ( Y t 1 , Y t 12 ) + e t
where f is a nonlinear function determined by the LSSVMs model and e t is the random errors. Thus, the hybrid forecast is
Y ^ t = ^ t + N t ^
Equations (12) and (13) can be identified as N t ^ , therefore the forecasted values can be achieved by summation of linear and nonlinear components. Figure 1 shows the functional flowchart of hybrid models.
In short, the proposed hybrid process methodology is divided into two parts. The ARIMA model is used to analyse the linear composition problem in Part I. Part II develops an LSSVM model to model the residuals from Part I. Because the ARIMA model in Part I cannot handle the nonlinear component of the data, the residuals of the linear model will include information about the nonlinearity. The LSSVM results can be utilised as forecasts of the ARIMA model’s error terms. The hybrid model defines various patterns by combining the distinct features and strengths of the ARIMA and LSSVM models. As a result, it is more effective to model linear and non-linear patterns separately with two different models and then re-hybridize the forecast results to improve overall modelling and forecasting performance.

2.5. Proposed Algorithm

  • Step 1: Three selected time series of COVID-19 cases datasets (1 October 2020–4 November 2022), namely daily new positive cases, daily new deaths cases, and daily new recovered cases, are generated in R programming Language.
  • Step 2: Each of the generated datasets is defined as { X 1 i = x 11 , x 12 , x 13 , ,   x n 1 } , { X 2 i = x 21 , x 22 , x 23 , ,   x 2 n } , and { X 3 i = x 31 , x 32 , x 33 , ,   x 3 n } for daily new positive cases, daily new deaths cases, and daily new recovered cases, respectively. Then, the best ARIMA (p, d, q) is selected after checking the autocorrelation function (ACF) plot of ARIMA (p, d, q) residuals. The best fitted value for daily new positive cases is ARIMA (2, 1, 2), while it is ARIMA (1, 1, 2) and ARIMA (0, 1, 1) for daily new fatalities cases and daily new recovered cases of COVID-19, respectively.
  • Step 3: The fitted value, Y t i = ( Y t 1 , Y t 2 , …, Y t m ) and the residuals ε t i = ( ε t 1 ,   ε t 2 ,   .   ,   ε t n ) .
  • Step 4: Combine the values in step 3 as a set of input variables to obtain the output Y t
  • Step 5: The ARIMA (p, d, q) is defined by the order of q. According to the information in step 4, Vector Machines is carried out to examine the residuals to obtain the output L t using R-programming Language.
  • Step 6: A fitted value of ARIMA with the hybridization of Vector Machines model is obtained for all sample data. Then, the residuals ε t is generated to obtain the forecasting result N t ^ .
  • Step 7: The framing data split randomly into training data and testing data for further Vector Machines modelling. Run the Vector Machines procedure using the “e1071” and “liquidSVM” package in R-Programming Language.
  • Step 8: The two modifiable parameters of the LSSVM technique (γ and σ) derived by objective function minimization such as mean square error (MSE). The grid-search method updates the parameters exponentially in the specified range using predetermined equidistant steps.
  • Step 9: Assume the split data as the processing data and the order q as in Step 5. Therefore, the combine forecast as in Equation (16): Y ^ t = ^ t + N t ^
  • Step 10: Estimate the model performance using the statistical measurement which are MSE, RMSE, MAE, and MAPE.

2.6. Forecasting Evaluation Criteria

In order to assess the overall performance of the proposed hybrid models, the one of a kind statistical measurements standard which accompanied by [15,16,36] including MAE (Mean Absolute Error), MAPE (Mean Absolute Percentage Error), MSE (Mean Squared Error), and RMSE (Root Mean Squared Error) are used.
MAE = 1 n t = 1 n | Y t ^ Y t |
MAPE = 1 n t = 1 n | Y t ^ Y t Y t | × 100
MSE = 1 n t = 1 n ( Y t ^ Y t ) 2
RMSE = 1 n t = 1 n ( Y t ^ Y t ) 2 = MSE
In time-series analysis, measurement tools such as Akaike’s information criterion (AIC) and the Bayesian information criterion (BIC) are commonly used to determine the appropriate length for distributed lag for the ARIMA model. As a result, model selection is based on the model with the lowest AIC and BIC values to provide measures of model performance, resulting in the selection of the best ARIMA model. Meanwhile, three parameters such as C are used as measurement tools to determine the best fitted model for LSSVMs models. Meanwhile, for the LSSVMs models, two parameters such as γ and σ are used as the measurement tools to determine the best fitted model.
Incorrect LSSVM model parameter selection can lead to over or underfitting of the training data. The parameter sets of the LSSVMs model with the lowest MSE value, as with the ARIMA model, will be selected for use in the best fitting model. As a result, for the hybrid models, the ARIMA first functioned as a pre-processor, filtering the linear pattern of datasets. The ARIMA model’s error term is then fed into the SVM in the hybrid models. LSSVMs were used to reduce the ARIMA error function.

3. Results and Discussion

3.1. Application of the Hybrid Model of COVID-19 in Malaysia

This section examined the proposed model’s performance in two ways: first, the performance of the proposed models compared to ARIMA, SVM, LSSVM, ARIMA–SVM, and ARIMA–LSSVM models; second, the percentage improvement of the proposed models compared to ARIMA and SVM models. Since the World Health Organization (WHO) declared COVID-19 to be a worldwide pandemic, the COVID-19 time-series datasets have been extensively studied. The predictive capability of the developed novel models was then compared using three well-known datasets of daily COVID-19 cases in Malaysia—daily new positive cases data, daily new fatalities cases data, and daily new recovered cases data—to demonstrate the performance of the proposed model in terms of accuracy, effectively, and accurately. All these data are reported from the 1 October 2020 to 4 November 2022 and retrieved from the COVIDNOW website at https://covidnow.moh.gov.my/, accessed on 10 January 2023.
The minimum value of new death, new cases, and new recovered cases in Table 1 is 0, 2600, and 1.8, respectively, while the maximum value of new cases, death, and recovered cases is 33,872.0, 592, and 33,406, respectively. Similarly, the mean and median for new cases, deaths, and recovered cases are 6322.7, 47.51, and 6415.5, respectively, where the parentheses indicate the median (3471, 11, 3447.0). The first quartile values for daily new cases, death cases, and recover cases are 1922, 4, and 1843, respectively. The number of daily new cases, deaths, and recoveries in the third quartile is 6824, 58, and 6775, respectively. Furthermore, the standard deviations for new cases, deaths, and recoveries are 7097.8, 81.12, and 7058.3 percentiles, respectively.
Furthermore, this section discusses the process of proposed models for both parts, i.e., Part I (Linear Modelling) and Part II (Nonlinear Modelling), using three well-known COVID-19 datasets, namely daily new positive cases, daily new deaths cases, and daily new recovered cases, to demonstrate the effectiveness of the proposed models. Both linear and nonlinear modelling, as well as the data used in this study, are carried out using R programming.
Part I (Linear Modelling)—ARIMA is used to generate the best ARIMA model for the daily new positive case dataset (2, 1, 2). ARIMA is the best fitting ARIMA model for the daily new death case dataset (1, 1, 2). Meanwhile, the best ARIMA model is reported as ARIMA in the case of the daily new recovered cases dataset (0, 1, 1). Table 2 summarizes the results of this ARIMA (p, d, q) model. Table 3 displays the estimates for all parameters. The p-values for all parameters are small, as shown in this table. As a result, for confirmed, recovered, and death cases, the models were statistically significant and could be used to forecast the future [37,38].
Part II (Nonlinear Modelling)—Based on the concepts of support vector machine design and the use of pruning algorithms in R-programming software, an optimal machine learning algorithm was created. For the daily new positive COVID-19 cases datasets, parameters γ = 264, σ = 0.008 show the smallest values of MSE i.e., 6,661,412 (see Table 4). Therefore, this parameters value was selected for use in the best-fitting model for the datasets of daily new positive COVID-19 cases. Whereas the smallest value of MSE is 250.887 and 21114252 (Table 4), with parameters γ = 877, σ   = 0.006 and γ = 334, σ   = 0.008 are selected as the best-fitting model for daily new death and daily new recovered cases of COVID-19, respectively.

3.1.1. New Positive Cases Data Forecasts

The daily new positive cases datasets series contains 765 data points and is recoded from 1 October 2020 to 4 November 2022 (see Figure 2). The number of daily new positive COVID-19 cases in Malaysia has increased significantly twice since July 2021, but has now dropped below 5000 new cases. However, it has continued and increased to a maximum of 33,406.00 around March–April 2022. This figure is expected to fall precipitously until 5 November 2022. The COVID-19 datasets have been extensively used with a wide range of linear and nonlinear time-series models, including ARIMA and machine learning methods [7,8,9,11,13,16,19,20,21,22,23,24,25,26]. The analysis of daily new positive cases of COVID-19 is critical as an indicator of the effectiveness of preventive measures that have been taken, are being taken, and will be taken by authorities to control the spread of this epidemic more effectively.
Therefore, a similar approach to that used by Aisyah et al. [15] is used to investigate the performance of the proposal models on daily new positive cases of COVID-19 datasets, where the dataset is divided into two samples, known as training sample and testing sample. According to Aisyah et al. [15] and Nurul Hila et al. [16], datasets should be divided into two parts to achieve the best results: 70–80% for training and the remaining 20–30% for testing [39,40]. The training data are used to assemble the models, while the testing data are used to evaluate the forecasting performances of the models based on statistical measurements. Thus, the daily new positive cases of the COVID-19 dataset are divided into two samples in this study: the training dataset and the test dataset. The training datasets contain 612 observations from day 1 to day 612, accounting for 80% of the datasets used exclusively to formulate from 1 October 2020 to 4 June 2022. In order to evaluate the forecasting performance of proposed models, the test sample datasets used approximately 153 observations from days 613–765 (20%) from the 5 June 2022 to the 4 November 2022
Table 5 displays the performance of the proposed model on the daily new positive COVID-19 case datasets. The proposed models produced results in terms of measurement error terms, namely MSE and MAE, which have smaller values of 10634.1142 and 46.54471, respectively. Similar results were obtained from the testing datasets, with MSE, MAPE, RMSE, and MAE values of 25478.114, 0.01547, 159.6182, and 75.6987, respectively. The findings are examined in greater detail using figures such as those shown in Figure 3a–e based on these numerical results. The estimated values for the proposed model (test sample) of daily new positive COVID-19 cases are shown in this figure. The proposed model line, as seen in this figure, closely matches the actual data. Figure 4a–e show the estimated values of our model for test sample data of ARIMA, SVM, LSSVM, ARIMA-SVM, and ARIMA–LSSVM models for COVID-19 cases. A comparison of the proposed model’s (ARIMA–LSSVM) lines for the test sample (Figure 4e) with the lines from the ARIMA, SVM, LSSVM, and ARIMA–SVM models (Figure 4a–d) clearly shows that the proposed model’s lines are somewhat close to the actual data. When we compared the performance of our proposal models to the performance of ARIMA, SVM, LSSVM, and ARIMA–SVM models, we discovered that our proposal models are efficient, accurate, and precise. In addition, the number of daily new positive COVID-19 cases is plotted, as shown in Figure 5. The daily new positive cases of COVID-19 for Malaysia are forecasted based on this figure for the next three weeks.
Based on Table 6, we examined the performance of the proposed models for the daily newly positive COVID-19 cases dataset by comparing the percentages of MSE, MAPE, RMSE, and MAE. The study hypothesis investigates the assumptions of the proposed hybrid model (ARIMA–LSSVM) approach to single ARIMA, SVM, LSSVM models, as well as hybrid ARIMA–LSSVM models. The proposed model outperformed the ARIMA–SVM model in MAE, MAPE, MSE, and RMSE, with improvements of 48.50%, 72.54%, 58.39%, and 35.50%, where the parentheses indicate an ARIMA, SVM, and LSSVM model, respectively, that results in the following: 80.96%, 89.80%, 91.48%, 70.81%; 80.61%, 89.97%, 90.72%, 69.54%; 63.20%, 79.59%, 69.31%, 44.60%. As a result of these findings (Table 5 and Table 6 and Figure 3, Figure 4 and Figure 5), it is possible to conclude that the proposed model produced greater accuracy and efficiency than ARIMA and SVM.

3.1.2. New Deaths Cases Data Forecasts

In addition to the Malaysian daily new positive COVID-19 cases datasets, the Malaysian daily new deaths cases datasets are taken into account and used to evaluate the performance of the proposed models. This dataset, like the daily new positive dataset and the daily new death case dataset, has a recording period of 1 October 2020 to 4 November 2022 (see Figure 6) and contains 765 data points divided into two samples. As the number of daily positive COVID-19 cases reported rises, so does the number of deaths, which now stands at around 600. The training dataset contains 612 observations (80%) from 1 October 2020 to 4 June 2022, and the test sample contains approximately 153 observations (20%) from 5 June 2022 to 4 November 2022 to evaluate the prediction performance of the proposed model.
A similar approach to the daily new positive cases of the COVID-19 dataset was used to study the performance of the proposed model on the daily new death cases of the COVID-19 dataset. The dataset was divided into two samples, namely training sample and testing sample. It accounts for approximately 80% of the daily new death cases in the COVID-19 dataset for the training sample (involving 612 observations with the period 1 October 2020 until 4 June 2022). The remaining 20% is for the test sample, which includes approximately 153 observations from 5 June 2022 to 4 November 2022.
As shown in Table 7, the performance of the proposed models using the daily new deaths datasets from COVID-19 is first characterized by statistical measurements such as MSE, MAPE, RMSE, and MAE. The results for the training data in this table show that the proposed model produces the smallest MSE and MAE values of 19.6422 and 1.03218, respectively, when compared to ARIMA, SVM, LSSVM and ARIMA–SVM. The same pattern can be seen in the test data, where all the statistical measures used have the lowest values when compared to the ARIMA, SVM, LSSVM, and ARIMA–SVM models.
The study then examines the estimated value of the suggested model for the COVID-19 case dataset for daily deaths, as shown in Figure 7a–e. This graph makes it abundantly clear that the proposed model line and the observed data are nearly identical. Additionally, Figure 8a–e each show the estimated values for the test sample for ARIMA, SVM, LSSVM, ARIMA–SVM and ARIMA–LSSVM the suggested models. Once more, it is obvious that when compared to ARIMA, SVM, LSSVM, and ARIMA–SVM models, the test sample lines for our suggested model (Figure 8e) are somewhat close to the actual data. This demonstrates that the outcomes of our suggested model are in line with prior findings and are more effective, accurate, and precise than those of ARIMA, SVM, LSSVM, and ARIMA–SVM models. The number of daily COVID-19 death cases is also plotted, just like in Figure 9. The daily new death cases of COVID-19 in Malaysia are anticipated to decline because of this number over the course of the following three weeks, indicating a downward trend.
As shown in Table 8, a similar method in the daily addition of positive COVID-19 case dataset is used to investigate the performance of the proposed model for the daily recorded death COVID-19 case dataset using percentages MSE, MAPE, RMSE, and MAE. Again, the percentage of improvement shows that our proposed model outperforms the ARIMA–LSSVM models for all statistical measures, with results of 5.07%, 1.82%, 3.80%, and 1.98%, respectively; there are also improvements (62.47%, 67.03%, 85.31%, 61.68%; 61.61%, 66.74%, 83.92%, 59.90%; 61.01%, 65.10%, 83.50%, 59.34%) for MAE, MAPE, MSE, and RMSE. The ARIMA, SVM, LSSVM, and ARIMA–SVM model results are shown in parentheses. The presented results (see Table 7 and Table 8 and Figure 7, Figure 8 and Figure 9) clearly show that our proposed model outperforms the ARIMA, SVM, LSSVM, and ARIMA–SVM models in terms of efficiency and accuracy.

3.1.3. New Recovered Cases Data Forecasts

The investigation to study the performance of the proposed model is continued with the dataset of new daily recovered cases of COVID-19 in Malaysia. Predicting Malaysia’s daily new recovered COVID-19 cases is just as important as the previous two datasets. The data used in this paper include daily observations from 1 October 2020 to 4 November 2022, for a total of 765 data points in the time series (Figure 10). The number of patients recovered from COVID-19 exhibits the same trend, with a significant increase twice. Beginning in July 2021, the number of recovered patients increases exponentially until it reaches over 22,500.00 in August 2021 (the time-series plot is shown in Figure 10) and then drops.
However, around March–April 2022, the number of recovered COVID-19 cases increased again to a maximum of 33,872.00, then decreased and showed a relatively stable movement after that. This dataset is also divided into two samples, namely the training dataset and the test dataset. The training dataset, which included 612 observations (80%) from 1 October 2020 to 4 June 2022, was used in the same way as the previous datasets to formulate the model. In contrast, the test sample uses approximately 153 observations (20%) for the period 5 June 2022–4 November 2022. Table 9 displayed the performance of the proposed model on the daily new recovered COVID-19 case datasets based on training and testing samples. The results in Table 9 clearly show that the proposed training sample model produces the smallest MSE and MAE values, with 47602.551 and 80.2214, respectively, when compared to the MSE and MAE models of the ARIMA, SVM, LSSVM, and ARIMA–LSSVM models. For the test sample, the same scenario as the training sample produced the smallest MSE, MAPE, RMSE, and MAE with values of 13004.11, 0.0125, 114.0351, and 54.14471, respectively, when compared to ARIMA, SVM, LSSVM, ARIMA–SVM and ARIMA–LSSVM models.
Figure 11a–e shows the estimated value of the dataset for daily new recovered COVID-19 cases for the test sample. Once more, this graph demonstrates how closely the predicted value from the proposed models seems to match the actual values. Figure 12a–e present an additional analysis of the outcomes of the proposed model. These plots (Figure 12a–e) show the predicted values for the test samples derived from ARIMA, SVM, LSAVM, ARIMA–SVM and ARIMA–LSSVM models. In these models, however, the proposed model is close to the true value because, as we shall see in Figure 11e, the proposed model dominates them. As shown in Figure 13, the number of daily new recovered COVID-19 cases is plotted. This figure makes it abundantly clear that the suggested model maintains the data’s original sharpness. The daily new recovered COVID-19 cases for Malaysia are predicted from this figure for the upcoming three weeks, and it suggests that these cases will rise in the days to come in Malaysia.
As shown in Table 10, further research was completed to determine how well the proposed models performed for the daily newly recovered COVID-19 case datasets for MSE, MAPE, RMSE, and MAE in terms of the percentage. When comparing the results of the proposed model to ARIMA, SVM, LSSVM, and ARIMA–SVM models, the results show a better improvement when looking at the percentage of improvement for statistical measurements such as MSE, MAPE, RMSE, and MAE, with results of 47.99%, 68.43%, 50.20%, and 9.42% improvement (86.02%, 91.99%, 95.21%, 78.11%; 85.43%, 91.69%, 94.57%, 76.71%; 81.06%, 88.91%, 91.32%, 70.54%). The results reported in the parentheses are the ARIMA, SVM, LSSVM, and ARIMA–SVM models. As a result, based on the findings, the proposed model has produced results that are more accurate and effective than those produced by ARIMA, SVM, LSSVM, and ARIMA–SVM models.

4. Conclusions

In conclusion, predicting the spread of COVID-19 with accuracy and efficiency is essential but frequently challenging for decision-makers, especially the front-line workers and health care authorities. Despite what might seem to be an endless spread of COVID-19, there have been numerous efforts to develop time-series models and ongoing research to enhance forecasting model efficacy. One of the most well-liked types of hybrid models that divide time series into linear and non-linear forms is the hybrid approach. In this study, a hybrid model that combines some linear and non-linear predictions is proposed. Utilizing three well-known COVID-19 datasets—daily new positive cases, daily new death cases, and daily new recovered cases—revealed that our proposed models were demonstrated as having the highest efficiency, accuracy, and precision. In comparison to ARIMA, SVM, LSSVM, and ARIMA–SVM models, the proposed model with cross-validation check based on MSE, RMSE, MAE, and MAPE makes the most accurate predictions. In terms of performance (the proposed models compared to ARIMA, SVM, LSSVM and ARIMA–SVM models) for both the training and testing datasets, the proposed models’ performance yields the smallest values of MSE, RMSE, MAE, and MAPE. This indicates that the proposed model’s predicted value is more closely aligned with the observed value. Therefore, our proposed models had a higher level of precision and could be suggested for COVID-19 forecasting. It can be concluded that the proposed model may be the most efficient and effective way to increase prediction accuracy performance, especially since it is important to anticipate and stop the spread of COVID-19 cases.

5. Limitations and Future Recommendation

In this research study, an attempt was made to predict the overall number of confirmed cases, fatalities, and recoveries of COVID-19 in Malaysia. Investigating SVM performance with various kernel functions and developing the best hyperparameters for the SVM forecasting model can help to increase the forecast’s accuracy in upcoming work. Since only one-step-ahead forecasting is considered in this paper, multi-step forecasts can be centralised in subsequent work. It has been demonstrated that multi-step forecasts can greatly increase the trading system’s realism [41,42]. Additionally, to improve the performance of the model in terms of efficiency and accuracy of dataset prediction, hybrid approaches such as bootstrap and double bootstrap methods [16,43,44] can be considered in the hybridization of ARIMA and SVM. Given the dearth of researchers using bootstrap in daily COVID-19 forecasting cases, it is a reliable method. Numerous studies have demonstrated that the bootstrap resampling method yields a more precise estimate [45]. Future studies should also consider (i) the clinical and behavioural aspects such as actions, cognition, and emotions and (ii) the possibility of the underreporting of cases and deaths, as well as delays in notifications, in order to avoid biased predictions, forecasts, and results.

Author Contributions

Conceptualization, A.A.K.A.H., W.I.A.W.M.N. and M.S.L.; Methodology, A.A.K.A.H., W.I.A.W.M.N., S.Z. and N.H.Z.; Software, A.A.K.A.H. and W.I.A.W.M.N.; Formal analysis, M.S.L.; Writing—original draft, W.I.A.W.M.N., M.S.L. and E.A.; Writing—review & editing, S.M.A.M., E.A., R.G. and M.T.A.; Supervision, M.S.L.; Project administration, M.S.L. and W.A.M.; Funding acquisition, M.S.L. and W.A.M. All authors have read and agreed to the published version of the manuscript.

Funding

The publication is partially sponsored by the Research Management Office, Universiti Malaysia Terengganu (UMT).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express his gratitude to the Research Management Centre, Universiti Malaysia Terengganu (UMT) and the Research Management Centre, Universiti Malaysia Perlis (UniMAP) for partially grant of the journal publication fee as well as to the editors and the referees for careful reading and for comments which greatly improved the paper.

Conflicts of Interest

The authors declare no conflict of interest. The authors identify and declare no personal circumstances or interests that may be perceived as inappropriately influencing the representation or interpretation of the reported research results. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Abdullah, M.T.; Lola, M.S.; Hisham, A.E.; Sabreena, S.; Nor Fazila, C.M.; Idham, K.; Dennis, C.Y.T. Framework of Measures for COVID-19 Pandemic in Malaysia: Threats, Initiatives and Opportunities. J. Sustain. Sci. Manag. 2022, 17, 8–18. [Google Scholar] [CrossRef]
  2. Ali, M.; Khan, D.M.; Aamir, M.; Khalil, U.; Khan, Z. Forecasting COVID-19 in Pakistan. PLoS ONE 2020, 15, e0242762. [Google Scholar] [CrossRef] [PubMed]
  3. WHO. Coronavirus Disease (COVID-19) in Malaysia. 2020. Available online: https://www.who.int/malaysia/emergencies/coronavirus-disease-(COVID-19)-in-Malaysia (accessed on 23 May 2020).
  4. KKM. COVID-19 Malaysia: Situasi Terkini 25 Oktober 2020. 2020. Available online: https://covid-19.moh.gov.my/terkini (accessed on 25 June 2022).
  5. Gecili, E.; Ziady, A.; Szczesniak, R.D. Forecasting COVID-19 confirmed cases, deaths and recoveries: Revisiting established time series modeling through novel applications for the USA and Italy. PLoS ONE 2021, 16, e0244173. [Google Scholar] [CrossRef]
  6. Awwad, F.A.; Mohamoud, M.A.; Abonazel, M.R. Estimating COVID-19 cases in Makkah region of Saudi Arabia: Space-time ARIMA modeling. PLoS ONE 2021, 16, e0250149. [Google Scholar] [CrossRef]
  7. Sahai, A.K.; Rath, N.; Sood, V.; Singh, M.P. ARIMA modelling & forecasting of COVID-19 in top five affected countries. Diabetes Metab. Syndr. Clin. Res. Rev. 2020, 14, 1419–1427. [Google Scholar] [CrossRef]
  8. Alzahrani, S.I.; Aljamaan, I.A.; Al-Fakih, E.A. Forecasting the Spread Of The COVID-19 Pandemic In Saudi Arabia Using ARIMA Prediction Model Under Current Public Health Interventions. J. Infect. Public Health. 2020, 13, 914–919. [Google Scholar] [CrossRef]
  9. Benvenuto, D.; Giovanetti, M.; Vassallo, L.; Angeletti, S.; Ciccozzi, M. Application of the ARIMA model on the COVID-2019 epidemic dataset. Data Brief 2020, 29, 105340. [Google Scholar] [CrossRef]
  10. Ceylan, Z. Estimation of COVID-19 prevalence in Italy, Spain, and France. Sci. Total Environ. 2020, 729, 138817. [Google Scholar] [CrossRef]
  11. Hernandez-Matamoros, A.; Fujita, H.; Hayashi, T.; Perez-Meana, H. Forecasting of COVID19 per regions using ARIMA models and polynomial functions. Appl. Soft Comput. 2020, 96, 106610. [Google Scholar] [CrossRef]
  12. Khan, F.M.; Gupta, R. ARIMA and NAR based prediction model for time series analysis of COVID-19 cases in India. J. Saf. Sci. Resil. 2020, 1, 12–18. [Google Scholar] [CrossRef]
  13. Kayode, O.; Fahimah, A.; Mustapha, R.; Jacques, D. Data Analysis and Forecasting of COVID-19 Pandemic in Kuwait Based on Daily Observation and Basic Reproduction Number Dynamics. Kuwait J. Sci. Special Issue 2021, 1–30. [Google Scholar] [CrossRef]
  14. Rahman, M.S.; Chowdhury, A.H.; Amrin, M. Accuracy comparison of ARIMA and XGBoost forecasting models in predicting the incidence of COVID-19 in Bangladesh. PLoS Glob. Public Health 2022, 2, e0000495. [Google Scholar] [CrossRef]
  15. Aisyah, W.I.W.M.N.; Muhamad Safiih, L.; Razak, Z.; Nurul Hila, Z.; Abd Aziz, K.A.H.; Elayaraja, A.; Nor Shairah, A.Z. Improved of Forecasting Sea Surface Temperature based on Hybrid ARIMA and Vector Machines Model. Malays. J. Fundam. Appl. Sci. 2021, 17, 609–620. [Google Scholar] [CrossRef]
  16. Nurul Hila, Z.; Muhamad Safiih, L.; Maman Abdurachman, D.; Fadhilah, Y.; Mohd Noor Afiq, R.; Aziz, D.; Yahaya, I.; Mohd Tajuddin, A. Improvement of Time Forecasting Models using A Novel Hybridization of Bootstrap and Double Bootstrap Artificial Neural Networks. Appl. Soft Comput. J. 2019, 84, 105676. [Google Scholar] [CrossRef]
  17. Lee, M.C. Using support vector machine with a hybrid feature selection method to the stock trend prediction. J. Expert Syst. Appl. 2009, 36, 10896–10904. [Google Scholar] [CrossRef]
  18. Vapnik, V.N. The Nature of Statistical Learning Theory, 1st ed.; Springer: New York, NY, USA, 1995. [Google Scholar]
  19. Sudheer, C.; Maheswaran, R.; Panigrahi, B.K.; Mathur, S. A hybrid SVM-PSO model for forecasting monthly streamflow. Neural Comput. Appl. 2013, 24, 1381–1389. [Google Scholar] [CrossRef]
  20. Chakraborty, T.; Chakraborty, A.K.; Biswas, M.; Banerjee, S.; Bhattacharya, S. Unemployment Rate Forecasting: A Hybrid Approach. Comput. Econ. 2020, 57, 183–201. [Google Scholar] [CrossRef]
  21. Zhang, G.P. Time series forecasting using a hybrid ARIMA and Neural Network Model. Neurocomputing 2003, 50, 159–175. [Google Scholar] [CrossRef]
  22. Terui, N.; Van Dijk, H. Combined forecasts from linear and nonlinear time series models. Int. J. Forecast. 2002, 18, 421–438. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, X.; Meng, M. A Hybrid Neural Network and ARIMA Model for Energy Consumption Forecasting. J. Comput. 2012, 7, 1184–1190. [Google Scholar] [CrossRef]
  24. Muhamad Safiih, L.; Nurul Hila, Z.; Mohd Tajuddin, A.; Vigneswary, P.; Mohd Noor Afiq, R.; Razak, Z.; Suffian, I.; Khalili, I. Improving the Performance of ANN-ARIMA Models for Predicting Water Quality in The Offshore Area of Kuala Terengganu, Terengganu, Malaysia. J. Sustain. Sci. Manag. 2018, 13, 27–37. [Google Scholar]
  25. Pai, P.F.; Lin, C.-S. A hybrid ARIMA and Support Vector Machines Model in Stock Price Forecasting. Int. J. Manag. Sci. 2005, 3, 497–505. [Google Scholar] [CrossRef]
  26. Lee, N.-U.; Shim, J.-S.; Ju, Y.-W.; Park, S.-C. Design and Implementation of the SARIMA–SVM time series analysis algorithm for the improvement of atmospheric environment forecast accuracy. Soft Comput. 2017, 22, 4275–4281. [Google Scholar] [CrossRef]
  27. Hao, Y.; Xu, T.; Hu, H.; Wang, P.; Bai, Y. Prediction and analysis of Corona Virus Disease 2019. PLoS ONE 2020, 15, e0239960. [Google Scholar] [CrossRef] [PubMed]
  28. Roy, S.; Ghosh, P. Factors affecting COVID-19 infected and death rates inform lockdown- related policymaking. PLoS ONE 2020, 15, e0241165. [Google Scholar] [CrossRef] [PubMed]
  29. Mahdavi, M.; Choubdar, H.; Zabeh, E.; Rieder, M.; Safavi-Naeini, S.; Jobbagy, Z.; Ghorbani, A.; Abedini, A.; Kiani, A.; Khanlarzadeh, V.; et al. A machine learning based exploration of COVID-19 mortality risk. PLoS ONE 2021, 16, e0252384. [Google Scholar] [CrossRef]
  30. Singhal, T. A Review of Coronavirus Disease-2019 (COVID-19). Indian J. Pediatr. 2020, 87, 281–286. [Google Scholar] [CrossRef] [Green Version]
  31. Qu, Z.; Li, Y.; Jiang, X.; Niu, C. An innovative ensemble model based on multiple neural networks and a novel heuristic optimization algorithm for COVID-19 forecasting. Expert Syst. Appl. 2023, 212, 118746. [Google Scholar] [CrossRef]
  32. Zivkovic, M.; Bacanin, N.; KVenkatachalam Nayyar, A.; Djordjevic, A.; Strumberger, I.; Al-Turjman, F. COVID-19 cases prediction by using hybrid machine learning and beetle antennae search approach. Sustain. Cities Soc. 2021, 66, 102669. [Google Scholar] [CrossRef]
  33. Melin, P.; Sánchez, D.; Castro, J.R.; Castillo, O. Design of Type-3 Fuzzy Systems and Ensemble Neural Networks for COVID-19 Time Series Prediction Using a Firefly Algorithm. Axioms 2022, 11, 410. [Google Scholar] [CrossRef]
  34. Sarah, M. The Future of Pandemics. News-Medical. Available online: https://www.news-medical.net/health/The-Future-of-Pandemics.aspx (accessed on 17 January 2022).
  35. Suykens, J.A.K.; Vandewalle, J. Least Squares Support Vector Machine Classifiers. Neural Process. Lett. 1999, 9, 293–300. [Google Scholar] [CrossRef]
  36. Naeem, M.; Yu, J.; Aamir, M.; Khan, S.A.; Adeleye, O.; Khan, Z. Comparative analysis of machine learning approaches to analyse and predict the COVID-19 outbreak. Peer J. Comput. Sci. 2021, 17, e746. [Google Scholar] [CrossRef] [PubMed]
  37. Qiang, X.; Aamir, M.; Naeem, M.; Ali, S.; Aslam, A.; Shao, Z. Analysis and Forecasting COVID-19 Outbreak in Pakistan Using Decomposition and Ensemble Model. Comput. Mater. Contin. 2021, 68, 842–856. [Google Scholar] [CrossRef]
  38. Adhikari, S.P.; Meng, S.; Wu, Y.-U.; Mao, Y.-P.; Ye, R.-X.; Wang, Q.-Z.; Sun, C.; Sylvia, S.; Rozelle, S.; Raat, H.; et al. Epidemiology, causes, clinical manifestation and diagnosis, prevention and control of coronavirus disease (COVID-19) during the early outbreak period: A scoping review. Infect. Dis. Poverty 2020, 9, 29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Ahmadini, A.A.H.; Naeem, M.; Aamir, M.; Dewan, R.; Alshqaq, S.S.A.; Mashwani, W.K. Analysis and Forecast of the Number of Deaths, Recovered Cases, and Confirmed Cases from COVID-19 for the Top Four Affected Countries Using Kalman Filter. Front. Phys. 2021, 9, 629320. [Google Scholar] [CrossRef]
  40. Alessa, A.A.; Alotaibie, T.M.; Elmoez, Z.; Alhamad, H.E. Impact of COVID-19 on Entrepreneurship and Consumer Behaviour: A Case Study in Saudi Arabia. J. Asian Financ. Econ. Bus. 2021, 8, 201–210. [Google Scholar] [CrossRef]
  41. Huck, N. Pairs trading and outranking: The multi-step-ahead forecasting case. Eur. J. Oper. Res. 2010, 207, 1702–1716. [Google Scholar] [CrossRef]
  42. Nisbet, R.; Elder JMiner, G. Chapter 11—Model Evaluation and Enhancement. In Handbook of Statistical Analysis and Data Mining Applications; Academic Press: Cambridge, MA, USA, 2018; pp. 215–233. [Google Scholar] [CrossRef]
  43. Nurul Hila, Z.; Muhamad Safiih, L. The Performance of BB-MCEWMA Model: Case Study on Sukuk Rantau Abang Capital Berhad, Malaysia. Int. J. Appl. Bus. Econ. Res. 2016, 14, 63–77. [Google Scholar]
  44. Nurul Hila, Z.; Muhamad Safiih, L.; Nur Shazrahanim, K. Modelling Moving Centreline Exponentially Weighted Moving Average (MCEWMA) with bootstrap approach: Case study on sukuk musyarakah of Rantau Abang Capital Berhad, Malaysia. Int. J. Appl. Bus. Econ. Res. 2016, 14, 621–638. [Google Scholar]
  45. Muhamad Safiih, L.; Nurul Hila, Z.; Mohd Noor Afiq, R.; Hizir, S. Double Bootstrap Control Chart for Monitoring SUKUK Volatility at Bursa Malaysia. J. Teknol. 2017, 79, 149–157. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flowchart process for hybrid ARIMA-LSSVM models.
Figure 1. Flowchart process for hybrid ARIMA-LSSVM models.
Diagnostics 13 01121 g001
Figure 2. Malaysian daily new positive COVID-19 cases (1 October 2020 to 4 November 2022).
Figure 2. Malaysian daily new positive COVID-19 cases (1 October 2020 to 4 November 2022).
Diagnostics 13 01121 g002
Figure 3. Results obtained from the proposed model for daily new positive COVID-19 cases dataset: (a) actual data vs. ARIMA model, (b) actual data vs. LSSVM models, (c) actual data vs. SVM model, (d) actual data vs. ARIMA–SVM models, (e) actual data vs. ARIMA–LSSVM models.
Figure 3. Results obtained from the proposed model for daily new positive COVID-19 cases dataset: (a) actual data vs. ARIMA model, (b) actual data vs. LSSVM models, (c) actual data vs. SVM model, (d) actual data vs. ARIMA–SVM models, (e) actual data vs. ARIMA–LSSVM models.
Diagnostics 13 01121 g003aDiagnostics 13 01121 g003b
Figure 4. Models’ prediction of daily new positive COVID-19 cases dataset (20% test sample): (a) actual data vs. ARIMA model, (b) actual data vs. SVM model, (c) actual data vs. LSSVM models, (d) actual data vs. ARIMA–SVM models, (e) actual data vs. ARIMA–LSSVM models.
Figure 4. Models’ prediction of daily new positive COVID-19 cases dataset (20% test sample): (a) actual data vs. ARIMA model, (b) actual data vs. SVM model, (c) actual data vs. LSSVM models, (d) actual data vs. ARIMA–SVM models, (e) actual data vs. ARIMA–LSSVM models.
Diagnostics 13 01121 g004aDiagnostics 13 01121 g004b
Figure 5. Actual and three weeks ahead forecasted values of ARIMA, SVM, LSSVM, ARIMA–SVM, and ARIMA–LSSVM models for new cases of COVID-19 of the 80% training and 20% testing set.
Figure 5. Actual and three weeks ahead forecasted values of ARIMA, SVM, LSSVM, ARIMA–SVM, and ARIMA–LSSVM models for new cases of COVID-19 of the 80% training and 20% testing set.
Diagnostics 13 01121 g005
Figure 6. Malaysian daily new deaths COVID-19 cases (1 October 2020 to 4 November 2022).
Figure 6. Malaysian daily new deaths COVID-19 cases (1 October 2020 to 4 November 2022).
Diagnostics 13 01121 g006
Figure 7. Results obtained from the proposed model for daily new death COVID-19 cases dataset: (a) actual data vs. ARIMA model, (b) actual data vs. LSSVM models, (c) actual data vs. SVM model, (d) actual data vs. ARIMA–SVM models, (e) actual data vs. ARIMA–LSSVM models.
Figure 7. Results obtained from the proposed model for daily new death COVID-19 cases dataset: (a) actual data vs. ARIMA model, (b) actual data vs. LSSVM models, (c) actual data vs. SVM model, (d) actual data vs. ARIMA–SVM models, (e) actual data vs. ARIMA–LSSVM models.
Diagnostics 13 01121 g007aDiagnostics 13 01121 g007b
Figure 8. Models’ prediction of daily new death COVID-19 cases dataset (20% test sample): (a) actual data vs. ARIMA model, (b) actual data vs. SVM model, (c) actual data vs. LSSVM models, (d) actual data vs. ARIMA–SVM models, (e) actual data vs. ARIMA–LSSVM models.
Figure 8. Models’ prediction of daily new death COVID-19 cases dataset (20% test sample): (a) actual data vs. ARIMA model, (b) actual data vs. SVM model, (c) actual data vs. LSSVM models, (d) actual data vs. ARIMA–SVM models, (e) actual data vs. ARIMA–LSSVM models.
Diagnostics 13 01121 g008aDiagnostics 13 01121 g008bDiagnostics 13 01121 g008c
Figure 9. Actual and three-weeks-ahead forecasted values of ARIMA, SVM, LSSVM, ARIMA–SVM and ARIMA–LSSVM models for daily new deaths COVID-19 cases of the 80% training and 20% testing set.
Figure 9. Actual and three-weeks-ahead forecasted values of ARIMA, SVM, LSSVM, ARIMA–SVM and ARIMA–LSSVM models for daily new deaths COVID-19 cases of the 80% training and 20% testing set.
Diagnostics 13 01121 g009
Figure 10. Malaysian daily new recovered COVID-19 cases (1 October 2020 to 4 November 2022).
Figure 10. Malaysian daily new recovered COVID-19 cases (1 October 2020 to 4 November 2022).
Diagnostics 13 01121 g010
Figure 11. Results obtained from the proposed model for daily new recovered COVID-19 cases dataset: (a) actual data vs. ARIMA model, (b) actual data vs. LSSVM models, (c) actual data vs. SVM model, (d) actual data vs. ARIMA–SVM models, (e) actual data vs. ARIMA–LSSVM models.
Figure 11. Results obtained from the proposed model for daily new recovered COVID-19 cases dataset: (a) actual data vs. ARIMA model, (b) actual data vs. LSSVM models, (c) actual data vs. SVM model, (d) actual data vs. ARIMA–SVM models, (e) actual data vs. ARIMA–LSSVM models.
Diagnostics 13 01121 g011aDiagnostics 13 01121 g011b
Figure 12. Models’ prediction of daily new recovered COVID-19 cases dataset (20% test sample): (a) actual data vs. ARIMA model, (b) actual data vs. SVM model, (c) actual data vs. LSSVM models, (d) actual data vs. ARIMA–SVM models, (e) actual data vs. ARIMA–LSSVM models.
Figure 12. Models’ prediction of daily new recovered COVID-19 cases dataset (20% test sample): (a) actual data vs. ARIMA model, (b) actual data vs. SVM model, (c) actual data vs. LSSVM models, (d) actual data vs. ARIMA–SVM models, (e) actual data vs. ARIMA–LSSVM models.
Diagnostics 13 01121 g012aDiagnostics 13 01121 g012b
Figure 13. Actual and three weeks ahead forecasted values of ARIMA, SVM, LS–SVM, ARIMA–SVM and ARIMA–LSSVM models for daily new recovered COVID-19 cases of the 80% training and 20% testing set.
Figure 13. Actual and three weeks ahead forecasted values of ARIMA, SVM, LS–SVM, ARIMA–SVM and ARIMA–LSSVM models for daily new recovered COVID-19 cases of the 80% training and 20% testing set.
Diagnostics 13 01121 g013
Table 1. Descriptive statistics of COVID-19 daily new cases, death, and recovered cases of Malaysia.
Table 1. Descriptive statistics of COVID-19 daily new cases, death, and recovered cases of Malaysia.
New CaseNew DeathNew Recovered
Min2.60000 × 10201.8
1st Qu1.9220 × 10341.8430 × 103
Median3.4710 × 103113.4470 × 103
Mean6.4155 × 10347.50986.3227 × 103
3rd Qu6.8240 × 103586.7750 × 103
Max3.3406 × 1045923.3872 × 104
SD7.0978 × 10381.12157.0583 × 103
Table 2. The best ARIMA (p, d, q) model selection.
Table 2. The best ARIMA (p, d, q) model selection.
COVID-19 Daily CasesARIMA (p, d, q)AICBIC
Daily New Positive Cases(2, 1, 2)12,564.5412,587.73
Daily New Deaths Cases(1, 1, 2)6930.126948.63
Daily New Recovered Cases(0, 1, 1)13,044.7413,054.01
Table 3. Parameter estimates of ARIMA models and their p-values.
Table 3. Parameter estimates of ARIMA models and their p-values.
Model ParametersEstimateZ-Statp-Value
New Case ARIMA (2, 1, 2)
θ 1 1.2408120.0852.2 × 10−16
θ 2 −0.9715−98.3202.2 × 10−16
φ 1 −1.2628−42.2252.2 × 10−16
φ 2 0.873848.1022.2 × 10−16
Recovered Case ARIMA (0, 1, 1) 2.2 × 10−16
φ 1 −0.3473−9.9532.2 × 10−16
Death Case ARIMA (1, 1, 2) 2.2 × 10−16
θ 1 0.859519.8522.2 × 10−16
φ 1 −1.6196−35.6512.2 × 10−16
φ 2 0.703920.4322.2 × 10−16
Table 4. LSSVMs Model Parameters for the daily new COVID-19 cases datasets.
Table 4. LSSVMs Model Parameters for the daily new COVID-19 cases datasets.
COVID-19 Daily CasesLSSVM ParameterMSE
γ = 11, σ = 0.00811,432,512
γ = 38, σ = 0.00810,235,488
Daily New Positive Casesγ = 74, σ = 0.0089,025,413
γ = 110, σ = 0.0088,014,123
γ = 264, σ = 0.0086,661,412
γ = 25, σ = 0.0061678.364
γ = 56, σ = 0.0061233.481
Daily New Deaths Casesγ = 277, σ = 0.006965.143
γ = 436, σ = 0.006554.368
γ = 877, σ = 0.006250.887
γ = 54, σ = 0.00828,412,113
γ = 89, σ = 0.00827,140,039
Daily New Recovered Casesγ = 125, σ = 0.00826,412,142
γ = 275, σ = 0.00823,032,256
γ = 334, σ = 0.00821,114,252
Table 5. Performance measures of the proposed model for daily new positive COVID-19 cases datasets.
Table 5. Performance measures of the proposed model for daily new positive COVID-19 cases datasets.
MODELSTRAINTEST
MSEMAEMSEMAPERMSEMAE
ARIMA929,843.169611.0274298,988.280.15167546.7982397.57
SVM8,355,184.4832001.644274,588.160.15421524.0116390.3848
LSSVM1084.1527739.538783,026.5500.07580288.1432205.6450
ARIMA–SVM42,552.713790.3484561,223.4740.05633247.4337146.9841
ARIMA–LSSVM10,634.114246.5447125,478.1140.01547159.618275.6987
Table 6. Percentage improvement of the proposed models with other forecasting models (the COVID-19 cases of daily new positive cases).
Table 6. Percentage improvement of the proposed models with other forecasting models (the COVID-19 cases of daily new positive cases).
ModelMAEMAPEMSERMSE
ARIMA80.959654989.8002241791.4785576270.80857252
SVM80.6092091789.9682251590.7213355469.53918577
LSSVM63.1896228979.5910290269.3132931644.60455773
ARIMA–SVM48.4987151772.536836558.3850566935.49051726
Table 7. Performance measures of the proposed model for daily new deaths COVID-19 cases datasets.
Table 7. Performance measures of the proposed model for daily new deaths COVID-19 cases datasets.
MODELSTRAINTEST
MSEMAEMSEMAPERMSEMAE
ARIMA697.99911.80836.067410.568382.463211.92791
SVM1409.1921.80065.389200.536872.321461.85605
LSSVM505.18111.43095.389200.536872.321461.85605
ARIMA–SVM49.44593.538120.926300.190880.963030.76230
ARIMA–LSSVM19.64221.032180.891140.187410.944000.72364
Table 8. Percentage improvement of the proposed models with other forecasting models (the COVID-19 cases of daily new death cases).
Table 8. Percentage improvement of the proposed models with other forecasting models (the COVID-19 cases of daily new death cases).
ModelMAEMAPEMSERMSE
ARIMA62.4650528367.0273408685.3126787261.67602437
SVM61.6059253966.7358892483.9233493459.90434808
LSSVM61.0118261965.0921079683.4643360859.33593514
ARIMA–LSSVM5.0714941621.817896063.7957465181.976054744
Table 9. Performance measures of the proposed model for daily new recovered COVID-19 cases datasets.
Table 9. Performance measures of the proposed model for daily new recovered COVID-19 cases datasets.
MODELSTRAINTEST
MSEMAEMSEMAPERMSEMAE
ARIMA1,802,678.36804.4378271,462.220.1560521.0203387.2768
SVM7,636,804.131890.917239,672.000.1504489.5630371.6573
LSSVM1,206,113.52723.9413149,871.530.1127387.1324285.9190
ARIMA–SVM99,205.699136.851926,108.020.0396161.5797104.1002
ARIMA–LSSVM47,602.55180.221413,004.110.0125114.035154.14471
Table 10. Percentage improvement of the proposed models with other forecasting models (the COVID-19 cases of daily new recovered cases).
Table 10. Percentage improvement of the proposed models with other forecasting models (the COVID-19 cases of daily new recovered cases).
ModelMAEMAPEMSERMSE
ARIMA86.0191186391.9871794995.2096059678.11311767
SVM85.4315494491.6888297994.5742055876.70675684
LSSVM81.0629199288.9086069291.3231619170.54364347
ARIMALSSVM47.9878905168.4343434350.1911290129.42485968
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

K Abdul Hamid, A.A.; Wan Mohamad Nawi, W.I.A.; Lola, M.S.; Mustafa, W.A.; Abdul Malik, S.M.; Zakaria, S.; Aruchunan, E.; Zainuddin, N.H.; Gobithaasan, R.U.; Abdullah, M.T. Improvement of Time Forecasting Models Using Machine Learning for Future Pandemic Applications Based on COVID-19 Data 2020–2022. Diagnostics 2023, 13, 1121. https://doi.org/10.3390/diagnostics13061121

AMA Style

K Abdul Hamid AA, Wan Mohamad Nawi WIA, Lola MS, Mustafa WA, Abdul Malik SM, Zakaria S, Aruchunan E, Zainuddin NH, Gobithaasan RU, Abdullah MT. Improvement of Time Forecasting Models Using Machine Learning for Future Pandemic Applications Based on COVID-19 Data 2020–2022. Diagnostics. 2023; 13(6):1121. https://doi.org/10.3390/diagnostics13061121

Chicago/Turabian Style

K Abdul Hamid, Abdul Aziz, Wan Imanul Aisyah Wan Mohamad Nawi, Muhamad Safiih Lola, Wan Azani Mustafa, Siti Madhihah Abdul Malik, Syerrina Zakaria, Elayaraja Aruchunan, Nurul Hila Zainuddin, R.U. Gobithaasan, and Mohd Tajuddin Abdullah. 2023. "Improvement of Time Forecasting Models Using Machine Learning for Future Pandemic Applications Based on COVID-19 Data 2020–2022" Diagnostics 13, no. 6: 1121. https://doi.org/10.3390/diagnostics13061121

APA Style

K Abdul Hamid, A. A., Wan Mohamad Nawi, W. I. A., Lola, M. S., Mustafa, W. A., Abdul Malik, S. M., Zakaria, S., Aruchunan, E., Zainuddin, N. H., Gobithaasan, R. U., & Abdullah, M. T. (2023). Improvement of Time Forecasting Models Using Machine Learning for Future Pandemic Applications Based on COVID-19 Data 2020–2022. Diagnostics, 13(6), 1121. https://doi.org/10.3390/diagnostics13061121

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop