Residuals perfectly symmetric about zero against fitted values

Consider a modelling a response $Y$ against two categorical variables (which can take $4\times 2=8$ possible combinations). We have 16 values for the response, with two values for every combination of categorical variables.

When plotting the studentised residuals against the fitted values (so we have two residuals for every one of the 8 distinct fitted $y$ values), I am getting perfect symmetry about $0$. I know symmetry about zero is expected as the studentised residuals have the $t$-distribution, but perfect symmetric seems rather odd and makes me uncertain about the linear model.

If this is indeed true, could someone post a sketch proof? I am guessing this is linked to all the variables here being categorical (so matrix $X$ in the model $Y = X\beta +\varepsilon$ only contains the entries $0$ or $1$) somehow leading to the the fitted values coming out to be the average of the two values for every separate combination of categorical variables.

• At each level you have one mean fitted to two values. That mean will be at the average of the observations, so the corresponding residuals must be equidistant and opposite in sign. – Glen_b Nov 30 '14 at 6:07