2个计数变量

Question

lapply功能我R中很环保，所以有可能是一个非常简单的解决办法是：

我要计算一个方阵的列向量之间的平均相关性：

x<-matrix(rnorm(10000),ncol=100) 
aux<-matrix(seq(1,10000)) 
loop<-sapply(aux,function(i,j) cov(x[,i],x[,j]) 
cor_x<-mean(loop) 

当评估sapply行时，我得到错误'下标越界'。我知道我可以通过脚本来做到这一点，但有什么办法可以在一行代码中实现这一点？

Answer 1

问题归因于aux。由于您有100列，因此变量aux必须从1到100。但是您的aux是沿行的x的序列，因此范围从1到10000。它将与下面的代码工作：

aux <- seq(1, 100) 
loop <- sapply(aux, function(i, j) cov(x[, i], x[, j])) 

之后，您可以计算平均协方差：

cor_x <- mean(loop) 

如果你想排除重复字段（例如，COV（X，Y）本质上是相同的以COV（Y，X）），你可以使用：

cor_x <- mean(loop[upper.tri(loop, diag = TRUE)]) 

如果你也想排除COV（X，X），即方差，你可以使用：

cor_x <- mean(loop[upper.tri(loop)]) 

Answer 2

不需要任何循环。只需使用mean(cov(x))，这非常有效。