satterthwaite

Author: strengejacke

NAMESPACE

@@ -114,6 +114,8 @@ S3method(ci,systemfit)
 S3method(ci,varest)
 S3method(ci,zerocount)
 S3method(ci,zeroinfl)
+S3method(ci_kenward,default)
+S3method(ci_kenward,glmmTMB)
 S3method(cluster_discrimination,cluster_analysis)
 S3method(cluster_discrimination,default)
 S3method(cluster_performance,dbscan)
@@ -138,7 +140,6 @@ S3method(display,parameters_pca_summary)
 S3method(display,parameters_sem)
 S3method(display,parameters_simulate)
 S3method(display,parameters_standardized)
-S3method(dof_satterthwaite,lmerMod)
 S3method(equivalence_test,MixMod)
 S3method(equivalence_test,estimate_contrasts)
 S3method(equivalence_test,estimate_means)

R/ci_satterthwaite.R

@@ -2,13 +2,18 @@
 #' @export
 ci_satterthwaite <- function(model, ci = 0.95, ...) {
   df_satter <- dof_satterthwaite(model)
+  if (inherits(model, "glmmTMB")) {
+    component <- "conditional"
+  } else {
+    component <- "all"
+  }
   out <- lapply(ci, function(i) {
     .ci_dof(
       model = model,
       ci = i,
       dof = df_satter,
       effects = "fixed",
-      component = "all",
+      component = component,
       method = "satterthwaite",
       ...
     )

R/dof_satterthwaite.R

@@ -1,17 +1,5 @@
 #' @rdname p_value_satterthwaite
 #' @export
 dof_satterthwaite <- function(model) {
-  UseMethod("dof_satterthwaite")
-}
-
-
-#' @export
-dof_satterthwaite.lmerMod <- function(model) {
-  insight::check_if_installed("lmerTest")
-
-  parameters <- insight::find_parameters(model, effects = "fixed", flatten = TRUE)
-  lmerTest_model <- lmerTest::as_lmerModLmerTest(model)
-  s <- summary(lmerTest_model)
-
-  stats::setNames(as.vector(s$coefficients[, 3]), parameters)
+  insight::get_df(model, "satterthwaite")
 }

man/p_value_kenward.Rd

@@ -29,6 +29,12 @@ mp2 <- model_parameters(model, digits = 5, ci_method = "s", effects = "fixed")
 mp3 <- model_parameters(model, digits = 5, ci_method = "kr", effects = "fixed")
 mp4 <- model_parameters(model, digits = 5, ci_method = "wald", effects = "fixed")
 mp5 <- model_parameters(model3, digits = 5, ci_method = "kr", effects = "fixed")
+mp6 <- model_parameters(
+  model3,
+  digits = 5,
+  ci_method = "satterthwaite",
+  effects = "fixed"
+)
 
 test_that("model_parameters, ci_method default (residual)", {
   expect_equal(
@@ -85,6 +91,8 @@ test_that("model_parameters, ci_method satterthwaite", {
     c(24.57489, 4.93385, -1.91805, -0.05477, -2.96368, -4.41987, -0.16904, -0.05117),
     tolerance = 1e-3
   )
+  expect_equal(mp2$SE, mp6$SE, tolerance = 1e-3)
+  expect_equal(mp2$df_error, mp6$df_error, tolerance = 1e-3)
 })
 
 test_that("model_parameters, ci_method kenward", {
@@ -114,13 +122,13 @@ test_that("model_parameters, ci_method kenward", {
 
   expect_warning(expect_warning(expect_warning(
     {
-      mp6 <- model_parameters(model4, digits = 5, ci_method = "kr", effects = "fixed")
+      mp7 <- model_parameters(model4, digits = 5, ci_method = "kr", effects = "fixed")
     },
     regex = "Model was not fitted"
   )))
 
-  expect_equal(mp5$SE, mp6$SE, tolerance = 1e-3)
-  expect_equal(mp5$df_error, mp6$df_error, tolerance = 1e-3)
+  expect_equal(mp5$SE, mp7$SE, tolerance = 1e-3)
+  expect_equal(mp5$df_error, mp7$df_error, tolerance = 1e-3)
 })
 
 test_that("model_parameters, ci_method wald (t)", {