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