\begin{Schunk}
\begin{Sinput}
 rmse <- vector(nsims, mode="numeric")
 for (i in 1:nsims){ rmse[i] <- res[[i]]$sigma}
 hist(rmse, breaks=20, probability=T, main="", xlab="RMSE (est. std. dev. of error term)")
 lines(density(rmse), lty=4)
 legend("topright",legend=c(paste("mean=",round(mean(rmse),3)),paste("std.dev.=",round(sd(rmse),3)) ))
\end{Sinput}
\end{Schunk}
