18/04/26 18:58:06.22 UTg/BRA4.net
# URLﾘﾝｸ(www.igaku-shoin.co.jp)
# 相対リスク＝1.82を論文記載の対照群のリスク（36人/10万人）
# 1.82（95％信頼区間0.79-4.20）
rr=1.82
r0=36/10^5
r1=rr*r0
nn=seq(10^4,10^5,by=100)
rr2p.e <- function(n){
Epi::twoby2(matrix(c(n*r1,n*r0,n-n*r1,n-n*r0),ncol=2),print=FALSE)$p.value
}
plot(nn,sapply(nn,rr2p.e),type='l') ; abline(h=0.05, lty=3)
n=uniroot(function(n) rr2p.e(n)-0.05, c(10^4,10^5))$root
mat=matrix(c(n*r1,n*r0,n-n*r1,n-n*r0),ncol=2)
Epi::twoby2(mat)