r - R中按组划分的spearman相关性

Question

你如何在 R 中按组计算 Spearman 相关性。我发现以下链接按组讨论 Pearson 相关性。但是当我尝试用 spearman 替换类型时，它不起作用。

https://stats.stackexchange.com/questions/4040/r-compute-correlation-by-group

Answer 1

对于基本的 R 解决方案，这个怎么样:

df <- data.frame(group = rep(c("G1", "G2"), each = 10),
                 var1 = rnorm(20),
                 var2 = rnorm(20))

r <- by(df, df$group, FUN = function(X) cor(X$var1, X$var2, method = "spearman"))
# df$group: G1
# [1] 0.4060606
# ------------------------------------------------------------ 
# df$group: G2
# [1] 0.1272727

然后，如果您想要 data.frame 形式的结果:

data.frame(group = dimnames(r)[[1]], corr = as.vector(r))
#   group      corr
# 1    G1 0.4060606
# 2    G2 0.1272727

编辑 : 如果您喜欢 plyr -基于解决方案，这里是一个:

library(plyr)
ddply(df, .(group), summarise, "corr" = cor(var1, var2, method = "spearman"))