### R code from vignette source 'cond-survival-sol.rnw'

###################################################
### code chunk number 1: cond-survival-sol.rnw:23-38
###################################################
prob <- matrix( c(
1,43.44,42.13,36.60,34.29,
2,22.80,19.11,24.00,21.86,
3,16.74,14.60,21.02,15.67,
4,13.84,10.62,15.59,13.54,
5,11.00,8.69,14.55,11.40,
6,10.13,7.36,9.95,11.17,
7,8.67,5.65,11.37,8.99,
8,7.97,5.51,8.69,8.55,
9,7.42,5.37,10.07,8.14,
10,7.75,5.94,5.16,7.26,
11,4.91,5.66,7.14,2.57,
12,6.72,5.42,6.06,5.63,
13,6.20,6.25,5.00,2.13), ncol=5, byrow=T )
prob


###################################################
### code chunk number 2: cond-survival-sol.rnw:42-45
###################################################
prob <- prob[,-1]
colnames(prob) <- c("Col.M","Col.F","Rect.M","Rect.F")
prob


###################################################
### code chunk number 3: cond-survival-sol.rnw:52-57
###################################################
prob <- prob/100
surv1 <- 1-prob[1,]
surv3 <- surv1 * (1-prob[2,]) * (1-prob[3,])
surv5 <- surv3 * (1-prob[4,]) * (1-prob[5,])
( DK.surv <- rbind( surv1, surv3, surv5 ) )


###################################################
### code chunk number 4: cond-survival-sol.rnw:60-61
###################################################
round( DK.surv*100, 1 )


###################################################
### code chunk number 5: cond-survival-sol.rnw:94-98
###################################################
cond3 <- surv3/surv1
cond5 <- surv5/surv1
DK.cond <- rbind( cond3, cond5 )
round( DK.cond*100, 1 )


###################################################
### code chunk number 6: cond-survival-sol.rnw:113-114
###################################################
round( (S.DK <- apply( 1-prob, 2, cumprod ) )*100, 1 )


###################################################
### code chunk number 7: cond-survival-sol.rnw:118-119
###################################################
sweep( as.array(S.DK), 2, S.DK[1,], "/" )


###################################################
### code chunk number 8: cond-survival-sol.rnw:124-131
###################################################
US.surv <- matrix( c(
1, 68, 69, 74, 74,
3, 46, 48, 48, 50,
5, 36, 39, 35, 39), ncol=5, byrow=TRUE )
US.surv <- US.surv[,-1]/100
colnames( US.surv ) <- colnames( DK.surv )
US.surv


###################################################
### code chunk number 9: cond-survival-sol.rnw:134-138
###################################################
acond3 <- US.surv[2,]/US.surv[1,]
acond5 <- US.surv[3,]/US.surv[1,]
US.cond <- rbind(acond3,acond5)
round( US.cond*100, 1 )


###################################################
### code chunk number 10: cond-survival-sol.rnw:158-160
###################################################
round( US.surv/DK.surv, 2 )
round( US.cond/DK.cond, 2 )