inst/unitTests/testing.R

test_findCDS <- function() {
  fcds <- findCDS("../extdata/unit_test.fasta", startCodon = c("ATG"), stopCodon = c("TAG", "TAA", "TGA"))
  checkTrue(all(as.data.frame(fcds) ==
      data.frame(seqnames = factor(c("ut_1", "ut_3"), levels = c("ut_1", "ut_2", "ut_3")), start = c(1,2), end = c(7,10), width = c(7,9), strand = factor(c("*", "*"), levels = c("+", "-", "*")), frame = c(0, 1))))
}

test_readRibodata <- function() {
  datadir <- system.file("extdata", package = "riboSeq")
  chlamyFasta <- paste(datadir, "/rsem_chlamy236_deNovo.transcripts.fa", sep = "")  
  fastaCDS <- findCDS(fastaFile = chlamyFasta,
                      startCodon = c("ATG"),
                      stopCodon = c("TAG", "TAA", "TGA"))
  rR <- readRibodata(paste(datadir, "/chlamy236_plus_deNovo_plusOnly_Index17_short", sep = ""), replicates = c(1))
  fCs <- frameCounting(rR, fastaCDS)
  checkTrue(all(frameShift(rC = fCs)[,2] == c(0,2,0,1)))
}

Try the riboSeqR package in your browser

Any scripts or data that you put into this service are public.

riboSeqR documentation built on Nov. 8, 2020, 8:23 p.m.