\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlstd{fit2}\hlkwb{<-}\hlkwd{lm}\hlstd{(Y}\hlopt{~}\hlstd{X1}\hlopt{+}\hlstd{Z}\hlopt{+}\hlstd{Z}\hlopt{:}\hlstd{X1)}
\hlstd{mu0.2}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit2,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,}\hlkwc{Z}\hlstd{=}\hlnum{0}\hlstd{)))}
\hlstd{mu1.2}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit2,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,}\hlkwc{Z}\hlstd{=}\hlnum{1}\hlstd{)))}
\hlkwd{c}\hlstd{(mu0.2,mu1.2,mu1.2}\hlopt{-}\hlstd{mu0.2)}
\end{alltt}
\begin{verbatim}
+ [1]  4.286084  2.744744 -1.541341
\end{verbatim}
\begin{alltt}
\hlstd{fit3}\hlkwb{<-}\hlkwd{lm}\hlstd{(Y}\hlopt{~}\hlstd{X2}\hlopt{+}\hlstd{Z}\hlopt{+}\hlstd{Z}\hlopt{:}\hlstd{X2)}
\hlstd{mu0.3}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit3,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,}\hlkwc{Z}\hlstd{=}\hlnum{0}\hlstd{)))}
\hlstd{mu1.3}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit3,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,}\hlkwc{Z}\hlstd{=}\hlnum{1}\hlstd{)))}
\hlkwd{c}\hlstd{(mu0.3,mu1.3,mu1.3}\hlopt{-}\hlstd{mu0.3)}
\end{alltt}
\begin{verbatim}
+ [1]  4.331889  2.647514 -1.684375
\end{verbatim}
\begin{alltt}
\hlstd{fit4}\hlkwb{<-}\hlkwd{lm}\hlstd{(Y}\hlopt{~}\hlstd{Z)}
\hlstd{mu0.4}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit4,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,}\hlkwc{Z}\hlstd{=}\hlnum{0}\hlstd{)))}
\hlstd{mu1.4}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit4,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,}\hlkwc{Z}\hlstd{=}\hlnum{1}\hlstd{)))}
\hlkwd{c}\hlstd{(mu0.4,mu1.4,mu1.4}\hlopt{-}\hlstd{mu0.4)}
\end{alltt}
\begin{verbatim}
+ [1]  4.288448  2.770900 -1.517548
\end{verbatim}
\end{kframe}
\end{knitrout}