Processing math: 100%
スキップしてメイン コンテンツに移動

感染症モデル入門 (5) 再生産数Rtとそれに基づくリスク管理


再生産数R_tを理解する。

1.再生産数とはなにか? SIRモデルから導出する。

離散的なSIR感染症モデルの2番目の式(2)を次のように変形する。

\begin{eqnarray}   {I_{t + 1}} &=& {I_t} + \beta {S_t}{I_t} - \gamma {I_t} \hfill \label{eqtn:1} \\    \Rightarrow \,\,\,\,\,\,{I_{t + 1}} - {I_t} &=& \left( {\beta {S_t} - \gamma } \right){I_t} \hfill \nonumber \\    \Rightarrow \,\,\,\,\,\,\frac{{{I_{t + 1}} - {I_t}}}{{{I_t}}} &=& \left( {\beta {S_t} - \gamma } \right) \end{eqnarray}
 最後の式1の左辺の「感染者数の伸び率」が0を超える(下回る)ことは、感染人数I_tが増加(減少)することを意味する。言い換えれば、右辺が、
\begin{gathered}   \beta {S_{t}} - \gamma  > 0 \,\,\,\,\,\,\Rightarrow \,\,\,\,\,\,{S_{t}} > \frac{\gamma }{\beta } \hfill \\ \end{gathered}
であれば、感染者数は増加する。両辺を\gamma/\betaで割ると, 明日の時点における再生産数を得ることができる。
\begin{equation} \boxed{\,\,{{\bar R}_t} \equiv {S_t}\frac{\beta }{\gamma } > 1\,\,} \end{equation}
 これを時点tにおける実行再生産数と呼ぶ。回復人数を示すR_tと区別するために、Rの上にバーを付けて再生産数\bar{R}_tとしている。この意味について以下で議論する。その前に感染症の「閾値定理」と「基本再生産数」を定義する。

2.基本再生産数\bar R_0の直感的解釈。

上の式で時点ゼロt=0に置き換えたものを感染症の「閾値定理」とよぶ。
\begin{equation} \boxed{\,\,{{\bar R}_0} \equiv {S_0}\frac{\beta }{\gamma } > 1\,\,} \end{equation}
 つまり、感染症が今後増加するかどうかは、この値が1を超えているかどうかの「閾値」を示すので「閾値定理」と呼ばれる。
 また上の式の右辺を全人口Nで置き換えたものを「基本再生算数」と呼ぶ。
\begin{equation} \boxed{\,\,{{\bar R}_0} \equiv N\frac{\beta }{\gamma }\,\,} \end{equation}
 これは感染症の流行が始まる時点ゼロ(t=0)で、それ以降何も対策を講じなければ、1人の人が、回復するか死亡するまでに、何人に感染させるかを示す式となっている。
稲葉寿 (2015) 「基本再生産数$\bar R_0$の数理」の表1より再掲

3.実行再生産数\bar R_tの直感的解釈

3.1 実行再生産数\bar R_t

実行再生産数\bar R_tを再度示すと次のようになる。
\begin{equation}    \bar R_t \equiv S_t\frac{\beta}{\gamma} \end{equation}
 あるいはまた、
\boxed{\,\,{{\bar R}_t} \equiv {S_t}\frac{\beta }{\gamma } = {S_t}\beta \boxed{\frac{1}{\gamma }}\,\,}
とかくことができる。ここで右辺の最後の項\frac{1}{\gamma}は、感染日数の確率分布に指数分布という厳しい仮定を設定すると「平均感染日数」とも解釈できる(この導出を理解することは必要ないが、関心のある人ために、初等の数理統計の知識の範囲での説明を数学付録に示す)。
 これから言えることをまとめてみよう。
1) 閾値あるいは臨界値: この値が時点tで1を超えていると感染者数が今後増加する。つまり、感染症は収束しないことを意味する。
2) 3つの要因: 右辺は1つの変数S_tと、2つのパラメータ(定数)、つまり感染率\betaと回復率\gammaからなりたっている。
3)パラメータの比: 2つのパラメータは比(分数)の形をとっている。つまり感染率\betaを回復率\gammaで割ったものが右辺にあらわれている。これは、1人の人が感染をし、その後回復する可能性の大きさを「比率」で表している。この比が大きくなればなるほど、感染は深刻化していくことを表していることは、直感的にも理解できよう。
 これを、具体的な人数で表現するために、その時の感染可能人数をS_tかけて、感染者がどのくらい増えるかを示す。数値例で考えてみよう。
