### R code from vignette source 'std-colon-sol.rnw'

###################################################
### code chunk number 1: std-colon-sol.rnw:28-37
###################################################
M <- matrix(
c(10,1157,46.0,0.9,22,1109,41.9,2.0,0.44
,76,809,32.0,9.4,68,786,29.7,8.6,1.09
,305,455,18.0,67,288,524,19.8,55,1.22
,201,102,4.0,196,354,229,8.6,155,1.27), nrow=4, byrow=T )
M <- data.frame(M)
names(M) <- c("mca","mpy","mp","mr",
              "fca","fpy","fp","fr","rr")
M


###################################################
### code chunk number 2: std-colon-sol.rnw:45-51
###################################################
rates <-
with( M, c( sum(mca)/sum(mpy)*100,
            sum(fca)/sum(fpy)*100 ) )
rates[3] <- rates[1]/rates[2]
names(rates) <- c("M rate","F rate","M/F RR")
round( rates, 2 )


###################################################
### code chunk number 3: std-colon-sol.rnw:55-62
###################################################
wm <- with( M, mpy/sum(mpy) )
rates <-
with( M, c( sum(mca/mpy*wm)*100,
            sum(fca/fpy*wm)*100 ) )
rates[3] <- rates[1]/rates[2]
names(rates) <- c("M rate","F rate","M/F RR")
round( rates, 2 )


###################################################
### code chunk number 4: std-colon-sol.rnw:68-68
###################################################



###################################################
### code chunk number 5: std-colon-sol.rnw:69-71
###################################################
WSP <- c(96,24,100,90,90,80,80,60,60,60,60,50,40,40,30,20,10,5,3,2)
WSP


###################################################
### code chunk number 6: std-colon-sol.rnw:77-89
###################################################
wt <- sum(WSP[1:8])
wt[2] <- sum(WSP[9:12])
wt[3] <- sum(WSP[13:16])
wt[4] <- sum(WSP[17:19])
wt <- wt/sum(wt)
wt
rates <-
with( M, c( sum(mca/mpy*wt)*100,
            sum(fca/fpy*wt)*100 ) )
rates[3] <- rates[1]/rates[2]
names(rates) <- c("M rate","F rate","M/F RR")
round( rates, 2 )


###################################################
### code chunk number 7: std-colon-sol.rnw:95-102
###################################################
wy <- c(35,20,20)
rates <-
with( M[1:3,], c( sum(mca/mpy*wy)/1000,
                  sum(fca/fpy*wy)/1000 ) )
rates[3] <- rates[1]/rates[2]
names(rates) <- c("M cum.rate","F cum.rate","M/F RR")
round( rates, 4 )


###################################################
### code chunk number 8: std-colon-sol.rnw:107-111
###################################################
rates <- 1 - exp( -rates )
rates[3] <- rates[1]/rates[2]
names(rates) <- c("M cum.risk","F cum.risk","M/F RR")
round( rates, 4 )