Add tests for estimate_grouplevel.coxme (#553)

* Add tests for estimate_grouplevel.coxme

* fix

* fix warning

Author: strengejacke

DESCRIPTION

@@ -1,7 +1,7 @@
 Type: Package
 Package: modelbased
 Title: Estimation of Model-Based Predictions, Contrasts and Means
-Version: 0.12.0.14
+Version: 0.12.0.15
 Authors@R:
     c(person(given = "Dominique",
              family = "Makowski",
@@ -58,6 +58,7 @@ Suggests:
     coda,
     collapse,
     correlation,
+    coxme,
     curl,
     discovr,
     easystats,
@@ -113,4 +114,4 @@ Roxygen: list(markdown = TRUE)
 Config/Needs/check: stan-dev/cmdstanr
 Config/Needs/website: easystats/easystatstemplate
 LazyData: true
-Remotes: easystats/insight
+Remotes: easystats/insight, easystats/parameters

NEWS.md

@@ -29,6 +29,8 @@
 * `estimate_expectation()` and `estimate_relation()` now support objects of
   class `htest`.
 
+* `estimate_grouplevel()` now supports models from package *coxme*.
+
 * Documentation of the `display()` method for *modelbased* objects has been
   added.
 

inst/WORDLIST

@@ -102,6 +102,7 @@ conditionspeed
 counterfactuals
 confounder
 confounders
+coxme
 d'être
 d’être
 dat

tests/testthat/test-estimate_grouplevel.R

@@ -140,3 +140,38 @@ test_that("estimate_grouplevel - Bayesian, rstanarm", {
   expect_identical(dim(out), c(3L, 4L))
   expect_named(out, c("Group", "Level", "Parameter", "Coefficient"))
 })
+
+
+skip_if_not_installed("withr")
+
+withr::with_environment(
+  new.env(),
+  test_that("estimate_grouplevel, coxme", {
+    skip_on_cran()
+    skip_if_not_installed("coxme")
+    skip_if_not_installed("survival")
+
+    Surv <- survival::Surv
+    rats <- survival::rats
+    lung <- survival::lung
+
+    set.seed(1234)
+    rats$grp <- sample(letters[1:3], nrow(rats), replace = TRUE)
+
+    data(eortc, package = "coxme")
+    d <<- coxme::eortc
+    d2 <<- rats
+    d3 <<- lung
+
+    m1 <- coxme::coxme(Surv(y, uncens) ~ trt + (1 | center), data = d)
+    m2 <- coxme::coxme(Surv(time, status) ~ ph.ecog + age + (1 | inst), d3)
+    m3 <- coxme::coxme(Surv(time, status) ~ rx + (1 + rx | litter) + (1 | grp), d2)
+
+    out <- estimate_grouplevel(m1)
+    expect_identical(dim(out), c(37L, 5L))
+    out <- estimate_grouplevel(m2)
+    expect_identical(dim(out), c(18L, 5L))
+    out <- estimate_grouplevel(m3)
+    expect_identical(dim(out), c(203L, 5L))
+  })
+)