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

Question

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

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

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

Answer 1

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

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

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

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

Answer 2

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

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

Do something then append that value do an empty vector