Nothing
# Setup -------------------------------------------------------------------
expected_statistics <- read_stats("../data/stats.json")
# stat() and stats() ------------------------------------------------------
test_that("BF stats works", {
statistics <- list()
lm1 <- lm(Fertility ~ ., data = swiss)
lm2 <- update(lm1, . ~ . - Examination)
BF10 <- 1 / exp((BIC(lm2) - BIC(lm1)) / 2)
BF_stats <- custom_stats(
method = "BF BIC method",
statistics = c(
custom_stat(name = "BF", value = BF10, subscript = "10"),
custom_stat(name = "BF", value = 1 / BF10, subscript = "01")
)
)
statistics <- add_stats(
list = statistics,
output = BF_stats,
notes = "Wagenmakers (2007) method for calculating Bayes factors"
)
statistics$BF_stats$statistics[[1]]
expected_statistics$BF_stats$statistics[[1]]
expect_equivalent(
statistics$BF_stats,
expected_statistics$BF_stats,
tolerance = 0.00001
)
})
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.