Estimated marginal means (+/-SE) in mixed-effect model for two or more fixed effects

After specifying a mixed-effect model (using lme4 package) with 2+ fixed factors, for example:

model = lmer(DV ~ facA * facB + (1|subject), data=DATA)

I would like to estimate the marginal means (+/-SE) for each cell of facA and/by facB adjusted for each other (optionally allowing their interaction), so I can use these values to create a plot in any software I wish. Thank you!