\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{Z}\hlopt{+}\hlstd{Z}\hlopt{:}\hlstd{X1}\hlopt{+}\hlstd{Z}\hlopt{:}\hlstd{X2)}
\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) 11.206903   0.233694  47.955492        0
+ X1          -3.216690   0.118417 -27.164058        0
+ X2           1.775659   0.059201  29.993931        0
+ Z            3.109948   0.285001  10.912047        0
+ X1:Z        -0.927392   0.145289  -6.383066        0
+ X2:Z         2.018438   0.072729  27.752731        0
\end{verbatim}
\begin{alltt}
\hlstd{mu0.1}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit1,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,X2,}\hlkwc{Z}\hlstd{=}\hlnum{0}\hlstd{)))}
\hlstd{mu1.1}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit1,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,X2,}\hlkwc{Z}\hlstd{=}\hlnum{1}\hlstd{)))}
\hlkwd{c}\hlstd{(mu0.1,mu1.1,mu1.1}\hlopt{-}\hlstd{mu0.1)}
\end{alltt}
\begin{verbatim}
+ [1]  4.435704  2.565313 -1.870391
\end{verbatim}
\end{kframe}
\end{knitrout}