计算 R 循环输出的平均值

Question

我需要计算 R for 循环 output 的均值。例如：

for (x in 1:3) {
  print(x)
  vector <- c(sample(1:2, 200, replace = TRUE))
  print(vector)
  vmean <- (mean(vector))
  print(vmean)
}
 

这给了我：

\[1\] 1

\[1\] 1.535

\[1\] 2

\[1\] 1.495

\[1\] 3

\[1\] 1.545

如何结合这三种方法来计算方法的平均值？

Answer 1

我们可以创建vmean初始化并连接（或创建特定长度的向量并分配）以创建means vector ，然后获取均值的means

vmean <- c()
for (x in 1:3) {
  print(x)
  vector <- c(sample(1:2, 200, replace = TRUE))
  #print(vector)
  vmean <- c(vmean, mean(vector))
  #print(vmean)
}

-输出

> vmean
[1] 1.520 1.495 1.460
> mean(vmean)
[1] 1.491667

---

或者使用带有colMeans的replicate

vmean <- colMeans(replicate(3, sample(1:2, 200, replace = TRUE)))
vmean
[1] 1.495 1.535 1.575
mean(vmean)

注意：值不同，因为没有设置种子