test_that("mutation models are validated", {
alleles = 1:3
afr = c(.2,.3,.5)
rate = 0.1
M = mutationModel("prop", alleles=alleles, afreq=afr, rate=rate)
expect_silent(validateMutationModel(M))
})
test_that("male and female models are equal by default", {
alleles = 1:3
afr = c(.2,.3,.5)
rate = 0.1
m = mutationMatrix("prop", alleles=alleles, afreq=afr, rate=rate)
M = mutationModel("prop", alleles=alleles, afreq=afr, rate=rate)
expect_identical(M$female, m)
expect_identical(M$male, m)
})
test_that("stepwise mutation model splits correctly in male/female", {
alleles = c(1, 1.5, 2,3)
M = mutationModel("step", alleles = alleles,
rate = list(male=0, female=0.1),
rate2 = list(male=0.2, female=0),
range = list(male=2, female=3))
matM = mutationMatrix("step", alleles = alleles,
rate = 0, rate2 = 0.2, range = 2)
matF = mutationMatrix("step", alleles = alleles,
rate = 0.1, rate2 = 0, range = 3)
expect_identical(M, mutationModel(list(male=matM, female=matF)))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.