tests/testthat/test_kmerize.R

In coregenomics/kmap: Generate mappable regions of the genome for a given K-mer length

context("Chopping views into k-mers")

test_that("kmerize returns BSgenomeViews-class", {
    expect_is(kmerize(views, kmer), "BSgenomeViews")
})
test_that("kmerize returns ranges no wider than the kmer size", {
    kmer_ <- 10
    expect_lte(max(width(kmerize(views, kmer_))), kmer_)
    kmer_ <- 100
    expect_lte(max(length(width(kmerize(views, kmer_)))), kmer_)
    kmer_ <- 1
    expect_lte(max(width(kmerize(views, kmer_))), kmer_)
})
test_that("kmerize throws error for invalid kmer size", {
    expect_error(kmerize(views, -1), "kmer")
    expect_error(kmerize(views, 0), "kmer")
    expect_error(kmerize(views, 10.5), "kmer")
})