tests/testthat/test-report.aov.R

In report: Automated Reporting of Results and Statistical Models

test_that("report.aov", {
  model <- anova(lm(Sepal.Width ~ Species, data = iris))
  r1 <- report(model)
  expect_equal(
    c(
      ncol(as.report_table(r1, summary = TRUE)),
      nrow(as.report_table(r1, summary = TRUE))
    ),
    c(7, 2),
    ignore_attr = FALSE
  )
  expect_equal(
    as.report_table(r1, summary = TRUE)$Mean_Square[1],
    5.6724,
    tolerance = 0.01
  )

  model <- aov(Sepal.Width ~ Species, data = iris)
  r2 <- report(model)
  expect_equal(
    c(
      ncol(as.report_table(r2, summary = TRUE)),
      nrow(as.report_table(r2, summary = TRUE))
    ),
    c(7, 2),
    ignore_attr = FALSE
  )
  expect_equal(
    as.report_table(r2, summary = TRUE)$Mean_Square[1],
    5.6724,
    tolerance = 0.01
  )

  r3 <- report(model)
  expect_equal(
    c(
      ncol(as.report_table(r3, summary = TRUE)),
      nrow(as.report_table(r3, summary = TRUE))
    ),
    c(7, 2),
    ignore_attr = FALSE
  )
  expect_equal(
    sum(as.report_table(r3, summary = TRUE)$Mean_Square),
    5.787854,
    tolerance = 0.01
  )

  data <- iris
  data$Cat1 <- rep_len(c("X", "X", "Y"), nrow(data))
  data$Cat2 <- rep_len(c("A", "B"), nrow(data))

  model <- aov(Sepal.Length ~ Species * Cat1 * Cat2, data = data)
  r4 <- report(model)
  expect_equal(
    c(
      ncol(as.report_table(r4, summary = TRUE)),
      nrow(as.report_table(r4, summary = TRUE))
    ),
    c(7, 8),
    ignore_attr = FALSE
  )
  expect_equal(
    as.report_table(r4, summary = TRUE)$Mean_Square[1],
    31.6060,
    tolerance = 0.01
  )

  model <- aov(Sepal.Length ~ Species * Cat1 + Error(Cat2), data = data)
  r5 <- report(model, verbose = FALSE)
  expect_warning(report(model), "non-finite")
  expect_equal(
    c(
      ncol(as.report_table(r5, summary = TRUE)),
      nrow(as.report_table(r5, summary = TRUE))
    ),
    c(8, 5),
    ignore_attr = FALSE
  )
  expect_equal(
    as.report_table(r5, summary = TRUE)$Mean_Square[1],
    31.60607,
    tolerance = 0.01
  )

  # snapshot tests -----

  set.seed(123)
  expect_snapshot(
    variant = "windows",
    suppressWarnings(report(anova(lm(Sepal.Width ~ Species, data = iris))))
  )

  set.seed(123)
  expect_snapshot(
    variant = "windows",
    suppressWarnings(report(anova(lm(
      wt ~ as.factor(am) * as.factor(cyl),
      data = mtcars
    ))))
  )

  set.seed(123)
  expect_snapshot(
    variant = "windows",
    suppressWarnings(report(aov(wt ~ cyl + Error(gear), data = mtcars)))
  )
})