4) 数値例:t日目の感染可能人数がS_t=100人で、感染率が1パーセント\beta=0.01であり、回復率が10パーセント\gamma=0.1であったとしよう。この2つの比は\tfrac{\beta}{\gamma}=\tfrac{0.01}{0.1}=0.1=10パーセントとなる。これに感染可能人数100人を掛けるとt日めの実行再生産数は10人となる。

5)再生産数の推移:前回示したExcelによるSRIモデルのシミュレーションで、以下にしめした図で、右側のJ列にゼロ日(初期時点)からの推移を数値で、また下向きの↙でその推移が図で示されている。この場合、再生産数は急速に減衰していることがわかる。


4. 再生産数に基づくリスク管理と経済政策

実行再生算数を再度見てみよう。
\begin{equation}    \bar R_t \equiv S_t\frac{\beta}{\gamma}=\mbox{感染可能人数}\frac{\mbox{感染率}}{\mbox{回復率}} \end{equation}
 この数値を1以下にしなければいけない。したがって、対策としては、
 1) 感染率\betaを小さくする。
 2) 回復率\gammaを高める。
 3) t日の「感染可能人数」をS_t小さくする。
の3つの政策があることになる。それぞれの意味を考えてみよう。

4.1 政策1: 感染率\betaを小さくする。

感染率を小さくすることは、1) 身体属性にかかわることと、2) 社会的な距離を保つこと、2つに分類されよう。
 最初の「身体属性に関わること」は感染症にたいする抵抗力をつけよということであった。一般的には健康状態を保つこと、特に既往症を持つ人、高齢の人は特に注意しようということが強調された。第2にの「社会的な距離(social distance)をたもつこと」、日本では三密をさける(Avoid the “3Cs”)ということであった、つまり、密閉・密集・密接を回避するとことが推奨された。これにより新型コロナウィルス感染の特徴である飛沫感染リスクを低くすることができることが強調された。また、マスクをする、石鹸を用いた手洗いをする、消毒剤による手洗いなどこのことに付随して推奨された。マスクをすることは、WHOは当初懐疑的であったようにだ。しかし、最近のドイツの労働経済学者の実証研究 Mitze, Kosfeld, Rode, and Wälde (2020)によれば、マスクの使用は感染者を40パーセントも減らすこうかがあることをたしかめている。マスク恐るべし。

4.2  回復率\gammaを高める。

回復率は、感染した人が健常者として回復し、通常の社会的生活に復帰できるような確率と考えることができる。これには、1) 医療資源(お医者さん、看護師などの人的資源と、病院や保健所、感染症の研究所などの物理的な資源)を充実することと、2)既に感染した人に対する特効薬の開発を上げることができよう。

4.3 「感染可能人数」をS_t小さくする。

「感染可能人数」をS_t小さくすることが、今回の新型コロナ危機で一番問題に成った点である。これを基本再生産数との関係で表現すると、
{\bar R_e} \equiv \left( {1 - \alpha } \right){\bar R_0}= \boxed{\left( {1 - \alpha } \right)N}\left( {\frac{\beta }{\gamma }} \right) < 1.0
となる。ここで\alphaは削減率である。\alphaは結果が1以下になるように決定されるのが普通である。
 実はSIR感染症モデルは、人に対する感染症だけでなく、家畜や野生の動物にたいする感染症の伝播もこうしたモデルで記述できる。動物の場合には、ウィルスが人間に感染すると考えられる場合、\alphaは殺処分率となる。
 オランダでは、毛皮やその油脂を目的とするミンクにおいて、新型コロナウィルスに感染したミンクが、人間の感染にも影響するとのことで、3.5万匹のミンクを殺処分とすることが決定された鳥インフルエンザや豚インフルエンザなどでも、人間には感染しないとわかっていても、その波及が農家に甚大な経済損失を及ぼすと考えられる場合には、殺処分が実行されている
 人間の場合にはそうしたことは勿論できない。その代わりに、厳格なロックダウン、あるいはそれよりも緩やかな外出制限をおこなうことにより、感染可能人数を、一時的ではあるが、削減することが行われている。

5. 接触率「8割」削減の意味

今回のサイン型コロナウィルス危機では、特別事態宣言を出すかどうかの判断において、接触率8割削減が提案され、それに沿うような形の行動が取られた。その場合に説明は実行再生算数と接触率との関係ではなく、よりわかりやすい下の図に示すような説明がおこなわれた。

