### R code from vignette source 'Trial-tolb-sol.rnw'

###################################################
### code chunk number 1: Trial-tolb-sol.rnw:5-12
###################################################
library( Epi )
options(digits=3)
D <- c( 30,  21)
n <- c(204, 215)
names( D ) <-
names( n ) <- c("Tolbutamide", "Placebo")
cbind( D, n )


###################################################
### code chunk number 2: Trial-tolb-sol.rnw:16-18
###################################################
Q <- D/n
Q


###################################################
### code chunk number 3: Trial-tolb-sol.rnw:25-30
###################################################
QR <- Q[1]/Q[2]
SE.lqr <- sqrt( 1/D[1]-1/n[1] + 1/D[2]-1/n[2] )
QR.95lo <- QR / exp(1.96*SE.lqr)
QR.95up <- QR * exp(1.96*SE.lqr)
cbind( QR, SE.lqr, QR.95lo, QR.95up)


###################################################
### code chunk number 4: Trial-tolb-sol.rnw:35-40
###################################################
QD <- Q[1] - Q[2]
SE.qd <- sqrt( sum( Q*(1-Q)/n ) )
QD.95low <- QD - 1.96*SE.qd
QD.95up <- QD + 1.96*SE.qd
cbind( QD, SE.qd, QD.95low, QD.95up)


###################################################
### code chunk number 5: Trial-tolb-sol.rnw:49-50
###################################################
twoby2( cbind( D, n-D ) )