Skip to main content

Prediction and analysis of COVID-19 daily new cases and cumulative cases: times series forecasting and machine learning models



COVID-19 poses a severe threat to global human health, especially the USA, Brazil, and India cases continue to increase dynamically, which has a far-reaching impact on people's health, social activities, and the local economic situation.


The study proposed the ARIMA, SARIMA and Prophet models to predict daily new cases and cumulative confirmed cases in the USA, Brazil and India over the next 30 days based on the COVID-19 new confirmed cases and cumulative confirmed cases data set(May 1, 2020, and November 30, 2021) published by the official WHO, Three models were implemented in the R 4.1.1 software with forecast and prophet package. The performance of different models was evaluated by using root mean square error (RMSE), mean absolute error (MAE) and mean absolute percentage error (MAPE).


Through the fitting and prediction of daily new case data, we reveal that the Prophet model has more advantages in the prediction of the COVID-19 of the USA, which could compose data components and capture periodic characteristics when the data changes significantly, while SARIMA is more likely to appear over-fitting in the USA. And the SARIMA model captured a seven-day period hidden in daily COVID-19 new cases from 3 countries. While in the prediction of new cumulative cases, the ARIMA model has a better ability to fit and predict the data with a positive growth trend in different countries(Brazil and India).


This study can shed light on understanding the outbreak trends and give an insight into the epidemiological control of these regions. Further, the prediction of the Prophet model showed sufficient accuracy in the daily COVID-19 new cases of the USA. The ARIMA model is suitable for predicting Brazil and India, which can help take precautions and policy formulation for this epidemic in other countries.

Peer Review reports


With the widespread of the new coronavirus, it has become a serious threat to the health of people worldwide. This new virus was later named Coronavirus disease 2019(COVID-19), a kind of respiratory infectious disease with lung inflammation [1,2,3].COVID-19 shows more special transmission characteristics than previous infectious diseases, which leads to its faster transmission speed, wider transmission range, higher transmission risk and rapid epidemic spread, posing a significant threat to global public life and security [4,5,6,7].

Since discovering COVID-19 cases in Wuhan, Hubei Province, in December 2019, the epidemic and its variants have spread rapidly throughout the world. So far, the COVID-19 pandemic has spread widely in 188 countries [5], and over 278 million cases and just under 5.4 million deaths have been reported globally 2 according to WHO. At present, the pandemic spread in most countries is still growing and has not been effectively controlled. As of December 28, 2021, the cumulative prevalence of new coronal pneumonia was most outstanding in the USA, followed by India Brazil. Due to the fastest spread of the new coronavirus, the cumulative number of confirmed cases and the daily number of new cases are still increasing in the above countries. Therefore, an analysis of the current cumulative and the new number of cases of COVID-19 has essential research implications for predicting its prevalence trends.

Different models have been used to predict COVID-19 prevalence and mortality rate in recent studies. For example, multiple linear regression [8], Artificial Neural Network [9], multilayer perceptron [10] grey prediction model [11], simulation model [12], Holt model [13], LSTM model [14], and support vector regression [15, 16]. However, the spread of epidemic disease is random and will be affected by many factors [17, 18]. A large number of studies all show that the effect is not best achieved if only a single prediction tool is utilized to predict trends. In addition, the above statistical model can predict the development trend of the epidemic in the medium and long term.

However, as time-series data, COVID-19 cases have some dynamic fluctuation trend in the various situation with epidemic prevention and control [19, 20], which is suitable for establishing a time series model for prediction, but the commonly used single time series analysis model is challenging to capture the nonlinear part of the COVID-19 data series.

Considering that there are many factors affecting the prevalence of COVID-19, these factors lead to the complex characteristics of nonlinear, random and periodic data [21,22,23], so it is necessary to establish a prediction model and compare the effectiveness of different models. Exploring the prediction effect of the time series and machine learning method on the future of the epidemic can find the potential infectious risk in advance [10], avoid the outbreak of the epidemic, and provide data support for the decision-making of the prevention and treatment of infectious diseases.

The Automatic Regressive Integrated Moving Average (ARIMA) model has some advantages in its simple structure and immediate applicability. The ARIMA model has been applied to the prediction and estimation of prevalent diseases, such as typhoid fever [24], tuberculosis [25], influenza [26] and COVID-19 [27, 28]. Since ARIMA methods do not contain much mathematics or statistics, but also are capable of correlating regulation with short-term changing trends in the time series. So, the model is more suitable for predicting the short-term epidemic diseases. The Prophet is an open-source model that can handle time-series data with the advantages of taking strong seasonal effects, missing data, outliers, and changes in trends [29, 30]. And it is currently useful for predictive analysis of COVID-19. What’s more, the SARIMA and Prophet models can be used to capture some periodic or seasonal changes, further find the nonlinear fluctuations of data, and improve the accuracy of prediction results [31].

