\begin{Schunk}
\begin{Sinput}
 tb1 <- vector(nsims, mode="numeric")
 for (i in 1:nsims){ tb1[i] <- coefficients(res[[i]])[2,3]}
 hist(tb1, breaks=20, probability=T, main="", xlab="1000 t ratios", xlim=c(-2,7))
 nnotsig <- length(which(tb1< 1.983))  ##qt(0.025, df=100)))
 lines(density(tb1), lty=4)
 rtb1 <- range(tb1)
 xseq <- seq(from=rtb1[1],to=rtb1[2],length.out=200)
 legend("topright", legend=c("Simulated t stats", paste("mean=",round(mean(tb1),3)),paste("std.dev.=",round(sd(tb1),3)) ))
 text(-1.4, 0.13, pos=4, labels=c(paste("    ",nnotsig,"\nsample estimates\n not significant \n at 0.05")))
 rtb1 <- range(tb1)
 xseq <- seq(from=-2,to=rtb1[2],length.out=200)
 lines(xseq, dt(xseq, df=98))
 text(-1.2, 0.3, "t dist")
\end{Sinput}
\end{Schunk}
