ゴンペルツ曲線とは何か？（14）

#緊急事態宣言の効果を定量的に検証する
緊急事態宣言が出されなかった 2020.7月〜9月の第2波と緊急事態宣言が出された 2021.3月〜6月の第4波を比較して定量的に効果を検証します。
※ノイズが少ない大阪のデータを用います。
##緊急事態宣言が出されなかった時の経過
2020.7.8~2020.9.7 の新規感染者数と同時期の梅田駅前の人出を示します。
![大阪第2波検証210620.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/6fa1ff0c-bfcb-bbac-d192-c6c052e70e43.jpeg)

```math 
\ y=513.6e^{-0.06346(x-28.25)-e^{-0.06346(x-28.25)}}
```
という二重指数関数に当てはまります。
上り坂と下り坂で人流に大きな変化がないことがわかります。

累計感染者数の実測値と、上式を積分した式は、ほぼ完全に一致します。
![大阪第2波累計検証210620.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/dbfbcae0-ad00-fc15-563b-a096ae86906e.jpeg)

```math 
\ Y=8093e^{-e^{-0.06346(x-28.25)}}
```
上図を対数化したものが以下の図です。
![大阪検証第2波対数210620.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/ba8431a4-0c10-6db5-60e8-9b022846e0af.jpeg)
感染者数が指数関数的に増加する場合、累計感染者数を対数化したグラフは上向きの直線になるはずです。
しかし 2020.7.8 の対数グラフの傾きを *A*<sub>0</sub> とすると日々 *e*<sup>-0.06346</sup> 倍に減衰しています（一定減衰仮定）。これはまさに**二重指数関数（ゴンペルツ曲線）**の特徴です。
減衰係数（0.06346）と 開始〜ピークまでの日数 には以下のような関係が成り立ちます。
[※ゴンペルツ曲線とは何か？（10）](https://qiita.com/BMJr/items/5a3fa10d72d107498604)

```math 
\ 減衰係数=\frac{2}{開始〜ピークまでの日数}
```

人流に変化がなくても収束するという現象は緊急事態宣言が出されなかった日本各地の第2波で見られました。
以下は 2020.7月〜9月の各地の新規感染者数（上段）と駅前の人出（下段）です。
![兵庫福岡愛知第2波210620.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/a8aa1d9e-c81f-a32e-3eae-6498588920d2.jpeg)
第二波の経過中、全く人流が変化しなかった兵庫県の累計感染者数と対数化データを以下に示します。
![兵庫第二波累計対数.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/037deabc-ad5c-a5bf-bc68-22707d5f45bd.jpeg)
大阪と同様、累計感染者数は上記の二重指数関数と一致します。対数グラフの傾きも大阪と同様に一定減衰仮定を満たします。
**実効再生産数 = 接触率 × 基本再生産数** という仮説が正しいとすると対数グラフの傾きが一定減衰するためには人流が以下の右図のように指数関数的に減少しなければいけません。
![兵庫人出実際と仮説.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/8247fa9b-f744-3b1d-de31-5780ee02566b.jpeg)
・・・しかし実際の人流（上図：左）には全く変化がありませんでした。
**実効再生産数 = 接触率 × 基本再生産数** という概念は指数関数を感染者数に当てはめるために考案されたものですが検証が必要であると思われます。

行動抑制の強度が一定でも上昇->下降の変化を来すのは、一定の重力加速度で上昇->下降する砲弾の軌跡に似ています。※砲弾の軌跡は以下のような二次関数となります。
![砲弾の軌跡210620.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/8e9d7d4d-693b-6bf4-9202-66cc401c0403.jpeg)
行動抑制の強度は重力加速度に似ています。
##緊急事態宣言が出された時の経過
上記をふまえて大阪の第4波（2021.3.13~6.13）を検証します。2021.4.25 に緊急事態宣言が出される前をオレンジ色、出された後を赤色にしています。
![大阪第4波新規検証210620.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/d3225529-3ea9-a1fa-69d1-eb8d861252a3.jpeg)
宣言前（2021.3.13~4.24）の新規感染者数から

```math 
\ y=3122e^{-0.05135(x-38.96)-e^{-0.05135(x-38.96)}}
```
という二重指数関数を導きました。
宣言の2週間後くらいから減少の角度が予測式に比べて急峻になっています。
累計感染者数は以下のようになります。
（※流行初期のノイズを減らすために 2021.3.23 からカウントしています）
![大阪第4波累計検証210620.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/62eceacc-96c1-fb21-3ac7-dbf619bdb82f.jpeg)

```math 
\ Y=60791e^{-e^{-0.05135(x-38.96)}}
```
から予測される 2021.3.23~6.13 の累計感染者数は 56930人ですが実測値は 53125人で 3805人少なくなっています。
累計感染者数を対数化すると以下のようになります。
![大阪第4累計対数検証210620.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/bc11873e-1c8b-2caa-e867-0b457c194221.jpeg)
2021.3.23 の対数グラフの傾きを *A*<sub>0</sub> とすると日々 *e*<sup>-0.05135</sup> 倍に減衰しています。
注目すべきは、この減衰が流行初期から一定の割合を保っていることです。
上図を見ると全く緊急事態宣言の効果がないように見えますが宣言発令2週間後の 2021.5.9 からの対数グラフを拡大すると以下のようになります。
![大阪第4対数検証拡大210620.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/d5a7f1d6-392e-9d61-39ba-91654da07dad.jpeg)
2021.5.9 以降で対数グラフの傾きが日々 *e*<sup>-0.07183</sup> 倍に減衰しています。
###もし 2021.3.23 に発令されたらどうなったか？
今回は発令がピークアウトの後だったので、宣言による累計感染者数の減少効果が 3805人と限定的なものになりました。
もし流行初期の 2021.3.23 に緊急事態宣言が発令されたらどうなったでしょうか？
発令2週間後の 2021.4.6 を境に減衰係数が 0.05135 から 0.07183 に変化した場合、
![大阪早期発令の効果210620.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/b8dccce8-7c36-7b50-94df-49a8e188b2bc.jpeg)
ピークは 2021.4.11 になり 2021.6.13 の累計感染者数は 32503人になります。2021.4.25 に発令する場合より約2万人（40%）減となります。
##まとめ
- 行動抑制の強度は、砲弾の軌跡における重力加速度に相当する。
- 人流が一定の場合、累計感染者数を対数表示した時の傾きは一定減衰仮定を満たす。
- 緊急事態宣言は減衰係数を増大させて収束を早める効果がある。
- ピークを過ぎた後に発令した場合、効果は限定的である。
- 流行初期に発令すると大きな効果が得られるため、発令は早ければ早いほど良い。しかし発令からピークアウトまでにタイムラグがあるので見栄えは悪い。（今回の第4波のようにピークアウトに合わせて発令した方が印象として「効いた感じ」が得られやすい）

・・・本稿が新型コロナ対策の一助になれば幸甚です。

[ゴンペルツ曲線とは何か？（1）](https://qiita.com/BMJr/items/24fca066926d471c651e)
[ゴンペルツ曲線とは何か？（2）](https://qiita.com/BMJr/items/1f9135edfb5e68b4f7da)
[ゴンペルツ曲線とは何か？（3）](https://qiita.com/BMJr/items/4e2782e8f1d3c1db6801)
[ゴンペルツ曲線とは何か？（4）](https://qiita.com/BMJr/items/f1665fcf61fc1e182311)
[ゴンペルツ曲線とは何か？（5）](https://qiita.com/BMJr/items/6e48af3017941589bd74)
[ゴンペルツ曲線とは何か？（6）](https://qiita.com/BMJr/items/60858ef7148732b04160)
[ゴンペルツ曲線とは何か？（7）](https://qiita.com/BMJr/items/854d7ea041ac426c3e79)
[ゴンペルツ曲線とは何か？（8）](https://qiita.com/BMJr/items/acad4ec612086d20ac79)
[ゴンペルツ曲線とは何か？（9）](https://qiita.com/BMJr/items/775ae42bb5e0cc6dc461)
[ゴンペルツ曲線とは何か？（10）](https://qiita.com/BMJr/items/5a3fa10d72d107498604)
[ゴンペルツ曲線とは何か？（11）](https://qiita.com/BMJr/items/2055c27c0cb5596816df)
[ゴンペルツ曲線とは何か？（12）](https://qiita.com/BMJr/items/0cf2f2d2801ef52cda60)
[ゴンペルツ曲線とは何か？（13）](https://qiita.com/BMJr/items/608222a2e39cb01f2f52)