コロナ感染推移 Excel近似曲線による予測試行

kensyou_jikenbo @kensyou_jikenbo

[メモ　コロナ将来予測] Excelの「近似曲線作成」と「延長予測」機能を試してみた。 対象は、急増で連日マスコミ報道の「報告日別集計」ではなく、原理的に検査数の影響を受けにくい「発症日別集計」を平均潜伏期間5日で換算した「推定感染日別集計」 約1か月の多項式近似で良く合っていると思われる。 pic.twitter.com/6fZMqo43di

拡大

kensyou_jikenbo @kensyou_jikenbo

これを曲線をExcelで30日間延長し添付。 なお「発症日別集計（推定感染日別も）」はデータ確定までに時間がかかるため、前掲近似は昨日時点確定済みの12月24日分までを対象。 グラフの形状だけでも見て頂くと、急カーブの上昇。 12月24日は711人だったが、1月7日予測は約3000人。 pic.twitter.com/KBTpug5WJM

拡大

kensyou_jikenbo @kensyou_jikenbo

但し前述のようにデータ確定に時間がかかるので、報告日別集計での確定は約2週間後とすると、1月7日のデータは1月21日頃になる。 また「発症数＜報告数」だが、単純計算で報告数を予測してみる。 「3000÷711⇒約4.2倍」なので1月7日報告数2447人を4.2倍で約1万人。 つまり1月21日は約1万人予測（東京)

kensyou_jikenbo @kensyou_jikenbo

仮に1月7日時点で潜在的に約1万人にも達していれば、それが報告に現れて来る過程で、宣言の効果が見えやすくなると思われ、判定も、しやすくなると推察。 まず今後2週間の推移が注目。

kensyou_jikenbo @kensyou_jikenbo

[メモ　西浦氏の計算] 同氏予測が各社で報道されているが1月5日NHK記事で計算方法を推測してみる www3.nhk.or.jp/news/special/c… ①1月5日：新規陽性者の報告数1278人（添付図左側） ②1月5日～2月28日：55日間 ③平均世代時間：５日 （東洋経済が西浦氏監修の計算で使用） toyokeizai.net/articles/-/351… （続く） pic.twitter.com/2HUWr6NpYq

拡大

kensyou_jikenbo @kensyou_jikenbo

④　②/③＝11世代（2月末まで） ⑤実行再生産数：1月5日～3月末「1.1で一定」と仮定（添付図） ⑥　⑤④から1.1^11≒2.85（1.1の11乗） ⑦1278人x2.85≒≒3600人（前ツィ図の2月末3500人予測と大体合う） ⇒これが予測の計算根拠とすると、余りにも簡易過ぎないか??? （続く） pic.twitter.com/N2IJbxvWgq

拡大

kensyou_jikenbo @kensyou_jikenbo

1月5日時点の実効再生産数を「3月末まで一定」は、比較用仮定としても現実味が欠け過ぎではないか。 但し西浦氏の主旨としては「対策実施した時の感染者数予測（添付）」の参考提示と警鐘が目的と推察可能。 しかし前記のように元々の予測手法に疑問有り。 数理疫学とはこのようなもの？ pic.twitter.com/Ot5yK7E2FA

拡大

kensyou_jikenbo @kensyou_jikenbo

[メモ　Excel近似曲線の延長による予測] 個人的にも予測してみた ■近似曲線の作成対象期間 ・開始11月20日：第3波の中で更に新たな増加の波が始まった頃 ・終了12月24日：人流などの環境が変わる可能性がある年末年始休みの前（発症日別集計で12月29日） pic.twitter.com/pCPp10lb9w

拡大

kensyou_jikenbo @kensyou_jikenbo

⇒近似曲線を2週間延長した結果添付 (1)1月7日予測：約3500人 一方で前掲西浦氏の「実効再生産数1.1固定」計算 (2)2月末予測：約3600人 (1)は「発症日別集計」からの予測で、毎日話題に上る「報告日別集計」からの(2)より無症状者の分が少ない。 それで時期が2か月近く違うのに(1)(2)が同程度。 pic.twitter.com/nFUetnpYQU

拡大

kensyou_jikenbo @kensyou_jikenbo

