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

![大阪全経過Grapher.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/9ce90ad2-63c3-9853-5656-67fb77f289dc.jpeg)
#大阪のコロナ発生数 全経過を1つの式で表す。
上図は Mac の Grapher で作成した 2020.7.7~2021.2.14 の大阪 COVID-19 新規発生数です。

```math 
\ y=544e^{-0.0667(x-30)-e^{-0.0667(x-30)}}+1087e^{-0.0417(x-148)-e^{-0.0417(x-148)}}+1278e^{-0.1(x-190)-e^{-0.1(x-190)}}
```
という1つの式で全経過を表しています。

エクセルを使う場合は以下のようにします。
A の列に日数（0~222）を入力し B の列に以下の数式を入力します。（コピペ可能）

=544*EXP(-0.0667*(A1-30)-EXP(-0.0667*(A1-30)))+1087*EXP(-0.0417*(A1-148)-EXP(-0.0417*(A1-148)))+1278*EXP(-0.1*(A1-190)-EXP(-0.1*(A1-190)))

|列A|列B|
|:--------:|:-------:|
|0|2.46822778|
|1|3.72121537|
|2|5.44014703|
|~中略~|~中略~|
|220|111.912116|
|221|104.413095|
|222|97.4651426|
列A と列B を選択し「挿入->散布図」から描画できます。
![大阪全経過エクセル.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/7e450ba6-5330-9328-a51d-8528604ad2ba.jpeg)
![大阪経過エクセル.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/b4ee7c56-ea38-4cd2-8e54-59764ba0c48c.jpeg)
##算出法の解説
[ゴンペルツ曲線とは何か？（8）](https://qiita.com/BMJr/items/acad4ec612086d20ac79)
[ゴンペルツ曲線とは何か？（10）](https://qiita.com/BMJr/items/5a3fa10d72d107498604)
にて「開始〜ピークまでの日数」と「ピーク時発生数」という2つのパラメータから流行の全体像を描けることを示しています。
![大阪開始日とピーク.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/b4b4b6b6-1e39-169d-5d80-12f6d7e525f0.jpeg)
2020. 7. 7 開始
30日後の 8. 6 ピーク 200人
2020.10.15 開始
48日後の 12. 2 ピーク 400人
2020.12.24 開始
20日後の 1.13 ピーク 470人（※前の波からの上乗せ分）
これら6個のデータを使います。
![大阪前の波と後ろの波.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/32c72328-2535-016c-7935-d8f0b7291450.jpeg)
2020.12.24 からの波を単独で計算すると上図の黒い波になり、前の波（緑）と足し合わせると青い波になります。
2020.7.7 を *x* = 0 として、それぞれの波を平行移動させます。元データと計算結果がきれいに fit しています。
![大阪全経過Final.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/9366a40f-49ba-679f-9118-0deb522cc240.jpeg)
下り坂の勾配が「開始〜ピークまでの日数」と「ピーク時発生数」によって規定されることがわかります。
###※追記（2021.2.20）
2020.3.25 を day 0 として 2021.2.19 (day 331）までの全経過を表すと下図のようになります。
![OsakaFullCourse.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/8a2c1369-06d1-5287-a7fa-c795cc9ffe6e.jpeg)

```math 
\ y=190e^{-0.1176(x-18)-e^{-0.1176(x-18)}}+544e^{-0.0667(x-134)-e^{-0.0667(x-134)}}+1087e^{-0.0417(x-252)-e^{-0.0417(x-252)}}+1278e^{-0.1(x-294)-e^{-0.1(x-294)}}
```
という1つの式で全経過が表せます。

エクセルを使う場合は以下のようにします。
A の列に日数（0~331）を入力し B の列に以下の数式を入力します。（コピペ可能）

=190*EXP(-0.1176*(A1-18)-EXP(-0.1176*(A1-18)))+544*EXP(-0.0667*(A1-134)-EXP(-0.0667*(A1-134)))+1087*EXP(-0.0417*(A1-252)-EXP(-0.0417*(A1-252)))+1278*EXP(-0.1*(A1-294)-EXP(-0.1*(A1-294)))

すべての波にゴンペルツ曲線が当てはまっていることがわかります。
##まとめ
ゴンペルツ曲線はピークの前後数日の範囲を超えると値がほぼ 0 になって無視できるため上記のような重ね合わせが可能となります。今後、いくつ新たな波が出現しても1つの式で表すことができます。
※大阪のデータを用いたのは、ゴンペルツ曲線がきれいに当てはまるためです。
###※追記 第4波の予測（2021.4.22）
2021.4.22 の時点で大阪の第4波の経過を予測しました。

- 大阪の第4波がゴンペルツ曲線に従う。
- 2021.3.13 が第4波の始点、2021.4.21 がピークである。

と仮定しました。
![大阪第4波.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/c85a2e9a-2852-45af-ddb7-01341d3ebb02.jpeg)
2021.3.13 を day 0 とすると

```math 
\ y=3262e^{-0.0513(x-39)-e^{-0.0513(x-39)}}
```
となり 2021.6月末くらいに 100人程度に落ち着くという予測になります。第4波の総感染者数は 63608人となります。
※上記は**終息までの最短コース**であり、新たな波が重なれば終息が後ろにずれ込みます。
###※追記　予測の検証（2021.6.29）
2021.4.25 に緊急事態宣言が発令された影響か、感染者数の収束が予測に比べて早くなりました。予測式との差分から緊急事態宣言の効果を定量化できます。
[ゴンペルツ曲線とは何か？（14）](https://qiita.com/BMJr/items/4757e087752500913f1b) では「緊急事態宣言の効果を定量的に検証する」と題して検証しました。

[ゴンペルツ曲線とは何か？（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)
[ゴンペルツ曲線とは何か？（12）](https://qiita.com/BMJr/items/0cf2f2d2801ef52cda60)
[ゴンペルツ曲線とは何か？（13）](https://qiita.com/BMJr/items/608222a2e39cb01f2f52)
[ゴンペルツ曲線とは何か？（14）](https://qiita.com/BMJr/items/4757e087752500913f1b)