緊急事態宣言から一ヵ月、効果は?数理モデルを用いて検証する。一ヵ月、効果は?数理モデルを用いて検証する。専門TwitterFacebook0はてブ0Pocket02020.05.07データの出典は奥村教授のgithub及びジャッグジャパン株式会社のcsvファイルより、表示されているデータは5/7時点になっています。慣れてきたものの、終わりの見えない感染拡大に不安を募らせながら生活している人もいらっしゃるかと思われます。そこで今回は、既存のデータをもとに、日本、そして今回は東京のデータも分析し、どれほど非常事態宣言、それに基づく自粛・休業要請などが効果を上げたのか、前回の立てた発症者数を推定する数理モデルと、新たに作った感染者数を推定する数理モデルの2つを用いて考察していきたいと思います。前回の記事をご覧になっていない方は、まず先にそちらをご覧になることを強くお勧めします。推定してみる:正規分布型と発症日分布型データの出典は奥村教授のgithubより、表示されているデータは4/27時点になっています。新型コロナウイルスの感染拡大に伴い、4月7日に全国に緊急事態宣言が出されてから早くも半月が経とうとしています。報道では毎日のように新規感染者数が...scienceboy.jp2020.04.28また最後にExcelでの実装データも掲載しています。DLも可能なのでぜひご覧ください。筆者は統計学者でも感染症専門家でもないただの文系学部1年です。感染症数理モデルに関しては北海道大学の西浦博教授などが取り組んでいらっしゃいます。以下に紹介させていただきます。潜伏期間の確率密度を利用して感染時刻をノンパラメトリック推定する逆計算をしています。次週末を目処に説明会見つきでRの推定コードを出せるかとは思います。前回紹介した正規分布モデルですが、問題点がありました。以下にその理由を述べます。確定日-発症日の分布を考慮していない前回紹介した正規分布型では潜伏期間の2~14日に分布しており、検査を受けて結果が返ってくるまでを2日と仮定して確定日-感染日を4~16日の間で正規分布をとるというモデルを作りました。問題点を解決し、より実際の感染者数像に近いモデルを作るために、今回は発症日分布型で用いた「確定日-発症日」確率密度に潜伏期間の推定確率密度(正規分布を用いる)を掛け合わせ、「確定日-感染日」の推定確率密度を導出し、これをパラメーターとして感染者数の推移の推定をしていきます。「確定日-感染日」の推定確率密度を導出するまずこれが「確定日-発症日」の確率密度です。出典:ジャッグジャパン株式会社のcsvファイル,確定日-発症日を確率密度関数に変換したもの次に「発症日-感染日」(潜伏期間)の推定確率密度です。今だ疫学的なデータが見つからないので、今回も前回の正規分布型と同様に正規分布に従うと仮定し、平均値8,標準偏差3の正規分布の確率密度関数で代用します。これを掛け合わせ、「確定日-感染日」の推定確率密度を得ます。これをパラメータにし、前回と同様の手法で実際の感染者数を推定します。推移・増加率を用いて、4/7に出された緊急事態宣言とそれに伴う行動変容がどれほどの効果を上げたか検証していきます。今回紹介した感染日分布型の値になります。推定される感染者数と発症者数の推移を見てみましょう。統計上、その逆の順に推測がしにくくなります(実測値から離れたトレンドは推測しにくいため)また、感染者数が4/11、発症者数4/20までと途中で終わっているのは、それ以降を推測するデータがまだ存在しないためです。見たところ、感染拡大及び発症拡大のピークは過ぎているように思われます。行動変容の効果を検証するには、感染者の増加率がどれほどのペースで下がっているかが重要になるので、増加率のデータを見てみましょう。縦棒は非常事態宣言が出された4/7を指しています。感染者についてですが、緊急事態宣言はそこまで大きな影響を与えていないように思われます。対して発症者数は顕著に減少しています。前回も述べた通り、行動変容がもっともはやく表れるのは感染者数のほうであり、発症者の顕著な増加率減少をもって行動変容は効果があったと断言することは早計といえます。発症者増加率は実測値の影響が大きく、感染拡大のトレンドを判断するにはやはり感染者増加率のほうが適していると考えられます。感染者増加率に関しても1以下を継続して下回り、減少傾向にあるのは好ましい結果だといえるでしょう。感染者や発症者の推移を近似すればもう十数日後には全国の新規陽性確定者も二桁で安定し始めると考えられますが、どうなるのでしょうか。東京も同様に感染拡大のピークは過ぎ去ったように感じられます。全国の推移と異なり、やや増加率の減少は遅い傾向にあるようです。増加率についても見てみましょう。増加率は発症者、感染者ともに全国よりもわずかに高い数値で推移しており、感染者推移が1以下で推移している点は好ましいものの、発症者推移では1に迫る日もあるなど、予断を許さない状態であるといえるでしょう。感染者推移では非常事態宣言の効果が全国よりも薄いと感じられるのも、東京の特徴です。個人的には東京のみが自粛・休業要請を延長して行うというシナリオも十分考えられると思います。まとめ今回は二つの数理モデルを用いて感染者数と発症者数が実際どのように推移しているのかを推測してみました。厚労省が5/6にPCR検査の基準を変更したことなどもあり、これからのデータに関しては既存のデータとリニアにつなぐことが難しくなると考えられますが、継続して検証を続けていきたいと思います。緊急事態宣言とそれに伴う行動変容が機能したかを検証することは、それをやめるうえで重要です。吉村知事が大阪モデルと称して自粛解除のための出口戦略を提示して話題になりましたが、出口戦略を考えるうえでは統計的な解析や疫学的な見地が不可欠であり、さもなければ予期せぬ事態の悪化を招くのみです。これらのモデルを製作することで、皆様が日々の報道に一喜一憂せず、根拠を持って冷静に事態を静観することができるだけでなく、これから提示される「出口戦略」が果たして本当に正しい政策なのかを見極めることもできるようになるのではないか、と思います。最後に数理モデルの実装データ(Excel)を置いておきます。参考までにご覧ください。気ままに生きる文系大学生。写真と絵と自作PCが好き。前回の定例会議で新規メンバーGoneの追加の話がありまして、本日が初記事になる予定でしたが、諸事情で代打の私が投稿する形になります。本記事はタイトル記載作品のネタバレを含みます。ADVの読みやすさ前置きしておくが僕はこの分野に精通している訳ではない。踏み切れなかった「気力」を持ち合わせていた程度だと考えて欲しい。ADVには文字数という制約がある。小説のように1面で何百文字を埋めることはできない。【Web漫画】ラフターラフターep.1ブログ毎日投稿のネタが尽きたわけではない。最近流行りのTwitter漫画。自分も描いてみたんですよね。バズらないかな〜とか邪な感情を抱きながら。タイトルは『ラフターラフター』。コメントをキャンセルメールアドレスが公開されることはありません。※サイト次回のコメントで使用するためブラウザーに自分の名前、メールアドレス、サイトを保存する。気ままに生きる文系大学生。写真と絵と自作PCが好き。defined.")}))}();タイトルとURLをコピーしました。