tests/testthat/test-sim_unlinked_explicit.R

context("sim_unlinked_explicit")

test_that("sim_unlinked_explicit", {
  phy  <- ape::read.tree(text = "(t1:10,(t3:2,t2:2):8);")

  sequences <- nodeSub::sim_unlinked(phy, node_time = 0.1,
                                     rate1 = 0.01,
                                     rate2 = 0.01)
  testthat::expect_true(class(sequences$alignment) == "phyDat")

  sequences_2 <- nodeSub::sim_unlinked_explicit(phy = phy,
                                                rootseq = sequences$root_seq,
                                                node_time = 0.1,
                                                rate1 = 0.01,
                                                rate2 = 0.01)

  d1 <- calc_dist(sequences$alignment, root_sequence  = sequences$root_seq)
  d2 <- calc_dist(sequences_2$alignment, root_sequence  = sequences$root_seq)

  testthat::expect_equal(sum(d1), sum(d2), tolerance = 0.1)
})

Try the nodeSub package in your browser

Any scripts or data that you put into this service are public.

nodeSub documentation built on Aug. 19, 2022, 5:23 p.m.