### R code from vignette source 'Prev-trial1-sol.rnw'

###################################################
### code chunk number 1: Prev-trial1-sol.rnw:9-11
###################################################
cases <- c(474, 402)
rates <- c(56.3, 47.5) # per 10000 years


###################################################
### code chunk number 2: Prev-trial1-sol.rnw:14-15
###################################################
names( rates ) <- c("BetaCarotene","Placebo")


###################################################
### code chunk number 3: Prev-trial1-sol.rnw:26-28
###################################################
pyears <- (cases/rates)*10000
pyears


###################################################
### code chunk number 4: Prev-trial1-sol.rnw:38-43
###################################################
ratio <- rates[1]/rates[2]
SE.logr <- sqrt(sum(1/cases))
ratio.95low <- ratio/exp(1.96*SE.logr)
ratio.95up  <- ratio*exp(1.96*SE.logr)
cbind(ratio, SE.logr, ratio.95low, ratio.95up)


###################################################
### code chunk number 5: Prev-trial1-sol.rnw:63-68
###################################################
diff <- rates[1] - rates[2]
SE.diff <- sqrt(sum(rates^2/cases))
diff.95low <- diff - 1.96*SE.diff
diff.95up <- diff + 1.96*SE.diff
cbind(diff, SE.diff, diff.95low, diff.95up)


###################################################
### code chunk number 6: Prev-trial1-sol.rnw:77-81
###################################################
Z <- diff/SE.diff
P <- 1 - pchisq( Z^2, 1 )
test.diff <- cbind(Z, P)
test.diff


###################################################
### code chunk number 7: Prev-trial1-sol.rnw:86-89
###################################################
Z <- log(ratio)/SE.logr
P <- 1 - pchisq( Z^2, 1 )
( test.ratio <- cbind(Z, P) )


###################################################
### code chunk number 8: Prev-trial1-sol.rnw:93-96
###################################################
tt <- rbind( test.diff, test.ratio )
rownames( tt ) <- c("diff","ratio")
tt