### R code from vignette source 'surv-rate.rnw'

###################################################
### code chunk number 1: surv-rate.rnw:487-497
###################################################
D <- 12
Y <- 1276.3/1000
m0 <- glm( D   ~ 1, offset=log(Y), family=poisson )
m1 <- glm( D/Y ~ 1, weights=Y, family=poisson )
m2 <- glm( D/Y ~ 1, weights=Y, family=poisson(link=identity) )
library( Epi )
round( rbind( ci.lin( m0, E=T )[,c(1,2,5:7)],
              ci.lin( m1, E=T )[,c(1,2,5:7)],
              ci.lin( m2 )[,c(1,2,NA,5:6)] ), 3 )
round( c( 1/sqrt(D), sqrt(D)/Y ) , 3 )