日経新聞電子版、2020/4/25 日経サイエンス20206月号新型コロナ感染症、接触削減「8割必要」モデルで算出より。この元資料は「新型コロナウイルス感染症対策専門家会議 「新型コロナウイルス感染症対策の状況分析・提言」 2020514
つまり、接触率を大きくすることにより、「緊急事態宣言の直前の全国で1日500人の新規感染者数が100人に未満に抑え込める「日数」に何日が必要であるか」という説明が行われた。接触率を8割削減することにより、図の赤の実線で示されているように、8割削減すれば2週間ていどで目標を達成することができるとされた。結果的には上の図の3案のなかでもっとも「厳格」な?8割削減案が実施された。

数学付録: 1/\gammaは平均感染期間を示す
\ref{eqtn:1}は感染人数を表す差分方程式、ここで、感染人数の減少、つまり、感染⇒回復に関する部分だけをとりだすと、I_t = I_{t - 1} - \gamma {I_{t - 1}}となる。これは、感染者がどのように減っていくかを示している。この差分方程式を初期時点から始まり、順次代入を繰り返していくと、
\begin{eqnarray}  \mbox{1日目} \qquad {I_1} &=& {I_0} - \gamma {I_0} = \left( {1 - \gamma } \right){I_0} \hfill \nonumber \\ \mbox{2日目} \qquad {I_2} &=& {I_1} - \gamma {I_1} = \left( {1 - \gamma } \right){I_1} = \left( {1 - \gamma } \right)\left( {1 - \gamma } \right){I_0} = {\left( {1 - \gamma } \right)^2}{I_0} \hfill \nonumber \\    \vdots  \hfill \nonumber \\ \mbox{t日目} \qquad {I_t} &=& {\left( {1 - \gamma } \right)^t}{I_0} \approx {e^{ - \gamma t}}{I_0} \nonumber \end{eqnarray}
となる。最後の近似式、{\left( {1 - \gamma } \right)^t}{I_0} \approx {e^{ - \gamma t}}{I_0}は、{e^x} \approx \left( {1 + x} \right)\,\,\, \Rightarrow \,{\left( {{e^x}} \right)^t} = \boxed{{e^{xt}} \approx {{\left( {1 + x} \right)}^t}}という関係を利用している。ここでeはネピア数。e^xx=0 の周りでテイラー展開をして0次と1次の項だけを考慮すると近似的に(1+x)となる。両辺をt乗すると最終結果を得る。この近似はもしが十分に小さい値を取る時に成立する。元の式で、 x\equiv \gamma、すなわち、回復率であるので、ほぼ成立すると言ってよいだろう。この式は
\frac{{{I_t}}}{{{I_0}}} = {e^{ - \gamma t}}
と書き直すことができる。この式の左辺は、ゼロ期とt期の感染者数の比率である変化率がどのように減少していくかをしめしているので、右辺は、不確実な回復期日数Tが指数分布に従う時に、回復日数Tt日以上である確率(生存確率=1-分布関数)を表している。つまり、
\Pr \left( {\tilde T > t} \right) = {e^{ - \gamma t}}
と書くことができる。指数分布の期待値は、
E\left[ {\tilde T} \right] = \frac{1}{\gamma }
と書くことができる。この結果から次のようなことが言える。
注意1: 指数分布の形状、従って、その平均や分散は、\gammaだけできまる。
注意2 : 指数分布の平均と分散は等しい(厳しい仮定であり現実には合わない、分散のほうが平均より大きい(過分散)、つまり、「平均回復期間<回復期間の分散」であることのほうが普通。
注意3 : 回復率の推定は、ちょっと難しいかもしれないが、COVID-19にかかった人の平均入院期間を計算できれば、回復率\gammaの推定は可能.

回復率\gammaはファイナンスの世界では、信用リスク分析におけるデフォルト強度として馴染み深い。このγの大きさを変えると、指数分布の形状がどのように変わるかを示している。この計算をおこなうExcelシートは、ダウンロード可能なSRIモデルのシートで異なるタブで計算できる。


指数分布の形状と平均・分散を決定するγ

参考文献

Mitze, T., Kosfeld, R., Rode, J., and Wälde, K. (2020). Face Masks Considerably Reduce COVID-19 Cases in Germany: A Synthetic Control Method Approach., IZA DP No. 13319.

Abaluck, J., Chevalier, J. A., Christakis, N. A., Forman, H. P., Kaplan, E. H., Ko, A., and Vermund, S. H. (2020). The case for universal cloth mask adoption and policies to increase supply of medical masks for health workers. Available at SSRN 3567438.

コメント

このブログの人気の投稿

「8割おじさん」の数理モデル

米国の週刊誌Newwweek誌の日本語版『ニューズウィーク日本版』の最新の 2020年6月9日号 では「日本モデル」と題した特集号で、北海道大学医学部の西浦教授が、「8割おじさん」の数理モデル」という論文を寄稿している。冊子体の本誌でも読めるし、また楽天マガジン・dマガジンなどのオンラインの雑誌読み放題でも全てのページが配信されている。 ニューズウィークとしては異例のアカデミックな論文である。数式こそ無いものの、変数名や学術用語が散見される異例の論文に成っている。同誌で経済や政治関係の記事に親しんでいる人にとっては、ちょっと難しい論文かもしれない。 しかし、非常に示唆に富んだ論文である。特に、「8割削減」の意味を丁重に説明している。最近、8割削減に関して、そこまでする必要がなかったという批判があちこちから出ているが、これをよんでから批判をすべきであろう。 疫学における感染症モデルは勿論、マクロ経済や資産価格決定モデル(CAPMやブラック=ショールズモデルなど)も、所詮は複雑な現実の抽象化に過ぎない。しかし、よいモデルというのは、抽象の革新をつかみ、リスク管理を行う場合、大きな間違いをしない基準となりうる。 また、この号では、”日本のコロナ対策は過剰だったか”という、西浦教授と國井修(世界エイズ・結核・マラリア対策基金)との対談記事も掲載されている。ともにお医者さんでありかつモデリングを研究している二人の言葉は、経済やファイナンスのモデルを用いて仕事をしている、アナリスト、エコノミスト、研究者にとっても学ぶことが大きいとおもう。 國井氏はこう述べている「・・・私達モデラーはリスク評価についてはアンダーアクト(控えめに言う)よりは、オーバーリアクト(大げさに言う)して話をすべきと、肝に命じながらやってきた・・・・」 このことが経済や金融市場のリスク分析に当てはまるかどうか、人それぞれ異なる意見をもっているかもしれない、しかし、1つの教訓だろう。 ただし、西浦先生の論文は、数式を使わないようにしたために、高校数学を理解している人にとっては、帰って難しくなっているきらいがある。むしろ、数式を使って説明したほうが良いかもしれない。ところが、多くの感染症に関する論文や本では、特に日本語で書かれた本や論文では連立微分方程式をつかって説明している。そうす...

感染症モデル入門 (3) Excelで感染症 SIRモデルをシミュレーションする

感染症流行をExcelでシミュレーションする(ファイルがダウンロード出来ます)。 1.準備:ファイルのダウンロードをする。 SIRモデルのわかりやすい(つまり最も簡単な)Excelシートを作成した。Excelファイルを用いたシミュレーションは以下の二通りの方法で実行できる。 1.1 「Google スプレッドシート」を用いてオンラインでSIRモデルのシミュレーションを行う 。  オンラインでシミュレーションをする場合には、Google スプレッドシートをつかいます。「 ここ 」をクリックしてください。「感染症数理モデル_入門と応用_森平_2020-05-29」というファイルがGoogle スプレッドシート上で開くはずです。次のような画面を見ることができます。 携帯などから利用していて、通信環境などが良くないときには、オフラインでの実行をおすすめします。>ファイル?オフラインで実行する、とすれば、自分の携帯にダウンロードしてあるGoogle スプレッドシート上でこのプログラムが実行できるはずです。 1.2 Excelファイルとしてダウンロードして、自分のPC上でシミュレーションをおこなう。  Excelシートとして使い時は、上で示したGoogle スプレッドシートのメニュー上で >ファイル>ダウンロード>MicroSoft Excel (.xlsx)  としてファイルをExcelプログラムしてダウンロードしてください。ダウンロードしたファイルを実行すると、次のような画面が現れるはずです。 2.初期値とパラメータ値を設定する。 2.1 Excelのシートで説明しよう。1行目は変数名とパラメー名前である。2行目の薄緑色のセルに具体的な数値を与える必要がある。ここでは次のような値を初期値とパラメータ値として設定している。 1)  感染可能人数の初期値:S_0=100人 2)  感染者人数の初期値:I_0=1人 3)  回復人数の初期値:R_0=0人 4)  感染率:\beta=0.01,つまり1パーセント 5)  回復率:\beta=0.1,つまり10パーセント (削減率については今回は計算結果に影響の無いようにしている) これらの...