Therefore, it is of great practical significance to predict the daily new cases and cumulative confirmed cases of COVID-19 all over the world. This study establishes the SARIMA, ARIMA and Prophet models to predict the daily new cases and cumulative confirmed cases of COVID-19 in the United States, Brazil and India in the next 30 days ( As shown in Fig. 1), and evaluate the prediction accuracy of the model to provide a further reference for the prediction and early warning of infectious diseases.

Fig. 1
figure 1

The proposed methodology of the COVID-19 forecasting model

Materials and methods

Data collection

This article is based on the official WHO website, and MS Excel 2019 was used to build a COVID-19 time-series database. To create a stable and effective ARIMA model, at least 30 observations are required [32]. Cumulative cases and daily confirmed cases from the three countries of the USA, India, and Brazil, as of May 1, 2020, through November 30, 2021, were selected for train data of the construction of disease prediction models and the cumulative cases and daily confirmed cases of those three countries during next 30 days (December 1, 2021, to December 30, 2021) will be forecasted by fitted models. A statistical description of this raw data is presented in Table 1. Forecast the model prediction performance for confirmed case data for the next month with 95% relative confidence intervals (December 1, 2021–December 30, 2021).

Table 1 Descriptive statistics on the cases of COVID-19 in USA, Brazil, and India

SARIMA and ARIMA model

ARIMA is a type of algorithm for the analysis and forecasting of time series data, namely the Box—Jenkin model, first proposed by Box and Jenkins in the 1970s [32]. The ARIMA (p, d, q) model is known as the differential autoregressive moving average model. Due to the seasonal feature of the raw data, the SARIMA model (seasonal autoregressive integrated moving average), as an extension of ARIMA, is also often used for time series forecasting after seasonal adjustment. Such model is to apply mathematical models to non-stationary time series after smoothing the data, which is used to estimate and extrapolate the state of something at some point in the future by analyzing the pattern of historical data and making future predictions based on that pattern and historical data from the past and the present [33]. The cumulative number of confirmed cases and daily new cases of COVID-19 is a random series with nonlinear or seasonal character, so the model can be considered suitable for forecasting. ARIMA simulates and estimates the state of something at some point in the future. The ARIMA model includes the following steps [34]: Step 1: Assessment of the model; Step 2: The model parameters were estimated; Step 3: Check the hypotheses of the model validation; Step 4: Modeling predictions.The structure of the ARIMA (p, d, q) model is Eq. (1).

$${y}_{t}={\varnothing }_{1}{y}_{t-1}+{\varnothing }_{2}{y}_{t-2}+\dots +{\varnothing }_{p}{y}_{t-p}+{e}_{t}-{\theta }_{1}{e}_{t-1}-{\theta }_{2}{e}_{t-2}-\dots -{\theta }_{q}{e}_{t-q}$$

In Eq. (1), ϕa(a = 1,2,…,p) and \(\theta\) b(b = 0,1,2,…,q) are parameters of the model. yt and ɛt represent the original value and arbitrary error at time step t. The arbitrary error represented by ɛt represents σ2 with zero mean and standard deviation. Taking the value q = 0 in Eq. (1) works as A.R. model with order p, and for p = 0, it becomes the M.A. model with q order. So (p, q) are both important factors to determine the ARIMA model.

The Prophet model

The Prophet is a powerful and fast open-source time series model developed by Facebook. which could well handle the impact of missing values and outliers in the time series on the prediction and is suitable for the prediction analysis of the COVID-19 epidemic[35,36,37]. They are combined in the following equation.

$$Y\left(t\right)=g\left(t\right)+s\left(t\right)+h\left(t\right)+{\varepsilon }_{t}$$

where \(Y\left(t\right)\) indicates the trend indicator data at time t; \(g\left(t\right)\) indicates the trend term and is the portion of the time series in which there is a non-cyclical trend of change; \(s\left(t\right)\) indicates the period term and is the portion of the time series that exhibits a periodicity of change; \(h\left(t\right)\) indicates a holiday term and is the portion of the sequence that is affected by holidays and, since data from this study do not have an effect of the holiday term in trend projections, this one was not considered; \({\varepsilon }_{t}\) It is an error term which accounts for any unusual changes not accommodated by the model. \({\varepsilon }_{t}\) denotes errors due to unusual changes.

