{'fraction_non_alphanumeric': 0.030102476515798463, 'fraction_numerical': 0.015554471147980968, 'mean_word_length': 4.793109540636042, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Covid‐19 cases prediction using SARIMAX Model by tuning hyperparameter through grid search cross‐validation approach

Abstract SARS‐Coronavirus was first detected in December 2019, later named COVID‐19, and declared a pandemic by the World Health Organization (WHO). As prediction models assist policymakers in making decisions based on expected outcomes. Existing models were only used to anticipate a smaller range of data resulting in irrelevant predictions. Our research focuses on predicting COVID‐19 confirmed, recovered, and deceased Indian cases for 20 days ahead. Tuning of hyperparameters is performed with a grid search cross‐validation approach. The dataset is collected from the Kaggle. Our forecast indicates that the count of confirmed and deceased cases is higher whereas, recovered cases prediction shows a decreasing trend. The R 2 Score achieved is 0.5112 and root‐mean‐square error (RMSE) is 1251 using optimized SARIMAX. Finally, Monte Carlo simulation has also been performed to justify the prediction accuracy as compared to other models such as linear, polynomial, prophet, and SARIMAX without grid search cross validation.

reliable trend predictions in this paper. The study was conducted in India since the both the population and density of the country are quite large, which often is the cause of rapid disease spread. However, the major concern for everybody is to figure out when would the COVID-19 pandemic end and would it return back. Despite the availability of various vaccines and well-established treatment methods, the COVID-19 pandemic continues to kill large number of people around the world. In different parts of the world, many COVID-19 variants have been discovered. Medical professionals are searching for new strategies, procedures, and technology to contain the virus and prevent further human loss of life.
Machine learning has been extensively used in forecasting or predicting future trends (Han & Chi, n.d.). The rate at which the coronavirus infects humans is quickly increasing, resulting in large number of infections (Andrea & Giuseppe, 2020). For getting an accurate and reliable forecast of increase, the autoregressive integrated moving average (ARIMA) model is widely used in time series prediction. In this work, seasonal ARIMA with eXogenous factors (SARIMAX) statistical model is castoff to achieve a 20-day forecast of aggregate COVID-19 cases for India. SAR-IMAX modelling is one of the high ranking modelling strategies for forecasting a time series. The SARIMAX models are always represented by a set of parameters, like (p, d, q) Â (P, D, Q, S). The auto-regression order is represented by p, the degree of trend variance is symbolized by d, and the (MA) order is symbolized by q. The value of an exogenous variable 'x' is determined outside the model and imposed on the model. In other words, variables that have an impact on a model yet are not influenced by it. P, D, Q and S are Seasonal ARIMA (SARIMA) combinations. Timeseries information is a collection of numerical numbers with a time component. There are two forms of time series data, namely stationary and non-stationary data. A stationary time series data have no associated forms with respect to time, whereas non-stationary time series data have certain patterns that are also mentioned as seasonality. The mean and variance of non-stationary data are not constant throughout. As a result, by computing the variance between two subsequent assumptions, non-stationary data can be turned into stationary data. Differencing is the term used for this technique. This method avoids seasonality and tendencies (Fanoodi et al., 2019). In order to create analytical models that are based on the time-series dataset, the following conventions are made (ArunKumar et al., 2021).
Absent of anomalies and outliers in time series data amounts to the following.
1. Univariate data (time-series dataset consists of only a single variable, as the SARIMA model regresses a variable with its own previous values) 2. Assumption of stationary data means variance and mean are constant throughout the time period.
3. Both Error terms and Model parameters are considered being constant w.r.t time.
Machine Learning is now a fast expanding and quickly developing field. It optimizes the performance of computers by programming them with data. Using training data or previous experiences, it learns to optimize the parameters for computer programs. It can also forecast the future using the data. Machine Learning can also assist in developing a mathematical model based on data statistics. Machine Learning's major goal is to learn from data feed without the intervention of humans. In other words, it learns from given data (experience) and produces the desired output by searching for trends/patterns in the data (Rao & Gudivada, 2018). The research in this paper consists of finding a better prediction model for everyday recovered cases, confirmed cases, and deceased cases in India. We have used an optimized SARIMAX model and tuned the hyperparameters using grid search cross validation to obtain better and optimal results. This provides an accurate COVID-19 case prediction in a short span of time. When it comes to implementing machine learning models for time series prediction, computational efficiency is crucial. Machine learning algorithms update model parameters automatically depending on data, but users are frequently required to set some additional parameters, called hyperparameters, which have a big impact on prediction accuracy.

| Motivation
The motivation behind predicting the COVID-19 cases is due to the exponential increase of cases of public healthcare fighting against COVID-19.
Proposed forecasting models may inspire and assist policymakers in making decisions based on expected outcomes. The Indian government has enacted a number of measures and prohibitions. The government's quick decisions have slowed the spread of COVID-19 to a great extent. However, the pandemic continues to spread despite these decisions. Indeed, the future spread predictions may be useful for future planning and controlling the COVID-19 spread. Furthermore, asymptomatic corona cases would have a significant role in controlling the disease's spread around the world. This has motivated us to incorporate existing examples in order to provide reliable trend predictions.

| Contribution of this research
The major contribution of the proposed research includes: 1. To anticipate virus confirmed, recovered and deceased cases in order to analyse how accurate the optimal SARIMAX model is in predicting and forecasting such events.

2.
To optimize the prediction model by tuning the hyperparameters and adopting a set of features that allows explaining the dependencies of the future COVID-19 cases.
3. To optimize the prediction results and advance the accuracy of the SARIMA model with the support of a grid search cross-validation approach.
4. To validate the proposed model with Monte Carlo Simulation and compare it with state-of-art models.
The proposed model is simple, consistent, and thorough. One of the strengths of the proposed model is the accuracy of the predictions it provides as opposed to the existing research works lacking accuracy. The findings of the proposed results are also compared to state-of-the-art models to evaluate model selection, fitting, and predicting accuracy. The majority of existing hyperparameters tweaking approaches are general in nature. The reason why the existing approaches are not accurate enough lies in the fact that they use only fewer range data to predict, which often is not very useful for predictions. As a result, most of the existing approaches are not predicting accurately due to sudden spikes in cases. It motivated us to aim for the most accurate predicting model using a grid search cross-validation approach for recovered, confirmed, and deceased cases. As a result, the proposed approach can be useful in disease control methods and in tracking the spread of COVID-19. In order to properly control the on-going pandemic, healthcare organizations and governments must be able to regulate and track the spread, which would alleviate people's fear and allow them to focus on preparing them for the next phase of the epidemic. It is well known that the epidemic has significantly reduced job and work opportunities for citizens and organizations, resulting in a financial burden on families. This study has the potential to assist governments in acting and making sound decisions, as well as planning for the future, in order to reduce public fear and lift up their morale to deal with the possibility of the next phases of the epidemic.
The research article is designed as follows. Section 1 shows the introduction. Section 2 presents an analysis of the existing methods. In particular, it analyses the prediction of the past survey on the COVID-19 dataset using machine learning techniques. Section 3 presents the methodology, dataset description, and model used in the work. Section 4 discusses the data validation. Finally, Sections 5 and 6 present the results and conclusion of the research in this article.

| LITERATURE REVIEW
COVID-19 is a communicable disease that has affected all walks of life and has severely curtailed crowded activities (Almutairi et al., 2021). As the coronavirus continued to persist, experts began to investigate COVID-19's future. There have been several types of research conducted for predicting different outcomes from the cases of SARS-CoV2 including the death, recovered and confirmed cases. Time series forecasting difficulties have been extensively researched in the literature, with COVID-19 forecasting emerging as a new topic. Forecasting models can be used to predict the disease's influence on the community, which can help to control the epidemic. In this section, a number of works have been listed. The focus of this review has been on the prediction in the cases which use the ARIMA model and its variants. For, several works have been considered., which has led to conclude that optimization of the model can be done with the variants of the ARIMA model to improve prediction accuracy.
Research in (Gecili et al., 2021) has provided four variations of time series prediction models for counting the infected, death, and recovery cases. The models were accessible in the prediction package in R language. Four models were selected for application to the publicly available dataset of COVID-19 for the Italy and the United States. The Holt and Trigonometric Exponential Smoothing State-Space model with Box-Cox transformation (TBATS) models was found to have fewer prediction errors and shorter prediction ranges as compared to the ARIMA and Spline models. When compared to the other models, the ARIMA model produced more consistent results. The Akaike information criterion (AIC) values of the TBATS and Spline models were identical and frequently lower than those of the ARIMA model, indicating that the TBATS and Spline models were better suited to the data. For the final time period, the mean absolute percentage error (MAPE) values for ARIMA, TBATS, Spline, and Holt were 6.1%, 8.4%, 6%, and 6.9%, respectively. The ARIMA model can also be adjusted to produce more exact results.
In Fanelli and Piazza (2020), the authors created the ARIMA(p, d, q) model and investigated the unique COVID-19 epidemiological pattern in the three most impacted nations of Europe, namely Spain, Italy, and France. The data used ranged from 21 February to 15 April 2020. The author had considered the model's several orders (p, d, and q) and chose the best performing order based on the lowest MAPE values for the three countries. The confirmed cases account for 10% to 20% of the total number of people who became infected. The apparent COVID-19 death rate in Italy was between 4% and 8%, but it was much lower in China (1% to 3%). Based on the estimation, 2500 ventilation units were considered to be a reasonable amount for the peak demand that health authorities in Italy should consider in their strategic planning. Further, this model can be used for the Indian dataset also.
In (Benvenuto et al., 2020), on the basis of Johns Hopkins epidemiology data, COVID-19 cases were forecasted using ARIMA. For the period between 20 January and 10 February 2020, the COVID-19 dataset is available. An auto regression (AR) model, an MA model, and an ARIMA model make up the SARIMA model . To stabilize the time-series log transformation, differences were preferred. The ARIMA model was examined using the autocorrelation function (ACF) and partial autocorrelation (PACF) correlograms. According to the findings, ARIMA (1,0,4) outperformed ARIMA (1,0,3) in forecasting COVID-19 incidence. Seasonality has little effect on the incidence or prevalence of COVID-19, according to ACF and PACF. Data on prevalence and incidence are forecasted with relative 95% confidence intervals. However, the trouble with this technique is that, while the number of confirmed occurrences increases, the frequency drops.
In (Singh et al., 2020), an amalgam model that is Wavelet-ARIMA is used for predicting the death cases that are majorly afflicted in five countries by COVID-19. France, Italy, Spain, the United Kingdom, and the United States are the five countries involved. Using discrete wavelet decomposition, the method divided the input dataset into component series, which were then separately subjected to an appropriate econometric model. COVID-19 death instances have increased significantly, according to the projection. The accuracy was investigated throughout the last 66 days. To build a better hybrid model, a mixed econometric model was merged with discrete wavelet decomposition of the dataset to more precisely forecast death occurrences. The Wavelet-ARIMA model produced an 80% better outcome for Italy, Spain, and the United Kingdom, and a 50% better outcome for France and the United States. Around half of the errors were reduced using the hybrid ARIMA model. When compared to the ARIMA model, the forecast obtained by the hybrid Wavelet-ARIMA model reduced errors by nearly 50%. In the case of Italy, Spain and the United Kingdom, the hybrid model outperformed the ARIMA model by more than 80%, whereas for France and the United States, it outperformed the ARIMA model by about 50%.
In (Hernandez-Matamoros et al., 2020), the authors have proposed an ARIMA model for each country with a lesser root-mean-square error (RMSE), to predict the evolution of the 145 countries that are scattered into 6 regions. The ARIMA model gives the possibility to predict the behaviour of a virus. The authors proposed a relation between the countries that lie in the same geographic areas and predicted the cases of COVID-19 with an RMSE mean value of 144.81. As a result, the model revealed a link between the anticipated error and the population of each 1 million people. The model was evaluated with 10% of the genuine data to determine the optimal parameter of the model. It was concluded that the ARIMA model can predict an average RMSE of 144.81.
In (Khan & Gupta, 2020), to forecast COVID-19 infection cases in India, a univariate time series model was presented. The data was collected between 31 January and 25 March 2020, and it was checked with data collected between 26 March and 4 April 2020. To test the accuracy of forecasted models, a nonlinear autoregressive neural network was built. For the next 50 days, the model anticipated COVID-19 cases. Based on data as of 4 April 2020, the results revealed an increase in the trend of true and expected numbers of cases of roughly 1500 per day. Based on the Bayesian Information Criteria (BIC) values and the overall highest R 2 values of 0.95, the optimal ARIMA (1,1,0) model was chosen. The Levenberg-Marquardt optimization training technique (LM) was used to optimize the NAR model architecture, which consisted of ten neurons and had the greatest R 2 value of 0.97.
In Duan and Zhang (2020) COVID-19 transmission was studied and anticipated, as well as the upper and lower bounds of the expected values, for a 7-day period from 27 April to 3 May 2020, using a 95% confidence level of new confirmed cases every day. The Wind database contains daily fresh confirmed cases in South Korea and Japan from 20 January 2020, until 26 April 2020. The model's orders were ARIMA (6,1,7) for Japan, which is the best fit, and ARIMA (2,1,3) for South Korea, which is also the best fit. As a result, the estimated ARIMA model did an excellent job of capturing the dependent structure of the daily new confirmed cases time series. Finally, the anticipated value, as well as the upper and lower bounds, were calculated using the ARIMA model. Limitations of the expected value for daily new confirmed cases under a 95% confidence level for the seven-day periods ending on 27 April 2020, and 3 May 2020, were stated.
In view of the foregoing discussion, the research in this article attempts to forecast the number of fatalities, recoveries, and confirmed cases using the SARIMAX model for the Indian dataset. This model with different orders generates a 20-day forecast of aggregate cases of  in the case of India as it is used for predicting the time series data. Previous research works were mainly focused on constructing methods to achieve an accurate and time-efficient model for predicting the spread of this pandemic situation. The major drawbacks of which included a lack of accuracy. As discussed, most of the existing approaches do not predict accurately due to sudden spikes in cases. It has motivated us to find the most accurate predicting model using a grid search cross-validation approach for confirmed, recovered, and deceased cases.

| Data
In our evaluation, the Indian dataset is collected from Kaggle (Kaggle, 2021). This is the world's largest dataset repository. The information gathered pertains to all COVID-19 cases reported in India between 1 February 2020 and 9 March 2021. Our dataset is in CSV format and contains many attributes related to COVID-19. The attributes include the date, timestamp, daily confirmed, total number of cases, everyday recovered, overall recovered, daily deceased, and overall deceased cases of India as shown in Table 1. Training and a test set were created from the dataset. Table 2, contains information about the experimental setup for the SARIMAX model.
The data were first evaluated and presented as a graph for confirmed instances (Figure 1), with the x-axis denoting days and the y-axis representing the total number of confirmed cases. The daily data of confirmed COVID-19 cases is depicted in Figure 1.

| SARIMAX model
Box and Jenkin first proposed the SARIMAX model in 1976 (Box & Jenkins, 1976). SARIMAX is a time-series data analysis statistical tool. To forecast the future values, the SARIMAX analyses the alterations between values in a time series. An SARIMAX model is made up of three elements.
The Auto regression, integrated regression, and MA are abbreviated as AR, I, and M, respectively. As in (Hernandez-Matamoros et al., 2020), every element is a parameter. Other models consist of the autoregressive (AR) model, the MA model, and the SARIMA model (Fattah et al., 2018 F I G U R E 1 Per day statistics of daily confirmed cases. effectively, the time series must be stationery and trend-free. The ADF (Augmented Dickey-Fuller) approach was utilized to determine whether or not the time-series data was stationary, including the log transformation and differences (Cheung & Lai, 1995).  SARIMAX model is given by Equation (1) as mentioned below, (Prabhakaran, 2020) where Y tÀ1 shows the first lag of the series, Y tÀ2 shows the second lag of the series and so on. β 1 shows the coefficient of first lag that the model examines, β 2 shows the coefficient of second lag that the model examines and so on. / shows the intercept term examined by the model, Y t is the output value which depends on its own lagged value and lagged predicted errors, € 1 and € tÀ1 states the errors of the equations and so on.

AR MODEL
Auto regressive model is one where the value of Y t depends only on its own lags, which is given by Equation (2) (Prabhakaran, 2020) F I G U R E 9 Prediction of daily confirmed cases using proposed SARIMAX with grid search CV model. where Y tÀ1 shows the first lag of the series, β 1 shows the coefficient of the first lag that the model examines and / shows the intercept term examined by the model.

MA MODEL.
Moving average is the model where the value Y t depends only on the lagged predicted errors. This is given by Equation (3) below (Prabhakaran, 2020).
where € 1 and € tÀ1 are the errors of Equations (4) and (5) as mentioned below a. Prediction for daily confirmed cases in India Figures 9-11 show the prediction graphs for the confirmed cases, where the orange colour graph shows the prediction curve and the blue colour curve shows the original data points. The prediction results were obtained through the optimal SARIMAX (4, 1, 5)Â(4, 1, 5, 9) by tuning the hyperparameters of the SARIMAX model to obtain the best order through grid search cv. Hyperparameter tuning is an optimization loop that is built on top of the machine learning model learning to categorize the collection of hyper parameters that results in the lowest validation set error. As a result, a validation set must be specified, as well as a loss (Bissuel, 2020). Table 4 shows the result of the SARIMAX optimized model for daily confirmed cases generated by the model. Figures 12-14, show the prediction graph for the recovered cases, where the orange colour graph shows the prediction curve and the blue colour curve shows the original data points. The prediction results were obtained through the optimal SAR-IMAX (5,1,1)Â(5,1,1,9). Table 5, shows the result of the SARIMAX optimized model for daily recovered cases generated by the model. Figures 15-17, show the prediction graph for the deceased cases, where the orange colour graph shows the prediction curve and the blue colour curve shows the original data points. The prediction results were obtained through the optimal SAR-IMAX (5,1,1)Â(5,1,1,9) through grid search cv by tuning the hyperparameters.

c. Prediction for daily deceased cases
F I G U R E 1 1 Prediction of daily confirmed cases (magnified view). Table 6 shows the result of SARIMAX optimized model for daily deceased cases generated by the model. d. ACF and PACF graphs for cases For the ARIMA model parameters, an ACF (auto-correlation function) graph and PACF (partial autocorrelation) correlogram were constructed. In a timeseries, the association between present observation with the past time-steps or lags observation.
Autocorrelation function (ACF), the plot of autocorrelation versus the lags is called autocorrelation plot. Hence, ACF illustrates the linear relation between the assumption at 't' time and observation at the past time 'tÀk ' (ArunKumar et al., 2021). ACF assists researchers in classifying knowledge linked to the previous concurrent finding. The partial ACF is used to determine the level of interface between observations made between two-time intervals of elimination. PACF aids in determining the degree of correctness of the current variables using its preceding values, while maintaining specified constant values (Makridakis et al., 2008). The ACF and PACF plots are shown in Table 7, The formula (Nayak et al., 2019), where k is the lag, provides the difference between Yt and Yt À k. Lag k autocorrelation indicate the association between the assumptions that are k time periods separated or aside.
Partial autocorrelation function (PACF), the intervening observations are focused while computing the association between the two assumptions at distinct times Y t represents the time series. So, PACF between the two observations, say Y t and Y tÀ2 can be shown in the equation below (ArunKumar et al., 2021).
F I G U R E 1 2 Prediction of daily recovered cases using the proposed SARIMAX with grid search CV model.
F I G U R E 1 3 Prediction of daily recovered cases using the proposed SARIMAX with grid search CV model showing the upper and lower confidence bound.
In Table 7 (see inline figures) the light blue colour indicates the values over which the autocorrelations are statistically significantly different from zero. In simple words, it describes how strongly the series' current value is connected to its previous values. Because ACF incorporates all of these factors while determining correlations, it is referred as a 'full auto-correlation plot'. We can predict a progressive decline in the ACF plot because the present has a good association with the previous lags. We expect PACF to drop sharply after close lags because these lags can record variation so well that we do not require previous lags to forecast present (Salvi, 2019) e. SARIMAX optimized model validation For hyperparameter selection in model, grid search approach is used. Minor changes in hyperparameter values can result in significantly different forecast results, selecting great hyperparameter values is crucial. Grid search is a strategy for determining the best hyperparameters for a model. Finding hyperparameters in training data, unlike parameters, is impossible. As a result, we develop a model for each combination of hyperparameters in order to determine the best hyperparameters. Cross-validation assesses how well a model generalizes to a new dataset. To acquire a good approximation of how well a predictive model performs, cross-validation is performed.
We have two datasets with this method: an independent dataset and a training dataset. We can divide a single dataset into two sets by partitioning it. GridSearchCV has two methods: 'fit' and 'score' (Ayuya, 2021). Figure 18 shows the proposed workflow for optimizing the SARIMAX model using grid search cv.
Algorithm for SARIMAX model using grid search approach:

| DATA VALIDATION
The disparities between the raw series observed and the anticipated values derived from the two techniques were compared to demonstrate the efficacy of the prediction strategy utilized in this investigation. Because the combined and chosen estimates for measuring bias and model accuracy as analytical approaches have been widely used, the R-Square and root-mean-square error (RMSE) were chosen as measurements (Christodoulos et al., 2011). The model is validated using the performance metrics as mentioned in Table 8,

| Monte Carlo simulation
Monte-Carlo sampling is a prominent tool for comparing the forecasting algorithms' performance. Multiple patches of a dataset are randomly picked in this procedure, and then tests of the forecasting systems' performance are run. It gives you the average of the error numbers you got for each data patch. The Monte-Carlo strategy's most remarkable feature is that it ensures an accurate comparison of forecasting systems and eliminates the possibility of skewed findings acquired by chance (Bokde et al., 2020).
Step 1: For each input random variable, define a distribution of possible inputs.
Step 2: Create random outputs from those distributions.
Step 3: Use that set of outputs to do a deterministic calculation.
Step 4: Combine the results of the separate computations to arrive at a final result.
Where, the parameter object indicates the output of prediction error function. Size shows the quantity of time series data used in the strategy. Iteration represents the number of iterations performed. The fval shows the flag to view predicted values in each iteration. Finally, figs shows the flag to display the plots for each iteration.
ÞR-Square score tells how close the forecasted value is to the regression line (Chordia & Pawar, 2021). The R-squared (R 2 ) score is a simple to compute and explain gauge of confidence (Lup on & Gaggin, 2015). It is the level to which a data point fits the linear regression, thereby indicating how well the regression line forecasts real values. The coefficient of determination gives you to compute how much variance the model's selfstanding variables have displayed. It gives us the model's goodness-of-fit, as well as a score that is always between 0 and 1 (Renaud & Victoria-Feser, 2010). Whereas the true value at certain ith assumption, is the forecasted value at ith assumptions, shows the mean of the overall observations, and is the overall count of observations Þ =n RMSE is defined as the variance or root of the residuals. It's a metric for how the residuals distribute around the best-fit line. It will be easy to read because the units are the same as the output units. This is often negative in nature, and a lower RMSE value enhances model performance. Whereas Y i shows the predicted value, X i shows the real value and n is the count of observations. (Moody, 2019; Willmott & Matsuura, 2005)

| RESULTS
Results obtained through the optimal SARIMAX model have shown a great improvement in the performance, tuning the hyperparameters through grid search cv accessing the best parameter for the SARIMAX model, which also improves the performance. Hence, the parameters (p, d, q)Â(P,D, F I G U R E 2 0 Residual error for predicting daily recovered cases using the SARIMAX optimized model.
F I G U R E 1 9 Residual error for predicting daily confirmed cases using the SARIMAX optimized model.
F I G U R E 2 1 Residual error for predicting daily deceased cases using the SARIMAX optimized model. Q,S) can be optimally obtained. The 'event' description and the data collection method, employed in the present case, determine the outcome of the prediction and estimation. Case definition and data gathering was preserved in real time for the future comparisons or viewpoints. In general, the fitted values and predicted values, produced using two methods (R 2 , RMSE), came quite close to the real incidence of COVID-19 disorders.
Although more data is needed to provide a more precise forecast, the virus's transmission looks to be significantly decreasing. Furthermore, despite the fact that the number of confirmed cases continues to rise, the frequency has decreased significantly. The number of cases will hit an all-time high if the virus does not evolve any new mutations. Figures 19-21 show the residual error plots for predicting daily confirmed, recovered and deceased cases.
The comparison has been made with other models showing the forecast for the 30-day time period, the metrics obtained are shown in Tables 9 using the proposed SARIMAX model, Table 5 shows the comparison of various state-of-art models with the proposed model. The comparison shows that our model is performing better in terms of minimum loss. The SARIMAX model with a grid search approach is producing better results as compared to SARIMAX without a grid search approach. This table shows the performance metrics of confirmed cases in India. The findings are correct, with a steady rise in the number of confirmed and deceased cases and a falling graph for recovered cases.
The Monte Carlo() function extends the prediction errors() module by allowing you to evaluate different approaches for randomly selected patches of the input time series Table 10 shows the results obtained through Monte Carlo using various models.

| CONCLUSION
In this paper, prediction of confirmed, recovered, and deceased Indian cases of COVID-19 obtained for 20 days in advance, until 29 March 2021, using an optimized SARIMAX statistical model with grid search cross validation approach. This study answers two of today's most pressing questions; namely when the COVID-19 epidemic would cease, and would there a chance of a second rebound if people resumed their normal lives.
Despite fast viral mutation and the structure of the dataset dependent on time and date, the research attempted to limit data variability by using just the dataset. Predicted results are more accurate with usage of the optimal SARIMAX model. Our forecast indicates that the count of confirmed as well as deceased cases are higher, whereas recovered cases prediction shows a decreasing trend. Proposed model can be extended to predict the cases using longer dataset. Furthermore, a hybrid model can be created incorporating proposed model. On the basics of our predictions, strategy administrators of healthcare should take necessary decisions at the right time in providing healthcare aids to the public as well as the agencies responsible for transporting equipment to hospitals. The prediction would also help the policymakers to frame their policies based on the pandemic situation. The result presented in this article emphasizes the significance of societal distancing and execution of various protective ways of COVID-19.