context("Feature Ignore Singletons")
test_that("we can add remove singletons to a models", {
model <- coal_model(10, 1) + feat_ignore_singletons()
expect_true(has_ign_singletons(model))
expect_false(has_ign_singletons(coal_model(10, 1)))
})
test_that("the function to remove singletons works", {
ssl <- list(create_test_segsites(), create_test_segsites())
expect_equal(remove_singletons(ssl),
list(create_test_segsites()[, -c(1, 3)],
create_test_segsites()[, -c(1, 3)]))
expect_equal(remove_singletons(list(create_empty_segsites())),
list(create_empty_segsites()))
})
test_that("singletons are removed from simulations", {
model <- coal_model(2, 1) +
feat_mutation(10) +
feat_ignore_singletons() +
sumstat_sfs()
expect_equal(simulate(model)$sfs, 0)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.