Question

    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次。

提前致谢

Answer 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())

如何在矩阵中保存变量索引以计算出现次数

1 个答案: