sampling.dist.1<-NULL
for(sample.count in 1:10000){
  set.seed(sample.count)
  sample.var.1<-var(rnorm(2,2.5,1.118))
  sampling.dist.1<-c(sampling.dist.1, sample.var.1)
}
mean(sampling.dist.1)
var(sampling.dist.1)

hist(sampling.dist.1,  freq=F, col="grey", xlab="", xlim=c(-1, 6), breaks=100)
par(new=T)
plot(density(sampling.dist.1), axes=F, main="", xlim=c(-1, 6), lwd=2, col="blue")