[過去ﾛｸﾞ] 数学 統計に詳しい人が語るコロナウイルス (1002ﾚｽ)
上下前次1-新
抽出解除 必死ﾁｪｯｶｰ(本家) (べ) 自ID ﾚｽ栞 あぼーん
このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索 歴削→次ｽﾚ 栞削→次ｽﾚ 過去ﾛｸﾞﾒﾆｭｰ
464(2): 2020/05/07(木)02:09 ID:VnQvkZ57(1/8) AAS
>>463
Wilsonの式を使ったんだろうね。
method x n mean lower upper
1 agresti-coull 3 312 0.0096 0.0019 0.029
2 asymptotic 3 312 0.0096 -0.0012 0.020
3 bayes 3 312 0.0112 0.0016 0.023
4 cloglog 3 312 0.0096 0.0027 0.026
省8
465(1): 2020/05/07(木)02:30 ID:VnQvkZ57(2/8) AAS
Wilsonの式での95％信頼区間幅を1％以内にしたいなら
サンプルサイズは38411が必要という計算になった。
Rでの算出プログラム
library(binom)
x=3
n=312
binom.wilson(x,n)
省11
467: 2020/05/07(木)18:13 ID:VnQvkZ57(3/8) AAS
>>466
95％信頼区間の下限境界は
事前分布をJeffreysで
> qbeta(0.95,0.5+50,0.5,lower=F)
[1] 0.9624989
事前分布を一様分布で
> qbeta(0.95,1+50,1,lower=F)
省1
468: 2020/05/07(木)18:52 ID:VnQvkZ57(4/8) AAS
特異度の95%の信頼区間の下限値を0.99にするのに必要なサンプルサイズは
事前分布を一様分布で
> uniroot(function(x) fn(x)-0.99, c(100,500))$root
[1] 297.0728
Jeffreysで
> uniroot(function(x) fn(x,0.5,0.5)-0.99, c(100,500))$root
[1] 190.8606
省8
469(1): 2020/05/07(木)20:13 ID:VnQvkZ57(5/8) AAS
外部ﾘﾝｸ:www.buzzfeed.com
のデータを使って、不明なものは一様分布（ベータ分布の形状母数(1,1)）に事前分布を設定してMCMCしてみる。
x=c(3,33)
n=c(312,1000)
m=50
N=length(n)
shape1=1
省12
471(1): 2020/05/07(木)22:36 ID:VnQvkZ57(6/8) AAS
>>470
同一キットじゃないから、ご指摘の通り。
しかも神戸大の方では陰性検体での確認はされていないから、神戸大方の陽性率が高いのは偽陽性を含む可能性もあるね。
大阪市大だけのデータでやってみると。
画像ﾘﾝｸ[png]:i.imgur.com
472: 2020/05/07(木)22:46 ID:VnQvkZ57(7/8) AAS
>>471
使える情報が少なすぎて信頼区間が広すぎ。
有病率:0-87%
キットの感度:0-70%
ってコイントス変わらんな。
473: 2020/05/07(木)22:48 ID:VnQvkZ57(8/8) AAS
結局のところ、断定的な結論は出せないねということを数字で確認しているだけｗ
上下前次1-新書関写板覧索設栞歴
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ
ぬこの手 ぬこTOP 0.054s