tests/testthat/test-FindPeak_FindPeakGene.R

In liuwd15/tomoda: Tomo-seq data analysis

context("Test finding peak genes.")

library(SummarizedExperiment)
data(zh.data)
zh <- createTomo(zh.data)

test_that("find peak",
          {
            expect_equal(findPeak(c(0:5, 5:0), threshold=1, length=4), c(3,10))
            expect_equal(findPeak(c(0:2, 2:0), threshold=1, length=4), c(0,0))
          })

test_that("find peak genes",
          {
            peak_genes <- findPeakGene(zh, nperm=0)
            expect_equal(ncol(peak_genes), 4)
            peak_genes <- findPeakGene(zh, nperm=1e4)
            expect_equal(ncol(peak_genes), 6)
            peak_genes <- findPeakGene(zh, length=ncol(zh))
            expect_true(is.null(peak_genes))
          })