\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlstd{fit.psr.2}\hlkwb{<-}\hlkwd{lm}\hlstd{(Y}\hlopt{~}\hlstd{Z}\hlopt{+}\hlstd{Z}\hlopt{:}\hlstd{X1}\hlopt{+}\hlstd{Z}\hlopt{:}\hlstd{X2}\hlopt{+}\hlstd{Z}\hlopt{:}\hlstd{X1}\hlopt{:}\hlstd{X2}\hlopt{+}\hlstd{ps.fit.c)}
\hlkwd{round}\hlstd{(}\hlkwd{coef}\hlstd{(}\hlkwd{summary}\hlstd{(fit.psr.2)),}\hlnum{6}\hlstd{)}
\end{alltt}
\begin{verbatim}
+              Estimate Std. Error   t value Pr(>|t|)
+ (Intercept)  3.085174   0.027244 113.24310        0
+ Z            3.964332   0.135035  29.35781        0
+ ps.fit.c     2.126006   0.046459  45.76071        0
+ Z:X1        -0.874477   0.083430 -10.48155        0
+ Z:X2         1.894645   0.042690  44.38134        0
+ Z:X1:X2      0.627663   0.039762  15.78543        0
\end{verbatim}
\begin{alltt}
\hlstd{mu0.psr}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit.psr.2,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,X2,}\hlkwc{Z}\hlstd{=}\hlnum{0}\hlstd{,ps.fit.c)))}
\hlstd{mu1.psr}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit.psr.2,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,X2,}\hlkwc{Z}\hlstd{=}\hlnum{1}\hlstd{,ps.fit.c)))}
\hlkwd{c}\hlstd{(mu0.psr,mu1.psr,mu1.psr}\hlopt{-}\hlstd{mu0.psr)}
\end{alltt}
\begin{verbatim}
+ [1]  4.439440  2.532145 -1.907295
\end{verbatim}
\end{kframe}
\end{knitrout}