R matches Psy output for contrast sigma (sqrt(MSE * sum (c^2/n)) when: for between contrasts: n = residual df for within contrasts: n = n subs