C3 & [<- ~ < > abs al as.vector be cbind coef d data.frame dnorm for lm matrix mean Mu nZ pnorm predict rep rmvn rnorm Sigma sigY sigZ solve sum t T zgrid C4 abline agest.samp apply c expression lines par plot quantile zest.dr.samp zest.samp zgrid ztrue C1 & [<- ^ ~ < > abs as.vector c cbind coef data.frame for length library lm matrix mean outer pnorm predict rep rmvn rnorm round sample seq set.seed sum T C2 & [<- < > abs box for hist length mean pnorm ps.fitted rep round seq sum T Y Z C2a be c lines Mu par plot sum zgrid ztrue A1 $ & [<- ~ < > al as.vector be cbind data.frame dnorm fitted for lm matrix mean Mu nZ predict rep rmvn rnorm sd Sigma sigY sigZ sum T zgrid A0 [<- ^ abs as.vector c cbind length library mean outer rep rmvn rnorm round sample seq set.seed sum T A4 $ & [<- ~ < > al as.vector be cbind dnorm for lm matrix Mu nZ predict rep rmvn rnorm sd Sigma sigY sigZ sum T zgrid A6 $ & [<- ^ ~ < > al as.vector be cbind d data.frame dnorm fitted for lm matrix mean Mu n nreps nZ predict rep rmvn rnorm sd Sigma sigY sigZ sum T X1 Z0 zgrid A8 $ [<- ~ al as.vector be cbind data.frame dnorm for lm matrix Mu nZ predict rmvn rnorm sd Sigma sigY sigZ zgrid A10 $ [<- ^ ~ al as.vector be cbind data.frame dnorm for lm matrix Mu nZ predict rmvn rnorm sd Sigma sigY sigZ Z0 zgrid