16/11/15 10:42:17.30 hp7RgOLv.net
SJW=function(A,B,C,D){ # Relative Risk cohort study
DR=matrix(c(A,B,C,D),2,byrow=TRUE)
fl=function(x)Epi::twoby2(DR,alpha=x,print=FALSE)$measures[1,2]
fu=function(x)Epi::twoby2(DR,alpha=x,print=FALSE)$measures[1,3]
le=100
xx=seq(0,1,length=le) ; xx=xx[c(-1,-le)]
y1=sapply(xx,fl)
y2=sapply(xx,fu)
xlim=c(min(y1),max(y2))
plot(y1,xx,log="x",xlim=xlim,type="l",lwd=2,xlab="Relative Risk(log scale)",ylab="p.value",main="P.value function(Cohort Study)")
lines(y2,xx,lwd=2)
abline(v=1,lty=3)
abline(h=0.05,lty=3,col="gray")
Epi::twoby2(DR)
return(DR)
}
Subjects CVD events % HR p value 95%CI
Placebo 2333 282 12.1
Empagliflozin 10mg 2345 243 10.4 0.85 0.07 0.72-1.01
Empagliflozin 25mg 2342 247 10.5 0.86 0.09 0.73-1.02
Empagliflozin 10mg+25mg 4687 490 10.5 0.86 0.04 0.74-0.99
dev.off()
layout(matrix(c(1,2,3,3),2,byrow=TRUE))
E10=SJW(243,2345-243,282,2333-282) ; chisq.test(E10,correct = FALSE) ; prop.test(c(243,282),c(2345,2333))
E25=SJW(247,2342-247,282,2333-282) ; chisq.test(E25,correct = FALSE) ; prop.test(c(247,282),c(2342,2333))
E1025=SJW(490,4687-490,282,2333-282) ; chisq.test(E1025,correct=FALSE) ; prop.test(c(490,282),c(4687,2333))