感度70%特異度90%で
有病率と陽性的中率・陰性的中率の関係をグラフにしてみた。
https://i.imgur.com/HGqPv2y.jpg
陽性的中率が0.8になるのは有病率が0.63のとき
そのRのコードはこれ。
rm(list=ls())
pr2pv <- function( # prevalence to predicative value
pr ,# prevalence
sn=0.7, # sensitibity=TP/(TP+FN)
sp=0.9) # specificity=TN/(TN+FP)
{
N=1 # polutaion million, billion,or any proper unit
si=pr*N # sick population
he=(1-pr)*N # healthy population
TP=si*sn
FN=si*(1-sn)
TN=he*sn
FP=he*(1-sn)
PPV=TP/(TP+FP)
NPV=TN/(TN+FN)
PV=c(PPV=PPV,NPV=NPV)
return(PV)
}
prev=seq(1e-7,1,length.out = 1000)
plot(prev,sapply(prev, function(x) pr2pv(x)['PPV']),bty='l',type='l',
ylab='predicative vale',xlab='prevalence(log)',main='sensitity=0.7,specificity=0.9',log='x',lwd=2)
lines(prev,sapply(prev, function(x) pr2pv(x)['NPV']),lty=3,lwd=2)
legend('center',bty='n',legend=c('Posivive Predicative Value','Negative Predicative Value'),lty=c(1,3),lwd=2)
abline(h=0.8,col='gray')
uniroot(function(x) pr2pv(x)['PPV']-0.8, c(0,1))
数学 統計に詳しい人が語るコロナウイルス
■ このスレッドは過去ログ倉庫に格納されています
2020/03/10(火) 11:16:04.35ID:H1fx2jVB
■ このスレッドは過去ログ倉庫に格納されています
ニュース
- 「いいの?前科ついちゃうよ」万引きした女子大学生から10万円を脅し取ったか 元コンビニ店長の男(53)逮捕 [煮卵★]
- 大谷翔平、育休でチームを離脱 球団が発表 第2子誕生へ…週末には復帰予定 長女誕生から1年 [(´?ω?`)知らんがな★]
- ランドセルにくぎ刺される「国に帰れ」など言われ、転校を余儀なくされた海外からの転校生 仙台市教育委員会が「いじめ重大事態」と認定 [煮卵★]
- 【芸能】立川志らく「ヤクザの事務所が隣近所にあれば、強盗だとか闇バイトが来ない」 暴力団の“利点”指摘に違和感 [冬月記者★]
- 【埼玉県警】国道で持ち運び可能なオービス盗まれる 速度取り締まり中 [nita★]
- 面識ない女子高生に覚醒剤注射し性的暴行、被告に懲役15年判決…福岡地裁小倉支部「身勝手極まりない犯行」 [少考さん★]
- 【地上波/DAZNほか】 FIFAワールドカップ2026 総合スレ★105【メキシコ/カナダ/アメリカ】
- 【地上波/DAZNほか】 FIFAワールドカップ2026 総合スレ★106【メキシコ/カナダ/アメリカ】
- 【MLB】ドジャース vs オリオールズ
- 函館競馬1回3日目
- 巨専】
- はません ★3
- 【高市】ナフサと石油がないは虚偽。レアアースも肥料もなく、気象異常で小麦もないのが確定。だから真実を話せと政府ヒアリングされる [784319933]
- 【悲報】大阪万博の日本製EVバスさん、不具合多発で『中国製EVバスが悪い』と歴史修正されて終わる [153736977]
- 【高市悲報】ヒゲの佐藤「東北大学に行きたかったが貧乏なので防衛大に行きました😤」過去の記事が発掘される [359965264]
- トリッカルって生活保護者を馬鹿にしてたってマジ?
- 円終了。外資が外貨を担保に円を借りて投資すると円安でノーレバでレバ3倍。更にファンドのレバで安全に30倍とかになる。だから円安高市 [784319933]
- 愛国保守「中国と戦うどォ😭」👈核600発に230万の軍隊持つレアアース,薬,肥料の輸出大国と日本が対立するメリット、不明… [941632843]