tests/testthat/test-get_cscg.R

In cubar: Codon Usage Bias Analysis

test_that("CSCg = 0 if all CSC are 0s", {
    ct <- get_codon_table()
    ct <- ct[ct$aa_code != "*", ]
    ct$csc <- rep(0, nrow(ct))
    cf <- matrix(sample.int(64), nrow = 1, ncol = 64)
    colnames(cf) <- sort(names(Biostrings::GENETIC_CODE))
    # works for a single sequence
    expect_equal(get_cscg(cf, ct), 0)
    # works for multiple sequences
    expect_equal(get_cscg(cf[c(1, 1), ], ct), c(0, 0))
})