tests/testthat/test-gendist_distribution.R

In HopkinsIDD/phylosamp: Sample Size Calculations for Molecular and Phylogenetic Studies

test_that("gendist_distribution returns linked_prob==1 when max_dist 0 and mean_gens_pdf[1]>0", {

    expect_equal(gendist_distribution(mut_rate = 1, mean_gens_pdf = c(0.2, 0.2, 0.2,
        0.2, 0.2), max_dist = 0)[, "linked_prob"][[1]], 1)

    expect_equal(gendist_distribution(mut_rate = 2, mean_gens_pdf = c(1), max_dist = 0)[,
        "linked_prob"][[1]], 1)
})

test_that("gendist_distribution fails when parameters invalid", {

    expect_error(gendist_distribution(mut_rate = -1, mean_gens_pdf = c(0.5, 0.5),
        max_link_gens = 1))

    expect_error(gendist_distribution(mut_rate = 1, mean_gens_pdf = 0))

    expect_error(gendist_distribution(mut_rate = 1, mean_gens_pdf = c(0.2, -0.2,
        0.6)))

    expect_error(gendist_distribution(mut_rate = 1, mean_gens_pdf = c(0.2, -0.2,
        0.6), max_gens = 0))

    expect_error(gendist_distribution(mut_rate = 1, mean_gens_pdf = c(0.2, -0.2,
        0.6), max_dist = -1))

    expect_error(gendist_distribution(mut_rate = 1, mean_gens_pdf = c(0.2, -0.2,
        0.6), max_link_gens = 0))
})