感度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
■ このスレッドは過去ログ倉庫に格納されています
ニュース
- 【速報】GMO不正アクセス 攻撃者が顧客の個人情報の一部をネットに公開「毎日正午に順次公開」★2 [Hitzeschleier★]
- 【速報】GMO、ランサム被害にあった顧客に新しいPC提供へ [Hitzeschleier★]
- 【兵庫】三田市の路上で17歳男性に約11分間わいせつ行為 85歳無職の男を逮捕 [煮卵★]
- 【声明】GMO「公開された情報にアクセス・拡散しないで」情報拡散者に対して法的措置へ [Hitzeschleier★]
- “万引きGメン”の男を私服警官が現行犯逮捕 勤務中のホームセンターで長靴盗んだか 容疑認める 広島 [煮卵★]
- サイバー攻撃ハッカー集団、GMOが27日12:00までに1億6000万払わなければ、顧客攻撃と声明 GMOは政府に防御要請、警察自衛隊へ ★2 [お断り★]
- トランプ「ディーゼルが高くて農家は不安よな トランプ、動きます」 [465483804]
- GMO熊谷社長、お気持ち表明「───今こそ、官民連携で日本を守ろう。」 [762037879]
- 辺野古転覆事故遺族が抗議団体と面談 「謝って、何もせず、批判が通り過ぎるのを待つ。ほとぼりが冷めたら元どおり。それは許されない」 [634830743]
- GMOハッカー集団、優しいwww「毎日お昼にウキウキウォッチング、個人情報を一人ずつ公開していく!」 [779857986]
- 楽天証券、不正アクセスで損害生じた場合50%補償 [237216734]
- 足クジータ…