1
    for (i in 1:10){

      z1<-rnorm(1000,0,1)
      z2<-rnorm(1000,0,1)
      z3<-rnorm(1000,0,1)
      z4<-rnorm(1000,0,1)
      z5<-rnorm(1000,0,1)
      z6<-rnorm(1000,0,1)
      z7<-rnorm(1000,0,1)

      X<-cbind(z1,z2,z3,z4,z5,z6)

      Z<-apply(X,2,mean)
       test[i]<-which(Z>0.01, arr.ind=TRUE)

    }

我有兴趣在这里将条件下的 Z 的索引保存在另一个变量中,以便在最后计算有多少索引已完全填充条件。

例如让第一个输出是 3 4 第二个输出 6 第三个 1 2 3 5

现在我想计算输出这 3 个输出出现了多少个、两个、三个等。

在所有这 3 个输出中,1 出现 1 次 2 出现 1 次 3 出现 2 次 4 出现 1 次 5 出现 1 次 6 出现 1 次。

提前致谢

4

1 回答 1

1

两种方式:

未矢量化

do.one <- function() 
  which(colMeans(replicate(7, rnorm(1000))) > .01)
replicate(10, do.one(), simplify=FALSE)

矢量化

do.one <- function() 
  colMeans(replicate(7, rnorm(1000))) > .01
replicate(10, do.one())
于 2013-05-28T22:05:57.660 回答