以上の比較から見えること ①「発症日別」グラフは、疫学調査の基本とされる「エピカーブ＝流行曲線」。それを延長しての予測も基本的な科学的手法。 つまり基本に忠実と思える予測結果が、西浦氏予測と大きく異なる。どなたか或る程度知識をお持ちの方に、ざっくりでも両方を評価して頂きたいところ。

kensyou_jikenbo @kensyou_jikenbo

②Excel生成の多項式近似グラフは、指数関数ではないが、急激な増加を予測している。 基本に忠実な予測ではあるが、このような急増が実際に起きるだろうか？⇒これは科学的というより個人的感想と推察になるので、今後実際のデータで照合継続予定。 以上

sosuke @nsosuke

東京の感染者数は gotoやめて激増し始め ピークは推定感染日で12/29-30あたり そこから順調に減少傾向で 昨日今日の陽性者の数字がようやく緊急事態宣言当日の影響の数字(感染推定から2週間程度) 何をやっているんだろうか

kensyou_jikenbo @kensyou_jikenbo

@nsosuke [参考　近似曲線による予測の試行結果] 12月24日までのデータで近似曲線を作成し、1月7日まで延長。 途中の「12月25日～30日」までは近似と実データが大体合致と見えます（＝近似は信憑性有りと推定） そして12月30日頃ピークですが、曲線との対比で「急峻なピークアウト」がより明瞭になります。 pic.twitter.com/38UFC0bq4e

拡大

拡大

kensyou_jikenbo @kensyou_jikenbo

@THoriyama @y_mizuno @sekimuyuutaroo @fxitizi @takavet1 @Dolphy_Cooder @tosiakic @carera964 @nsosuke 補足３： 「藤原かずえ」氏の解析も有り（図抜粋添付） ”陽性者増加は国民の気の緩みではなく気象のせいです” agora-web.jp/archives/20498… ⇒これと当方近似曲線（再掲）の比較をしたら、どうなるか。 他の予測方法等も含め、数字に強い人などの解析期待（当方は文系技術者ｗのため数字に弱い） pic.twitter.com/gvgUsx18Ms

拡大

拡大

責務有太郎 @sekimuyuutaroo

@kensyou_jikenbo @THoriyama @y_mizuno @fxitizi @takavet1 @Dolphy_Cooder @tosiakic @carera964 @nsosuke 東京五輪は『ワクチンなしの大会』と考えないといけない jiji.com/jc/article?k=2… 中でも巨額の契約を結ぶ米ＮＢＣテレビはＩＯＣに大 きな影響力を持ち、昨年３月の１年延期決定にも意向 が反映された。 ジュクジュクと横ばいor感染拡大のさなかっぽいかな😭

kensyou_jikenbo @kensyou_jikenbo

@sekimuyuutaroo @THoriyama @y_mizuno @fxitizi @takavet1 @Dolphy_Cooder @tosiakic @carera964 @nsosuke 改めて考えてみると、政府の判断タイミングで一番早そうなのは（再延期する場合に）「緊急事態宣言の再延期決定」時かも知れません。 それまでにも何らかの事情で決定が有るかも知れませんが、基本的には少なくとも2月までは引っ張ると思われ、そうなると早いのは2月7日の前。 注目したいと思います

責務有太郎 @sekimuyuutaroo

@kensyou_jikenbo @THoriyama @y_mizuno @fxitizi @takavet1 @Dolphy_Cooder @tosiakic @carera964 @nsosuke フロリダ州が横槍😱 フロリダ州の財務CFOが五輪“代替開催”立候補する 手紙をIOC会長へ　ツイッターは大炎上 chunichi.co.jp/article/191704 フロリダ州　感染者163万9914人　死者2万5000人超 東京都　　　感染者　9万5534人　死者809人 両方を比較　感染者 17倍　 　　　死者 31倍 アメリカって能天気😭

責務有太郎 @sekimuyuutaroo

