test_that("class methods", {
lik <- likelihood("normal", 0, 1)
expect_equal(names(lik), c("family", "parameters", "likelihood_function"))
pri <- prior("normal", 0, 1)
expect_equal(names(pri), c("family", "parameters", "prior_function"))
pred <- lik * pri
expect_equal(
names(pred),
c(
"integral",
"marginal_function",
"evidence_function",
"posterior_function",
"conditional_function",
"weighted_likelihood_function",
"prediction_function"
)
)
# normal distribution
lik <- likelihood("normal", 0, 1)
pri <- prior("normal", 0, 1)
prod1 <- lik * pri
prod2 <- pri * lik
expect_equal(prod1$integral, prod2$integral)
pr1 <- prior("cauchy", 0, 1, c(-Inf, Inf))
pr2 <- make_prior(new("cauchy"), 0, 1, c(-Inf, Inf))
expect_equal(pr1, pr2)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.