tests/testthat/test-weighting.R

In brendanf/spoar: SIMD Partial Order Alignment in R

sequences <- c(
    "CATAAAAGAACGTAGGTCGCCCGTCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGAT",
    "ATAAAGGCAGTCGCTCTGTAAGCTGTCGATTCACCGGAAAGATGGCGTTACCACGTAAAGTGATAATGATTAT",
    "ATCAAAGAACGTGTAGCCTGTCCGTAATCTAGCGCATTTCACACGAGACCCGCGTAATGGG",
    "CGTAAATAGGTAATGATTATCATTACATATCACAACTAGGGCCGTATTAATCATGATATCATCA",
    "GTCGCTAGAGGCATCGTGAGTCGCTTCCGTACCGCAAGGATGACGAGTCACTTAAAGTGATAAT",
    "CCGTAACCTTCATCGGATCACCGGAAAGGACCCGTAAATAGACCTGATTATCATCTACAT"
)

weights <- seq(1, 21, 4)

repsequences <- rep(sequences, weights)

test_that("weighting works", {
    expect_equal(
        spoaConsensus(sequences, w = weights),
        spoaConsensus(repsequences)
    )
    expect_equal(
        rep(spoaAlign(sequences, w = weights), weights),
        spoaAlign(repsequences)
    )
})