- 17 名前：１３２人目の素数さん mailto:sage [2020/03/11(水) 06:30:05.15 ID:hVKkfTiV.net]
- >>14
ご指摘ありがとうございます。
プログラムにバグがありました。
>10は撤回します。
正しくは
https://i.imgur.com/jO5diiG.jpg
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*sp
FP=he*(1-sp)
PPV=TP/(TP+FP)
NPV=TN/(TN+FN)
PV=c(PPV=PPV,NPV=NPV)
return(PV)
}
ご指摘のとおり、有病率36.36%のときに感度0.7，特異度0.9で陽性的中率が0.8になりました。
|
 |
