実効再生産数Rtは「1人の感染者が平均で何人に感染させるかを表した量」ではない

直観的であるが故に広まってしまった定義の誤り、および西浦博シミュレーションのデタラメを指摘する。
4
MAKIRINTARO @MAKIRIN1230

平均感染期間を平均発症間隔と等しいとしてしまったSIRモデルでの新規陽性者数の推移予想はこんな感じ。再生産数が1.4程度でもかなりかなり過大評価になってしまうことがわかる。ここでは実効再生産数を用いたモデルの結果比較するために、感染する可能性がある人は無数にいるとしている。 pic.twitter.com/YqMB6ODq2B

2021-01-09 22:32:33
拡大
拡大
拡大
拡大
MAKIRINTARO @MAKIRIN1230

西浦さんが平均発症間隔をいくらにしているかは発表されていないが、5日程度にしていることは今回の計算との比較で分かる。ただし再生産数が1.1程度の場合、平均感染期間と平均発症間隔を混同しているかどうかは混同していても差が小さいのでわからない。 twitter.com/MAKIRIN1230/st… pic.twitter.com/03rg52r8qW

2021-01-09 22:42:52
拡大
MAKIRINTARO @MAKIRIN1230

それにしても再生産数が1.1から1.2に変わるだけで、それが2ヶ月続くと、かなり新規陽性者数が増えてしまう。 twitter.com/MAKIRIN1230/st…

2021-01-09 22:47:24
MAKIRINTARO @MAKIRIN1230

またこのような出鱈目なシミュレーションになっていなければいいんだけれど、その可能性は低い。 twitter.com/MAKIRIN1230/st…

2021-01-09 22:51:53
MAKIRINTARO @MAKIRIN1230

これらのシミュレーションでは介入前の倍加時間が2.2日になっているんだけど、彼らは本気で倍加時間が2.2日と考ええていたんだろうか? そうだとしたらヤバすぎる。 pic.twitter.com/W7MGCSLS0O

2020-05-01 21:35:13
uncorrelated @uncorrelated

@MAKIRIN1230 公開されたソースコードは確認されましたか? github.com/contactmodel/C…

2021-01-14 06:11:30
MAKIRINTARO @MAKIRIN1230

実効再生産数が1.3から減少していった場合の報告日別陽性者数の過去7日間の移動平均の推移はこんな感じ。 pic.twitter.com/sYaRBFtdHH

2021-01-10 08:19:39
拡大
拡大
拡大
拡大
MAKIRINTARO @MAKIRIN1230

1ヶ月間で報告日別陽性者数の過去7日間の移動平均を500人程度に減らすためには、実効再生産数を0.5ぐらいまで下げる必要があることがわかる。 twitter.com/MAKIRIN1230/st… pic.twitter.com/Cjw8i1qEZV

2021-01-10 08:24:23
拡大
MAKIRINTARO @MAKIRIN1230

この場合でも、その1ヶ月間で約65000人の新規陽性者が出てしまう。 twitter.com/MAKIRIN1230/st…

2021-01-10 08:32:33
MAKIRINTARO @MAKIRIN1230

それにしても、このシミュレーション結果、0.65 scenario では1/10頃にはピークを迎えているようなので、それでは年末に遡って対策を強化したことになるけど。相変わらずどんだけデタラメなシミュレーションしてるんだろう。 pic.twitter.com/nDfrMYjJCd

2021-01-10 09:17:17
拡大
MAKIRINTARO @MAKIRIN1230

前回の緊急事態宣言では、実効再生産数は0.5ぐらいまでしか下がらなかった。またそれが1.5から0.5まで下がるのに約1ヶ月かかっている。つまり、この1ヶ月で報告日別陽性者数の過去7日間の移動平均を500人以下にすることはまず無理ということである。 twitter.com/MAKIRIN1230/st… pic.twitter.com/zNRQ03scmx

2021-01-10 10:37:37
拡大
MAKIRINTARO @MAKIRIN1230

10日間で陽性者数が1000人から2000人になることが異常なのではなく、分科会の専門家の能力の低さが異常なのである。

2021-01-11 08:37:08
MAKIRINTARO @MAKIRIN1230

実効再生産数の間違った解釈が広まったのも日本の専門かの能力の酷さのあらわれのひとつ。 「発症間隔の間に」という文言を抜いて世間に広めた責任は大きい。

2021-01-11 08:44:15
MAKIRINTARO @MAKIRIN1230

