计算两个列表之间的皮尔逊系数

Question

我有两个不同的列表（a和b），其中包含626257个向量，每个向量包含44个数字条目。一个列表包含样本数据，另一个列表用作参考。现在，我要计算两个列表的所有条目之间的皮尔逊相关性。我将值存储在变量（r）中。

不幸的是，“ r”的输出确实仅包含“ NA”条目。

这里的代码生成两个伪列表。

a = replicate(626257,rep(10,44),simplify = FALSE)
b = replicate(626257,rep(3,44),simplify = FALSE)

这里是计算相关性的代码。

r = lapply(seq_along(a), function(ind)cor(a[[ind]], b[[ind]]))
View(r)

Answer 1

您可以使用mapply轻松完成此任务。正如AkselA指出的那样，您首先需要模拟一些差异的数据（例如，使用rnorm(44)而不是rep, 10, 44)）。见下文：

a <- replicate(626257,rnorm(44),simplify = FALSE)
b <- replicate(626257,rnorm(44),simplify = FALSE)
r <- mapply(cor, a, b)