Prophet uses the Fourier series to forecast the seasonality effects, and the seasonality models are specified as the periodic functions of t[38, 39]. The arbitrary smoothing of seasonal effects with a scaling time variable using Fourier series is represented as:

$$s\left(t\right)={\sum }_{n=1}^{\infty }{a}_{n}\mathrm{cos}\frac{2n\pi t}{p}+{b}_{n}\mathrm{sin}\frac{2n\pi t}{P}$$

where P is the period and, for a given value of N, to fit the seasonality model, the parameters a1, a2,…, an and b1, b2,…, bn need to be estimated.

Analytical tools and model evaluation

ACF and PACF test

The ACF is a complete autocorrelation function that provides us with the autocorrelation value for any sequence with lag values. In brief, it describes the degree of correlation between the current value of that sequence and its past value. PACF is a partial autocorrelation function. Rather than finding correlations of lags like ACF with the current, it finds correlations of the residuals with the next lag value. An ACF shows the linear relationship between the observations at time t and previous observations at time t − n. The ACF and PACF for a given time series X can be defined as:

$$\mathrm{PACF}\left({X}_{t},{X}_{t-2}\right)=\frac{Covariance\left({X}_{t},{X}_{t-2}/{X}_{t-1}\right)}{\sqrt{Variance\left({X}_{t}/{X}_{t-1}\right)}\sqrt{Variance\left({X}_{t-2}/{X}_{t-1}\right) }}$$

where in the ACF plot, n is the lag (or difference between \({X}_{t}\) and \({X}_{t-n}\)); in the PACF plot between observed values \({X}_{t} and {X}_{t-2}\), n = 2.

Performance indices

Three indexes were employed in accessing model fitting and forecasting efficiency: namely Root Mean Square Error (RMSE), Mean Absolute Error (MAE), and Mean Absolute Percentage Error (MAPE), and were applied to test the predictive accuracy of the developed models. Lower RMSE, MAE, and MAPE values indicate a better data fit. The formulations of these criteria are expressed Eqs. (6)–(8), respectively [40]. A logarithmic approach may be necessary to make the time series stationary after differencing. This approach takes the log value of each point, followed by differencing. Bayesian information criterion (BIC) is a class of information criteria to measure the goodness of fit of a statistical model. It builds on the concept of entropy and can weigh the complexity of the estimated model against the goodness of fit of this model to the data. This information helps assess the model's parameters and how well the model performed. In this study, to prevent the excessive model complexity caused by the excessive model accuracy. Therefore, the function sets the lower value.

$$\begin{array}{c}RMSE=\sqrt{\frac{SSE}{n}} =\sqrt{\frac{\sum_{i=1}^{n}{\left({Y}_{i}-{\overline{Y} }_{i}\right)}^{2}}{n}}\end{array}$$
$$\mathrm{MAE}=\frac{1}{n}\sum_{i=1}^{n}({Y}_{i}-{\overline{Y} }_{i})$$
$$\begin{array}{c}MAPE=\frac{100}{n}\times \sum_{i=1}^{n}\left|\frac{\left({Y}_{i}-{\overline{Y} }_{i}\right)}{{Y}_{i}}\right|\end{array}$$
$$\mathrm{BIC}=-2\mathrm{log}L\left(\widehat{\theta }\right)+n\mathrm{log}N$$

In Eq. (6), (7), and (8), where \({\mathrm{Y}}_{i}\) is the actual expected output, \({\overline{Y} }_{i}\) Is the model's prediction, i = 1…n and n is the number of observations. In Eq. (9),\(\mathrm{log}L\left(\widehat{\theta }\right)\) is the likelihood function, N is the number of observations, and n is the number of model parameters.

Data analysis

Since the new confirmed cases of COVID-19 has periodically or Seasonal characteristics. The SARIMA model and Prophet model were used to predict next 30 days COVID daily new cases and comfirmed cases data, and the Prophet and SARIMA model were constructed for the prediction of daily cumulative cases. The three models are used for the forecast and simulations of this study based on R 4.1.1 software with forecast and prophet package. Before applying the prediction model, we use logarithmic conversion to process the original data to make the time series more stable and weaken the collinearity of the model, so as to improve the accuracy of prediction.Due to the periodicity of daily new cases, the seasonal components are eliminated,Considering that the daily number of new cases in COVID-19 has the characteristics of periodicity and seasonality,hence, the ARIMA and Prophet model are constructed for the cumulative confirmed case data and, in addition, the SARIMA and prophet model are applied for the daily new confirmed cases.