実効再生産数とは「(平均)発症間隔の間に新規感染者数が何倍になるかを表した量」で、「1人の感染者が平均で何人に感染させるかを表した量」ではない。 twitter.com/MAKIRIN1230/st…

2021-01-11 08:58:20
MAKIRINTARO @MAKIRIN1230

(平均)発症間隔ははっきりとはわかっていないのだから、それを用いる実効再生産数ではなく、期間を明確に例えば1週間として「1週間の間に新規感染者数(陽性者数)が何倍なったかを表した量」を用いた方がよい。 twitter.com/MAKIRIN1230/st… pic.twitter.com/CdL8bog5oe

2021-01-11 09:08:41
拡大
MAKIRINTARO @MAKIRIN1230

実効再生産数をRで一定とし、発症間隔の平均をT、その分散を無視すると、時刻tの新規感染者数A(t)は A(t)=A(0)×R^(t/T) と表せ、SIRモデルにおいて基本再生産数をR0で一定とし、平均感染期間をT0とすると、集団の人数が十分多いなら、時刻tの新規感染者数A0(t)は A0(t)=A(0)×e^(R0-1)t/T0 と表せる。

2021-01-11 11:15:54
MAKIRINTARO @MAKIRIN1230

T0=Tと誤って同一視した場合、RとR0がそれぞれ 1 に十分近いならそれらを同一視しても 、A(t)とA0(t)は近似的に等しいとみなせるが、それ以外は全く別ものである。

2021-01-11 11:24:17
MAKIRINTARO @MAKIRIN1230

この違いを認識できていない例として、この西浦さんの見当はずれのシミュレーションがある。 twitter.com/MAKIRIN1230/st…

2021-01-11 11:27:46
MAKIRINTARO @MAKIRIN1230

計算例で比較してみると、こんな感じ。 twitter.com/MAKIRIN1230/st…

2021-01-11 11:34:57
MAKIRINTARO @MAKIRIN1230

実効再生産数についての出鱈目な説明と西浦さんの出鱈目なシミュレーションの関係は「鶏が先か、卵が先か」と同じなんだろう。

2021-01-11 12:08:19
Hiroshi Makita Ph.D. 誰が日本のコロナ禍を悪化させたのか?扶桑社8/18発売中 @BB45_Colorado

@MAKIRIN1230 やっぱりあの人おかしいよね。 さすがにデータが揃ってきた(データの欠陥もはっきりした)今頃あれでは。

2021-01-11 15:23:29
MAKIRINTARO @MAKIRIN1230

@BB45_Colorado 実効再生産数を「一人の感染者が平均で何人感染させるかを表した量」という出鱈目を言っている時点で専門家としてアウトです。さらに発症間隔と感染期間を混同するだけでなく、別の形で定義されている再生産数を混同してシミュレーションしているのも酷いです。

2021-01-11 15:32:14
MAKIRINTARO @MAKIRIN1230

@BB45_Colorado 最近公開されたシミュレーション結果も自粛の影響がすぐ出るものになっていて、びっくりです。 twitter.com/MAKIRIN1230/st…

2021-01-11 15:35:15
石野竜三🌟語り芝居プロジェクト @Ryuzou_Ishino

色々なデータとそこから導かれた数字を見て推察するのが好きなのですが、数字の提議があやふやになっている事があってヤバいと思ったりします😵。学生の頃とは違って勉強をしても頭に入りきらなくなってる気がするなぁ😅。 twitter.com/MAKIRIN1230/st…

2021-01-11 15:40:11
Hiroshi Makita Ph.D. 誰が日本のコロナ禍を悪化させたのか?扶桑社8/18発売中 @BB45_Colorado

これ、的確に144文字以内で西浦氏のペテンを指摘しており、惚れ惚れする。 数値と用語の意味をきっちり理解しているからできること。 正直、僕には無理。 twitter.com/MAKIRIN1230/st…

2021-01-11 18:10:31
MAKIRINTARO @MAKIRIN1230

この東洋経済ONLINEのサイトでの実効再生産数の説明文は出鱈目だけど、簡易的な計算式は適切です。 平均世代時間と平均発症間隔は等しいと考えられるので、実効再生産数を(7日/平均発症間隔)乗したものが、7日間で新規陽性者数が何倍になっているのかを表しています。 pic.twitter.com/QSYUC70Ozb

2021-01-11 18:50:57
拡大