\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlstd{mu.x0}\hlkwb{<-} \hlstd{Xb} \hlopt{%*%} \hlstd{be}
\hlstd{mu.x1}\hlkwb{<-} \hlstd{Xb} \hlopt{%*%} \hlstd{be} \hlopt{+} \hlstd{Xp} \hlopt{%*%} \hlstd{psi}
\hlstd{mu.hat.a0}\hlkwb{<-}\hlkwd{sum}\hlstd{(W0}\hlopt{*}\hlstd{(Y}\hlopt{-}\hlstd{mu.x0))}\hlopt{+}\hlkwd{mean}\hlstd{(mu.x0)}
\hlstd{mu.hat.a1}\hlkwb{<-}\hlkwd{sum}\hlstd{(W1}\hlopt{*}\hlstd{(Y}\hlopt{-}\hlstd{mu.x1))}\hlopt{+}\hlkwd{mean}\hlstd{(mu.x1)}
\hlkwd{c}\hlstd{(mu.hat.a0,mu.hat.a1,mu.hat.a1}\hlopt{-}\hlstd{mu.hat.a0)}
\end{alltt}
\begin{verbatim}
+ [1]  4.455653  2.539759 -1.915894
\end{verbatim}
\begin{alltt}
\hlstd{mu.tilde.a0}\hlkwb{<-}\hlkwd{mean}\hlstd{(w0}\hlopt{*}\hlstd{(Y}\hlopt{-}\hlstd{mu.x0))}\hlopt{+}\hlkwd{mean}\hlstd{(mu.x0)}
\hlstd{mu.tilde.a1}\hlkwb{<-}\hlkwd{mean}\hlstd{(w1}\hlopt{*}\hlstd{(Y}\hlopt{-}\hlstd{mu.x1))}\hlopt{+}\hlkwd{mean}\hlstd{(mu.x1)}
\hlkwd{c}\hlstd{(mu.tilde.a0,mu.tilde.a1,mu.tilde.a1}\hlopt{-}\hlstd{mu.tilde.a0)}
\end{alltt}
\begin{verbatim}
+ [1]  4.455630  2.539749 -1.915881
\end{verbatim}
\end{kframe}
\end{knitrout}