tests/testthat/test_get_bam_bed.R

In rujinwang/SCOPE: A normalization and copy number estimation method for single-cell DNA sequencing

context("getbambed")
library(SCOPE)
library(WGSmapp)

test_that("File preparation works", {
  bamfolder <- system.file("extdata", package = "WGSmapp")
  bamFile <- list.files(bamfolder, pattern = '*.dedup.bam$')
  bamdir <- file.path(bamfolder, bamFile)
  sampname_raw <- sapply(strsplit(bamFile, ".", fixed = TRUE), "[", 1)
  bambedObj <- get_bam_bed(bamdir = bamdir, sampname = sampname_raw)
  bamdir <- bambedObj$bamdir
  sampname_raw <- bambedObj$sampname
  ref_raw <- bambedObj$ref

  expect_equal(length(bamFile), length(sampname_raw))
  expect_gte(length(as.character(unique(seqnames(ref_raw)))), 22)
})