f1 <-
function() {
nSamples = 1000
nReplicates = 100
data = rpois(nSamples, lambda = 20)
for (i in 2:nSamples) {
data = rbind(data, rpois(nSamples, lambda = 20))
}
apply(data, 1, mean)