\begin{Schunk}
\begin{Sinput}
 (trunc.mean <- mean(trunc.hat[2, ]))
\end{Sinput}
\begin{Soutput}
[1] 94.00859
\end{Soutput}
\begin{Sinput}
 ## Empirical standard error of the mean is:
 (trunc.mean.sd <- sd(trunc.hat[2, ]))
\end{Sinput}
\begin{Soutput}
[1] 0.7054973
\end{Soutput}
\begin{Sinput}
 ## Mean of within sample estimates of stardard error:
 (trunc.se.mean <- mean(trunc.hat[4, ]))
\end{Sinput}
\begin{Soutput}
[1] 0.7127951
\end{Soutput}
\begin{Sinput}
 ## Theory-based std.err based using parameters (ignoring trunctation)
 stderr.theory <- sigma/sqrt(N)
 stderr.theory
\end{Sinput}
\begin{Soutput}
[1] 1.224745
\end{Soutput}
\end{Schunk}
