QQ plot blank in check model for glmmTMB with tweedie distribution

DESCRIPTION

@@ -1,7 +1,7 @@
 Type: Package
 Package: performance
 Title: Assessment of Regression Models Performance
-Version: 0.10.9.2
+Version: 0.10.9.3
 Authors@R: 
     c(person(given = "Daniel",
              family = "Lüdecke",

NAMESPACE

@@ -457,6 +457,7 @@ S3method(r2_mcfadden,poissonirr)
 S3method(r2_mcfadden,poissonmfx)
 S3method(r2_mcfadden,polr)
 S3method(r2_mcfadden,probitmfx)
+S3method(r2_mcfadden,serp)
 S3method(r2_mcfadden,truncreg)
 S3method(r2_mcfadden,vglm)
 S3method(r2_mckelvey,default)

NEWS.md

@@ -4,6 +4,10 @@
 
 * Rudimentary support for models of class `serp` from package *serp*.
 
+## General
+
+* Improved error messages for `check_model()` when QQ-plots cannot be created.
+
 ## Bug fixes
 
 * Fixed issue in `check_normality()` for t-tests.

R/check_model.R

@@ -14,7 +14,7 @@
 #' @param check Character vector, indicating which checks for should be performed
 #'   and plotted. May be one or more of `"all"`, `"vif"`, `"qq"`, `"normality"`,
 #'   `"linearity"`, `"ncv"`, `"homogeneity"`, `"outliers"`, `"reqq"`, `"pp_check"`,
-#'   `"binned_residuals"` or `"overdispersion"`, Not that not all check apply
+#'   `"binned_residuals"` or `"overdispersion"`. Note that not all check apply
 #'   to all type of models (see 'Details'). `"reqq"` is a QQ-plot for random
 #'   effects and only available for mixed models. `"ncv"` is an alias for
 #'   `"linearity"`, and checks for non-constant variance, i.e. for
@@ -49,7 +49,7 @@
 #'
 #' @details For Bayesian models from packages **rstanarm** or **brms**,
 #' models will be "converted" to their frequentist counterpart, using
 #' [`bayestestR::bayesian_as_frequentist`](https://easystats.github.io/bayestestR/reference/convert_bayesian_as_frequentist.html).
 #' A more advanced model-check for Bayesian models will be implemented at a
 #' later stage.
 #'
@@ -77,7 +77,7 @@
 #' plots are helpful to check model assumptions, they do not necessarily indicate
 #' so-called "lack of fit", e.g. missed non-linear relationships or interactions.
 #' Thus, it is always recommended to also look at
 #' [effect plots, including partial residuals](https://strengejacke.github.io/ggeffects/articles/introduction_partial_residuals.html).
 #'
 #' @section Homogeneity of Variance:
 #' This plot checks the assumption of equal variance (homoscedasticity). The

R/check_model_diagnostics.R

@@ -51,12 +51,23 @@
     }
   }
 
-  if (is.null(res_)) {
+  if (is.null(res_) || all(is.na(res_))) {
     if (verbose) {
+      if (is.null(model_info$family)) {
+        fam <- "model"
+      } else {
+        fam <- paste0("`", model_info$family, "`")
+      }
       insight::format_alert(
-        sprintf(
-          "QQ plot could not be created. Cannot extract residuals from objects of class `%s`.",
-          class(model)[1]
+        paste(
+          sprintf(
+            "QQ plot could not be created. Cannot extract residuals from objects of class `%s`.",
+            class(model)[1]
+          ),
+          sprintf(
+            "Maybe the model class or the %s family does not support the computation of (deviance) residuals?",
+            fam
+          )
         )
       )
     }

tests/testthat/test-check_model.R

@@ -51,3 +51,24 @@ test_that("`check_model()` works for quantreg", {
   x <- check_model(qm, verbose = FALSE)
   expect_s3_class(x, "check_model")
 })
+
+test_that("`check_model()` warnings for tweedie", {
+  skip_if_not_installed("glmmTMB")
+  skip_if_not_installed("lme4")
+  data(sleepstudy, package = "lme4")
+  set.seed(123)
+  d <- sleepstudy[sample.int(50), ]
+  m <- suppressWarnings(glmmTMB::glmmTMB(Reaction ~ Days,
+    data = d,
+    family = glmmTMB::tweedie
+  ))
+  expect_message(
+    expect_message(
+      expect_message(
+        check_model(m, iterations = 1, verbose = TRUE),
+        regex = "Not enough model terms"
+      ),
+      regex = "QQ plot could not"
+    )
+  )
+})