頑張るネムネコ、あのシミュレーションに挑む!! [お前らに質問]
頑張ったネムネコ
昨日、紹介したシミュレーションの結果はあまりにひどいので、ネムネコは、データから、ロトカ・ヴォルテラの方程式
の係数a、b、c、dの求める手法を新たに開発しーー新たに開発なんて大袈裟な。それに、ネムネコが知らないだけで、既に誰かが同じ方法を提案しているに違いない!!――、その手法で得られた係数をもとにシミュレーションしてみたにゃ。
端からお前らに期待なんかしていなかったけれど、
そもそも知らないから、お前らに「どうやったら求められるのか」と質問したんだケロ。
そして、いつものように、誰の手も借りることなく、みずから解決!!
試行錯誤的にいろいろな数値を試し、「あ〜でもない、こうでもない」と、迷いたくないにゃ。
結果は、こちら↓。
このシミュレーションのほうが、紹介したシミュレーションの結果よりも、ずーと、ずっと、データをよく再現していると思うにゃ。
ネムネコの圧勝だケロ。
ネムネコが考えるに、ここに示した結果は、間違いなく、最適な解の一つ!だケロよ。
――最適とは「最も適す」の意味。なのに、最適な解の一つとは、これ如何に!!――
それはそれとして、
これで課題のサメの問題を解くことができる!!
そして、例によって、自画自賛のこの動画を埋め込むのであった。
ベイズの定理と疾病検査の信頼性についてのよもやま話 [高校の統計]
ベイズの定理と疾病検査の信頼性についてのよもやま話
日本でも新型コロナウイルスが流行しているので、「疾病検査の信頼性」に関する問題を一つ。
問題
ある感染症に感染しているかどうか調べる検査法がある。この検査法では、感染者に対して95%が陽性の結果を受け、また、、非感染者の5%に陽性の結果が出る。国民の人口の0.5%が感染症に感染しているとするとき、この検査によって陽性判定を受けたヒトが真の感染者である確率を求めよ。
ベイズの定理を用いて解くこともできるけれど、条件付き確率などの説明をしないといけないので、ベイズの定理を用いない、泥臭い方法でこの問題を解くことにする。
【解答例】
1000人検査するとすると、1000人のうちの0.5%、つまり、
1000×0.005=5人
がこの感染症に罹っていて、残りの
1000−5=995人
が非感染者の人数。
感染者5人のうち、この検査で陽性判定の受ける人数は
5×0.95=4.75人
非感染者995人のうち、この検査で陽性判定を受けるヒトの数は
995×0.05=49.75人
したがって、この検査で陽性の結果になるヒトの数は
4.75+49.75=54.5人
陽性判定を受けた54.5人のうちで実際に感染症に罹っているヒトの数は4.75人だから、答えは
4.75÷54.5≒0.087=8.7%
(解答終)
【ベイズの定理を用いた解答】
感染者である事象をA、非感染者である事象をB、検査で陽性が出る事象をEとする。
問題の条件より、
よって、ベイズの定理より求めるべき確率は
(解答終)
この結果は衝撃的だろう。
なんたって、
検査で陽性判定を受けたヒトの中で本当に感染症に感染しているヒトの割合はたったの9%で、残りの91%は感染していないのに陽性判定を受けた偽陽性のヒトなんだから。
このパラドクス的な結果が起きた理由は、全人口に占める感染者の割合、確率が0,5%と低いからなんですがね。
テレビのワイドショーのコメンテータが「国民の不安を払拭するために、希望するヒトは誰でも新型コロナウイルスの検査を受けれるようにするべきだ」と声高に主張しているけれど、そんなことをしたら、新型コロナウイルスに感染していないのに陽性判定を受ける偽陽性のヒトが山のように出て、結果、医療崩壊してしまうにゃ。
新型コロナウイルスのPCR検査の感度(新型コロナウイルスに感染しているヒトが検査を受けて陽性の結果が出る確率)は70%程度って話があるようだにゃ。
チョット小耳に挟んだ話なのでこの話の真偽はわからないけれど、
この話が本当だとしたら、
希望するヒトは誰でもPCR検査を受けれるようにしたら、もっととんでもない結果が出るケロよ。
仮に、
この検査の特異度(感染していないヒトが検査で陰性になる確率)も同じように70%――非感染者の30%が検査によって陽性判定を受ける!!――とし、日本人の0.1%(日本人の12万人以上)が新型コロナウイルスに感染しているとすると、この検査で陽性判定を受けたヒトの中に占める感染者の割合は
感度、特異度がともに90%だとしても、
日本人の0.01%である約1万2千人のヒトが新型コロナウイルスに感染しているとし、感度、特異度がともに99%であったとしても、
となり、検査によって陽性判定を受けたヒトの中でウイルスに感染しているヒトの割合は1%に満たず、陽性判定を受ける約99%が新型コロナウイルスに感染していない偽陽性のヒトってことになる。そして、この条件だと、日本人すべてを検査したら、感染していないのに陽性判定を受ける偽陽性のヒトが間違いなく100万人を突破する!!のであった。
つ・ま・り、
この手の検査は、たとえその検査がどんなに優れたものであったとしても、感染が(強く)疑われるヒトを対象に行わないと、検査自体が無意味なものになってしまうってことを示している。
参考までに、感度、特異度が70%(0.7)のとき、検査対象となる集団の感染率を横軸に、この検査で陽性判定を受けたヒトの中に占める感染者の割合を縦軸にとり、この関係を示すにゃ。
この図を見ると、検査対象となる集団の感染率が10%(0.1)のとき、検査によって陽性判定を受けるヒトの中に占める感染者の割合は約20%(0.2)、感染率が20%(0.2)のとき37%(0.37)程度であることがわかる。
感度、特異度が90%(0.9)のときは次のようになる。
調査対象の10%(0.1)が感染しているときに、ようやく、陽性判定を受けたヒトの中に占める感染者の割合が50%(0.5)に到達する。そして、クルーズ船ダイアモンド・プリンセス号のように感染率が20%近くでも、陽性判定者に占める感染者の割合は70%程度。裏を返せば、陽性判定を受けたヒトの中で感染していない偽陽性の人の割合が30%(0.3)もあったてこと。
ひょっとしたら、ネムネコが小耳に挟んだ「新型コロナウイルス検査の精度は70%」という話はこのことを言っているのかもしれない。
さらに、感度、特異度が99%(0.99)の場合も示すにゃ。
テレビなどのマスメディア、特にワイドショー受けしないからかもしれないけれど、感染率と偽陽性の人の関係について言及する(感染症の)専門家って、ほとんど、テレビのワイドショーなんかに登場しないよね。
この手の話は、医師国家試験の問題として出題されるらしいから、感染症の専門家ならば誰しも知っている基礎的な知識のはずなのに、マスコミに忖度し、この手の話をしないのだろうか。
YouTubeで、この問題について触れているお医者さんの動画があったので、紹介するにゃ。
高校の数学だと、
条件付き確率は
で定義し、したがって、
そして、
であるとき、事象AとEは独立と習ったかもしれない。
という見慣れない記号が出てきているけれど、これは、高校数学で習った条件付き確率のことでにゃ。見慣れない記号に幻惑されてはいけない。
韓国発入国制限140ヶ国…康京和外相「過度の制限避けるべき」 YouTubeより [ひとこと言わねば]
大阪市、前倒しで4月から給食無償化 産経 [ひとこと言わねば]
大阪市、前倒しで4月から給食無償化https://t.co/nF73cQq78V
— 産経ニュース (@Sankei_news) March 17, 2020
新型コロナウイルスの影響で経済的な負担が大きい子育て世代を支援。松井一郎市長「所得が減る人が非常に増えてくると思う。ひとり親家庭は収入激減もある。皆さんに安心して子育てしてもらいたい」
韓国、ドル資金不足の対応策を18日発表へ=政府関係者 ロイター [ひとこと言わねば]
韓国、ドル資金不足の対応策を18日発表へ=政府関係者 https://t.co/XKT7ELa33L
— ロイター (@ReutersJapan) March 17, 2020
マレーシアのイスラム教集会で新型コロナ大規模感染 産経 [ひとこと言わねば]
マレーシアのイスラム教集会で新型コロナ大規模感染https://t.co/S2z2bVZT04
— 産経ニュース (@Sankei_news) March 17, 2020
イスラム教の大規模集会の参加者から感染が相次いで確認されている。
16日までに400人弱が感染。既に参加者は東南アジアに拡散しており、大規模な宗教行事がクラスターとなっている実態が浮かぶ。#新型コロナウイルス
中国本土の新型コロナ新規感染者が増加、入国者の感染が原因 ロイター [ひとこと言わねば]
中国本土の新型コロナ新規感染者が増加、入国者の感染が原因 https://t.co/vNYDaEvkYl
— ロイター (@ReutersJapan) March 17, 2020
というわけで、本当に、中国が新型コロナウイルス対策で成功しているのか、実に疑わしいにゃ。
疑い深いネムネコは、地方の農村部には、新型コロナウイルス感染者が少なからぬ人数いると踏んでいる。この少なからぬヒトたちが千人単位なのか、万人単位なのかはわからないけれど…。
政府 米国便で搭乗前に体温測定へ 38度以上は渡航認めない方向 NHK [ひとこと言わねば]
政府 米国便で搭乗前に体温測定へ 38度以上は渡航認めない方向 #nhk_news https://t.co/qZmvAFA1RZ
— NHKニュース (@nhk_news) March 17, 2020
首相や閣僚 新型ウイルス経済対策検討で有識者から意見聴取へ NHK [ひとこと言わねば]
首相や閣僚 新型ウイルス経済対策検討で有識者から意見聴取へ #nhk_news https://t.co/3DCzl8WIAh
— NHKニュース (@nhk_news) March 17, 2020
専門家会議「無症状の帰国者も2週間不要不急の外出控えて」 NHK [ひとこと言わねば]
専門家会議「無症状の帰国者も2週間不要不急の外出控えて」 #nhk_news https://t.co/hZ8ZuZ2Hp9
— NHKニュース (@nhk_news) March 17, 2020
日銀 3兆4000億円のドル資金を供給 各国中央銀行と協調 NHK [ひとこと言わねば]
日銀 3兆4000億円のドル資金を供給 各国中央銀行と協調 #nhk_news https://t.co/EZx77j3J3q
— NHKニュース (@nhk_news) March 17, 2020
山陽新幹線の利用者 去年同期比54%減 新型ウイルスの影響 NHK [ひとこと言わねば]
山陽新幹線の利用者 去年同期比54%減 新型ウイルスの影響 #nhk_news https://t.co/VytsuMTCQ8
— NHKニュース (@nhk_news) March 17, 2020
新型コロナウイルス イラン支援に25億円拠出へ 外相 NHK [ひとこと言わねば]
新型コロナウイルス イラン支援に25億円拠出へ 外相 #nhk_news https://t.co/cji5nuEESu
— NHKニュース (@nhk_news) March 17, 2020
「08年の金融危機より深刻」 大胆な経済対策に含み=文大統領 聯合ニュース [ひとこと言わねば]
「08年の金融危機より深刻」 大胆な経済対策に含み=文大統領 https://t.co/MpggcO8Hsw
— yonhapnews (@yonhapjp) March 17, 2020
トランプ大統領「感染拡大 7月か8月までという見方も」 NHK [ひとこと言わねば]
トランプ大統領「感染拡大 7月か8月までという見方も」 #nhk_news https://t.co/sdHPULF0A7
— NHKニュース (@nhk_news) March 16, 2020
ナイル川クルーズ船ツアー 参加者など感染確認26人に 厚労省 NHK [ひとこと言わねば]
ナイル川クルーズ船ツアー 参加者など感染確認26人に 厚労省 #nhk_news https://t.co/7BV5x5oay3
— NHKニュース (@nhk_news) March 16, 2020
欧州の渡航中止勧告地域 滞在歴ある外国人も入国拒否へ 政府 NHK [ひとこと言わねば]
欧州の渡航中止勧告地域 滞在歴ある外国人も入国拒否へ 政府 #nhk_news https://t.co/nACZ3ibgyn
— NHKニュース (@nhk_news) March 16, 2020
新型ウイルス 世界各国で対策強化 入国禁止や外出大幅制限も #nhk_news https://t.co/r3dBAdR61G
— NHKニュース (@nhk_news) March 17, 2020
中小企業から資金繰りの相談急増 政府 金融機関に支援要請 NHK [ひとこと言わねば]
中小企業から資金繰りの相談急増 政府 金融機関に支援要請 #nhk_news https://t.co/Tp9vqMju8v
— NHKニュース (@nhk_news) March 16, 2020
迫る年度末 企業の資金繰り支え経済悪化回避が課題 政府・日銀 #nhk_news https://t.co/V4k0PrK89Y
— NHKニュース (@nhk_news) March 16, 2020
「経済対策は前例にとらわれず考える」経済再生相 #nhk_news https://t.co/KdE2XhDORs
— NHKニュース (@nhk_news) March 17, 2020
NY株式市場 3000ドル超の値下がり 過去最大 NHK [ひとこと言わねば]
【速報 JUST IN 】NY株式市場 3000ドル超の値下がり 過去最大 #nhk_news https://t.co/wKto9VVs46
— NHKニュース (@nhk_news) March 16, 2020
株価 600円超下落 ことし最安値を更新 #nhk_news https://t.co/cliRCgcmRg
— NHKニュース (@nhk_news) March 17, 2020
原油価格 大幅値下がり 30ドル割れ 世界経済に懸念 #nhk_news https://t.co/SESMLYguyf
— NHKニュース (@nhk_news) March 16, 2020
フィリピン、株式市場の取引を無期限停止 新型コロナ対策で https://t.co/kV1BwtXwjK
— ロイター (@ReutersJapan) March 17, 2020
新型コロナウイルス 中国以外の感染者数が中国上回る WHO NHK [ひとこと言わねば]
新型コロナウイルス 中国以外の感染者数が中国上回る WHO #nhk_news https://t.co/T01KtKa0My
— NHKニュース (@nhk_news) March 16, 2020
アイルランドの新型コロナ感染者、3月末には1.5万人と予想 https://t.co/BgIO3sfeYS
— ロイター (@ReutersJapan) March 17, 2020