# takes a series of samples from a series of sample sizes and calculates the mean of the sample variability and the variance of the sample variability
population.size = 1000000
population.mean = 300
population.std = 50
plot.delay = 0.0
maxsample.size = 100
samplespersample.size = 50
meansample.var = 0
varsample.var = 0
sample.var = 0
population.data = rnorm(population.size,population.mean,population.std)
population.mean = mean(population.data)
par(mfrow=c(2, 1))
par(mar=c(1,1,1,1))
# plot more samples
for (plot.counter in 2:maxsample.size)
{
for (sample.counter in 1:samplespersample.size)
{
sample.data = sample(population.data,plot.counter,replace = TRUE)
sample.var[sample.counter] = var(sample.data)
}
meansample.var[plot.counter-1] = mean(sample.var)
varsample.var[plot.counter-1] = var(sample.var)
}
barplot(meansample.var)
title(main='Mean of Sample Variances')
barplot(varsample.var)
title(main='Variance of Sample Variances')