COVID-19 cases description

This paper aims to perform a statistical, observational and predictive analysis of a cumulative confirmed case and new confirmed cases dataset from COVID-19 in three countries (USA, Brazil, and India) by ARIMA, SARIMA and prophet model. The COVID-19 data set is divided into 579 days of training samples and the next 30-day of prediction samples (Figs.2, 3). All statistical procedures were performed on the transformed COVID-19 data.

Fig. 2
figure 2

Daily variation of cumulative confirmed cases in USA, Brazil and India from May 1, 2020, to November 30, 2021

Fig. 3
figure 3

Comparison of daily new confirmed cases in USA, Brazil and India from May 1, 2020, to November 30, 2021

The cumulative number of confirmed cases in the USA is much higher than that in India and Brazil, showing a higher trend than India, with Brazil showing the slowest trend (Fig. 2). The analysis of new cases in three countries shows mainly three peaks of daily new cases in the USA, and they last for a long time. There is a marked peak in daily changes in new cases in India, while the overall trend in new cases in Brazil is relatively stable (Fig. 3).

The increasing trend of cumulative confirmed cases and new confirmed case data could provide more evidence to help citizens understand the time-series variation of COVID-19 in various countries. In addition, the daily new confirmed cases are different from the cumulative confirmed cases with a positive growth trend, which has unstable periodic fluctuations that increase or decrease and may occur seasonally. So we will specifically deal with this part of the period characteristics through the prophet and SARIMA models to facilitate the follow-up prediction.

Model construction and training

In constructing the ARIMA and SARIMA models, the autocorrelation coefficient ACF and the partial autocorrelation coefficient PACF are analyzed for the smooth time series, respectively. In addition, the indicator with the lowest BIC value was used as the best parameter for the ARIMA model. The optimal combination of model parameters for each indicator is shown in Table 2.

Table 2 Determine the best parameters of the SRIMA and ARIMA model

In the prophet model, because the data may have hidden periodicity, such as daily new case data, we mainly adjust the daily and seasonal cycle parameters and growth trend function to predict better. The data and relative and absolute errors of model fitting or prediction are shown in (Additional file 1).

Model fitting of COVID-19 cumulative cases and daily new cases

Based on the ARIMA and Prophet methods, we trained and studied the statistical data of the cumulative confirmed cases of COVID-19 in the USA, India and Brazil. The learning results are shown in Fig. 4.

Fig. 4
figure 4

ARIMA and Prophet model fitting for new cumulative cases in USA, Brazil and India from May 1, 2020, to November 30, 2021

In general, when the data fluctuation is slight, these two mos' training and learning results well captured the development trend of confirmed cases with good performance. The fitting results for the cumulative cases of the COVID-19 in the USA show that the average absolute error(MAE) of the ARIMA model is 1847.39, and its mean average relative error(MRE) is 0.00003 is lower than the prophet modle (0.02057). In India, the fitting average absolute error of the ARIMA model is 52.3197 and the mean average relative error is lower than 0.05, and the mean average relative error of Prophet model is also lower than 0.0005 (as shown in Additional file 1). In the training results in Brazil, the average absolute errors are 94.22 (ARIMA) and 2370.45 (Prophet), respectively, and the mean of average relative errors of ARIMA is lower than 0.0005, indicating that the fitting effect is better than that of the Prophet model (0.00001). This result verified that the ARIMA method has a particular advantage in learning the training set of cumulative confirmed cases with slow positive growth.

Similarly, we used the SARIMA method and the Prophet method to learn the data of daily new confirmed cases, and the results are shown in Fig. 5. For this kind of fluctuating data, these two methods are better to capture the fluctuation and trend of the data. Both performed better, especially in studying new cases in India (compared to the learning results in the USA and Brazil).

Fig. 5
figure 5

SARIMA and Prophet model fitting for daily new cases in USA, Brazil and India from May 1, 2020, to November 30, 2021

The MAE value of the Prophet model in the USA, India, Brazil is 7118.96, 2331.66 and 4819.42, which is lower than the ARIMA model (7877.08, 2847.32, 5661.42).

This indicates that the effect of the Prophet model is better than that of the SARIMA model in training data. And it is not difficult to find that the Prophet model has more advantages in learning the training time series data with large fluctuations through horizontal comparisons, such as the daily new case.

Model prediction of COVID-19 cumulative cases and daily new cases

