R，在 R 中运行实验 100000 次

Question

作为实验，我在 r 中生成了以下值:

>  q<- matrix(rbinom(10*1000,1,0.5),ncol=10)
>  for( i in 1:10){
+  a=a+q[1,i]
+ }
> a<-a/10
> 
> k<-sample(1:1000,1)
> b<-0
> for(i in 1:10){
+ b=b+q[k,i]
+ }
> b<-b/10
> 
> c<-rowMeans(q)
> c<-min(q)

所以，我想要的是重复实验100000次。例如，生成矩阵 q 100000 次，与 a、b、c 相同。最后，我想要一个向量，它返回 100000 次迭代的结果 a 的平均值、b 的平均值和 c 的平均值。

添加信息: 我已经按照你的建议做了:但它给了我以下消息:

> means<-rowMeans(replicate(100000,wrapper))
Error in rowMeans(replicate(1e+05, wrapper)) : 
  'x' must be an array of at least two dimensions

上一版本:这只是一个 () 的问题，我的意思是，它的工作原理如下:

> means<-rowMeans(replicate(100000,wrapper()))

Answer 1

您正在寻找replicate函数:

# Wrap your code above into a function that returns a, b, and c
wrapper <- function() {
  # Your code goes here
  return(c(a,b,c))
}

means <- rowMeans(replicate(100000, wrapper()))
names(means) <- c("a", "b", "c")