this should be

easystats · Aug 8, 2024 · c710673 · c710673
1 parent 4ecd1d0
commit c710673
Showing 1 changed file with 7 additions and 4 deletions.
diff --git a/R/demean.R b/R/demean.R
@@ -448,10 +448,13 @@ degroup <- function(x,
     })
     # create de-meaned variables by subtracting the group mean from each individual value
     person_means_list <- lapply(
-      seq_along(select),
-      function(i) dat[[select[i]]] - group_means_list[[i]][[length(by)]]
-      # select,
-      # function(i) dat[[i]] - stats::ave(dat[[i]], dat[[by[length(by)]]], FUN = gm_fun)
+      # seq_along(select),
+      # function(i) dat[[select[i]]] - group_means_list[[i]][[length(by)]]
+      select,
+      function(i) {
+        dat$higher_levels <- do.call(paste, c(dat[by], list(sep = "_")))
+        dat[[i]] - stats::ave(dat[[i]], dat$higher_levels, FUN = gm_fun)
+      }
     )
   } else {
     # cross-classified design: by > 1