The ARIMA and the Prophet model were used to predict the cumulative confirmed cases of COVID-19 in the next 30 days in three countries, as shown in Fig. 6. Obviously, the prediction trendy for the cumulative cases in the ARIMA model in Brazil and India is more similar to the actual data, which can be in line with the actual trend with a smaller relative error than the Prophet model (Additional file 1).

Fig. 6
figure 6

The prediction of the ARIMA and Prophet model for cumulative confirmed cases in the USA, Brazil and India from December 1, 2021, to December 30, 2021

While in the prediction of new cumulative cases and daily new cases the USA(Fig. 7), the Prophet model is closer to the actual change, unlike the trend in its prediction of daily cumulative COVID-19 confirmed cases and new cases in India and Brazil.

Fig. 7
figure 7

The prediction of the ARIMA and Prophet model for daily new cases in the USA, Brazil and India from December 1, 2021, to December 30, 2021

Accuracy evaluation

It can be seen from Table 3 that the values of RMSE and MAE in the Prophet model for predicting daily new cases and cumulative cases in the USA are lower than those in the SARIMA model. The training model’s R2 is similar to the R2 value of the SARIMA model, which is basically between 0.8 and 1.0, indicating that the prediction accuracy of the prophet model is higher and can be applied to the actual prediction of COVID-19. In the training dataset for daily new cases, the SARIMA model was over-fitting with poor generalization ability, whereas the Prophet model was fitted the cases and captured the seasonality hidden.

Table 3 Accuracy evaluation of ARIMA, SARIMA and Prophet on fitting and forecasting COVID-19 in USA, Brazil and India

But the RMSE, MAE and MAPE values of the cumulative cases were mostly higher in the Prophet model than the ARIMA model in India and Brazil's analysis, which revealed that the ARIMA model has a better ability to fit and predict the COVID-19 cumulative data with a positive growth trend. It should be noted that the fitting and predicted R2 value of the model for daily cumulative confirmed cases is close to 1 because its actual value is too large (tens of millions of cases).


With the recent outbreak of a new strain from Omicron, the incidence of a COVID-19 outbreak has reached a new height, increasing the burden on health systems in various countries. But according to the global confirmed case data, the Omicron infection cases have not peaked. According to USA media and health experts, the epidemic situation in the USA maybe even more acute in February. Omicron will continue to ravage through the states of the USA, with a proliferation of fatalities, a severe loss of hospital personnel, a shortage of labor, and several other social and economic confusion in an orderly fashion. The Brazilian Ministry of Health recently announced that multiple cases of Omicron Ba had been diagnosed in the country Cases of BA.2 subtype infection in Omicron. The new epidemic outbreak makes infection control for new strains of coronal pneumonia a tough challenge. Therefore, the construction of the prediction model of the epidemic trend of COVID-19 in most countries provides an essential reference for the prevention and control of the epidemic situation of COVID-19 and provides a decision-making basis for the prevention and control of the COVID-19 epidemic situation in the world.

In the model fitting and prediction stage, in comparing the prediction models between new cases and cumulative cases in different countries, the performance of the two models is better, but there are some differences in various country epidemics.

First of all, the prophet model has a unique advantage in evaluating the daily COVID-19 new cases and has higher accuracy for epidemic data with large fluctuations. This can be explained by The prophet model has better performance than ARIMA in predicting daily new cases and cumulative confirmed cases in the USA, so we can make a preliminary conclusion that when the data fluctuates greatly, its periodic characteristics become more obvious so that the prophet algorithm for decomposition calculation can capture it. This conclusion shows that the Prophet model is based on time series decomposition and machine learning fitting, which can automatically estimate long-term predictions for a wide range of forecasting problems, including marking out possible issues and getting predicted results in a faster time for further investigation by medical staff.

Secondly, in our current work, we have proved a seven-day period in daily new case data from different countries when fitting the SARIMA model. Previous studies used some regression techniques to eliminate time changes in effects and Poisson regression, binomial regression models, normalization, random forests, and other models[41,42,43,44,45,46]. However, the effect of these models on capturing the accuracy of the seasonal characteristics of events is not stable. In contrast, the analysis results of the SARIMA model are consistent with the actual data trend of COVID-19.

In addition, in this study, the prediction effect of the ARIMA and SARIMA model for India and Brazil with the accurate trend is better than the prophet model, indicating that this model may be more suitable for countries with less fluctuation of epidemic situation. It may be because the daily new cases of COVID-19 are greatly affected by prevention and control strategies and intervention policies leading to the number of cases changes fluctuating significantly, which is not easy to be well-identified by most models.

