test_that("No side effects when called empty", {
test_that("Side effect by show, width, or both", {
julia_progress(show = TRUE, width = "auto")

formula = weight ~ 1 + Diet + (1 + Diet | Chick),
data = subset(transform(ChickWeight, Diet = as.integer(Chick)), Chick %in% 1:5 & Time <= 20),

null_statistics <- permute_timewise_statistics(spec, nsim = 100)
out <- walk_threshold_steps(empirical_statistics, null_statistics, steps = 1:3)

empirical_ts <- compute_timewise_statistics(spec, statistic = "t")
empirical_chisqs <- compute_timewise_statistics(spec, statistic = "chisq")
test_that("tidy dims of empirical stats", {
null_ts <- permute_timewise_statistics(spec, statistic = "t")
null_chisqs <- permute_timewise_statistics(spec, statistic = "chisq")
expect_equal(dim(tidy(null_ts)), c(nrow(tidy(empirical_ts)) * 100, 4))
expect_equal(dim(tidy(null_chisqs)), c(nrow(tidy(empirical_chisqs)) * 100, 4))