bug fix for survey::svyglm() models with replicate weights

fix #240

NAMESPACE

@@ -73,6 +73,7 @@ S3method(model_get_terms,glmmTMB)
 S3method(model_get_terms,model_fit)
 S3method(model_get_weights,default)
 S3method(model_get_weights,model_fit)
+S3method(model_get_weights,svrepglm)
 S3method(model_get_weights,svyglm)
 S3method(model_get_xlevels,brmsfit)
 S3method(model_get_xlevels,default)

NEWS.md

@@ -8,6 +8,7 @@
 
 - fix in `supported_models`
 - bug fix when using `tidy_parameters()` for mixed models (#238)
+- bug fix for `survey::svyglm()` models with replicate weights (#240)
 
 # broom.helpers 1.14.0
 

R/model_get_weights.R

@@ -3,6 +3,10 @@
 #' This function does not cover `lavaan` models (`NULL` is returned).
 #'
 #' @param model a model object
+#' @note
+#' For class `svrepglm` objects (glm on a survey object with replicate weights),
+#' it will return the original sampling weights of the data, not the replicate
+#' weights.
 #' @export
 #' @family model_helpers
 #' @examples
@@ -69,6 +73,12 @@ model_get_weights.svyglm <- function(model) {
   stats::weights(model$survey.design)
 }
 
+#' @export
+#' @rdname model_get_weights
+model_get_weights.svrepglm <- function(model) {
+  model$survey.design$pweights
+}
+
 #' @export
 #' @rdname model_get_weights
 model_get_weights.model_fit <- function(model) {

tests/testthat/test-tidy_plus_plus.R

@@ -315,6 +315,17 @@ test_that("tidy_plus_plus() works with survey::svyglm", {
     res <- mod %>% tidy_plus_plus(),
     NA
   )
+
+  df_rep <- survey::as.svrepdesign(df)
+  mod_rep <- survey::svyglm(
+    response ~ age + grade * trt,
+    df_rep,
+    family = quasibinomial
+  )
+  expect_error(
+    res <- mod_rep %>% tidy_plus_plus(),
+    NA
+  )
 })
 
 test_that("tidy_plus_plus() works with survey::svycoxph", {