set.seed(3) lambda <- 0.2 num_sim <- 1000 sample_size <- 40 sim <- matrix(rexp(num_sim*sample_size, rate=lambda), num_sim, sample_size) row_means <- rowMeans(sim)