R中的多个模拟

时间:2019-01-29 07:03:39

标签: r for-loop indexing seq

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

public function subCategories()
{
    return $this->hasMany(SubCategory::class, 'id', 'cat_id');
}

1 个答案:

答案 0 :(得分: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)

希望这会有所帮助