@kensyou_jikenbo @THoriyama @y_mizuno @fxitizi @takavet1 @Dolphy_Cooder @tosiakic @carera964 @nsosuke 東京五輪『ワクチンなしの大会』 五輪で「医療スタッフ1万人に依頼」　橋本五輪相が説明 asahi.com/articles/ASP1V… 1万人程度の方に依頼をしてスタッフ確保を図っている 地域医療に支障が生じてしまってはならない コロナ下で、「医療スタッフ1万人」って😱 地域医療に支障がないように出来るの？

責務有太郎 @sekimuyuutaroo

@kensyou_jikenbo @THoriyama @y_mizuno @fxitizi @takavet1 @Dolphy_Cooder @tosiakic @carera964 @nsosuke 「ホワイトリスト」から脱落・・・ ＥＵ、入域許可国リストから日本除外へ －新型コロナ感染の増加受け news.yahoo.co.jp/articles/bb07b… 部外秘の協議だとして匿名で話した同関係者による と、日本で新型コロナウイルスの感染件数が増加し ていることが理由。 五輪大丈夫？😭

kensyou_jikenbo @kensyou_jikenbo

@sekimuyuutaroo @THoriyama @y_mizuno @fxitizi @takavet1 @Dolphy_Cooder @tosiakic @carera964 @nsosuke 専門家代表がこれですから、もはや何でもありでしょう news.yahoo.co.jp/articles/acbc6… ＜尾身会長「緊急事態宣言の効果が恐らく今週末あるいは来週初めに分かると思う。そこを見てこれからの対策でそれが解除の方法、あるいは解除の時期等に重要な影響を及ぼすと思っております」＞ →国会で、まるっきり虚偽

MIZUNO Yoshiyuki　水野義之 @y_mizuno

@kensyou_jikenbo @sekimuyuutaroo @THoriyama @fxitizi @takavet1 @Dolphy_Cooder @tosiakic @carera964 @nsosuke 私も、尾身会長と同じように、見ておりますが。系統的に減らすのは、そんなに簡単には行かないと思うので。

kensyou_jikenbo @kensyou_jikenbo

@y_mizuno @sekimuyuutaroo @THoriyama @fxitizi @takavet1 @Dolphy_Cooder @tosiakic @carera964 @nsosuke 水野先生からレス頂けると非常に有難いです。 というのも記事見せて頂きました rad-it21.com/ai/mizuno20180… そして近代科学の方法論に思いを巡らせながら、今の事態を考えて書いたのが先のツィです。 ピークアウトの「明確な事実」を基にしないのでは、第一の「明晰の規則」から反しているという見方です pic.twitter.com/HdObVJUXWs

拡大

らじらじ @THoriyama

@kensyou_jikenbo @y_mizuno @sekimuyuutaroo @fxitizi @takavet1 @Dolphy_Cooder @tosiakic @carera964 @nsosuke 緊急事態宣言がピークアウト確認後に発令されるのは前回と同じですね。 問題は、今回は高止まり(発症数だけで500付近に留まりそう)する懸念があるので、その際に何を持って『効果』と強弁なさるのか、まずはお手並み拝見と言った所ですね（その予防線で濃厚接触者を絞り始めたとか？）

kensyou_jikenbo @kensyou_jikenbo

@THoriyama @y_mizuno @sekimuyuutaroo @fxitizi @takavet1 @Dolphy_Cooder @tosiakic @carera964 @nsosuke ＜ピークアウト確認後…前回と同じ＞ →「約10日後に発令」という点も同様ですね。 ただ昨年は分かってないことだらけでしたが、今回は約1年経過。 例えば我々でも集計の違いを認識。しかも今回「報告日別」でも明確。 これを隠蔽するのか、そこまで保身と無責任な人たちなのか、どうも読めません。 pic.twitter.com/AUTUSnVknk

拡大

kensyou_jikenbo @kensyou_jikenbo

@THoriyama @y_mizuno @sekimuyuutaroo @fxitizi @takavet1 @Dolphy_Cooder @tosiakic @carera964 @nsosuke 補足： ＜何を持って『効果』と強弁＞ →この点がポイントになるでしょうね。 「宣言の効果」を語る場合に、まず前掲図の「宣言前の大幅減少（起点の取り方によっては約半分）」の説明が必要。 しかしそれを表面化させると「えっ！宣言って何だったのか？」となりますｗ やはり誤魔化すのでしょうかね