tests/testthat/test-kendall.R

In YunlongJiao/kernrank: Kernel Methods for Rank Data

context("Kendall kernel")
library(kernrank)

n.rep <- 50
n.min <- 5
n.max <- 10

# Errors ------------------------------------------------------------------

n <- sample(n.min:n.max, 1)
x <- rnorm(n); x[1] <- NA
expect_error(kendall_total(x, x))

n <- sample(n.min:n.max, 1)
x <- rnorm(n); x <- as.matrix(x)
expect_error(kendall_total(x, x))

# cor vs cor.fk -----------------------------------------------------------

for (i in seq(n.rep)) {
  n <- sample(n.min:n.max, 1)
  x <- rnorm(n)
  y <- rnorm(n)
  expect_equal(cor(x, y, method = "kendall"), kendall_total(x, y))
  expect_equal(cor(x, y, method = "kendall"), pcaPP::cor.fk(x, y))
}

for (i in seq(n.rep)) {
  n <- max(sample(n.min:n.max, 1), 6)
  x <- rep(1:n.min, length.out=n)[sample(n)]
  y <- rep(1:n.min, length.out=n)[sample(n)]
  expect_equal(cor(x, y, method = "kendall"), kendall_total(x, y))
  expect_equal(cor(x, y, method = "kendall"), pcaPP::cor.fk(x, y))
}