Abstract
In this article, we analyze the growth pattern of Covid-19 pandemic in India from March 4th to July 11th using regression analysis (exponential and polynomial), auto-regressive integrated moving averages (ARIMA) model as well as exponential smoothing and Holt-Winters smoothing models. We found that the growth of Covid-19 cases follows a power regime of (t2,t,..) after the exponential growth. We found the optimal change points from where the Covid-19 cases shift their course of growth from exponential to quadratic and then from quadratic to linear. After that, we saw a sudden spike in the course of the spread of Covid-19 and the growth moved from linear to quadratic and then to quartic, which is alarming. We have also found the best fitted regression models using the various criteria such as significant p-values, coefficients of determination and ANOVA etc. Further, we search the best fitting ARIMA model for the data using the AIC (Akaike Information Criterion) and provide the forecast of Covid-19 cases for future days. It was noticed that the ARIMA model fits better the Covid-19 cases for small regions. ARIMA (5, 2, 5) and ARIMA (5, 2, 3) are the best possible models for modeling Covid-19 cases for March 4th to July 10th and June 1th to July 10th, respectively.
1. Introduction
The Covid-19 pandemic has created a lot of havoc in the world. It is caused by a virus called SARS-CoV-2, which comes from the family of coronaviruses and is believed to be originated from the unhygienic wet seafood market in Wuhan, China but it has now infected around 215 countries of the world. With more than 13.2 million people affected around the world and more 575,000 deaths (As of July 14th, 2020), it has forced people to stay in their homes and has caused huge devastation in the world economy. (Singh & Jadaun [1], [2], Gupta et al. [3]).
In India, the first case of Covid-19 was reported on 30th January, which is linked to the Wuhan city of China (as the patient has travel history to the city). On 4th March, India saw a sudden hike in the number of cases and since then, the numbers are increasing day by day. As of 14th July, India has more than 908,000 cases with more than 23,000 deaths and is world’s 3rd most infected country. ([4]).
Since the outbreak of the pandemic, scientists across the world have been indulged in the studies regarding the spread of the virus. Lin et al. [5] suggested the use of the SEIR (Susceptible-Exposed-Infectious-Removed) model for the spread in China and studied the importance of government-implemented restrictions on containing the infection. As the disease grew further, Ivorra et al. [6] suggested a θ-SEIHRD model that took into account various special features of the disease. It also included asymptomatic cases into account (around 51%) in order to forecast the total cases in China (around 168500). Giordano et al. [7] also suggested an extended SIR model called SIDHARTHE model for cases in Italy which was more customized for Covid-19 in order to effectively model the course of the pandemic to help plan a better control strategy.
Petropoulos and Makridakis [8] suggested the use of exponential smoothing method to model the trend of the virus, globally. Kumar et al. [9] gave a review on the various aspects of modern technology used to fight against COVID-19 crisis.
Apart from the epidemiological models, various data-oriented models were also suggested in order to model the cases and predict future cases for various disease outbreaks from time-to-time. Various time-series models were also suggested in order to model the cases and predict future cases. ARIMA and Seasonal ARIMA models are widely used by researchers in order to model and predict the cases of various outbreaks. In 2005, Earnest et al. [10] conducted a research to model and predict the cases of SARS in Singapore and predict the hospital supplies needed using this model. Gaudart et al. [11] modelled malaria incidence in the Savannah area of Mali using ARIMA. Zhang [12] compared Seasonal ARIMA model with three other time series models to compare Typhoid fever incidence in China. Polwiang [13] also used this model to determine the time-series pattern of Dengue fever in Bangkok.
For Covid-19 as well, various researchers tried to model the cases through ARIMA. Ceylan [14] suggested the use of Auto-Regressive Integrated Moving Average (ARIMA) model to develop and predict the epidemiological trend of Covid-19 for better allocation of resources and proper containment of the virus in Italy, Spain and France. Chintalapudi et al. [15] suggested its use for predicting the number of cases and deaths post 60-days lockdown in Italy. Fanelli and Piazza [16] analyzed the dynamics of Covid-19 in China, Italy and France using iterative time-lag maps. It further used SIRD model to model and predict the cases and deaths in these countries. Zhang et al. [17] developed a segmented poisson model to analyze the daily new cases of six countries in order to find a peak point in the cases.
Since the spread of the virus started to grow in India, various measures were taken by the Indian Government in order to contain it. A nationwide lockdown was announced on March 25th to April 14th, which was later extended to May 3rd. The whole country was divided into containment zones (where large number of cases were observed from a relatively smaller region), red zones (districts where risk of transmission was high and had higher doubling rates), green zones (districts with no confirmed case from last 21 days) and orange zones (which didn’t fall into the above three zones). After the further extension of the lockdown till May 17th, various economic activities were allowed to start (with high surveillance) in areas of less transmission. Further, the lockdown was extended to May 31st and some more economic activities have been allowed as per the transmission rates, which are the rates at which infectious cases cause new cases in the population, i.e. the rate of spread of the disease. This was further extended to June 8th, with very less rules and specially the states were given the responsibility of setting the lockdown rules. The air and rail transport became open for general public. Post June 8th, we see that the restrictions are nominal with even shopping malls and religious places open for general public. Now, the responsibility of imposing restrictions lies with the respective state governments.
On the other hand, Indian scientists and researchers are also working on addressing the issues arising from the pandemic, including production of PPE kits and tests kits as well as studying the behaviour of spread of the disease and other aspects of management. Various mathematical and statistical methods have been used for predicting the possible spread of Covid-19. The classical epidemiological models (SIR, SEIR, SIQR etc.) suggested the increasing trend of the virus and predicted the peaks of the pandemic. Early researches showed the pandemic to reach its peak by mid-May. They also showed that the basic reproduction number (R0) and the doubling rates are lower in India, with comparison to European nations and USA. A tree-based model was proposed by Arti and Bhatnagar [18] and Bhatnagar [19] in order to study and predict the trends. They suggest that lockdown and social-distancing in India has played a significant role to control the infection rates. But now, as the lockdown restrictions are minimal, the cases in India are growing at an alarming high rate. Chatterjee et al. [20] suggests growth of the pandemic through power law and its saturation at the later stages. Due to the complexities in the epidemic models of Covid-19, various researchers have been focusing on the data in order to forecast the future cases. Chatterjee et al. [20, 21] and Ziff & Ziff [22] suggest that after exponential growth, the total count follows a power regime of t3, t2, t and before flattening out, where ‘t’ refers to time. It can therefore be realized that there is an urgent need to model and forecast the growth of Covid-19 in India as the virus is in the growing stage here.
In India, the most affected states are Maharashtra with over 2,60,000 cases (as of 14th July, 2020), Tamil Nadu (around 1,42,000 cases), Delhi (around 1,13,000 cases) and Gujarat (around 42,000 cases). The greatest number of cases per million have been seen in the national capital of Delhi (5740 cases per million). (Refer [23] for population estimates). Many states and union territories like, Kerala, Karnataka, Andaman and Nicobar Islands, Daman and Diu etc. which had recovered from majority of the cases have experienced a second wave of infections. This might be attributed to decreased travel restrictions and minimal lockdown measures. In their research, Singh & Jadaun [1] studied the significance of lockdown in India and suggested that the new Covid-19 cases would stop by the end of August in India with around 350,000 total cases. While some states may see an early stopping of new cases such as Telangana (mid-June), Uttar Pradesh and West Bengal (July-end) etc., the badly affected states of Maharashtra, Tamil Nadu and Gujarat will achieve this by August-end.
Since a proven vaccine and medication is yet to be developed by the researchers then in such a scenario, modeling the present situation and forecasting the future outcome becomes crucially important in order to utilize our resources in the most optimal way. Therefore, the article aims to study the growth curve of Covid-19 cases in India and forecast its future course. Since the disease is still in its growing age and very dynamic in nature, no model can guarantee for perfect validity for future. We therefore, need to develop the understanding of the present situation of the pandemic.
In this article, we first study the growth curve using regression methods (exponential, linear and polynomial etc.) and propose an optimal model for fitting the cases till July 11th. Further, we propose the use of time-series models for forecasting the future observations on Covid-19 cases. Here, we reach the best-fitted ARIMA model for forecasting the Covid-19 cases. We also compare these results with Exponential Smoothing (Holt-Winters) model. This study will help us to understand the course of spread of SARS-CoV-2 in India better and help the government and the people to optimally use the resources available to them. As this study presents the evaluation of the pandemic, it may also help the practitioners and researchers to build further models/analysis with other co-factors.
2. Statistical Methodologies
In this section, we briefly present the statistical techniques used for analyzing the Covid-19 cases in India. Here, we used usual regression (exponential, polynomial), times series (ARIMA) and exponential smoothing models.
2.1 Exponential-Polynomial Regression
Regression is a statistical technique that attempts to estimate the strength and nature of relationship between a dependent variable and a series of independent variables. Regression analyses may be linear and non-linear. A regression is called linear when it is linear in parameters e.g. y = β0 + β1t+ ∈and y = β0 + β 1t+ β 2t2 + β 3t3+∈, ∈∼ N(0, α2), where y is response variable, t denotes the indepenet variable, β0 is the intercept and other βs are known as slopes.
A non-linear regression is a regression when it is non-linear in its parameters e.g.. In the beginning of the spread of a disease, we see that the new cases are directly proportional to the existing infected cases and may be represented by , where k is the proportionality constant. Solving this differential equation, we get that, at the beginning of a pandemic, Thus, at the beginning of a disease, the growth curve of the cases grows exponentially.
As the disease spreads in a region, governments start to take action and people start becoming conscious about the disease. Thus, after some time, the disease starts to follow a polynomial growth rather than continuing to grow exponentially.
In order to fit an exponential regression to our data, we linearize the equation by taking the natural logarithm of the equation and convert it to a linear regression in first order.
We estimate the parameters of a linear regression of order p as following-
Let the model of linear regression of order with ∈∼ N(0, σ2)and i = 1,2,..,N. Let represent the residual sum of square (RSS).
We get the best estimates of these coefficients using the ordinary least squares (OLS) method. We will use this technique of the OLS in order to estimate the coefficients of our proposed model. (Refer Montgomery et al. [24])
Since we know that the growth curve of the disease changes after some time point, exponential to polynomial, we propose to use the following joint regression model with change point µ, where we take and p is the order of the polynomial regression model and t stands for the time (an independent variable).
During the analysis, we found that a suitable choice of f2 (t) is a quadratic or a cubic model. Once the order of the polynomial is kept fixed, an optimum value of the change point can be obtained by minimizing the residuals/errors. We can obtain the OLS estimates of the parameters of the model (1) as given below.
The least square estimates (LSEs) of the parameters, Θ = {θ 1, θ 2,µ, β 0, β 1, β 2, β 3,……, βp}can be obtained by minimizing the residual sum of squares (RSS) as given by- where, and are the estimates value of yi from the exponential and polynomial regression models, respectively and is the size of the data set.
The estimates of, Θ = {θ 1, θ 2,µ, β 0, β 1, β 2, β 3,……, βp}can be obtained by using the OLS method. We suggest to use the following algorithm while µ is kept fixed.
In order to find the optimal value of µ, i.e. the turning point between the exponential and polynomial growth, we will use the technique of minimizing the residual sum squares in section 3.
We will use MAPE (Mean Absolute Percentage Error) in order to evaluate the performance of the model, as given by] where, yt is the observed value at time point t and ŷt is an estimate of yt.
In order to make the results easy to interpret, we will also use Accuracy (%) measure defined by
2.2 ARIMA Model
The Auto-Regressive Integrated Moving Averages method gauges the strength of one dependent variable relative to other changing variables. It is one of the most used time-series models in diverse fields of data analysis as it takes into account the changing of trends, periodic changes as well as random disturbances in the time-series data. It is used for both better understanding of the data as well as forecasting, see Brockwell et al. [25].
Autoregressive models (AR) is effectively merged with the Moving Averages models (MA) to formulate a useful time-series model, ARIMA model. The Autoregression (AR) element of the model shows a changing variable that regresses on its own prior values and the Moving Average (MA) element incorporates the dependency between an observation and a residual error from a moving average model applied to prior observations. However, this model can only be applied to stationary data. Since many real-life datasets consist of an element of non-stationarity, in order to model such datasets, ARIMA model was developed. This model is open for non-stationary data as the Integrated (I) factor of the model represents the differencing of raw observations to allow the time-series to become stationary.
Here, we may refer the reader to follow Box et al. [26] and Box et al. [27] for more details on ARIMA model, estimation and its application.
The general forms of AR(p)and MA (q) models can be respectively represented as the following equations: where Øs and θs are autoregressive and moving averages parameters, respectively, Yt represents value of time-series at time point t, εt represents the random disturbance at time point t and is assumed to be independently and identically distributed (i.i.d.) with mean 0 and variance σ 2.
The ARMA (p,q) model can be represented as- where α is an intercept.
The differenced stationary time-series can be modelled as an ARMA model in order to use ARIMA model on the time-series data, see Ceylan [14], He & Tao [28] and Manikandan et al. [29]. The ARIMA model is generally denoted as ARIMA (p,d,q), where p is the order of auto-regression, d is the degree of difference and q is the order of moving average.
The first step to model the time-series by ARIMA is to determine the time-series data for stationarity. The Augmented Dickey-Fuller (ADF) Test may be applied to determine if the time series after differencing is stationary or not. The ADF test is applied to test the null hypothesis for the presence of a unit root (which indicates non-stationarity of the series).
In order to deduce the ARIMA (p,d,q) model, we can proceed as follows:
We have the ARMA (p’,q) represented as follows (as per Eqn. 5) It can be equivalently written as- where L is the lag operator, such that-La(Yt) =Yt-a. Now, assume that the polynomial (1 − has a unit root (i.e. a factor of (1 − L)of multiplicity d. Then, Eqn. (6) can be re-written as- where,p ′= p −d.
This can be generalized as- This defines an ARIMA (p,d,q) process with drift .
The second step is to plot the graphs of the Autocorrelation function (ACF) and the Partial Autocorrelation Function (PACF) to determine the most-likely values of p and q.
The final step is to obtain the optimal values of p, d and q by using the AIC (Akaike Information Criterion), for more details see https://en.wikipedia.org/wiki/Akaike_information_criterion. These information criteria may be used for selecting the best fitted models. Lower the values of criteria, higher will be its relative quality. The AIC is given by- where K=number of model parameters, ℓ = maximized value of log − likeli hood function.
2.3 Exponential Smoothing
Exponential smoothing is one of the simple techniques to model time-series data where the past observations are assigned weights that are exponentially decreasing over time. We propose the following models, for modelling of Covid-19 cases (see Holt [30] and Winters [31]).
For single exponential smoothing, let the raw observations be denoted by {yt} and {st}denote the best estimate of trend at time t. Then, s0= y0, st= αyt+(1− α)(st−1) where α ε (0,1) denotes the data smoothing factor.
For double exponential (Holt-Winters) smoothing, let the raw observations be denoted by {yt}, smoothened values by {st}, and {bt} denotes the best estimate of trend at time t. Then, where α ∈ (0,1) denotes the data smoothing factor and β ∈ (0,1) denotes the trend smoothing factor. For the forecast at t =(N + m) days, FN+m is calculated by
3. Analysis of Covid-19 cases in India
For this study, we have used the data available at GitHub, provided by Centre for Systems Science and Engineering (CSSE) at John Hopkins University (see [32]). For this study, we use R software. (see R Core Team [33]).
3.1 Exponential-Polynomial Regressions
We know that at the beginning of the spread of the disease in India, the growth was exponential and after some time, it was shifted to polynomial. We first obtain optimum turning point of the growth, i.e. when did the growth rate of the disease shifted to polynomial regime from the exponential. We consider both quadratic and cubic regression model for second part of the data. We will also discuss the types of polynomial growth (with their equations) in India.
In order to find the turning point of the growth curve, we follow the Algorithm 1, given in the previous section. Using that, we evaluate the RSS for all the days (from March 4th) and find the date on which it is minimum. The change points of growth curve for cubic and quadratic regressions are presented in Figure 1 depending upon the size of the data set. From Figure 1, we can confirm that the growth rate of Covid-19 cases was exponential till April 5th and then after it follows the polynomial growth regime while we use the Covid-19 cases till May 2nd.
We call the region of exponential growth in India as Region I. The coefficients of the model are presented in Table 2.
We see that after the exponential regime (till April 5th), the growth curve follows a polynomial growth till May 2nd. After this, we again see a change in the behavior of the growth curve. In Tables 3, 4, 5 and 6, we try to model these growth curves through regression analysis.
Having evaluated the coefficients for various models (i.e. linear, quadratic and cubic) as well as the important statistics (i.e. R2 values, p-values of the models as well as individual coefficients and F-statistic), we will select the best fitting models. In order to select the best fitting models for Region II (April 6th to May 2nd); III (May 3rd to May 15th), IV (May 16th to May 31st) and V (June 1st to July 11th), we have the following steps. We select that model which has high R2 values, significant p-value, high F-statistic and where the p-values of all the variables are significant.
We see for Region II, from Table 3 that the linear model is having a relatively lower F-statistic and R2 values in comparison to the Quadratic and Cubic models. So, we eliminate the possibility of linear fitting. Further, we see that the p-values, F-statistics and the R2 values are quite significant in both Quadratic as well as the Cubic models. But, if we look at the individual p-values of the coefficients, we see that the individual p-values are not significant for the Cubic model. On the other hand, the individual p-values are significant for the Quadratic model. Thus, we can conclude that the Quadratic model is the best fitting model for Region II (April 6th to May 2nd).
For Region III, from Table 4, that all the three models have high F-statistic values, high p-values and high R2 values. But we notice that the coefficient individual p-values are not significant in both Quadratic and Cubic models. Thus, we conclude that the Linear model is the best fitting model for Region III (May 3rd to May 15th).
For Region IV, from Table 5, we see that the R2 values for all the models are very high. All the models also have significant p-values. The F-statistic of both Quadratic and Cubic models are also high. But, the coefficient individual p-values are not significant in the Cubic model. Thus, we conclude that the Quadratic model is the best fitting model for Region IV (May 16th to May 31st).
For Region V, from Table 6, we see that the R2 values of all the models are very high (quadratic, cubic, quartic and quintic models have exceptionally high). All the models also have significant p-values. The F-statistic of Quadratic, Cubic, Quartic and Quintic models is high. F-statistic value of Quartic model is the highest. The coefficient individual p-values of Quartic model are also significant. Thus, we conclude that the Quartic model is the best fitting model for Region V (June 1st to July 11th)
Note: For Region V, due to spike in the cases, we also checked the fitting of exponential regression in this region, see Table 7. The RSE for Exponential Regression is 4178 and MAPE (%) is 0.85%. Both of these values are quite larger than those of Quartic model. (Refer Table 9 for RSE and MAPE values of Quartic model in Region V). Thus, we conclude that Quartic model is the best-fitting model for Region V (1st June to 11th June).
All the ANOVA (Table 8) for Region II, III, IV and V suggest significant p-values for its coefficients and suggest that the models fit well the respective regions.
Thus, according to our study, the growth of the virus was exponentially increasing from March 4th to April 5th. Then after, the covid 19 cases grew by following a quadratic rate from April 6th to May 2nd. After May 3rd, we experienced a linear growth. But after May 15th to May 31st, we experienced a sudden rise in the rate of growth of the pandemic and have seen quadratic growth again. Further, for the period of June 1st to July 11th, we see experienced a quartic (4-degree polynomial) growth, which is very alarming. (See Table 9 for best fitted regression models). Figure 2 shows the best fitted regression models to the daily cumulative cases of Covid-19 in India from March 4th to July 11th.
3.2 Time series Models fitting
We use cumulative confirmed cases from March 4th to July 10th for fitting time series models and data from July 11th to July 14th is used for model validation.
First, we check the stationarity of the data using the ADF Test. ADF test statistic is 6.3915 with p-value 0.99, which indicates that the growth of Covid-19 cases is not stationary. The ARIMA models may be useful over the ARMA models. In order to make time series stationary, we may apply difference operator. In Figure 6, we can see that 2nd order differencing make the series almost stationary. Thus, we can take d =2. The ACF and PACF plots are shown in Figure 3. ACF plot indicates that there is high autocorrelation (more than 0.5) up to time lag 15 days.
We then obtain the optimal ARIMA parameters (p, d, q) by using the AIC. As suggested by Box et al. [26, 27], we take various possible combinations of p =1,2,…,5, q =1,2,…,5 and d =1,2 and compute the AIC.
According to the AIC, the ARIMA (5, 2, 5) is the best fitted model for the Covid-19 cases, India (see Table 10). Estimates of ARIMA (5, 2, 5) parameters and MAPE are shown in Table 11.
Estimates of the Holt-Winters exponential smoothing and exponential smoothing models are given in Table 12. According to the MAPE and accuracy measures, the ARIMA (5, 2, 5) is a better model than the Holt-Winters exponential smoothing and usual exponential smoothing models. From this, we can conclude that the ARIMA model gives the best fit for the cases of Covid-19, followed by Holt-Winters model. The forecasting values along with 95% confidence intervals are shown in Table 13.
We now present ARIMA and Holt-Winters models for 1st June to 10th July, last segment of the study.
We have seen that our time-series data is non-stationary and thus, we select the most optimal values of (p,d,q), which has the least AIC. ACF and PACF plots are given in Figure 5. According to AIC, ARIMA(5, 2, 3) is the best-fitting model for the time-series data from June 1st to July 10th, with AIC=634.18. Estimates of ARIMA (5, 2, 3) model with the corresponding MAPE and Accuracy are given in Table 14.
Estimates of the Holt-Winters exponential smoothing and exponential smoothing models are given in Table 15. According to the MAPE and accuracy measures, the ARIMA (5, 2, 3) is a better model than the Holt-Winters exponential smoothing and usual exponential smoothing models. From this, we can conclude that the ARIMA model also gives the best fit for the cases of Covid-19 from 1st June to 10th July, followed by Holt-Winters model. The forecasting values along with 95% confidence intervals are shown in Table 16 and Figure 4. We observe that the ARIMA model captures the trend well and estimates the cumulative cases properly.
4. Conclusions
From the regression analysis, we conclude that the spread of Covid-19 disease grew exponentially from March 3rd to April 5th. Further, from April 6th to May 2nd, the cases followed a quadratic regression. From May 3rd to May 15th, we see a linear growth of the pandemic with average daily cases of 3584. After May 15th to May 31st, we again saw a spike in the cases that lead to a quadratic growth of the pandemic. And, from June 1st to July 11th, we saw a major spike in the growth of the pandemic as it has followed quartic growth.
Verma et al. [21] showed the four stages of the epidemic, S1: exponential, S2: power law, S3: linear and S4: flat. We saw that the course of Covid-19 in India followed this regime till May 15th. But after the linear trend from May 3rd to May 15th, the spread has again reached the quadratic growth and from June 1st to July 11th, India is witnessing a quartic growth. This might be attributed to the relaxation of lockdown measures in the country. Though it was much likely that the cases would start to reduce post-linear stage growth as the total cases may start to follow a square root equation, i.e.. And this might lead to reduction in the daily number of cases (as ,) leading to flattening of the curve. But, due to reduced restrictions, we see a reverse trend, which might be alarming and we suggest the imposition of strict lockdown in order to reverse this trend of pandemic growth. If we continue to open our economy in this way, we might go back to the exponential growth of the pandemic and this would lead to huge destruction to human lives and cause a greater impact on our economy.
We also observe that some cities have been the hotspots of the disease, such as Delhi (more than 1,31,000 cases), Mumbai (more than 94,000 cases), Chennai (more than 78,000 cases), Thane (more than 63,000 cases) etc. as on 14th June, 2020. While the other states and cities have seen a slower growth of the pandemic, these cities have seen explosive growths. Due to the opening of air and rail transport in the country, the virus is likely to spread in the other regions as well as people from these cities (especially metro cities) are travelling to different states. Thus, it is highly advisable that the country should go back to its lockdown phase until we see reduction in trend.
In time series analysis, we conclude that the ARIMA (5, 2, 5) is the best fitting model for the cases of Covid-19 from 4th March to 10th July with an accuracy of 97.38%. The basic exponential smoothing is not very accurate for our case but we see that the Holt-Winters model is around 97.11% accurate. Both ARIMA (5, 2, 5) and Holt-Winters models suggest a rise in the number of cases in the coming days. We observed that both the ARIMA and Holt-Winters models capture the data well and the actual data from 11th July validates the forecasts well as they lie in the predicted confidence intervals. But, while validating the model, the actual values are always near to the Upper Confidence Limits, it might be possible that in further days, our model might underestimate the cases. This might be possible because of the changing trend of the growth of the pandemic in India.
Thus, we used segmented time-series models and took a data from 1st June to 10th July to build separate ARIMA and Holt-Winters models. We concluded that ARIMA (5, 2, 3) is the best fitting model for Covid-19 cases in the given time period with an accuracy of 99.86%. The basic exponential smoothing is not very accurate or this case as well but, the Holt-Winters model is around 99.78% accurate. We also observe that the ARIMA and Holt-Winters models capture the data well and the actual data from 11th July validates the forecasts and lies near to the estimates.
We may also conclude that the cases of Covid-19 will rise in the coming days and the situation may turn alarming if proper measures are not followed. Since the economic activities have started in the country, people need to be more careful while going out. And explosion of the pandemic in the whole country can cause a serious damage to human lives, healthcare system as well as the economy of the country. Thus, there is an urgent need of imposing strict lockdown measures to curb the growth of the pandemic. We must also learn to lead our lives by following all the precautions even if the lockdown restrictions are relaxed and the economic activities are resumed.
Data Availability
The data is available at GitHub, provided by John Hopkins University. 15.https://github.com/CSSEGISandData/COVID-19/blob/master/csse_covid_19_data/csse_covid_19_time_series/time_series_covid19_confirmed_global.csv
Acknowledgments
Dr. Vikas Kumar Sharma greatly acknowledges the financial support from Science and Engineering Research Board, Department of Science & Technology, Govt. of India, under the scheme Early Career Research Award (file no.: ECR/2017/002416).