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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.