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

ドイツのメルケル首相は語る、再生算数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パーセント （削減率については今回は計算結果に影響の無いようにしている） これらの...
コメントを投稿
続きを読む

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

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

感染症モデル入門 (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$としている。この意味について以下で議論する。その前に感染症...
コメントを投稿
続きを読む