Coronavirus outbreak (COVID- 19) was first appeared in the Wuhan city, Hubei Province, in China. The virus was identified in the first half of January 2020 [1] . The epidemiological features of the disease are still unknown, and the number of total cases and the number of total deaths varies day by day. When the rapid spread and serious consequences of the disease were observed, precautions were taken and positive cases began to be recorded after the second half of January. Because the number of total cases and the number of total deaths are used for examining the course of the outbreak, modelling of these indicators are the important issue. The model results are valuable for determining the appropriate preventions. In the literature, it was observed that some models were used for short-term forecasts considering the previous data. It is possible to reduce the number of total cases and total deaths by taking necessary precautions. At this point, it is very important to make future estimates.
In this study, it was aimed to examine the course of outbreak in China with alternative statistical models by taking into the number of total cases and the number of total deaths for the day it became clear until April 20 and also we want to select the best model for defining with model selection criteria. For modelling Weibull, Negative Exponential, Von Bertalanffy, Janoscheck, Lundqvist-Korf and Sloboda models were used and AR, MA, ARMA, Holt, Brown and Damped models were used for time series.
The data was obtained between January 22 and April 21, 2020 in China records. The number of total cases and the number of total deaths were used for the calculations. The reasons for choosing China for modeling are that, it is the first country to fight the outbreak, therefore it is possible to observe the natural course of the outbreak and it is the country that struggles for the longest time.
The number of total cases and the total deaths of coronavirus disease in the world were recorded daily by following the web address https://www.worldometers.info/coronavirus/.
For coronavirus outbreak, the change in the number of total cases and the number of total deaths over time was analyzed separately. In the analysis of time-total case and time-total deaths, six different nonlinear models were used including Weibull, Negative exponential, Von Bertalanffy, Janoscheck, Lundqvist-Korf and Sloboda models (Table 1 ) [2] . For time series analysis, Box-Jenkins and exponential smoothing methods were used (Table 2 and  Table 3 ).
The models are examined in two subsections: first group models are Growth Curves and second group models are Time Series. The growth models descriptive equations are given in Table 1 [3] [4] [5] [6] [7] [8] . 
In Table 1 , Y t is the observed dependent variable as the number of total cases and the number of total death and t is the independent variable. In our models, t is the day. The A term is the is the asymptotic limit of the number of total cases and the number of total deaths as time goes to infinity, B is the proportion of the number of total cases to the number of total deaths. It is the proportion of the number of total cases to the number of total deaths, obtained after the estimated case / deaths with the initial value of time (days), to the highest number of total cases / deaths. The k term is the proportion of the maximum increase rate to the highest number of cases or death. The ࢽ , c, d are the changing points that occurs when the change in the estimated increase rate goes from increase to decrease. 
is the p. degree of autoregressive series [9] . MA (q) refers to the moving average model of order q. In this series,
is the White noise series [10] . ARMA (p, q)) model is expressed by both AR (p) and MA (q) processes [11] . 
. CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review)
The copyright holder for this preprint this version posted May 1, 2020. 
is the trend predicted value, p is the number of forecasting periods and ‫ݕ‬ ത ௧ ା is the forecasting value after p period [12] . In Damped trend method, if 0<<1, the trend is damped, if ߮ = 1, the equations become identical to the Holt's Linear Trend method. Tashman and Kruk (1996) defined that there may be value in allocating 1 , if applied in series with a strong tendency, with exponential trend [13] . The Brown's Single Parameter Linear Exponential Smoothing Model is more suitable, if there is an increasing or decreasing trend in the time series. In this model, the initial equations ‫ݕ‬ ௧ ଵ and ‫ݕ‬ ௧ ଶ are obtained by single exponential smoothing and double exponential smoothing, respectively [14] . For the estimation of post m process, the equation is given in the below [15] .
The exponential smoothing method is the method in which the estimates are constantly updated, taking into account the recent changes in the data [16] . In these methods, the weighted average of past period values is calculated and taken as the estimated value of future periods.
Estimation accuracy of the applied methods were evaluated with R 2 , Pseudo R 2 and BIC. Bayesian information criterion (BIC) was developed by Gideon E. Schwarz (1978) , who gave a Bayesian argument for adopting it [17] .
is the error variance . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted May 1, 2020. .
The data used in the study were the observations recorded between January 22 and April 21.
The parameters estimates that estimated from 91 daily number of cases with nonlinear mathematical models between 22 January-15 April 2020 in China were presented in Table 4 . Coefficient of determination (R 2 ), Pseudo R 2 and Mean Square Error (MSE) statistics were used to compare models. When Table 4 was examined, it was seen that the R 2 and MSE values of Weibull, and Janoscheck models were equal. The MSE of Sloboda model was slightly smaller than two model but R 2 was equal. The Sloboda model can be considered the most suitable model, through smaller MSE value, larger Pseudo R 2 value. The Weibull and Janoscheck models can also be chosen as an alternative model. 2559449.5 0.997
The prediction curves of growth models are given in Figure 1 . Box-Jenkins methods and exponential smoothing methods were used among time series models for total number of cases. Autocorrelation (ACF) and partial autocorrelation (PACF) graphs of the series were examined. When the ACF and PACF graphs are examined in Figure  2 , the first degree difference had been taken since the series were not stationary at the level. But the stationary assumption is not provided yet. The difference from the second degree was taken and the series became stationary. According to the ACF and PACF charts, the series quickly approached zero after the first delay in the ACF graph. In this case, since p = 0, d = 2 and q = 1, it was modeled by the integrated first degree moving averages method. In other words, the most suitable time series method was the ARIMA (0,2,1) model. In addition, exponential smoothing methods were used and the model performances were given in Table  5 .
. CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted May 1, 2020. 
The performance of the model was given in Table 6 and it was seen that model's estimations were successful like growth models. The parameter estimates of the ARIMA (0,2,1) model were given in Table 7 . The ARIMA (0,2,1) model in this study can be written as in the below.
The forecasts from ARIMA (0,2,1) model for 15 day were given in Table 5 . . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted May 1, 2020. As seen in Table 5 , the number of total cases continues increasingly, albeit at low speed. The number of total cases is forecasting to be 83311 on the 105th day of the outbreak, that is, on 5 May 2020.
Observed and predicted values of the total cases were shown in Figure 3 . The prediction and forecasting results of growth models for the number of total deaths
The parameters and the selection criteria of the growth models for total deaths were given in Table 8 . When Table 8 was examined regarding the number of total deaths in China, the most suitable models were Lundqvist-Korf and Sloboda growth models, respectively. The R 2 values of these models were found highest as 0.956 and also MSE values were of these are lower than . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted May 1, 2020. . the others. Lundqvist-Korf model can be considered as the most suitable model, since it had smaller MSE. The prediction curves of growth models were given in Figure 4 . The most suitable time series model was found as Holt linear trend exponential smoothing model among time series models for the number of deaths. The compliance of the models were given in Table 9 and it was seen that the predictions were successful as growth models. The parameter estimation of the Holt linear trend exponential smoothing model was presented in Table 10 . The observed and predicted values were given in Figure 5 . The forecasts of the number of total deaths for exponential smoothing model by using Holt linear trend exponential smoothing model for 15 day were given in Table 11 . The rate of increase in the number of deaths in China was decreasing and it was predicted that it will be between 3343-3355 in the period between April 16 and May 1, with a slight increase (Table  11 ). . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint this version posted May 1, 2020. The Holt linear trend exponential smoothing curve for the exponential smoothing model was given in Figure 5 . 
The best time series model was ARIMA (0,2,1) model for the number of cases. In a study, ARIMA model was used on the daily prevalence data of COVID-2019 from January 20, 2020 to February 10, 2020 and ARIMA (1,2,0) and ARIMA (1,0,4) models were obtained [18] . Logistic, Bertalanffy and Gompertz models were used to estimate the number of cases and deaths on COVID-19 disease in different regions in China before by Jia and et al. According to their study, the Logistics model was reported to be better than the others by the R 2 criterion conducted extensive research in quasi-experimental analysis method in various provinces in China and investigated the relationship between population and number of outbreak cases [19, 20] . In their study, they found that the correlation coefficients of the relationship between the population and the number of cases differed among regions. They have observed that the number of cases was high in regions with high population and there was a high correlation between them. They stated that factors such as immigrants, tourism and mobility plays an important role in this. Also the authors determined the number of cases with the epidemic growth model Fan et al. [19] . Roosa et al. (2020) analyzed the number of cases in some regions in China with generalized logistic growth model (GLM), Richards Model and Sub-Epidemic Model in a short-term (10 days). They found that the number of cases will increase. They estimated that the case increase (GLM) in the Guangdong and Zhejiang regions would be lower by using the Richards models and the rate would be higher by using the sub-epidemic model [21] .
In a study on the risk of infection of COVID-19 detected in a passenger ship in China in February 2020, it was observed that the risk of infection in those with close contact was higher than those with no close contact. The estimated of number of cases was obtained by back calculation method [22] . Relative Error (RMSRE) and R 2 . They found that the best method for modeling and estimating the number of total cases was the FPASSA-ANFIS method [23] .
Kuniya (2020) estimated the outbreak peak of coronavirus disease in Japan with the SEIR compartmental model [24] . In another study, the number of reproduction of the Wuhan novel coronavirus 2019-nCoV was estimated with the susceptible-exposed-infected-removed (SEIR) compartment model [25] .There are studies on coronavirus disease by different researchers using various statistical methods. Yuan [26] [27] [28] [29] . In our study, we compared the time series analysis using Weibull, Negative Exponential, Von Bertalanffy, Janoscheck, Lundqvist-Korf and Sloboda models, which are different from the methods used in previous studies. According to the literature, there is no publication about nonlinear mathematical models used in our study on coronavirus outbreak before.
The parameter estimates of the Weibull, Janoschek, Sloboda and Lundqvist-Korf were close to each other in the analysis of nonlinear growth models regarding the number of cases and deaths in China. Their R 2 and MKE statistics were similar. The parameter estimation and compliance statistics of the Negative Exponential and von Bertalanffy models, which are among the nonlinear models, differed in both the number of cases and the number of deaths. When nonlinear models are examined, R 2 value, which is used as a criterion in comparing models, was obtained from the lowest Negative Exponential model both in the number of total cases and the number of total deaths. The lowest MSE was obtained from Sloboda model for the number of total cases and from Lundqvist-Korf model for the number of total deaths.
In terms of time series analysis, the number of total cases and the number of deaths are modeled differently. The number of total cases were modelled with ARIMA(0,2,1) that is a moving averages method, while the number of total deaths were modeled by the Holt linear trend model, which is a exponential smoothing method. According to the estimation results, we estimate that the number of total cases and deaths will be increase and this will be a big danger. For this result, all necessary precautions must be taken against the danger.
. CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review)
The copyright holder for this preprint this version posted May 1, 2020. . https://doi.org/10.1101/2020.04. 26.20080465 doi: medRxiv preprint 