I want a marginal model, ideally fit via GLS. But the normality of residuals doesn't hold. It isn't much skewed, I don't want any transformations. It's just non-normal in shape. Yet still reporting means for that case is meaningful. The GEE estimation doesn't worry about that. Would you advice me to report rather GLS with sandwich SE, or just switch to GEE with sandwich ones?