Fix to 'nuisance' code (#503)

rvlenth · Jul 18, 2024 · 909821d · 909821d
1 parent 9732845
commit 909821d
Show file tree

Hide file tree

Showing 3 changed files with 4 additions and 3 deletions.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,8 +1,8 @@
 Package: emmeans
 Type: Package
 Title: Estimated Margin2l Means, aka Least-Squares Means
-Version: 1.10.3-090003
-Date: 2024-07-08
+Version: 1.10.3-090004
+Date: 2024-07-18
 Authors@R: c(person("Russell V.", "Lenth", role = c("aut", "cre", "cph"), 
     email = "russell-lenth@uiowa.edu"),
     person("Ben", "Bolker", role = "ctb"),

diff --git a/NEWS.md b/NEWS.md
@@ -10,6 +10,7 @@ title: "NEWS for the emmeans package"
     was previously not allowed, and we added an `"allow.na.levs"` option
     (defaults to `TRUE`) just in case we broke anything that used to work.
   * Better default contrasts in `qdrg()` (#501)
+  * Bug fix for nuisance factors when we have a multivariate response (#503)
 
 
 ## emmeans 1.10.3

diff --git a/R/ref-grid.R b/R/ref-grid.R
@@ -1126,7 +1126,7 @@ ref_grid <- function(object, at, cov.reduce = mean, cov.keep = get_emm_option("c
     basis$misc$avgd.over = paste(length(nuis), "nuisance factors")
     RA = rep(ra, k)
     basis$X = X[RA == ".main.grid.", , drop = FALSE]
-    basis$grid = grid[RA == ".main.grid.", , drop = FALSE]
+    basis$grid = grid[ra == ".main.grid.", , drop = FALSE]
     non.nuis = setdiff(names(ref.levels), info$nuis)
     basis$ref.levels = ref.levels[non.nuis]
     basis