\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlstd{fit1}\hlkwb{<-}\hlkwd{lm}\hlstd{(Y}\hlopt{~}\hlstd{X1}\hlopt{+}\hlstd{X2}\hlopt{+}\hlstd{Z1}\hlopt{+}\hlstd{Z2);}\hlkwd{round}\hlstd{(}\hlkwd{coef}\hlstd{(}\hlkwd{summary}\hlstd{(fit1)),}\hlnum{6}\hlstd{)}
\end{alltt}
\begin{verbatim}
+             Estimate Std. Error   t value Pr(>|t|)
+ (Intercept) 0.213573   0.078020  2.737426 0.006203
+ X1          0.910646   0.043347 21.008092 0.000000
+ X2          1.037746   0.030916 33.566295 0.000000
+ Z1          0.777611   0.087386  8.898532 0.000000
+ Z2          0.978058   0.062915 15.545821 0.000000
\end{verbatim}
\begin{alltt}
\hlstd{mu.00.hat}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit1,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,X2,}\hlkwc{Z1}\hlstd{=}\hlnum{0}\hlstd{,}\hlkwc{Z2}\hlstd{=}\hlnum{0}\hlstd{)))}
\hlstd{mu.10.hat}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit1,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,X2,}\hlkwc{Z1}\hlstd{=}\hlnum{1}\hlstd{,}\hlkwc{Z2}\hlstd{=}\hlnum{0}\hlstd{)))}
\hlstd{mu.01.hat}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit1,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,X2,}\hlkwc{Z1}\hlstd{=}\hlnum{0}\hlstd{,}\hlkwc{Z2}\hlstd{=}\hlnum{1}\hlstd{)))}
\hlstd{mu.11.hat}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit1,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,X2,}\hlkwc{Z1}\hlstd{=}\hlnum{1}\hlstd{,}\hlkwc{Z2}\hlstd{=}\hlnum{1}\hlstd{)))}
\hlkwd{c}\hlstd{(mu.00.hat,mu.10.hat,mu.01.hat,mu.11.hat)}
\end{alltt}
\begin{verbatim}
+ [1] 3.446219 4.223831 4.424277 5.201888
\end{verbatim}
\end{kframe}
\end{knitrout}