Question

我有以下数据框：

set.seed(1)
y <- data.frame(a1 = rnorm(5) , b1 = rnorm(5), c1 = rnorm(5),  a2 = rnorm(5), b2 = rnorm(5), c2 = rnorm(5))

我想获得列对的相关性： cor（a1，a2），cor（b1，b2），cor（c1，c2）

我尝试了以下内容，但NA显示为输出：

apply(y,2,function(x) cor(x[1],x[3]))

我想得到相当于

的结果

cor(y[,1],y[,4])
cor(y[,2],y[,5])
cor(y[,3],y[,6])

在我的实际数据框中，我有更多的列对。

有什么想法吗？

感谢您的支持。

Answer 1

num.vars <- length(y)
var1 <- head(names(y), num.vars / 2)
var2 <- tail(names(y), num.vars / 2)
mapply(cor, y[var1], y[var2])
#         a1         b1         c1 
#  0.2491625 -0.5313192  0.5594564

Answer 2

在名称上使用变量正则表达式的另一种方法。如果变量名称是任意顺序，这也适用。

nn <- 
unique(sub('([0-9]+)','',names(y )))

sapply(nn,function(x){
    xy = y[,grep(x,names(y))]
    cor(xy[,1],xy[,2])})
         a          b          c 
-0.7615458  0.5683647  0.5594564

计算R中数据框中列对的相关性

2 个答案: