18/12/11 12:05:49.81 k3SoZWT3.net
白色光とNBI拡大でのピロリ感染正診率
# pyroli陽性 白色光26/33 NBI拡大30/33
r1=26;r2=30;n1=n2=33
# pyroli陰性 白色光12/23 NBI拡大19/23
r1=12;r2=19;n1=n2=23
χ二乗検定 p-value = 0.05914
フィッシャーテストp-value = 0.05736でもどちらも有意差でないな。
prop.test(c(r1,r2),c(n1,n2))
Fisher.test(c(r1,r2),c(n1,n2))
ベイズでやる(JAGSでMCMC）と
ピロリ菌陰性のときの正診率の差の95%CIが０未満になるから有意差ありになる。
lower upper
-0.52247237 -0.02749815
js=as.matrix(jags4prop(r1,r2,n1,n2))
plotPost(js[,1],showMode=T)
plotPost(js[,2],showMode=T)
plotPost(js[,1]-js[,2],compVal=0,col=sample(colours(),1))
HDInterval::hdi(js[,1]-js[,2])
HDInterval::hdi(js[,1]/js[,2])