Automating the analysis of variance of orthogonal designs