Combine or subset emmGrid objects
## S3 method for class 'emmGrid' rbind(..., deparse.level = 1, adjust = "bonferroni") ## S3 method for class 'emmGrid' e1 + e2 ## S3 method for class 'emm_list' rbind(..., which, adjust = "bonferroni") ## S3 method for class 'emmGrid' x[i, adjust, drop.levels = TRUE, ...]
... |
In |
deparse.level |
(required but not used) |
adjust |
Character value passed to |
e1 |
An |
e2 |
Another |
which |
Integer vector of subset of elements to use; if missing, all are combined |
x |
An |
i |
Integer vector of indexes |
drop.levels |
Logical value. If |
A revised object of class emmGrid
The result of e1 + e2 is the same as rbind(e1, e2)
The rbind method for emm_list objects simply combines
the emmGrid objects comprising the first element of ....
rbind throws an error if there are incompatibilities in
the objects' coefficients, covariance structures, etc. But they
are allowed to have different factors; a missing level '.'
is added to factors as needed.
warp.lm <- lm(breaks ~ wool * tension, data = warpbreaks) warp.rg <- ref_grid(warp.lm) # Do all pairwise comparisons within rows or within columns, # all considered as one faily of tests: w.t <- pairs(emmeans(warp.rg, ~ wool | tension)) t.w <- pairs(emmeans(warp.rg, ~ tension | wool)) rbind(w.t, t.w, adjust = "mvt") update(w.t + t.w, adjust = "fdr") ## same as above except for adjustment ### Working with 'emm_list' objects mod <- lm(conc ~ source + factor(percent), data = pigs) all <- emmeans(mod, list(src = pairwise ~ source, pct = consec ~ percent)) rbind(all, which = c(2, 4), adjust = "mvt") # Show only 3 of the 6 cases summary(warp.rg[c(2, 4, 5)])
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.