结果为R`cor（）`，但带有`Kendall's W`

Question

我有一个数据框，想为每对变量计算肯德尔的W。

此函数kendall(df, correct = TRUE)计算整个数据帧的W。我想要一个更像cor(df, method=c("kendall"))的表，该表比较每对变量，但是该函数使用Kendall的tau而不是Kendall的W。

我想进行以下计算：

install.packages("irr")
library(irr)
df<-iris
kendall(df[,1:4], correct=TRUE)

通过这种方式：

cor(df[,1:4])

Answer 1

我不确定我是否理解正确，下面的代码有帮助吗？

my.kendall <- function(df) {
  func<-Vectorize(function(i,j){kendall(df[,c(i,j)])$value})
  outer(X=1:ncol(df), Y=1:ncol(df), FUN=func)
}