16/12/01 03:11:40.08 EFWMKZr0.net
有病 pN 無病(1-p)N
検査陽性 a TP se b FP 1-sp
検査陰性 c FN 1-se d TN sp
v=1/a + 1/b + 1/c + 1/d
se=a/(a+c)
sp=d/(b+d)
pLH=TP/FP=se/(1-sp)
nLH=FN/TN=(1-se)/sp
OR=pLH/nLH
S=(se/(1-se)) / (sp/(1-sp))
# log(OR) = alpha + beta*log(S) + error
y=log(OR)
x=log(S)
res.lm = lm(y~x,weights=1/v)
summary(res.lm)