« ハイライト十条店でカラフルジャンボチキンカツ ご飯大盛りを食す。Paypayが使えるようになってた。 | トップページ | 大阪王将で焼きそばと餃子を食す。 »

2020年5月 3日 (日)

コロナウイルス、中国、日本、韓国、アメリカでの感染者数を指数関数とロジスティック関数でフィッティングした。(7/5更新)そりゃ連日東京で100人越えすればそうなるよな、、、という増加。アメリカは本格的に指数関数に近づいている。韓国も全然終息しない(一時終息とかいってなのはなんだったのか)

(緊急事態宣言延長で、以前に書いたものが長くなったので分割しました)

まずは日本のPCR陽性者数は厚生労働省が毎日報告してます。

https://www.mhlw.go.jp/stf/seisakunitsuite/bunya/0000164708_00001.html

中国のはWikiPediaに詳しい。

https://en.wikipedia.org/wiki/2019%E2%80%9320_coronavirus_pandemic_in_mainland_China

さて、感染症の微分方程式として

SIRモデル

SIRモデル(感染症の微分方程式、コロナウイルスなども)の計算をカシオの高精度計算サイトkeisan.casio.jpにUP! 4段4次のルンゲクッタ法を使用。

SEIRモデル

SEIRモデル(コロナウイルス感染の微分方程式その2)をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP!ルンゲクッタ法を使用。

の計算ができる自作式をUPしましたが、データのフィッティングには

もっと簡単なSISモデルから出てくるロジスティック曲線を使います。

https://en.wikipedia.org/wiki/Compartmental_models_in_epidemiology

https://ja.wikipedia.org/wiki/%E3%83%AD%E3%82%B8%E3%82%B9%E3%83%86%E3%82%A3%E3%83%83%E3%82%AF%E6%96%B9%E7%A8%8B%E5%BC%8F

ではまずは中国本土の5/2までのデータ。

うーん、4月くらいからなんか人工的に見えるのは気のせいだろうか。

Corona_china_0502

日本は?

まずリニアスケール。17000人くらいが上限のロジスティック関数でよくフィッティングできる。

Corona_japan_linear_0502

ログスケールにすると、なんか3月末とそれ以降で関数形が大きく変わって見える。

5/7に8万人の患者、という予想には至らないようでそこは緊急事態宣言の効果は出てるんでしょう。

Corona_japan_log_0502

 

---

5/7更新

緊急事態宣言から1か月、8万人よりはずっと少ない。検査数も増えてないとはいえ、これはこれで成果ではないかと。

リニアスケール:

Corona_japan_linear_0507

ログスケール:

Corona_japan_log_0507

そして中国。狙ったように毎日1人とか2人とか増やしている。

逆に怪しいというのがわからんのかな(インフルエンザだってもっと患者いるだろう)。

Corona_china_linear_0507

 

5/17更新

本日から韓国も追加。こう比較するとまだ日本は一定になってない感じ。

Corona_japan_log_0517 Corona_china_log_0517 Corona_korea_log_0517

5/25追記

本日で緊急事態宣言が全国で解除。グラフではだいぶ収束はしているように見える(逆に韓国ちょっとずつ上がってる)。

Corona_japan_log_0525 Corona_korea_log_0525

中国はまあどこまで信用していいのか、、、ほとんど増えてない。

Corona_china_log_0525

 

5/30更新

今日からアメリカも追加。どの国も増えてるが、特にアメリカまったく減ってない感じがする。

これはリニアスケールのほうが分かりやすいので今日はそちらで。

Corona_korea_linear_0530 Corona_china_linear_0530 Corona_usa_linear_0530 Corona_japan_linear_0530

6/7更新

中国以外(まあどこまで正しいデータを出しているかは不明ですが)はまだまだ増えてますね。

USはひどすぎるが、韓国はなんで終息したみたいな雰囲気を醸し出したのか、、、

Corona_china_linear_0607 Corona_us_linear_0607 Corona_japan_linear_0607 Corona_korea_linear_0607

 

6/14更新:

東京で集団感染があったということで日本も増えてる、が韓国がめっちゃ増えてるとか中国も市場で感染とか、アメリカも全く減る気配がないなどまだまだ終息の見込みがない?

Corona_korea_linear_0613 Corona_japan_linear_0614 Corona_china_linear_0613 Corona_us_linear_0613

6/22更新:

どこも終息してないが、特にアメリカと韓国は増加の傾きが増えているようにみえる。第二波か、、、

Corona_korea_linear_0622 Corona_japan_linear_0622 Corona_usa_linear_0622 Corona_china_linear_0622

 

6/28更新:

どこも増加の速度が上がっているが、特にアメリカは指数関数に近づいているように見える、、、やばい。

Corona_korea_linear_0628 Corona_japan_linear_0628 Corona_china_linear_0628 Corona_usa_linear_0628

 

7/5更新:そりゃ連日東京で100人越えすればそうなるよな、、、という増加。

     アメリカは本格的に指数関数に近づいている。

Corona_korea_linear_0705 Corona_japan_linear_0705 Corona_us_linear_0705 Corona_china_linear_0705

« ハイライト十条店でカラフルジャンボチキンカツ ご飯大盛りを食す。Paypayが使えるようになってた。 | トップページ | 大阪王将で焼きそばと餃子を食す。 »

日記・コラム・つぶやき」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

« ハイライト十条店でカラフルジャンボチキンカツ ご飯大盛りを食す。Paypayが使えるようになってた。 | トップページ | 大阪王将で焼きそばと餃子を食す。 »

最近の記事

最近のコメント

2025年2月
            1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28  
フォト
無料ブログはココログ