问题描述
我编写了一个函数,它从带有异常值的正态分布中抽取随机变量。但是,我现在想取这些数据的平均值并将其附加到一个空向量 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)))