tests/testthat/test.calculate_max_correlation.r

context("calculate_max_correlation")

library(magrittr)
#library(testthat)

rm(list = ls())
set.seed(0)
numcol <- 1:10
nacol <- numcol
nacol[4:7] <- NA
rndcol <- rnorm(10)
faccol <- c("q", "w", "e", "r", "t", "y", "u", "i", "o", "p") %>% as.factor()

#### normal max ####

# self correlation
actual <-
  data.frame(
    a = numcol,
    b = numcol) %>%
  calculate_max_correlation(c("a", "b"))
expect_equal(actual, 1)

# self correlation (v2)
actual <-
  data.frame(
    a = numcol,
    b = faccol) %>%
  calculate_max_correlation(c("a", "b"))
expect_equal(actual, 0.2242, tolerance = 0.001)

# normal correlation
actual <-
  data.frame(
    a = numcol,
    b = rndcol) %>%
  calculate_max_correlation(c("a", "b"))
expect_equal(actual, 0.0633, tolerance = 0.001)

# pair-wise correlation
actual <-
  data.frame(
    a = numcol,
    b = rndcol,
    c = numcol) %>%
  calculate_max_correlation(c("a", "b", "c"))
expect_equal(actual, 1)

# column filtered correlation
actual <-
  data.frame(
    a = numcol,
    b = numcol,
    c = rndcol) %>%
  calculate_max_correlation(c("a", "c"))
expect_equal(actual, 0.0633, tolerance = 0.001)

# column filtered correlation
actual <-
  data.frame(
    a = numcol,
    b = nacol) %>%
  calculate_max_correlation(c("a", "b"))
expect_equal(actual, 1)
markanewman/mndredge documentation built on May 9, 2019, 5:52 a.m.