### R code from vignette source 'ATCB-vitE-sol.rnw'

###################################################
### code chunk number 1: ATCB-vitE-sol.rnw:36-42
###################################################
rate <- c(11.6,17.8)
D <- c(99,151)
names(rate) <- names(D) <- c("VitE","Plc")
D / (rate/10000)
cr5 <- 1 - exp( -rate/10000 * 5 )
round( cr5*100, 2 )


###################################################
### code chunk number 2: ATCB-vitE-sol.rnw:55-57
###################################################
rate[1]/rate[2]
cr5[1]/cr5[2]


###################################################
### code chunk number 3: ATCB-vitE-sol.rnw:64-66
###################################################
rate[1]-rate[2]
round( (cr5[1]-cr5[2])*100, 2 )


###################################################
### code chunk number 4: ATCB-vitE-sol.rnw:76-78
###################################################
(rate["Plc"]-rate["VitE"])/rate["Plc"]
(cr5["Plc"]-cr5["VitE"])/cr5["Plc"]