tests/testthat/test_residuals.R

In inbo/inlatools: Diagnostic Tools for INLA Models

context("distribution_check")
test_that("checks the model properties", {
  ds <- generate_data()
  model <- INLA::inla(
    poisson ~ f(group_id, model = "iid"),
    family = "poisson",
    data = ds,
    control.predictor = list(compute = TRUE)
  )
  expect_error(
    residuals(model, type = "deviance"),
    "only Pearson residuals are available"
  )
  expect_is(resid <- residuals(model), "numeric")
  expect_identical(length(resid), nrow(ds))

  model <- INLA::inla(
    poisson ~ f(group_id, model = "iid"),
    family = "gaussian",
    data = ds,
    control.predictor = list(compute = TRUE)
  )
  expect_error(
    residuals(model),
    "`gaussian` distribution not .* handled by `residuals"
  )

  selected <- sample(nrow(ds), ceiling(0.5 * nrow(ds)))
  ds$poisson[selected] <- NA
  ds$zipoisson[-selected] <- NA
  model <- INLA::inla(
    cbind(poisson, zipoisson) ~ f(group_id, model = "iid"),
    family = c("poisson", "poisson"),
    data = ds
  )
  expect_error(residuals(model), "Only single responses are handled")
})