mle.stat <- NULL
postscript("MLEdist.eps")
par(mfrow=c(2,2))
for(i in 1:1000) {
  x <- rnorm(10)
  mle.stat <- cbind(mle.stat,mean(x)*as.numeric(mean(x)>0))
}
hist(mle.stat,xlab="MLE",main="n=10")
abline(v=mean(mle.stat),col=2,lwd=2)
for(i in 1:1000) {
  x <- rnorm(100)
  mle.stat <- cbind(mle.stat,mean(x)*as.numeric(mean(x)>0))
}
hist(mle.stat,xlab="MLE",main="n=100")
abline(v=mean(mle.stat),col=2,lwd=2)
for(i in 1:1000) {
  x <- rnorm(500)
  mle.stat <- cbind(mle.stat,mean(x)*as.numeric(mean(x)>0))
}
hist(mle.stat,xlab="MLE",main="n=500")
abline(v=mean(mle.stat),col=2,lwd=2)
for(i in 1:1000) {
  x <- rnorm(1000)
  mle.stat <- cbind(mle.stat,mean(x)*as.numeric(mean(x)>0))
}
hist(mle.stat,xlab="MLE",main="n=1000")
abline(v=mean(mle.stat),col=2,lwd=2)
dev.off()