\begin{Schunk}
\begin{Sinput}
 cof <- sapply(reshet, coefficients)
 cofhomo <- sapply(reshomo, coefficients)
 par(mfrow=c(1,2))
 hhet <- hist(cof[2,], prob=T, main="Heteroskedastic Data", xlab="Slope Est.")
 lines(density(cof[2,]))
 hist(cofhomo[2,], prob=T, xlim=range(hhet$breaks), main="Homoskedastic Data", xlab="Slope Est.")
 lines(density(cofhomo[2,]))
\end{Sinput}
\end{Schunk}
