27

要找到两个矩阵 X 和 Y 的逐行相关性,输出应该具有 X 的第 1 行和 Y 的第 1 行的相关值,...,因此总共有 10 个值(因为有 10 行):

X <- matrix(rnorm(2000), nrow=10)
Y <- matrix(rnorm(2000), nrow=10)

sapply(1:10, function(row) cor(X[row,], Y[row,]))

现在,我应该如何将此函数应用于两个列表(每个列表包含大约 50 个数据帧)?

考虑列表 A 具有数据框 $1、$2、$3... 等等,列表 B 具有相似数量的数据框 $1、$2、$3。因此,该函数应应用于列表中的其他数据帧listA$1,listB$1listA$2,listB$2以此类推。listA$1最后,在比较 1 (和listB$1) 以及其他情况下,我将有十个值。

这可以使用“lapply”来完成吗?

4

1 回答 1

44

你似乎在寻找mapply. 这是一个例子:

listA <- list(matrix(rnorm(2000), nrow=10),
              matrix(rnorm(2000), nrow=10))
listB <- list(matrix(rnorm(2000), nrow=10),
              matrix(rnorm(2000), nrow=10))
mapply(function(X,Y) {
  sapply(1:10, function(row) cor(X[row,], Y[row,]))
  }, X=listA, Y=listB)
于 2013-09-25T10:52:05.340 回答