test_that("expit and logit are inverses", { expect_equal(logit(expit(c(0.3, 0.4))), c(0.3, 0.4)) }) test_that("real_to_simplex and simplex_to_real are inverses", { expect_equal(simplex_to_real(real_to_simplex(c(1, 3))), c(1, 3)) }) test_that("offspring_gf_2 and f1_gf_dr are same for tetraploids", { expect_equal( offspring_gf_2(alpha = 1/12, xi1 = 1/3, xi2 = 1/3, p1 = 2, p2 = 1), f1_gf_dr(alpha = 1/12, g1 = 2, g2 = 1, ploidy = 4) ) expect_equal( pvec_tet_2(alpha = 1/12, xi = 1/3, ell = 0), gamfreq_dr(alpha = 1/12, g = 0, ploidy = 4, log_p = FALSE) ) expect_equal( pvec_tet_2(alpha = 1/12, xi = 1/3, ell = 1), gamfreq_dr(alpha = 1/12, g = 1, ploidy = 4, log_p = FALSE) ) expect_equal( pvec_tet_2(alpha = 1/12, xi = 1/3, ell = 2), gamfreq_dr(alpha = 1/12, g = 2, ploidy = 4, log_p = FALSE) ) expect_equal( pvec_tet_2(alpha = 1/12, xi = 1/3, ell = 3), gamfreq_dr(alpha = 1/12, g = 3, ploidy = 4, log_p = FALSE) ) expect_equal( pvec_tet_2(alpha = 1/12, xi = 1/3, ell = 4), gamfreq_dr(alpha = 1/12, g = 4, ploidy = 4, log_p = FALSE) ) })