set.seed(123) m <- 2000 means <- 165 for (i in 1:m) {means[i]= mean(rnorm(i, 165, 5.5)) } plot(means, type = 'l', xlab = 'Numbers of samples', ylim = c(160,170)) abline(h=165, col='red')