tests/testthat/test-get_gc.R

test_that("return correct results", {
  seqs <- Biostrings::DNAStringSet(c('ATGATT', 'ATC', 'ATA', 'GCG'))
  cf <- count_codons(seqs)
  expect_equal(get_gc(cf), c(1/6, 1/3, 0, 1))
})

Try the cubar package in your browser

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

cubar documentation built on April 3, 2025, 8:58 p.m.