tests/testthat/testWindowVec.R

In bhklab/CREAM: Clustering of Genomic Regions Analysis Method

library(CREAM)

context("WindowVec")

test_that("Specify window size for each order of COREs", {
    InputData   <- read.table("A549_Chr21.bed", sep="\t")
    colnames(InputData) <- c("chr", "start", "end")
    MinLength <- 1000
    if(nrow(InputData) < MinLength){
        stop(paste( "Number of functional regions is less than ", MinLength, ".", sep = "", collapse = ""))
    }
    peakNumMin <- 2
    WScutoff <- 1.5
    WindowVecFinal <- WindowVec(InputData, peakNumMin, WScutoff)
    expect_equal_to_reference(WindowVecFinal,
                "windowVecFinal.rds")
})