做一些事情然后附加该值做一个空向量

我编写了一个函数，它从带有异常值的正态分布中抽取随机变量。但是，我现在想取这些数据的平均值并将其附加到一个空向量 1000 次（所以基本上是使用我的函数具有 1000 种不同平均值的向量）这是我到目前为止的尝试。

M = 1000
for (i in 1:M){
  means = numeric()
  m = mean(rnormout(100,1,10))
  means[i] = m[i]
}

但是，我的均值向量中似乎只有 1000 个 NA，有什么建议吗？谢谢！

有几种方法可以做到这一点。如果您想使用循环，诀窍是在启动循环之前定义一个空向量，并在其中存储值。这是一个使用 rnorm 从标准法线生成值的示例，因为我没有您的 rnormout 函数的代码：

M <- 1000
means <- vector()
for (i in 1:M){
  means[i] = mean(rnorm(100,1))
}

或者，最好根本不使用 for 循环。您可以使用 map_dbl 在一行中替换循环：

library(tidyverse)
M <- 1000
means <- map_dbl(1:M,~ mean(rnorm(100,1)))

我们可以在 sapply 中使用 base R（使用 R 4.1.0）

M <- 1000
sapply(seq_len(M),\(x) mean(rnorm(100,1)))