General characteristics of COVID-19 in China
Wuhan, Hubei Province shut down outward traffic beginning January 23, 2020, followed by the rest of Hubei Province. To find high-risk areas caused by imported cases, we drew a heatmap of the migration out of Hubei on January 22, 2020 (Fig. 1a), which indicated that people mainly migrated to Henan, Hunan, Chongqing, Jiangxi, Guangdong, Anhui, Sichuan, Jiangsu, Zhejiang, Beijing and Shanghai. A heatmap of the cumulative confirmed cases in Chinese provinces from January 22 to March 4, 2020 highlights similar provinces (Fig. 1b). Hubei Province was the location of the concentrated COVID-19 outbreak, followed by its neighbors (Henan, Anhui, Jiangxi, Hunan and Chongqing) and some economically developed and densely populated provinces (Guangdong, Zhejiang, Jiangsu, Shandong, Sichuan, Shanghai and Beijing). Thus, Sichuan, Guangdong, Beijing, Shandong, Chongqing, Zhejiang, Jiangxi, Anhui, Jiangsu, Hunan, Shanghai and Henan were selected as high-risk areas with imported cases for further analysis. In addition, since Hubei Province had the most severe epidemic, we also analyzed national data excluding Hubei Province to present the average epidemic in other provinces.
Table 1 shows the peak number of confirmed COVID-19 cases, the corresponding peak date and the cumulative number of confirmed cases in China. Fig. 2 shows a time series of confirmed COVID-19 cases in the identified provinces. The confirmed COVID-19 cases in Hubei Province and nationwide showed a rapid increase before February 4, followed by a decline, and gradually stabilized after February 18, 2020. In high-risk provinces with imported cases, the peak of confirmed cases was around January 30, 2020 in Sichuan, Guangdong, Zhejiang and Shanghai, and around February 2, 2020 in Beijing, Chongqing, Jiangxi, Anhui, Jiangsu, Hunan and Henan.
Table 1 Peak number of confirmed COVID-19 cases, corresponding peak date and cumulative confirmed cases in Hubei Province, China and twelve high-risk provinces from January 22 to March 4, 2020.
Area
|
Peak confirmed cases
|
Peak date
|
Cumulative cases
|
China
|
3893
|
2/4/2020
|
81047
|
Hubei
|
3156
|
2/4/2020
|
67990
|
China except Hubei
|
890
|
2/3/2020
|
13057
|
Sichuan
|
36
|
1/30/2020
|
538
|
Guangdong
|
127
|
1/31/2020
|
1325
|
Beijing
|
32
|
2/2/2020
|
411
|
Shandong
|
45
|
2/5/2020
|
757
|
Chongqing
|
38
|
2/2/2020
|
571
|
Zhejiang
|
132
|
1/29/2020
|
1209
|
Jiangxi
|
85
|
2/3/2020
|
937
|
Anhui
|
72
|
2/3/2020
|
991
|
Jiangsu
|
37
|
2/3/2020
|
631
|
Hunan
|
74
|
2/1/2020
|
1017
|
Shanghai
|
27
|
1/30/2020
|
329
|
Henan
|
109
|
2/3/2020
|
1273
|
Two outliers occurred in China and Hubei Province on February 12 and 13, as the National Health Commission of the PRC revised the definition of COVID-19 confirmed cases in Hubei Province on February 12, adding “clinical case” to “confirmed case,” and left the other provinces unchanged [22]. Another outlier was found in Shandong Province on February 20, corresponding to an outbreak at a prison with 200 confirmed cases [23]. The overall trend of confirmed cases in the other provinces increased first and then decreased.
Impact evaluation of emergency response
We fitted the growth curves at two different periods to assess the impact of the emergency response implemented in each province. Fig. 3 shows the growth curves of each area. The coefficients of the logistic growth curve models in two periods are referred in Supplementary Table S2 and S3. The fitted cumulative confirmed cases were close to the actual observed cases, and the R2 of all models was above 0.95.
The average growth rates of the two periods in Hubei Province, China and twelve high-risk provinces are presented in Table 2 and Fig. 4. The average growth rate decreased by 44.42% nationally and by 32.5% outside Hubei Province. The average growth rate in each province decreased significantly after the emergency response. The average growth rate in the twelve high-risk areas decreased by 29.8%, which was lower than that outside Hubei Province. Before the emergency response, the provinces with the highest average growth rates were ranked from highest to lowest as follows: Hunan, Hubei, Zhejiang, Shandong, Jiangxi, Jiangsu, Guangdong, Sichuan, Anhui, Henan, Chongqing, Beijing and Shanghai. Hubei, Shandong, Zhejiang, Jiangxi and Hunan had growth rates higher than the national average. After the emergency response, the average growth rate of each province from highest to lowest was Zhejiang, Hunan, Anhui, Shanghai, Jiangxi, Jiangsu, Hunan, Guangdong, Hubei, Chongqing, Beijing, Sichuan and Shandong. The growth rates of Guangdong, Zhejiang, Jiangxi, Anhui, Jiangsu, Hunan, Shanghai and Henan were higher than the national average.
Table 2 Comparison of the average growth rates before and after the emergency response in China, Hubei Province and twelve high-risk provinces.
Area
|
r_1a
|
r_2b
|
Percentage decrease
|
China
|
0.565
|
0.314
|
0.444
|
Hubei
|
0.614
|
0.328
|
0.466
|
China except Hubei
|
0.508
|
0.343
|
0.325
|
Sichuan
|
0.475
|
0.279
|
0.413
|
Guangdong
|
0.498
|
0.370
|
0.257
|
Beijing
|
0.443
|
0.308
|
0.305
|
Shandong
|
0.584
|
0.179
|
0.693
|
Chongqing
|
0.450
|
0.313
|
0.304
|
Zhejiang
|
0.603
|
0.435
|
0.279
|
Jiangxi
|
0.576
|
0.397
|
0.311
|
Anhui
|
0.469
|
0.417
|
0.111
|
Jiangsu
|
0.509
|
0.393
|
0.228
|
Hunan
|
0.625
|
0.418
|
0.331
|
Shanghai
|
0.440
|
0.402
|
0.086
|
Henan
|
0.468
|
0.393
|
0.160
|
Average of 12 high-risk areas
|
0.512
|
0.359
|
0.298
|
a: r_1: Average growth rate before the emergency response.
b: r_2: Average growth rate after the emergency response.
Prediction capacity evaluation of logistic growth curve models
We used cumulative confirmed case data, from January 22 to February 4, 2020, to simulate a short-term dynamic prediction. Table 3 shows the MAE and MAPE of the logistic growth curve model in each province. Fig. 5 shows the 1-step dynamic prediction of the logistic growth curve model in Hubei Province, China and twelve high-risk provinces. The 1-step dynamic prediction outperformed the rest, with a MAPE of 1.16%-5.45% in different areas. Except for the models for China, Hubei and Shandong provinces, which were affected by the three outliers mentioned above, the models showed predictions close to the observations.
Table 3 MAE and MAPE of the logistic growth curve model in Hubei Province, China and twelve high-risk provinces.
Area
|
MAE
|
MAPE(%)
|
1 out-of-sample
|
3 out-of-sample
|
7 out-of-sample
|
1 out-of-sample
|
3 out-of-sample
|
7 out-of-sample
|
China
|
1322.3
|
2170.57
|
2285.19
|
3.54
|
5.02
|
13
|
Hubei
|
1392.29
|
2472.08
|
2290.85
|
4.05
|
6.28
|
14.04
|
China except Hubei
|
1780.99
|
1781.32
|
1805.11
|
3
|
3.97
|
6.55
|
Sichuan
|
2.83
|
7.7
|
11.1
|
4.58
|
5.88
|
8.77
|
Guangdong
|
4.81
|
9.62
|
16.37
|
2.64
|
3.3
|
4.99
|
Beijing
|
2.41
|
3.39
|
3.77
|
3.43
|
3.97
|
5.78
|
Shandong
|
17.6
|
23.49
|
27.35
|
5.45
|
7.71
|
12.95
|
Chongqing
|
2.17
|
4.03
|
6.04
|
2.7
|
3.33
|
4.73
|
Zhejiang
|
5.24
|
13.98
|
17.22
|
3.4
|
4.27
|
6.52
|
Jiangxi
|
3.92
|
9.04
|
13.7
|
1.72
|
2.57
|
4.83
|
Anhui
|
3.59
|
9.33
|
15.49
|
1.16
|
1.82
|
4.3
|
Jiangsu
|
2.78
|
9.69
|
13.33
|
2.71
|
3.94
|
7.2
|
Hunan
|
4.15
|
10.34
|
14.66
|
2.03
|
2.86
|
5.07
|
Shanghai
|
1.75
|
3.13
|
4.28
|
4.48
|
5.28
|
7.5
|
Henan
|
4.69
|
8.73
|
14.2
|
1.68
|
2.35
|
4.28
|