Estimated Marginal Means, aka Least-Squares Means
Obtain estimated marginal means (EMMs) for many linear, generalized
linear, and mixed models. Compute contrasts or linear functions of EMMs,
trends, and comparisons of slopes. Plots and other displays.
Least-squares means are discussed, and the term "estimated marginal means"
is suggested, in Searle, Speed, and Milliken (1980) Population marginal means
in the linear model: An alternative to least squares means, The American
Statistician 34(4), 216-221 <doi:10.1080/00031305.1980.10483031>.
R package emmeans: Estimated marginal means
Note: emmeans is a continuation of the package lsmeans. The latter will eventually be retired.
Estimated marginal means (EMMs, previously known as least-squares means in the context of traditional regression models) are derived by using a model to make predictions over a regular grid of predictor combinations (called a reference grid). These predictions may possibly be averaged (typically with equal weights) over one or more of the predictors. Such marginally-averaged predictions are useful for describing the results of fitting a model, particularly in presenting the effects of factors. The emmeans package can easily produce these results, as well as various graphs of them (interaction-style plots and side-by-side intervals).
Estimation and testing of pairwise comparisons of EMMs, and several other types of contrasts, are provided. There is also a
cldmethod for display of grouping symbols.
Two-way support of the
glhtfunction in the multcomp package.
For models where continuous predictors interact with factors, the package's
emtrendsfunction works in terms of a reference grid of predicted slopes of trend lines for each factor combination.
Vignettes are provided on various aspects of EMMs and using the package. See the CRAN page
The package incorporates support for many types of models, including standard models fitted using
glm, and relatives, various mixed models, GEEs, survival models, count models, ordinal responses, zero-inflated models, and others. Provisions for some models include special modes for accessing different types of predictions; for example, with zero-inflated models, one may opt for the estimated response including zeros, just the linear predictor, or the zero model. For details, see
vignette("models", package = "emmeans")
Various Bayesian models (carBayes, MCMCglmm, MCMCpack) are supported by way of creating a posterior sample of least-squares means or contrasts thereof, which may then be examined using tools such as in the coda package.
Package developers may provide emmeans support for their models by writing
vignette("extending", package = "emmeans")
Versions and installation
CRAN The latest CRAN version may be found at https://CRAN.R-project.org/package=emmeans. Also at that site, formatted versions of this package's vignettes may be viewed.
Github To install the latest development version from Github, install the newest version (definitely 2.0 or higher) of the devtools package; then run
devtools::install_github("rvlenth/emmeans", dependencies = TRUE, build_opts = "") ### To install without vignettes (faster): devtools::install_github("rvlenth/emmeans")
Note: If you are a Windows user, you should also first download and
install the latest version of
For the latest release notes on this development version, see the NEWS file
Functions in emmeans
|as.emmGrid||Convert to and from emmGrid objects|
|add_grouping||Add a grouping factor|
|CLD||Extract and display information on all pairwise comparisons of estimated marginal means.|
|MOats||Oats data in multivariate form|
|emm_options||Set or change emmeans options|
|emm_list||The emm_list class|
|auto.noise||Auto Pollution Filter Noise|
|contrast||Contrasts and linear functions of EMMs|
|emmeans-package||Estimated marginal means (aka Least-squares means)|
|emmip||Interaction-style plots for estimated marginal means|
|extending-emmeans||Support functions for model extensions|
|emmeans||Estimated marginal means (Least-squares means)|
|emmobj||Construct an emmGrid object from scratch|
|emtrends||Estimated marginal means of linear trends|
|emmGrid-class||The emmGrid class|
|pigs||Effects of dietary protein on free plasma leucine concentration in pigs|
|str.emmGrid||Miscellaneous methods for emmGrid objects|
|oranges||Sales of oranges|
|lsmeans||Wrappers for alternative naming of EMMs|
|update.emmGrid||Update an emmGrid object|
|joint_tests||Compute joint tests of the terms in a model|
|emm||Support for multcomp::glht|
|hpd.summary||Summarize an emmGrid from a Bayesian model|
|plot.emmGrid||Plot an emmGrid or summary_emm object|
|rbind.emmGrid||Combine or subset emmGrid objects|
|ref_grid||Create a reference grid from a fitted model|
|as.mcmc.emmGrid||Support for MCMC-based estimation|
|xtable.emmGrid||Using xtable for EMMs|
|models||Models supported in emmeans|
|summary.emmGrid||Summaries, predictions, intervals, and tests for emmGrid objects|
|regrid||Reconstruct a reference grid with a new transformation or posterior sample|
|pwpp||Pairwise P-value plot|
|qdrg||Quick and dirty reference grid|
Vignettes of emmeans
Last month downloads
|License||GPL-2 | GPL-3|
|Packaged||2019-09-11 22:32:18 UTC; rlenth|
|Date/Publication||2019-09-12 05:10:03 UTC|
|suggests||bayesplot , biglm , brms , car , coda (>= 0.17) , ggplot2 , knitr , lattice , lme4 , lmerTest (>= 2.0.32) , MASS , mediation , mgcv , multcomp , multcompView , MuMIn , nlme , ordinal (>= 2014.11-12) , pbkrtest (>= 0.4-1) , rmarkdown , rsm , scales , testthat|
|enhances||CARBayes , coxme , gee , geepack , MCMCglmm , MCMCpack , mice , nnet , pscl , rstanarm , sommer , survival|
|imports||estimability (>= 1.3) , graphics , methods , mvtnorm , numDeriv , plyr , stats , utils , xtable (>= 1.8-2)|
|depends||R (>= 3.2)|
|Contributors||Henrik Singmann, Maxime Hervé, Jonathon Love, Paul Buerkner|
Include our badge in your README