\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{4}\hlstd{,}\hlnum{4}\hlstd{,}\hlnum{3}\hlstd{,}\hlnum{0}\hlstd{))}
\hlstd{nvec}\hlkwb{<-}\hlkwd{as.character}\hlstd{(}\hlkwd{c}\hlstd{(}\hlnum{0}\hlstd{,}\hlkwd{rep}\hlstd{(}\hlnum{1}\hlopt{:}\hlnum{5}\hlstd{,}\hlkwc{each}\hlstd{=}\hlnum{2}\hlstd{)))}
\hlstd{nvec}\hlkwb{<-}\hlkwd{paste}\hlstd{(}\hlstr{'M'}\hlstd{,nvec,}\hlkwc{sep}\hlstd{=}\hlstr{''}\hlstd{)}
\hlstd{nvec[}\hlkwd{c}\hlstd{(}\hlnum{2}\hlstd{,}\hlnum{4}\hlstd{,}\hlnum{6}\hlstd{,}\hlnum{8}\hlstd{,}\hlnum{10}\hlstd{)]}\hlkwb{<-}\hlkwd{paste}\hlstd{(nvec[}\hlkwd{c}\hlstd{(}\hlnum{2}\hlstd{,}\hlnum{4}\hlstd{,}\hlnum{6}\hlstd{,}\hlnum{8}\hlstd{,}\hlnum{10}\hlstd{)],}\hlstr{'a'}\hlstd{,}\hlkwc{sep}\hlstd{=}\hlstr{''}\hlstd{)}
\hlstd{nvec[}\hlkwd{c}\hlstd{(}\hlnum{3}\hlstd{,}\hlnum{5}\hlstd{,}\hlnum{7}\hlstd{,}\hlnum{9}\hlstd{,}\hlnum{11}\hlstd{)]}\hlkwb{<-}\hlkwd{paste}\hlstd{(nvec[}\hlkwd{c}\hlstd{(}\hlnum{3}\hlstd{,}\hlnum{5}\hlstd{,}\hlnum{7}\hlstd{,}\hlnum{9}\hlstd{,}\hlnum{11}\hlstd{)],}\hlstr{'b'}\hlstd{,}\hlkwc{sep}\hlstd{=}\hlstr{''}\hlstd{)}
\hlkwd{boxplot}\hlstd{(ests2,}\hlkwc{ylim}\hlstd{=}\hlkwd{range}\hlstd{(}\hlopt{-}\hlnum{3.5}\hlstd{,}\hlnum{3.5}\hlstd{),}\hlkwc{pch}\hlstd{=}\hlnum{19}\hlstd{,}\hlkwc{cex}\hlstd{=}\hlnum{0.5}\hlstd{,}\hlkwc{names}\hlstd{=nvec)}
\hlkwd{title}\hlstd{(}\hlstr{'IPW and AIPW analyses: estimated PS model'}\hlstd{)}
\hlkwd{abline}\hlstd{(}\hlkwc{h}\hlstd{=}\hlnum{1}\hlstd{,}\hlkwc{lty}\hlstd{=}\hlnum{2}\hlstd{,}\hlkwc{col}\hlstd{=}\hlstr{'red'}\hlstd{)}
\end{alltt}
\end{kframe}

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

}


\begin{kframe}\begin{alltt}
\hlcom{#Bias}
\hlstd{b}\hlkwb{<-}\hlkwd{sqrt}\hlstd{(n)}\hlopt{*}\hlkwd{apply}\hlstd{(ests2}\hlopt{-}\hlnum{1}\hlstd{,}\hlnum{2}\hlstd{,mean)}
\hlcom{#Variance}
\hlstd{v}\hlkwb{<-}\hlstd{n}\hlopt{*}\hlkwd{apply}\hlstd{(ests2,}\hlnum{2}\hlstd{,var)}
\hlcom{#MSE}
\hlstd{mse}\hlkwb{<-}\hlstd{b}\hlopt{^}\hlnum{2}\hlopt{+}\hlstd{v}
\hlstd{res2}\hlkwb{<-}\hlkwd{rbind}\hlstd{(b,v,mse)}
\hlkwd{colnames}\hlstd{(res2)}\hlkwb{<-}\hlstd{nvec}
\hlkwd{rownames}\hlstd{(res2)}\hlkwb{<-}\hlkwd{c}\hlstd{(}\hlstr{'Bias'}\hlstd{,}\hlstr{'Var.'}\hlstd{,}\hlstr{'MSE'}\hlstd{)}
\hlkwd{round}\hlstd{(res2,}\hlnum{2}\hlstd{)}
\end{alltt}
\begin{verbatim}
+      M0    M1a   M1b   M2a   M2b    M3a   M3b   M4a   M4b   M5a   M5b
+ Bias  0  -0.41  0.15 -0.15 -0.14  -0.28  0.15  0.00  0.02  6.19  6.56
+ Var.  0 571.46 84.49 80.33 65.25 303.00 84.49 19.42 18.93 55.76 28.30
+ MSE   0 571.63 84.51 80.35 65.27 303.08 84.51 19.42 18.93 94.01 71.31
\end{verbatim}
\end{kframe}
\end{knitrout}