\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlcom{#Monte Carlo study}
\hlstd{n}\hlkwb{<-}\hlnum{1000}
\hlstd{nreps}\hlkwb{<-}\hlnum{200}
\hlstd{d}\hlkwb{<-}\hlnum{0.05}
\hlstd{zstab.samp}\hlkwb{<-}\hlkwd{matrix}\hlstd{(}\hlnum{0}\hlstd{,}\hlkwc{nrow}\hlstd{=nreps,}\hlkwc{ncol}\hlstd{=nZ)}
\hlkwa{for}\hlstd{(irep} \hlkwa{in} \hlnum{1}\hlopt{:}\hlstd{nreps)\{}
    \hlstd{X}\hlkwb{<-}\hlkwd{rmvn}\hlstd{(n,}\hlkwc{mu}\hlstd{=Mu,Sigma)}
    \hlstd{Xal}\hlkwb{<-}\hlkwd{cbind}\hlstd{(}\hlnum{1}\hlstd{,X)}
    \hlstd{ps.true}\hlkwb{<-}\hlstd{Xal} \hlopt{%*%} \hlstd{al}
    \hlstd{Z}\hlkwb{<-}\hlkwd{rnorm}\hlstd{(n,ps.true,sigZ)}
    \hlstd{Xm}\hlkwb{<-}\hlkwd{cbind}\hlstd{(}\hlnum{1}\hlstd{,Z,X)}
    \hlstd{muval}\hlkwb{<-}\hlkwd{as.vector}\hlstd{(Xm} \hlopt{%*%} \hlstd{be)}
    \hlstd{Y}\hlkwb{<-}\hlkwd{rnorm}\hlstd{(n,muval,sigY)}

    \hlstd{fnum}\hlkwb{<-}\hlkwd{lm}\hlstd{(Z}\hlopt{~}\hlnum{1}\hlstd{)}
    \hlstd{fden}\hlkwb{<-}\hlkwd{lm}\hlstd{(Z}\hlopt{~}\hlstd{X)}
    \hlstd{ratio}\hlkwb{<-}\hlkwd{dnorm}\hlstd{(Z,}\hlkwd{predict}\hlstd{(fnum),}\hlkwd{sd}\hlstd{(fnum}\hlopt{$}\hlstd{residuals))}\hlopt{/}\hlkwd{dnorm}\hlstd{(Z,}\hlkwd{predict}\hlstd{(fden),}\hlkwd{sd}\hlstd{(fden}\hlopt{$}\hlstd{residuals))}

    \hlstd{fitw}\hlkwb{<-}\hlkwd{lm}\hlstd{(Y}\hlopt{~}\hlstd{Z,}\hlkwc{weight}\hlstd{=ratio)}
    \hlstd{zstab.samp[irep,]}\hlkwb{<-}\hlkwd{predict}\hlstd{(fitw,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(}\hlkwc{Z}\hlstd{=zgrid))}
\hlstd{\}}
\end{alltt}
\end{kframe}
\end{knitrout}