I want to calculate estimated group mean scores in a 2×2 Gaussian regression after obtaining the regression coefficients. Here is toy data. 100 observations each of region – a and b – and sex – m and f. I have designed the scores so there is a 5-point difference on average between regions a and b but no difference between m and f.