Research Article

Spatiotemporal Study of COVID-19 in Fars Province, Iran, October-November 2020: Establishment of Early Warning System

Table 5

The frequency 55 369 forecasted cases during 1 October–14 November 2020 accompanied by the results of time series analysis in the 26 cities of Fars province, Iran.

Time series analysis results
CityForecasted casesModelResidual analysis tests ()Order of min AICC YWAICBIC
n%Ljung-boxDiff sign pointsRank test statisticMcLeod-LiTurning pointsJarque-bera

Shiraz1489526.9MA (8)0.2300.300.690.0300.100<0.001082808303
Lar27624.99SARMA12 (2, 2)0.4000.90.490.0400.560<0.001010981101
Jahrom17163.10SARMA12 (3, 4)0.7700.730.37<0.0010.010<0.001012211225
Fasa17723.20ARMA (1, 1)0.4600.730.12<0.0010.090<0.001013711371
Kazerun16613AR (1)0.4400.300.040.0200.220<0.001012671266
Marvdasht22154SARMA12 (1, 1)0.6100.90.040.0100.560<0.001012471247
Darab14402.60ARMA (1, 1)0.9600.020.800.0400.180<0.001013601361
Abade16613SARMA12 (1, 1)0.9600.420.600.0400.180<0.001013601361
Zarindasht18833.40AR (1)0.8600.030.850.0100.9500.030013571356
Firoozabad14952.70SARMA12 (1, 2)0.0600.080.210.0300.010<0.001013911389
Lamerd14402.60SARMA12 (2, 2)0.1900.210.34<0.0010.00020.030013641368
Mamasani14952.70SARMA12 (1, 2)0.5800.040.320.0400.490<0.001013811380
Ghirkarzin14402.60ARIMA (3, 1)0.9800.730.510.0200.4500.010013641366
Safashahr14402.60ARIMA (1, 1)0.5000.450.760.00100.3500.020013481348
Sepidan19933.60SARMA12 (3, 1)0.1300.300.920.0300.9500.040013441342
Stahban14402.60ARMA (1, 1)0.9800.20.81<0.0010.180<0.001013911392
Eghlid14402.60MA (1)0.1400.730.170.0200.670<0.001013911388
Farashband14402.60SARMA12 (2, 2)0.9200.730.40<0.0010.0400.001013691370
Pasargad16613MA (1)0.0600.040.470.0100.4900.050013341332
Neireez14402.60AR (1)0.2300.450.340.0020.4500.010013641363
Arsenjan14402.60ARMA (2, 1)0.0800.420.140.0300.830<0.001013721367
Rostam14402.60AR (1)0.0900.070.20<0.001<0.001<0.001013211321
Khonj14402.60SARMA12 (2, 1)0.6800.060.05<0.0010.830<0.001012951291
Kharame14402.60SARMA12 (1, 3)0.4500.430.230.0010.0010.010013561357
Bavanat14402.60MA (1)0.2200.910.02<0.0010.310<0.001013081306
Sarvestan14402.60SARMA12 (2, 2)0.5400.90.960.0100.5900.130012601258

AR, autoregressive; MA, moving average; ARMA, autoregressive moving average; ARIMA, autoregressive integrated moving average; SARMA, seasonal autoregressive moving average. It is noteworthy to mention that if four out of the six residual analysis test statistics are statistically significant, it is enough to say that the model is a good fit; based on residual analysis McLeod- Li test statistic needs to be significant and the other five tests should be greater than or equal to 0.05; in addition, the order of min AICCYW, which assesses the mean of white noise residual, needs to be zero. The lower the AIC/BIC scores, the better the model fits.