Nothing
context("SumStat Pi")
test_that("calculation of nuc. div. is correct", {
seg_sites <- create_segsites(matrix(c(0, 0, 0, 1, 1,
0, 1, 1, 0, 0,
0, 1, 1, 1, 1), 5, 3), 1:3 / 5)
expect_equal(calc_nucleotide_div(list(seg_sites), 1:5), 1.6)
expect_equal(calc_nucleotide_div(list(seg_sites,
create_empty_segsites()), 1:5),
c(1.6, 0.0))
expect_equal(calc_nucleotide_div(list(seg_sites), 1:3), 4 / 3)
expect_equal(calc_nucleotide_div(list(seg_sites), c(1, 3, 5)), 2)
seg_sites <- list(create_segsites(matrix(c(0, 0, 0, 0,
0, 0, 1, 1,
0, 1, 1, 0,
1, 0, 1, 0), 4, 4, TRUE), 1:4 / 4))
expect_equal(calc_nucleotide_div(seg_sites, 1:4), 2)
})
test_that("nuc. div. statistik works", {
seg_sites <- create_segsites(matrix(c(0, 0, 0, 1, 1,
0, 1, 1, 0, 0,
0, 1, 1, 1, 1), 5, 3), 1:3 / 5)
model <- coal_model(5, 1)
pi <- sumstat_nucleotide_div("pi", 1)
expect_equal(pi$calculate(list(seg_sites), NULL, NULL, model), 1.6)
expect_error(sumstat_nucleotide_div("pi", 1:2))
expect_error(sumstat_nucleotide_div("pi", 2)$calculate(list(seg_sites),
NULL, NULL, model))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.