### R code from vignette source 'C:/Bendix/teach/NSCE/2022/pracs/surv-tongue-sol.rnw'

###################################################
### code chunk number 1: surv-tongue-sol.rnw:6-9
###################################################
N <- c(130,78,45,33,25,19,12)
D <- c(45,24,5,2,1,0,0)
L <- c(7,9,7,6,5,7,6)


###################################################
### code chunk number 2: surv-tongue-sol.rnw:14-20
###################################################
res <- data.frame( N=N, D=D, L=L,
                   eff.den =               N-L/2,
                  pr.death =            D/(N-L/2),
                   pr.surv =          1-D/(N-L/2),
                  cum.surv = cumprod( 1-D/(N-L/2) ) )
round( res, 3 )


###################################################
### code chunk number 3: curv
###################################################
plot( 0:7, c(1,res$cum.surv), pch=16, type="b", ylim=0:1,
      ylab="Survival", xlab="Time since diagnosis" )
abline( h=c(1:3/4) )