tests/testthat/test-get_gc3s.R

test_that("return correct results", {
    seqs <- Biostrings::DNAStringSet(c('ATGATT', 'ATC', 'ATA', 'GCGGCA'))
    cf <- count_codons(seqs)
    expect_equal(get_gc3s(cf), c(0, 1, 0, 0.5))
})

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.