16/11/19 22:31:01.03 +gOnxegu.net
# Exact法(Clopper-Pearson)信頼区間のconf.levelを変更した信頼区間を返す
# RR:Risk Ratio, l0,u0: conf.level=cl0での信頼区間をcl1の区間に変更
ConvExactCI=function(RR,l0,u0,cl0=0.95,cl1=0.90){ # 95%→90%に
L0 = function(n,RR,cl0) 1/(1+(n-n*RR+1)/(n*RR*qf((1-cl0)/2, 2*n*RR,2*(n-n*RR+1))))
fl = function(x,u)L0(x,RR,cl0)-u
(nL=uniroot(fl,u=l0,c(1,10^4))$root)
(rL=RR*nL)
(L1=L0(nL,RR,cl1))
U0 = function(n,RR,cl0) 1/(1+(n-n*RR)/((n*RR+1)*qf(1-(1-cl0)/2, 2*(n*RR+1), 2*(n-n*RR))))
fu = function(x,u)U0(x,RR,cl0)-u
(nU=uniroot(fu,u=u0,c(1,10^4))$root)
(nU=RR*nU)
(U1=U0(nU,RR,cl1))
CI=c(L1,U1,cl1)
names(CI)=c("lower","upper","conf.level")
return(round(CI,4))
}