16/11/15 13:11:45.79 hp7RgOLv.net
# p値関数 vcd::oddsratio, confintで信頼区間境界からグラフ化
CIfunc=function(A,B,C,D){ # Oddsratio 症例対照study
DO=matrix(c(A,B,C,D),2,byrow=TRUE)
oDO=vcd::oddsratio(DO,log=FALSE)
# confint(oDO,level=0.95)
f1=function(x)confint(oDO,level=1-x)[1]
f2=function(x)confint(oDO,level=1-x)[2]
xx=seq(0,1,by=0.01) ; xx=xx[-1]
y1=sapply(xx,f1)
y2=sapply(xx,f2)
xlim=c(min(y1),max(y2))
plot(y1,xx,xlim=xlim,log="x",type="l",lwd=2,xlab="Odds Ratio(log scale)",ylab="p.value",main="P.value function(Case-Control)")
lines(y2,xx,lwd=2)
abline(v=1,lty=3)
abline(h=0.05,lty=3,col="gray")
Epi::twoby2(DO)
return(DO)
}
DO=CIfunc(4,386,4,1250)
chisq.test(DO,correct=FALSE)