tests/testthat/test_matches_old.R

In jpalowitch/cbce: Correlation Bi-Community Extraction method

source("sim_eQTL_network.R")
source("cbce_old.R")
library(rlist)
library(pipeR)

context("Test ChiSq")

set.seed(123345)

comms <- function(exres, dx) {
  exres %>>% 
    list.filter(exists("StableComm") && min(sum(initial_set <= dx), sum(initial_set > dx)) > 1) %>>%
    list.map(StableComm)
}
areIdentical <- function(res1, res2, dx) {
 #Check extraction orders are the same
 ord1 <- res1$report[[4]]
 ord2 <- res2$report[[4]]
 
 identical(ord1, ord2) && 
 {
   s1 <- comms(res1$extract_res[ord1], dx)
   s2 <- comms(res2$extract_res[ord2], dx)
   all(list.zip(s1 = s1,s2 = s2) %>>% list.mapv(setequal(s1, s2)))
 }
}

test_that("ChiSq gives the same extraction", {
out <- 
capture.output({  
  sim <- sim_eQTL_network(make_param_list(cmin=3, cmax=5, b=1))
  res1 <- cbce(sim$X, sim$Y)
  res2 <- cbceNW_c(sim$X, sim$Y, twoSided = TRUE)
  expect_true(areIdentical(res1, res2, sim$dx))
  
  sim <- sim_eQTL_network(make_param_list(cmin=3, cmax=5, b=1))
  set.seed(12345)
  res1 <- cbce(sim$X, sim$Y, exhaustive = TRUE)
  set.seed(12345)
  res2 <- cbceNW_c(sim$X, sim$Y, twoSided = TRUE, exhaustive=TRUE)
  expect_true(areIdentical(res1, res2, sim$dx))
  
  sim <- sim_eQTL_network(make_param_list(cmin=3, cmax=5, b=1))
  res1 <- cbce(sim$X, sim$Y, start_nodes = 1:sim$dx)
  res2 <- cbceNW_c(sim$X, sim$Y, twoSided = TRUE, start_nodes = 1:sim$dx)
  expect_true(areIdentical(res1, res2, sim$dx))
})
})