test_that("lumping works for custum matrix", {
m1 = mutationMatrix("custom", matrix = matrix(1/3, ncol=3, nrow=3), alleles=1:3)
lumped1 = lumpedMatrix(m1, lump = 2:3)
expect_equal(dim(lumped1), c(2,2))
expect_equal(colnames(lumped1), c("1", "lump"))
expect_equal(attr(lumped1, 'lumpedAlleles'), c("2", "3"))
})
test_that("wrong input is caught when lumping matrix", {
mat = mutationMatrix("eq", alleles = 1:4, rate = 0.1)
expect_error(lumpedMatrix(mat, lump = 2:3, afreq = 1),
"Frequency vector does not match the number of alleles")
expect_error(lumpedMatrix(mat, lump = 2:3, afreq = 1:4),
"Allele frequencies do not sum to 1")
mod = mutationModel("eq", alleles = 1:4, rate = 0.1)
expect_error(lumpedModel(mod, lump = 2:3, afreq = 1),
"Frequency vector does not match the number of alleles")
expect_error(lumpedModel(mod, lump = 2:3, afreq = 1:4),
"Allele frequencies do not sum to 1")
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.