Good code. Thanks! But something funny happened on the way to the ANOVA: When my M matrix has Fisher-z-transformed data, I only get ANOVA stats for the 2nd factor and nothing for the 1st factor (running a 2x2 model). When the data are inverse-transformed the ANOVA works fine and I get stats for both factors. Any thoughts?