R中的多个模拟

Question

我正在尝试从R的指数分布中运行1000个模拟。我希望获得模拟列表及其均值向量。 但是，当我运行代码时，列表的前999个元素为NULL，因此前999个均值是NA。 我究竟做错了什么？

nosim <- 1000
n <- 40
r<- 0.2
means<- vector()
sims <- list()
set.seed(4993)
for(i in nosim){
  sims[[i]] = rexp(n, r)
}
means = sapply(sims,mean)

Answer 1

我认为问题在于您的变量nosim是标量而不是序列/数组。

您应该改用seq （或更快的seq.int ）（如seq(1,nosim) ）：

nosim <- 1000
n <- 40
r<- 0.2
means<- vector()   # not really necessary
sims <- list()
set.seed(4993)
for(i in seq(1,nosim)){    # or use `for(i in seq.int(nosim) ){`
  sims[[i]] = rexp(n, r)
}
means = sapply(sims,mean)

希望这可以帮助