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

#実測値は Gompertz 曲線に当てはまるか？　その２
##累計発生数の検証
[ゴンペルツ曲線とは何か？（3）](https://qiita.com/BMJr/items/4e2782e8f1d3c1db6801)
[ゴンペルツ曲線とは何か？（6）](https://qiita.com/BMJr/items/60858ef7148732b04160)
で累計発生数がゴンペルツ曲線に当てはまるかを検証する方法を示しました。
累計発生数を対数化した値を $\ F(t)$ とし、任意の時間間隔を *s* とします。
対数化された Gompertz 曲線では、以下の式が成り立ちます。

```math
\ F(t+s)=a+bF(t)
```
$\ F(t+s)$ は切片 *a*、傾き *b* の $\ F(t)$ の一次関数となります。
##新規発生数の検証
本稿では累計発生数ではなく**新規発生数**がゴンペルツ曲線に当てはまるかを大まかに推測する方法を示します。
[Winsor CP. The Gompertz Curve as a Growth Curve. Proc Natl Acad Sci U S A. 1932;18:1-8.](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1076153/pdf/pnas01729-0009.pdf)
にゴンペルツ曲線の変曲点に関する記載があり
「累計発生数が**最大値の 1/*****e*** に達した時が変曲点となる」と述べられています。
（ちなみにロジスティック曲線では最大値の 1/2 に達した時が変曲点になります。）

```math
\ G(t)=81301\times0.0738^{e^{-0.1503(t-8)}}
```
で表せるゴンペルツ曲線があるとします。
以下に累計発生数 *G*(*t*) と1日毎の差分から導いた新規発生数を図示します。
![gompertz_inflection.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/8db03366-ab54-2012-bc49-12a2011b8961.jpeg)
累計発生数が最大値の 1/*e* に達した時に新規発生数がピークに達していることがわかります。
観測開始からピークまでの発生数 (*a*) は総数の 1/*e* = 0.3679
ピークから最大値までの発生数 (*b*) は総数の 1 – 1/*e* = 0.6321
*a* : *b* = 1 : 1.72 ≒ 1 : √3 となります。
「1 : √3」といえば・・・三角定規の辺の比率ですね。
45°、45°、90°の三角定規の辺の長さは 1 : 1 : √2
30°、60°、90°の三角定規の辺の長さは 1 : 2 : √3
![twin_triangle.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/2c7d1390-f15f-4ccb-a165-a18858c94f44.jpeg)
2つの三角定規を上図のように組み合わせると底辺が 1 : √3 に分割できます。
###2枚の三角定規の法則
開始からピークまでの発生数：ピークから終焉までの発生数 ≒ 1 : √3 であることを利用して上記の新規発生数を2枚の三角定規に合わせてみましょう。
![twin_triangles_2.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/0c7dd330-5984-6e4a-be78-8dccdd947e93.jpeg)
・・・新規発生グラフの縦横比を調整するときれいに2枚の三角定規に fit します。
以下は日本の COVID-19 第1波の新規発生数です。
![Japan_1st.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/d9cc0ada-ed23-c36e-3be1-c64febf31311.jpeg)
縦横比を調整して2枚の三角定規に当てはめると
![Japan_1st_triangle.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/3a3a4fa5-f4e7-7ac7-82f5-33f023a2f04a.jpeg)
第2波はどうでしょうか
![Japan_2nd.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/edf1e0f5-1341-8b08-c9a6-9ff93e73e524.jpeg)
縦横比を調整して2枚の三角定規に当てはめると
![Japan_2nd_triangle.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/22367490-fcd3-a338-2b7d-db8a96fae5c2.jpeg)
割とよく fit します。
米国の第2波も縦横比を調整すると・・・
![US_2nd.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/84a9d703-8a01-e9eb-f4e4-4053e8124f4f.jpeg)
イタリアとドイツの第1波も同様に縦横比を調整すると・・・
![Italy_German.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/59c13a69-3bdf-fc47-618b-1898569a4249.jpeg)
よく fit しています。
※上記はゴンペルツ曲線が当てはまりそうなケースを恣意的に選んでおり、「全ての国の第○波が2枚の三角定規に当てはまる」とは限りませんのでご注意ください。
###※追記（2020.12.3）
2020.12.3 時点で日本の第3波に当てはめてみました。
![Japan_3rd_wave.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/193431/7d6ef52d-114f-24d3-19fd-07a9708ce4d8.jpeg)

- 11月の終わりがピークである。
- 第3波がゴンペルツ曲線に従う。
- 人々の行動変容の強度が一定である。

と仮定すると第3波の収束は 2021.1.20 頃となる見込みです。
流行開始からピークまでの発生数は総数の 1/*e* = 0.3679
ピークから収束までの発生数は総数の 1 – 1/*e* = 0.6321
なので、現時点で病床占有率が37%を超える地域は収束までに病床が満杯になる可能性があります。
###※追記（2021.3.8）
Qiita の読者には高校生の方もいらっしゃるので、勉強のために「ゴンペルツ曲線は累計発生数が**最大値の 1/*****e*** に達した時が変曲点となる」ことを以下に証明します。
累計発生数を

```math
\ G(t)=G_{max}e^{-e^{-kt}}
```
とします。*G*<sub>*max*</sub> は最大値、*k* は減衰係数です。
*G*(*t*) を2回微分したもの *d*<sup>2</sup>*G*/*dt*<sup>2</sup> が 0 の時、*G*(*t*) は変曲点となります。
まずは *G*(*t*) を微分してみましょう。

```math
\begin{align}
\ G(t)&=G_{max}e^{-e^{-kt}}\\
\ U=-e^{-kt}&,　V=-kt　とすると\\
\ G(t)&=G_{max}e^U\\
\ U&=-e^V\\
\frac{dG}{dt}&=\frac{dG}{dU}\frac{dU}{dV}\frac{dV}{dt}\\
\frac{dG}{dU}&=G_{max}e^U\\
\frac{dU}{dV}&=-e^V\\
\frac{dV}{dt}&=-k\\
\frac{dG}{dt}&=G_{max}e^U(-e^V)(-k)\\
\ &=kG_{max}e^Ue^V\\
\ &=kG_{max}e^{-e^{-kt}}e^{-kt}\\
\ &=kG(t)e^{-kt}\\
\end{align}\\
```
*dG*/*dt* をもう1回微分します。
(*F*(*t*)*G*(*t*)) ' = *F* '(*t*)*G*(*t*) + *F*(*t*)*G* '(*t*) という公式を使います。

```math
\begin{align}
\frac{dG}{dt}&=kG(t)e^{-kt}\\
\frac{d^2G}{dt^2}&=kG'(t)e^{-kt}+kG(t)(-ke^{-kt})\\
\ &=k^2G(t)e^{-2kt}-k^2G(t)e^{-kt}\\
\ &=k^2G(t)e^{-kt}(e^{-kt}-1)\\
\frac{d^2G}{dt^2}&=0　の時\\
\ e^{-kt}-1&=0\\
\ e^{-kt}&=1\\
\therefore t&=0\\
\ G(0)&=G_{max}e^{-e^0}\\
\ &=G_{max}e^{-1}\\
\ &=\frac{G_{max}}{e}\\
\end{align}\\
```

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