allow to directly pass models

R/check_residuals.R

@@ -29,16 +29,18 @@ check_residuals <- function(x, ...) {
   UseMethod("check_residuals")
 }
 
+#' @rdname check_residuals
 #' @export
-check_residuals.default <- function(x, ...) {
-  insight::format_error("`check_residuals()` only works with objects returned by `simulate_residuals()` or `DHARMa::simulateResiduals()`.") # nolint
+check_residuals.default <- function(x, alternative = c("two.sided", "less", "greater"), ...) {
+  if (insight::is_model(x)) {
+    check_residuals(simulate_residuals(x, ...), alternative = alternative)
+  } else {
+    insight::format_error("`check_residuals()` only works with objects supported by `simulate_residuals()` or `DHARMa::simulateResiduals()`.") # nolint
+  }
 }
 
-#' @rdname check_residuals
 #' @export
-check_residuals.performance_simres <- function(x,
-                                               alternative = c("two.sided", "less", "greater"),
-                                               ...) {
+check_residuals.performance_simres <- function(x, alternative = c("two.sided", "less", "greater"), ...) {
   alternative <- match.arg(alternative)
   ts <- suppressWarnings(
     stats::ks.test(

tests/testthat/test-check_residuals.R

@@ -11,4 +11,14 @@ test_that("check_singularity, lme4", {
     capture.output(print(out)),
     "Warning: Non-uniformity of simulated residuals detected (p = 0.019)."
   )
+  skip_if_not_installed("MASS")
+  set.seed(3)
+  mu <- rpois(500, lambda = 3)
+  x <- rnorm(500, mu, mu * 3)
+  x <- ceiling(x)
+  x <- pmax(x, 0)
+  quine.nb1 <- MASS::glm.nb(x ~ mu)
+  set.seed(123)
+  result <- check_residuals(quine.nb1)
+  expect_equal(result, 0.000665414, tolerance = 1e-3)
 })