population.size = 1000000 population.mean = 300 population.std = 50 plot.delay = 0.0 maxsample.size = 100 sample.var = 0 population.data = rnorm(population.size,population.mean,population.std) population.mean = mean(population.data) par(mar=c(1,1,1,1)) # plot more samples for (plot.counter in 2:maxsample.size) { sample.data = sample(population.data,plot.counter,replace = TRUE) sample.var[plot.counter-1] = var(sample.data) } barplot(sample.var)