# set sample size and number of samples set.seed(23456789) n <- 11 reps <- 1000 # perform random sampling samples <- replicate(reps, rnorm(n,10,2)) # 11 x 1000 sample matrix # compute sample variances sample_var <- rep(NA, reps) for (i in 1:reps) { sample_var[i] <- var(samples[,i]) } # check that 'sample_var' is a vector is.vector(sample_var) (mean(sample_var)) (var(sample_var)) hist(sample_var, breaks=40, prob=T, main=paste( "samples of size 11" ),col="steelblue") par(new=T) plot(density(sample_var), xlab="", axes=F, main="", col="black")