Nothing
test_that("get_nll agrees with other package for truncated normal likelihood", {
cat("\n\n", as.character(Sys.time()), "--------------- test_theory.R ---------------\n\n")
ours <- get_nll(mu = 1.2, tau = 0.23, yi = 0.25, sei = 0.2, tcrit = 1.96)
check <- -log(truncnorm::dtruncnorm(x = 0.25,
mean = 1.2,
sd = sqrt(0.23 ^ 2 + 0.2 ^ 2),
a = -Inf,
b = 0.2 * 1.96))
expect_equal(ours, check)
})
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.