\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlkwd{par}\hlstd{(}\hlkwc{mar}\hlstd{=}\hlkwd{c}\hlstd{(}\hlnum{2}\hlstd{,}\hlnum{2}\hlstd{,}\hlnum{0}\hlstd{,}\hlnum{0}\hlstd{))}
\hlkwd{boxplot}\hlstd{(Y}\hlopt{~}\hlstd{Z)}
\hlkwd{points}\hlstd{(}\hlnum{1}\hlopt{:}\hlnum{2}\hlstd{,}\hlkwd{c}\hlstd{(mu0,mu1),}\hlkwc{col}\hlstd{=}\hlstr{'red'}\hlstd{,}\hlkwc{pch}\hlstd{=}\hlnum{15}\hlstd{)}
\end{alltt}
\end{kframe}

{\centering \includegraphics[width=\maxwidth]{figure/CA3a-1} 

}


\begin{kframe}\begin{alltt}
\hlkwd{c}\hlstd{(}\hlkwd{mean}\hlstd{(Y[Z}\hlopt{==}\hlnum{0}\hlstd{]),}\hlkwd{mean}\hlstd{(Y[Z}\hlopt{==}\hlnum{1}\hlstd{]),}\hlkwd{mean}\hlstd{(Y[Z}\hlopt{==}\hlnum{1}\hlstd{])}\hlopt{-}\hlkwd{mean}\hlstd{(Y[Z}\hlopt{==}\hlnum{0}\hlstd{]))}
\end{alltt}
\begin{verbatim}
+ [1]  4.288448  2.770900 -1.517548
\end{verbatim}
\end{kframe}
\end{knitrout}