Estimation of group means when adjusting for covariates in generalized linear models