estimate_means function

Estimate Marginal Means (Model-based average at each factor level)