スキップしてメイン コンテンツに移動

ドイツのメルケル首相は語る、再生算数R_tを！

2020年4月15日 メルケル首相は語る:基本再生産数R0について


　2020年4月15日ドイツ政府は長く続いたロックダウンの段階的解除を決定した。しかし、メルケル首相は、気を緩めるとどうなるかを「ロベルト・コッホ研究所；Robert Koch Institut」の感染症モデルに基づいてドイツ国民に次のように説明をした。

「ロベルト・コッホ研究所：Robert Koch Institut（英文）」は、結核のワクチンを開発したことで有名なロベルト・コッホによって1891年（明治24年）に王立プロイセン感染症研究所として設立された。ちなみに、慶應義塾大学医学部の創立者でもある北里柴三郎は、1885年（明治18年）ドイツのベルリン大学に当初は2年の予定で、最終的には6年間留学した。そこで、コッホ博士に師事し、破傷風菌の発見をはじめとする様々な業績をあげた。

https://twitter.com/i/status/1250791059371646977

写真はその時のドイツ語での国民への呼びかけである。クリックすれば、その時の映像を見ることができる。

映像の下に字幕が流れ、英語で翻訳をしたものがながれている。基本再生産数に関する部分の概要を日本語で下に示そう。

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
（感染者数の）カーブが平坦になりました。私たちの健康システムに過大な負担をかけないように、カーブをフラットに維持する必要があります。


モデルからの観測値ですが、今の再生産数（R0）は1です。これは平均値ですが、１人が他の１人に感染させる数値です。これが

1人から1.1人になれば、10月までに医療崩壊が生じます。

1.2人になれば、つまりいまより20%多くの人が感染するようになると7月には医療システムの崩壊に、

1.3人になると6月に崩壊が生じます。
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

ここで注目すべきことは、赤字で示したように、メルケル首相は、再生産数R0に言及をしその意味をわかりやすく説明している。この演説はドイツ国民の多大の共感をよんだようである。

「ドイツ日本研究所」が5月27日に行った、ドイツと日本のCOVID-19への政策対応の比較検討をおこなったオンラインフォーラム「National Approaches to Systemic Risk Germany and Japan under the COVID-19 Crisis」でドイツ側を代表して報告をおこなったリスクの研究者であるOrtwin Renn (scientific director at the Institute for Advanced Sustainability Studies (IASS) in Potsdam)氏は、メルケル首相は、コッホ研究所と新型コロナ感染症リスクに関してコッホ研究所と高度な連携のもとで政策決定を行ったと述べていた(発表時に用いられたパワーポイント資料はここから入手可能である）。


また、メルケル首相は、東ドイツ崩壊まで、東ドイツの科学アカデミーで理論物理学の研究を行っていた。感染症のモデルは、このブログの「感染症モデル入門（1)」でも説明するように、非線形の連立微分方程式を用いる。その最も有名なモデルの1つである「SRIモデル」は、科学における触媒の役割を記述モデルがその基礎になっている。メルケル首相は、1986年に博士論文を提出して博士号 (Dr. rer. nat.) を取得。そのご物理学者として分析化学にの研究をおこなったとのことである。あくまでも推測であるが、彼女にとっては、感染症モデルの理解は、アドバイザーの説明を聞いたのであろうが、容易であったに違いない。再生産数R_t の説明も、彼女の政治家としての説得力と彼女の物理学の素養をもった分析科学者としての知識の決勝であるかもしれない。








ラベル:

コメント

コメントを投稿

ここに感想をお願いいたします。

このブログの人気の投稿

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

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

感染症モデル入門： (２) SIRモデル（離散型）の場合。易しいです

感染症モデル入門 (2) SIRモデル（離散型）の場合。易しい！ １．SIR感染症モデル：一日ごとの人数変化  時間間隔が１日という離散的な時間変化を考えると、前回示した微分方程式での記述は、次のような差分方程式に変換できる。 \begin{eqnarray}   \mbox{感染可能人数}  \qquad  {S_{t + 1}} &=& {S_t} - \beta {S_t}{I_t} \hfill \\   \mbox{感染者数} \qquad  {I_{t + 1}} &=&  {I_t} + \beta {S_t}{I_t} - \gamma {I_t} \hfill \\   \mbox{回復人数} \qquad {R_{t + 1}} &=&  {R_t} + \gamma {I_t} \end{eqnarray} これは連続的な時間の変化$dt$を離散的な時間変化$\Delta t$でおきかえたものである。$dt$が「1秒間隔」であったのに対し、$\Delta t$は1日単位と考えてみよう。実際のCOVID-19のデータは1日単位で発表されるので、そのように考えたほうが実際に合っている。 　時間経過を1日と考えたのであるから、感染可能人数を$dS_t \approx \Delta S_t$、感染者数を$dI_t \approx \Delta I_t$, 回復者数を$dR_t \approx \Delta R_t$と表現する。これらも1日あたりの人数の変化である。 　そうすると、例えば、感染可能人数を示す微分方程式は$\left( {\frac{{\Delta {S_t}}}{{\Delta t}} = \frac{{\Delta {S_t}}}{1} = \Delta {S_t} = {S_{t + 1}} - {S_t}} \right) =  - \beta {S_t}{I_t}$と変換できる。つまり${S_{t + 1}} - {S_t} =  - \beta {S_{t}}{I_{t}}$を得る。...
コメントを投稿
続きを読む

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

再生産数$R_t$を理解する。 １．再生産数とはなにか？　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} 　最後の式１の左辺の「感染者数の伸び率」が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$としている。この意味について以下で議論する。その前に感染症...
コメントを投稿
続きを読む