cross tableなデータをbootさせてみた

set.seed(101)
m<-100
replicate<-numeric(m)
for(i in 1:m)replicate[i]<-mean(fisher.test(matrix(boot.table(a),ncol=2))$p.value)
replicate
[1] 8.225812e-04 2.117310e-08 4.130529e-08 4.141086e-09
.......
[97] 2.097987e-05 4.882027e-05 3.701652e-04 1.034722e-05

こんな感じでRelative riskとかOdds比の95%CIを直接出せるようにしておこう。
boot.tableはgllmに含まれているcross tableをbootさせてくれるboot関数。
‥来週、時間取れればいいけど‥忙しいんだろうなぁ‥orz