tests/testthat/test-dnastring2dist.R

In kullrich/MSA2dist: MSA2dist calculates pairwise distances between all sequences of a DNAStringSet or a AAStringSet using a custom score matrix and conducts codon based analysis

data(hiv)

test_that("compareCodons()", {
    myscore <- iupacMatrix()
    myscore[1,1] <- -1
    h <- hiv |> dnastring2dist(score=myscore)
    expect_true(h$sitesUsed[1,1] == 150)
    h <- hiv |> dnastring2dist(model="IUPAC")
    expect_true(h$sitesUsed[1,1] == 273)
    h <- hiv |> dnastring2dist(model="K80")
    expect_true(h$sitesUsed[1,1] == 273)
    mask1 <- IRanges::IRanges(start=c(1,61,121), end=c(30,90,150))
    h <- hiv |> dnastring2dist(model="IUPAC", mask=mask1)
    expect_true(h$sitesUsed[1,1] == 183)
    region1 <- IRanges::IRanges(start=c(1,139), end=c(75,225))
    h <- hiv |> dnastring2dist(model="IUPAC", region=region1)
    expect_true(h$sitesUsed[1,1] == 162)
    h <- hiv |> dnastring2dist(model="IUPAC", mask=mask1, region=region1)
    expect_true(h$sitesUsed[1,1] == 105)
    myscore[1,4] <- 0.5
    h <- hiv |> dnastring2dist(score=myscore)
    expect_true(h$distSTRING[1,2] == h$distSTRING[2,1])
    h <- hiv |> dnastring2dist(score=myscore, symmetric=FALSE)
    expect_true(h$distSTRING[1,2] != h$distSTRING[2,1])
})