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

###################################################
### code chunk number 1: Lifetable-sol.rnw:27-36
###################################################
N <- c(234, 207, 169, 129, 105,  85,  73,  62,  49, 34)
D <- c(24, 27, 31, 17, 7, 6, 5, 3, 2, 4)
L <- c( 3, 11, 9, 7,13, 6, 6, 10, 13, 6)
res <- cbind( N, D, L, N-L/2,
                    D/(N-L/2),
                  1-D/(N-L/2),
          cumprod(1-D/(N-L/2)) )
colnames(res)[4:7] <- c("eff.N","mort","surv","Surv")
round( res, 5 )


###################################################
### code chunk number 2: surv
###################################################
S <- cumprod(1-D/(N-L/2))
plot( 0:10, c(1,S), pch=16, type="b", ylim=0:1,
      ylab="Survival", xlab="Time since diagnosis" )