test-residuals.R
In afex: Analysis of Factorial Experiments

context("residuals: check if it works")

data(obk.long, package = "afex")
between <- aov_car(value ~ treatment*gender + Error(id), data = obk.long, 
                   fun_aggregate = mean)
mixed <- aov_car(value ~ treatment * gender + Error(id/(phase*hour)), data = obk.long)
within <- aov_car(value ~ 1 + Error(id/(phase*hour)), data = obk.long)

## between data with correct order
obk2 <- aggregate(value ~ gender + treatment + id , data = obk.long, FUN = mean)
between2 <- aov_car(value ~ treatment*gender + Error(id), data = obk2)
between2lm <- lm(value ~ treatment*gender, data = obk2)

between3 <- aov_car(value ~ treatment*gender + Error(id), 
                    data = obk2[rev(seq_len(nrow(obk2))),])


## within data with correct order
obk3 <- obk.long[with(obk.long, order(id, phase, hour)), ]
within2 <- aov_car(value ~ 1 + Error(id/(phase*hour)), data = obk3)

test_that("Residuals", {
  expect_message(residuals(within), "Data was changed")
  expect_message(residuals(mixed), "Data was changed")
  expect_message(residuals(between), "Data was changed")
  expect_message(residuals(between3), "Data was changed")
  
  expect_message(residuals(between2), regexp = NA)
  expect_message(residuals(within2), regexp = NA)
  
  expect_equal(residuals(between2), residuals(between2lm))
  
  expect_is(suppressWarnings(residuals(within)),"numeric")
  expect_is(suppressWarnings(residuals(mixed)),"numeric")
  expect_is(suppressWarnings(residuals(between)),"numeric")
  
  expect_is(residuals(within, append = TRUE),"data.frame")
  expect_is(residuals(mixed, append = TRUE),"data.frame")
  expect_is(residuals(between, append = TRUE),"data.frame")
})


test_that("Fitted", {
  expect_message(fitted(within))
  expect_message(fitted(mixed))
  expect_message(fitted(between))
  
  expect_message(fitted(between2), regexp = NA)
  expect_message(fitted(within2), regexp = NA)
  
  expect_equal(fitted(between2), fitted(between2lm))
  
  expect_is(suppressWarnings(fitted(within)),"numeric")
  expect_is(suppressWarnings(fitted(mixed)),"numeric")
  expect_is(suppressWarnings(fitted(between)),"numeric")
  
  expect_is(fitted(within, append = TRUE),"data.frame")
  expect_is(fitted(mixed, append = TRUE),"data.frame")
  expect_is(fitted(between, append = TRUE),"data.frame")
})