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')