# creates a population and then draws a series of samples of increasing size and plots the variance population.size = 1000000 population.mean = 300 population.std = 100 sampleSize = 1000 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:sampleSize) { sample.data = sample(population.data,plot.counter,replace = TRUE) sample.var[plot.counter-1] = var(sample.data) } barplot(sample.var)