tests/testthat/test.get.scc.R

In hicrep: Measuring the reproducibility of Hi-C data

library(hicrep)
context("get.scc")

test_that("get.scc() returns a list of four elements", {
  
    data(HiCR1)
    data(HiCR2)
    processed <- prep(HiCR1, HiCR2, 1000000, 0, 5000000)
    scc.out = get.scc(processed, 1000000, 5000000)
  
    expect_is(scc.out, "list")
    expect_equal(length(scc.out), 4)
    expect_false(is.null(scc.out$corr))
    expect_false(is.null(scc.out$wei))
    expect_gt(scc.out$std, 0)
    expect_lt(scc.out$scc, 1)
    expect_gt(scc.out$scc, -1)
})