C1 [<- ^ abs as.vector c cbind exp for library log matrix mean outer rbinom rep return rmvn rnorm sample set.seed sum T C1a eX mean sum Y Z C2 [<- ^ abline apply as.character boxplot c colnames<- ests mean n par paste range rbind rep round rownames<- sqrt title var C3 [<- ~ al as.vector be binomial cbind expit fitted for glm matrix mean Mu n rbinom rmvn rnorm set.seed sig Sigma sum C4 [<- ^ abline apply as.character boxplot c colnames<- ests2 mean n par paste range rbind rep round rownames<- sqrt title var C5 [<- ~ al as.vector be binomial cbind data.frame expit fitted for glm lm matrix mean Mu n predict rbinom rmvn rnorm set.seed sig Sigma sum C6 [<- ^ abline apply as.character boxplot c colnames<- ests3 mean n par paste range rbind rep round rownames<- sqrt title var