mattansb · Oct 20, 2026 · Oct 13, 2026 · Oct 13, 2026 · Oct 13, 2026 · Oct 13, 2026
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,7 +1,7 @@
 Type: Package
 Package: parameters
 Title: Processing of Model Parameters
-Version: 0.30.0.10
+Version: 0.30.0.11
 Authors@R:
     c(person(given = "Daniel",
              family = "Lüdecke",
@@ -241,4 +241,4 @@ Config/testthat/parallel: false
 Config/Needs/website: easystats/easystatstemplate
 Config/Needs/check: stan-dev/cmdstanr
 Config/rcmdcheck/ignore-inconsequential-notes: true
-Config/roxygen2/version: 8.1.0
+Config/roxygen2/version: 8.1.1
diff --git a/NAMESPACE b/NAMESPACE
@@ -59,6 +59,7 @@ S3method(ci,list)
 S3method(ci,lm)
 S3method(ci,lm_robust)
 S3method(ci,lme)
+S3method(ci,lmerMod)
 S3method(ci,lmodel2)
 S3method(ci,logistf)
 S3method(ci,logitmfx)
@@ -289,6 +290,7 @@ S3method(model_parameters,lcmm)
 S3method(model_parameters,list)
 S3method(model_parameters,lm_robust)
 S3method(model_parameters,lme)
+S3method(model_parameters,lmerMod)
 S3method(model_parameters,lmodel2)
 S3method(model_parameters,logistf)
 S3method(model_parameters,logitmfx)
@@ -501,6 +503,7 @@ S3method(p_value,list)
 S3method(p_value,lm)
 S3method(p_value,lm_robust)
 S3method(p_value,lme)
+S3method(p_value,lmerMod)
 S3method(p_value,lmodel2)
 S3method(p_value,logistf)
 S3method(p_value,logitmfx)
@@ -852,6 +855,7 @@ S3method(standard_error,lcmm)
 S3method(standard_error,list)
 S3method(standard_error,lm_robust)
 S3method(standard_error,lme)
+S3method(standard_error,lmerMod)
 S3method(standard_error,lmodel2)
 S3method(standard_error,logistf)
 S3method(standard_error,logitmfx)

diff --git a/NEWS.md b/NEWS.md
@@ -1,5 +1,14 @@
 # parameters (devel)
 
+## Breaking Changes
+
+* `compare_parameters()`'s `ci_method` now defaults to `NULL`, so the value is
+  potentially set by the models' unique methods.
+
+* `model_parameters()`, `ci()`, `standard_error()`, and `p_value()` now default
+  to the Satterthwaite approximation for `lmerModLmerTest` models and suggest
+  doing so for `lmerMod` models (#989).
+
 ## Changes
 
 * `model_parameters()` for models from *blavaan* with `standardize` no longer

diff --git a/R/1_model_parameters.R b/R/1_model_parameters.R
@@ -196,6 +196,9 @@
 #' - Applies to *linear mixed models*. CIs computed using the
 #'   Wald method (SE and a *t-distribution with Satterthwaite df*); p-values
 #'   computed using the Wald method with a *t-distribution with Satterthwaite df*.
+#'   This is the default for `lmerModLmerTest` models when the corresponding
+#'   method argument (`ci_method` or `method`) is `NULL` (and
+#'   `bootstrap = FALSE`).
 #'
 #' `"kenward"`
 #' - Applies to *linear mixed models*. CIs computed using the Wald

diff --git a/R/ci_generic.R b/R/ci_generic.R
@@ -35,13 +35,13 @@
    )
  )

  if (method == "ml1") {
    # nolint
    return(ci_ml1(model, ci = ci))
   } else if (method == "betwithin") {
     return(ci_betwithin(model, ci = ci))
   } else if (method == "satterthwaite") {
-    return(ci_satterthwaite(model, ci = ci))
+    return(ci_satterthwaite(model, ci = ci, vcov = vcov, vcov_args = vcov_args, ...))
   } else if (method %in% c("kenward", "kr")) {
     return(ci_kenward(model, ci = ci))
   }

diff --git a/R/compare_parameters.R b/R/compare_parameters.R
@@ -69,23 +69,25 @@
 #' compare_parameters(list("linear model" = m1, "logistic reg." = m2))
 #' }
 #' @export
-compare_parameters <- function(...,
-                               ci = 0.95,
-                               effects = "fixed",
-                               component = "conditional",
-                               standardize = NULL,
-                               exponentiate = FALSE,
-                               ci_method = "wald",
-                               p_adjust = NULL,
-                               select = NULL,
-                               column_names = NULL,
-                               pretty_names = TRUE,
-                               coefficient_names = NULL,
-                               keep = NULL,
-                               drop = NULL,
-                               include_reference = FALSE,
-                               groups = NULL,
-                               verbose = TRUE) {
+compare_parameters <- function(
+  ...,
+  ci = 0.95,
+  effects = "fixed",
+  component = "conditional",
+  standardize = NULL,
+  exponentiate = FALSE,
+  ci_method = NULL,
+  p_adjust = NULL,
+  select = NULL,
+  column_names = NULL,
+  pretty_names = TRUE,
+  coefficient_names = NULL,
+  keep = NULL,
+  drop = NULL,
+  include_reference = FALSE,
+  groups = NULL,
+  verbose = TRUE
+) {
   models <- list(...)
 
   if (length(models) == 1) {
@@ -117,14 +119,23 @@ compare_parameters <- function(...,
     }
   }
 
-  supported_models <- vapply(models, function(i) {
-    insight::is_model_supported(i) || inherits(i, "lavaan") || inherits(i, "parameters_model")
-  }, TRUE)
+  supported_models <- vapply(
+    models,
+    function(i) {
+      insight::is_model_supported(i) ||
+        inherits(i, "lavaan") ||
+        inherits(i, "parameters_model")
+    },
+    TRUE
+  )
 
   if (!all(supported_models)) {
     if (verbose) {
       insight::format_alert(
-        sprintf("Following objects are not supported: %s", toString(model_names[!supported_models])),
+        sprintf(
+          "Following objects are not supported: %s",
+          toString(model_names[!supported_models])
+        ),
         "Dropping unsupported models now."
       )
     }
@@ -169,7 +180,15 @@ compare_parameters <- function(...,
       dat <- model
     } else {
       # set default-ci_type for Bayesian models
-      if (insight::is_bayesian_model(model, exclude = c("bmerMod", "bayesx", "blmerMod", "bglmerMod")) && !ci_method %in% c("hdi", "quantile", "ci", "eti", "si", "bci", "bcai")) { # nolint
+      if (
+        insight::is_bayesian_model(
+          model,
+          exclude = c("bmerMod", "bayesx", "blmerMod", "bglmerMod")
+        ) &&
+          (is.null(ci_method) ||
+            !ci_method %in% c("hdi", "quantile", "ci", "eti", "si", "bci", "bcai"))
+      ) {
+        # nolint
         ci_method_tmp <- "eti"
       } else {
         ci_method_tmp <- ci_method
@@ -217,7 +236,10 @@ compare_parameters <- function(...,
 
     # add zi-suffix to parameter names
     if (any(dat$Component == "zero_inflated")) {
-      dat$Parameter[dat$Component == "zero_inflated"] <- paste0(dat$Parameter[dat$Component == "zero_inflated"], " (zi)")
+      dat$Parameter[dat$Component == "zero_inflated"] <- paste0(
+        dat$Parameter[dat$Component == "zero_inflated"],
+        " (zi)"
+      )
     }
 
     # add suffix
@@ -235,15 +257,28 @@ compare_parameters <- function(...,
   names(object_attributes) <- model_names
 
   # merge all data frames
-  all_models <- suppressWarnings(Reduce(function(x, y) {
-    merge(x, y, all = TRUE, sort = FALSE, by = c("Parameter", "Component", "Effects", "Group"))
-  }, m))
+  all_models <- suppressWarnings(Reduce(
+    function(x, y) {
+      merge(
+        x,
+        y,
+        all = TRUE,
+        sort = FALSE,
+        by = c("Parameter", "Component", "Effects", "Group")
+      )
+    },
+    m
+  ))
 
   # find columns with model numbers and create new variable "params_order",
   # which is pasted together of all model-column indices. Take lowest index of
   # all model-column indices, which then indicates order of parameters/rows.
   model_cols <- which(startsWith(colnames(all_models), "model"))
-  params_order <- as.numeric(substr(gsub("NA", "", do.call(paste0, all_models[model_cols]), fixed = TRUE), 0, 1))
+  params_order <- as.numeric(substr(
+    gsub("NA", "", do.call(paste0, all_models[model_cols]), fixed = TRUE),
+    0,
+    1
+  ))
 
   all_models <- all_models[order(params_order), ]
   all_models[model_cols] <- NULL
@@ -253,11 +288,20 @@ compare_parameters <- function(...,
     all_models$Group <- NULL
   }
 
-  attr(all_models, "model_names") <- gsub("\"", "", unlist(lapply(model_names, insight::safe_deparse)), fixed = TRUE)
+  attr(all_models, "model_names") <- gsub(
+    "\"",
+    "",
+    unlist(lapply(model_names, insight::safe_deparse)),
+    fixed = TRUE
+  )
   attr(all_models, "output_style") <- select
   attr(all_models, "all_attributes") <- object_attributes
   attr(all_models, "parameter_groups") <- groups
-  class(all_models) <- c("compare_parameters", "see_compare_parameters", unique(class(all_models)))
+  class(all_models) <- c(
+    "compare_parameters",
+    "see_compare_parameters",
+    unique(class(all_models))
+  )
 
   all_models
 }
@@ -270,11 +314,12 @@ compare_models <- compare_parameters
 
 # helper ----------------------------
 
-
 .set_pretty_names <- function(x, pretty_names) {
   # check if pretty names should be replaced by value labels
   # (if we have labelled data)
-  if (isTRUE(getOption("parameters_labels", FALSE)) || identical(pretty_names, "labels")) {
+  if (
+    isTRUE(getOption("parameters_labels", FALSE)) || identical(pretty_names, "labels")
+  ) {
     attr(x, "pretty_names") <- attr(x, "pretty_labels", exact = TRUE)
     pretty_names <- TRUE
   }
@@ -288,15 +333,19 @@ compare_models <- compare_parameters
       match_pretty_names <- match(names(att$pretty_names), x$Parameter)
       match_pretty_names <- match_pretty_names[!is.na(match_pretty_names)]
       if (length(match_pretty_names)) {
-        x$Parameter[match_pretty_names] <- att$pretty_names[x$Parameter[match_pretty_names]]
+        x$Parameter[match_pretty_names] <- att$pretty_names[x$Parameter[
+          match_pretty_names
+        ]]
       }
     } else {
       match_pretty_names <- att$pretty_names[x$Parameter]
       if (anyNA(match_pretty_names)) {
         match_pretty_names <- match(names(att$pretty_names), x$Parameter)
         match_pretty_names <- match_pretty_names[!is.na(match_pretty_names)]
         if (length(match_pretty_names)) {
-          x$Parameter[match_pretty_names] <- att$pretty_names[x$Parameter[match_pretty_names]]
+          x$Parameter[match_pretty_names] <- att$pretty_names[x$Parameter[
+            match_pretty_names
+          ]]
         }
       } else {
         x$Parameter <- att$pretty_names[x$Parameter]

diff --git a/R/extract_parameters.R b/R/extract_parameters.R
@@ -1,7 +1,7 @@
 # generic function ------------------------------------------------------

 #' @keywords internal
 .extract_parameters_generic <- function(
  model,
  ci,
  component,
@@ -522,7 +522,8 @@
         ci = ci,
         vcov = vcov,
         vcov_args = vcov_args,
-        verbose = verbose
+        verbose = verbose,
+        method = ci_method
       )
       fun_args <- c(fun_args, dots)
       ci_df <- suppressMessages(do.call("ci", fun_args))
@@ -567,7 +568,13 @@
 
   # p value
   if (!is.null(vcov)) {
-    fun_args <- list(model, vcov = vcov, vcov_args = vcov_args, verbose = verbose)
+    fun_args <- list(
+      model,
+      vcov = vcov,
+      vcov_args = vcov_args,
+      verbose = verbose,
+      method = ci_method
+    )
     fun_args <- c(fun_args, dots)
     parameters <- merge(
       parameters,
@@ -588,7 +595,7 @@
   } else {
     parameters <- merge(
       parameters,
-      p_value(model, dof = dof, effects = "fixed"),
+      p_value(model, dof = dof, method = ci_method, effects = "fixed"),
       by = "Parameter",
       sort = FALSE
     )

diff --git a/R/format.R b/R/format.R
@@ -3,7 +3,7 @@
 #' @inheritParams print.parameters_model
 #' @rdname print.parameters_model
 #' @export
 format.parameters_model <- function(
  x,
  pretty_names = TRUE,
  split_components = TRUE,
@@ -259,7 +259,7 @@
 #' @rdname print.compare_parameters
 #' @inheritParams print.parameters_model
 #' @export
 format.compare_parameters <- function(
  x,
  split_components = TRUE,
  select = NULL,
@@ -736,7 +736,7 @@


 # footer: type of uncertainty interval
 .print_footer_cimethod <- function(x) {
  if (!isTRUE(getOption("parameters_cimethod", TRUE))) {
    return()
  }
@@ -869,7 +869,8 @@
   show_re_msg <- (identical(model_class, "glmmTMB") &&
     (!string_method %in% c("Wald z-", "Wald normal") ||
       !ci_method %in% c("wald", "normal"))) ||
-    ((identical(model_class, "lmerMod") || identical(model_class, "glmerMod")) &&
+    (length(model_class) == 1L &&
+      model_class %in% c("lmerMod", "lmerModLmerTest", "glmerMod") &&
       !ci_method %in% c("wald", "normal", "profile", "boot"))
 
   if (
@@ -954,5 +955,5 @@
      "Some coefficients seem to be rather large, which may indicate issues with (quasi) complete separation. Consider using bias-corrected or penalized regression models."
    )
  }
  return(NULL)
 }