close easystats/see#103

NEWS.md

@@ -22,6 +22,8 @@
 * `test_performance()` now calls `test_lrt()` or `test_wald()` instead of
   `test_vuong()` when package *CompQuadForm* is missing.
 
+* `check_collinearity()` better supports GAM models.
+
 ## Bug fixes
 
 * `model_performance()` for models from *lme4* did not report AICc when

R/check_collinearity.R

@@ -280,7 +280,7 @@ check_collinearity.zerocount <- function(x,
   assign <- .term_assignments(x, component, verbose = verbose)
 
   # any assignment found?
-  if (is.null(assign)) {
+  if (is.null(assign) || all(is.na(assign))) {
     if (verbose) {
       warning(insight::format_message(sprintf("Could not extract model terms for the %s component of the model.", component), call. = FALSE))
     }
@@ -442,8 +442,14 @@ check_collinearity.zerocount <- function(x,
     }
   }))
 
+  if (insight::is_gam_model(x)) {
+    model_params <- as.vector(unlist(unlist(insight::find_parameters(x)[c(component, "smooth_terms")])))
+  } else {
+    model_params <- insight::find_parameters(x)[[component]]
+  }
+
   as.numeric(names(parms)[match(
-    insight::clean_names(insight::find_parameters(x)[[component]]),
+    insight::clean_names(model_params),
     parms
   )])
 }

R/check_model_diagnostics.R

@@ -18,7 +18,7 @@
 
 
 .diag_qq <- function(model) {
-  if (inherits(model, c("lme", "lmerMod", "merMod", "glmmTMB"))) {
+  if (inherits(model, c("lme", "lmerMod", "merMod", "glmmTMB", "gam"))) {
     res_ <- sort(stats::residuals(model), na.last = NA)
   } else if (inherits(model, "glm")) {
     res_ <- sort(stats::rstandard(model, type = "pearson"), na.last = NA)
@@ -213,6 +213,8 @@
     {
       if (inherits(model, "merMod")) {
         stats::residuals(model, scaled = TRUE)
+      } else if (inherits(model, "gam")) {
+        stats::residuals(model, type = "scaled.pearson")
       } else if (inherits(model, c("glmmTMB", "MixMod"))) {
         sigma <- if (faminfo$is_mixed) {
           sqrt(insight::get_variance_residual(model))