context("Optimization of Exact_cond_midP_linear_rank_tests_2xc()")
fill_nchoosek_old <- function(c0, npj) {
npj_choose_x1j <- matrix(0, c0, max(npj) + 1)
for (j in 1:c0) {
for (x1j in 0:npj[j]) {
npj_choose_x1j[j, x1j + 1] <- choose(npj[j], x1j)
}
}
return(npj_choose_x1j)
}
test_that("New fill_nchoosek produces the same output as the old one", {
for (iter in seq_len(10)) {
c0 <- sample(1:5, 1L)
npj <- rpois(c0, lambda = 10)
expect_equal(fill_nchoosek(c0, npj), fill_nchoosek_old(c0, npj))
}
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.