Nothing
# Setup -------------------------------------------------------------------
statistics <- list()
# stats() and stat() ------------------------------------------------------
lm1 <- lm(Fertility ~ ., data = swiss)
lm2 <- update(lm1, . ~ . - Examination)
BF10 <- 1 / exp((BIC(lm2) - BIC(lm1)) / 2)
BF_stats <- stats(
method = "BF BIC method",
statistics = c(
stat(name = "BF", value = BF10, subscript = "10"),
stat(name = "BF", value = 1 / BF10, subscript = "01")
)
)
statistics <- add_stats(
list = statistics,
output = BF_stats,
notes = "Wagenmakers (2007) method for calculating Bayes factors"
)
stats
# tidy_stats_to_data_frame() ----------------------------------------------
df <- tidy_stats_to_data_frame(statistics)
# write_stats() -----------------------------------------------------------
write_test_stats(statistics, "tests/data/stats.json")
# Cleanup -----------------------------------------------------------------
rm(df, statistics, BF10, BF_stats, lm1, lm2)
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.