### R code from vignette source 'rates-sol.rnw'

###################################################
### code chunk number 1: rates-sol.rnw:8-10
###################################################
cases <- c(6, 53)
pyears <- c(0.89, 18.19)


###################################################
### code chunk number 2: rates-sol.rnw:16-19
###################################################
rates <- cases/pyears
names(rates) <- c("Kuwait", "Egypt")
rates


###################################################
### code chunk number 3: rates-sol.rnw:25-29
###################################################
SE.r <- rates / sqrt(cases)
CL.low <- rates - 1.96*SE.r
CL.up <- rates + 1.96*SE.r
cbind(rates, SE.r, CL.low, CL.up)


###################################################
### code chunk number 4: rates-sol.rnw:36-40
###################################################
SE.logr <- sqrt(1/cases)
CL.low <- rates/exp(1.96*SE.logr)
CL.up <- rates*exp(1.96*SE.logr)
cbind(rates, SE.logr, CL.low, CL.up)