20/03/11 06:30:05.15 hVKkfTiV.net
>>14
ご指摘ありがとうございます。
プログラムにバグがありました。
>10は撤回します。
正しくは
URLﾘﾝｸ(i.imgur.com)
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になりました。