R `cor ()` результаты в стиле, но с `Kendall's W` - PullRequest
0 голосов
/ 03 июля 2018

У меня есть датафрейм, и я хочу вычислить W Кендалла для каждой пары переменных.

Эта функция kendall(df, correct = TRUE) вычисляет W для всего кадра данных. Я бы хотел таблицу, похожую на cor(df, method=c("kendall")), которая сравнивает каждую пару переменных, однако эта функция использует тау Кендалла, а не W Кендалла.

Я хотел бы сделать этот расчет:

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

Таким образом:

cor(df[,1:4])

1 Ответ

0 голосов
/ 03 июля 2018

Я не уверен, правильно ли я понимаю, полезен ли приведенный ниже код?

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)
}
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...