\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlstd{W}\hlkwb{<-}\hlstd{(ps.1}\hlopt{^}\hlstd{Z1}\hlopt{*}\hlstd{(}\hlnum{1}\hlopt{-}\hlstd{ps.1)}\hlopt{^}\hlstd{(}\hlnum{1}\hlopt{-}\hlstd{Z1))}\hlopt{*}\hlstd{(ps.2}\hlopt{^}\hlstd{Z2}\hlopt{*}\hlstd{(}\hlnum{1}\hlopt{-}\hlstd{ps.2)}\hlopt{^}\hlstd{(}\hlnum{1}\hlopt{-}\hlstd{Z2))}
\hlstd{W00}\hlkwb{<-}\hlstd{(}\hlnum{1}\hlopt{-}\hlstd{Z1)}\hlopt{*}\hlstd{(}\hlnum{1}\hlopt{-}\hlstd{Z2)}\hlopt{/}\hlstd{W}
\hlstd{W10}\hlkwb{<-}\hlstd{(Z1)}\hlopt{*}\hlstd{(}\hlnum{1}\hlopt{-}\hlstd{Z2)}\hlopt{/}\hlstd{W}
\hlstd{W01}\hlkwb{<-}\hlstd{(}\hlnum{1}\hlopt{-}\hlstd{Z1)}\hlopt{*}\hlstd{(Z2)}\hlopt{/}\hlstd{W}
\hlstd{W11}\hlkwb{<-}\hlstd{(Z1)}\hlopt{*}\hlstd{(Z2)}\hlopt{/}\hlstd{W}
\hlstd{mu.00.tilde}\hlkwb{<-}\hlkwd{mean}\hlstd{(W00}\hlopt{*}\hlstd{Y)}
\hlstd{mu.10.tilde}\hlkwb{<-}\hlkwd{mean}\hlstd{(W10}\hlopt{*}\hlstd{Y)}
\hlstd{mu.01.tilde}\hlkwb{<-}\hlkwd{mean}\hlstd{(W01}\hlopt{*}\hlstd{Y)}
\hlstd{mu.11.tilde}\hlkwb{<-}\hlkwd{mean}\hlstd{(W11}\hlopt{*}\hlstd{Y)}
\hlkwd{c}\hlstd{(mu.00.tilde,mu.10.tilde,mu.01.tilde,mu.11.tilde)}
\end{alltt}
\begin{verbatim}
+ [1] 2.115634 4.925918 3.323802 5.960125
\end{verbatim}
\end{kframe}
\end{knitrout}