21/01/17 19:15:17.92 GiUUK9x4.net
有症者の方の数字もadjusted odds ratioと乖離した。
> r1=249
> n1=14173
> r2=15
> n2=6614
>
> k=1e6
> neg=rbeta(k,1+r1,1+n1-r1) # negative history
> pos=rbeta(k,1+r2,1+n2-r2) # positive history
> neg.o=neg/(1-neg)
> pos.o=pos/(1-pos)
> or=pos.o/neg.o
> quantile(or, c(0.025,0.975))
2.5% 97.5%
0.07616186 0.21299638
> summary(or)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.03219 0.11050 0.13240 0.13550 0.15710 0.37640
>