Nothing
set.seed(342)
a <- rexp(10, .1)
M <- 3000
test_that("Dirichlet random generation and fitting works", {
expect_error(rdirichlet(-10, a))
expect_error(rdirichlet(10, -a))
expect_silent(rd <- rdirichlet(M, a))
expect_is(rd, "matrix")
expect_true(all(dim(rd) == c(M, length(a))))
expect_true(all(rd >= 0))
expect_silent(est <- fit_dirichlet(rd, abstol = .001))
expect_length(est$alpha, length(a))
expect_equal(est$alpha, a, tol = .1)
expect_equal(est$sum, sum(a), tol = .1)
})
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.