What's more, the cumulative case data shows a linear change, which can be predicted by some models like ARIMA. It is also why the ARIMA model, as a traditional time series tool, effectively predicts cumulative confirmed cases in three countries. The results further show that the epidemic situation of COVID-19 has a rapidly changing trend, which may be related to the incubation period of the disease and the prevention and control measures taken by the country [references], which is helpful to provide the basis for epidemic prevention and control for health management departments.


This study introduced the epidemic situation of COVID-19 in the USA, Brazil, and India. The persistent trend and scope of the epidemic were estimated using ARIMA, SARIMA and Prophet models. And found that the Prophet model is more suitable for daily new cases with large fluctuations and has its unique advantages compared with the commonly used ARIMA model. Therefore, we further prove that the three models for the number of cases of COVID-19 can be used to accurately predict the development of infectious diseases and find their potential periodicity. The current work can provide governments regarding make emergency macroeconomic strategies and the allocation of medical resources, regulate social production activities, even give a reference for national economic development, and offer more basis and value for the prediction of COVID-19 epidemics.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available in the [WHO] repository,



Coronavirus disease 2019


Automatic regressive integrated moving average


Autocorrelation function


Partial autocorrelation function


Root mean square error


Mean absolute error


Mean absolute percentage error


  1. Wang L, Wang Y, Ye D, Liu Q. Erratum to “A review of the 2019 Novel Coronavirus (COVID-19) based on current evidence” [International Journal of Antimicrobial Agents 55/6 (2020) 105948]. Int J Antimicrob Agents. 2020;56(3): 106137.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Fauci AS, Lane HC, Redfield RR. Covid-19—navigating the uncharted. Mass Medical Soc. 2020;382:1268–9.

    CAS  Google Scholar 

  3. Wu Z, McGoogan JM. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72 314 cases from the Chinese Center for Disease Control and Prevention. JAMA. 2020;323(13):1239–42.

    Article  CAS  PubMed  Google Scholar 

  4. Bavel JJV, Baicker K, Boggio PS, Capraro V, Cichocka A, Cikara M, Crockett MJ, Crum AJ, Douglas KM, Druckman JN. Using social and behavioural science to support COVID-19 pandemic response. Nat Hum Behav. 2020;4(5):460–71.

    Article  PubMed  Google Scholar 

  5. Peeling RW, Olliaro PL, Boeras DI, Fongwen N. Scaling up COVID-19 rapid antigen tests: promises and challenges. Lancet Infect Dis. 2021;21(9):e290–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Cheng LL, Guan WJ, Duan CY, Zhang NF, Lei CL, Hu Y, Chen AL, Li SY, Zhuo C, Deng XL, et al. Effect of recombinant human granulocyte colony-stimulating factor for patients with coronavirus disease 2019 (COVID-19) and lymphopenia: a randomized clinical trial. JAMA Intern Med. 2021;181(1):71–8.

    Article  CAS  PubMed  Google Scholar 

  7. Covid C, Team R, Covid C, Team R, Bialek S, Bowen V, Chow N, Curns A, Gierke R, Hall AJM et al: Geographic differences in COVID-19 cases, deaths, and incidence—United States. Mmwr 2020, 69(15):465.

  8. Labrague LJ, De Los Santos JAA. COVID-19 anxiety among front-line nurses: predictive role of organisational support, personal resilience and social support. J Nurs Manag. 2020;28(7):1653–61.

    Article  PubMed  Google Scholar 

  9. Roustaei N, Allahyari E. COVID-19 screening: use of an artificial neural network. BMJ Support Palliat Care. 2021;31(12):003391.

    Google Scholar 

  10. Oshinubi K, Amakor A, Peter OJ, Rachdi M, Demongeot J. Approach to COVID-19 time series data using deep learning and spectral analysis methods. Aims Bioeng. 2022;9(1):1–21.

    Article  CAS  Google Scholar 

  11. Zhang J, Jiang Z. A new grey quadratic polynomial model and its application in the COVID-19 in China. Sci Rep. 2021;11(1):1–27.

    CAS  Google Scholar 

  12. Yang W, Cowling BJ, Lau EH, Shaman J. Forecasting influenza epidemics in Hong Kong. PLoS Comput Biol. 2015;11(7): e1004383.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Martinez EZ, Aragon DC, Nunes AA. Short-term forecasting of daily COVID-19 cases in Brazil by using the Holt’s model. Rev Soc Bras Med Trop. 2020;53: e20200283.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Bhimala KR, Patra GK, Mopuri R, Mutheneni SR. Prediction of COVID-19 cases using the weather integrated deep learning approach for India. Transbound Emerg Dis. 2022;69(3):1349–63.

    Article  CAS  PubMed  Google Scholar 

  15. Hazarika BB, Gupta D. Modelling and forecasting of COVID-19 spread using wavelet-coupled random vector functional link networks. Appl Soft Comput. 2020;96:

    Article  PubMed  PubMed Central  Google Scholar 

  16. Gupta AK, Singh V, Mathur P, Travieso-Gonzalez CM. Prediction of COVID-19 pandemic measuring criteria using support vector machine, prophet and linear regression models in Indian scenario. J Interdiscip Math. 2021;24(1):89–108.

    Article  Google Scholar 

  17. Oshinubi K, Rachdi M, Demongeot J. Analysis of reproduction number R(0) of COVID-19 using current health expenditure as gross domestic product percentage (CHE/GDP) across countries. Healthcare (Basel, Switzerland). 2021;9(10):1247.

    Google Scholar 

  18. Dorantes-Acosta E, Avila-Montiel D, Klunder-Klunder M, Juarez-Villegas L, Marquez-Gonzalez H. Survival and complications in pediatric patients with cancer and COVID-19: a meta-analysis. Front Oncol. 2020;10: 608282.

    Article  PubMed  Google Scholar 

  19. Kupferschmidt K. After Omicron, some scientists foresee “a period of quiet.” Science. 2022;375(6579):366–7.

    Article  CAS  PubMed  Google Scholar 

  20. Normile D. As normalcy returns, can China keep COVID-19 at bay? Science. 2020;368(6486):18–9.

    Article  CAS  PubMed  Google Scholar 

  21. Oshinubi K, Rachdi M, Demongeot J. Analysis of reproduction number R0 of COVID-19 using current health expenditure as gross domestic product percentage (CHE/GDP) across countries. Healthcare (Basel). 2021;9(10):1247.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Tan AX, Hinman JA, Abdel Magid HS, Nelson LM, Odden MC. Association between income inequality and county-level COVID-19 cases and deaths in the US. JAMA Netw Open. 2021;4(5): e218799.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Zheng Z, Peng F, Xu B, Zhao J, Liu H, Peng J, Li Q, Jiang C, Zhou Y, Liu SJ. Risk factors of critical & mortal COVID-19 cases: A systematic literature review and meta-analysis. J Infect. 2020;81(2):e16–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Zhang X, Liu Y, Yang M, Zhang T, Young AA, Li X. Comparative study of four time series methods in forecasting typhoid fever incidence in China. PLoS ONE. 2013;8(5): e63116.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Zheng YL, Zhang LP, Zhang XL, Wang K, Zheng YJ. Forecast model analysis for the morbidity of tuberculosis in Xinjiang, China. PLoS ONE. 2015;10(3): e0116832.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. He Z, Tao H. Epidemiology and ARIMA model of positive-rate of influenza viruses among children in Wuhan, China: a nine-year retrospective study. Int J Infect Dis. 2018;74:61–70.

    CAS  PubMed  Google Scholar 

  27. Ceylan Z. Estimation of COVID-19 prevalence in Italy, Spain, and France. Sci Total Environ. 2020;729: 138817.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Demongeot J, Oshinubi K, Rachdi M, Hobbad L, Alahiane M, Iggui S, Gaudart J. Ouassou IJJMCS: The application of ARIMA model to analyze COVID-19 incidence pattern in several countries. J Math Comput Sci. 2021;12:Article ID 10.

    Google Scholar 

  29. Anastassopoulou C, Russo L, Tsakris A, Siettos C. Data-based analysis, modelling and forecasting of the COVID-19 outbreak. PLoS ONE. 2020;15(3): e0230405.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Peter OJ, Shaikh AS, Ibrahim MO, Nisar KS, Baleanu D, Khan I, Abioye AI. Analysis and dynamics of fractional order mathematical model of COVID-19 in Nigeria Using Atangana-Baleanu Operator. Comput Mater Continua. 2021;2:1823–48.

    Article  Google Scholar 

  31. Abioye AI, Umoh MD, Peter OJ, Edogbanya HO, Oguntolu FA, Kayode O, Amadiegwu SJ. Forecasting of COVID-19 pandemic in Nigeria using real statistical data. Commun Math Biol Neurosci. 2021;2021:Article ID 2.

    Google Scholar 

  32. Bartholomew DJ, Box GEP, Jenkins GM. Time Series Analysis Forecasting and Control. Operat Res Quarter. 1971;22(2):199.

    Article  Google Scholar 

  33. Demongeot J, Oshinubi K, Rachdi M, Hobbad L, Alahiane M, Iggui S, Gaudart J, Ouassou I. The application of ARIMA model to analyze COVID-19 incidence pattern in several countries. J Math Comput Sci. 2021;12:Article ID.

    Google Scholar 

  34. Schaffer AL, Dobbins TA, Pearson SA. Interrupted time series analysis using autoregressive integrated moving average (ARIMA) models: a guide for evaluating large-scale health interventions. BMC Med Res Methodol. 2021;21(1):58.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Taylor SJ, Letham B. Forecasting at scale. Am Stat. 2018;72(1):37–45.

    Article  Google Scholar 

  36. Fang WX, Lan PC, Lin WR, Chang HC, Chang HY, Wang YH. Combine facebook prophet and LSTM with BPNN forecasting financial markets: the morgan Taiwan index. 2019 International Symposium on Intelligent Signal Processing and Communication Systems (ISPACS). 2019: 1-2.

  37. Hasan MN, Haider N, Stigler FL, Khan RA, McCoy D, Zumla A, Kock RA, Uddin MJ. The global case-fatality rate of COVID-19 has been declining since May 2020. Am J Trop Med Hyg. 2021;104(6):2176–84.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Zhao N, Liu Y, Vanos JK, Cao G. Day-of-week and seasonal patterns of PM2 5. concentrations over the United States: time-series analyses using the Prophet procedure. Atmos Environ. 2018;192:116–27.

    Article  CAS  Google Scholar 

  39. Papacharalampous GA, Tyralis H. Evaluation of random forests and Prophet for daily streamflow forecasting. Adv Geosci. 2018;45:201–8.

    Article  Google Scholar 

  40. Becerra M, Jerez A, Aballay B, et al. Forecasting emergency admissions due to respiratory diseases in high variability scenarios using time series: a case study in Chile. Sci Total Environ. 2020;706: 134978.

    Article  CAS  PubMed  Google Scholar 

  41. Dhungana HN, Ghimire S. Commentary: statistical modeling for the prediction of infectious disease dissemination with special reference to COVID-19 spread. Front Public Health. 2021;9: 735857.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Benz LS, Lopez MJ. Estimating the change in soccer’s home advantage during the Covid-19 pandemic using bivariate Poisson regression. Adv Stat Anal. 2021;27(7):1–28.

    Google Scholar 

  43. Murakami D, Matsui T. Improved log-Gaussian approximation for over-dispersed Poisson regression: application to spatial analysis of COVID-19. PLoS ONE. 2022;17(1): e0260836.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Lu JQ, Lu JY, Wang W, Liu Y, Buczek A, Fleysher R, Hoogenboom WS, Zhu W, Hou W, Rodriguez CJ, et al. Clinical predictors of acute cardiac injury and normalization of troponin after hospital discharge from COVID-19. EBioMedicine. 2022;76: 103821.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Melin P, Sánchez D, Monica JC, Castillo O. Optimization using the firefly algorithm of ensemble neural networks with type-2 fuzzy integration for COVID-19 time series prediction. Soft comput. 2021;13(1):1-38.

    Google Scholar 

  46. Fernandes FT, de Oliveira TA, Teixeira CE, Batista AFM, Dalla Costa G, Chiavegatto Filho ADP. A multipurpose machine learning approach to predict COVID-19 negative prognosis in Sao Paulo, Brazil. Sci Rep. 2021;11(1):3343.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We thank the medical staff for their outstanding contribution and sacrifice during the COVID-19 pandemic.


This study was funded by the grants from the National Natural Science Foundation of China (12031010) Special Grant for the Prevention and Control of Infectious Diseases (2018ZX10713003). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations



Each author is expected to have made substantial contributions to the conception and design of the work; Yanding Wang, Ding Wang and Zehui Yan: the acquisition, analysis; DW, ZL, LZ, XG and MY: interpretation of data; YW, DW and ZY: the creation of software used in the work; YW, WZ, YW, ZY and DW have drafted the work and substantively revised it. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Wenyi Zhang or Yong Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

The results of predictive model.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, Y., Yan, Z., Wang, D. et al. Prediction and analysis of COVID-19 daily new cases and cumulative cases: times series forecasting and machine learning models. BMC Infect Dis 22, 495 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Prophet
  • COVID-19
  • Epidemiological dynamics prediction