\begin{Schunk}
\begin{Sinput}
  grabTofB <- function(){
      x <- rpois(1000, lambda=17)
      y <- 2 + .223 * x + 15 * rnorm(1000)
      seB <- coef(summary(lm(y~x)))[2,3]
  }
  rep1000<- replicate(1000,grabTofB())
  m1000 <- round(mean(rep1000),2)
  sd1000 <- round(sd(rep1000),2)
  h1000 <- hist(rep1000, prob= TRUE, main="Distribution of T Statistics for Slope", xlab="T stat for slope estimate")
  legend("topright",legend=c(paste("mean=",m1000),paste("sd=",sd1000)))
\end{Sinput}
\end{Schunk}
