closing report as work has moved over to patch #10123

Hi Muhali,

I didn't write this function nor have I been involved in maintaining the Statistics package. I came across this bug submission searching the internet and I had a stab fixing it. The changes include:

- fixing the error degrees of freedom bug when modeltype != 'full'
- fixing the bug that required replication in all combinations of factors
- fixing the order of the factors; now they match the input and are consistent with Matlab
- listing the error term sum sq., df and meanSS at the bottom of the ANOVA table
- adding a row at the bottom of the table with the total sum sq. and df
- tabulating p-values in APA style
- fixing the format of the Sum Sq, MeanSS and Fval so the table isn't ugly when numbers become large
- changed slightly the interaction formatting in the table to match Matlab better
- fixed minor issue with Factor listing in Table title

Note that the function still requires numeric or logical IDs for the group matrix (whereas Matlab supports cell arrays of strings etc.)

I ran some tests including a 3-way ANOVA test but please feel free to test it yourself.

I made the changes on top of a recent patch #10123 at: https://savannah.gnu.org/patch/index.php?10123

Hopefully someone involved in the statistics package will be able to update the package with this patch!

Best

Andy

I just noticed the doc mentioning that the function only works for the full model. Default is the linear model, however, which is the one most people would probably choose.

Because of that, all but the full model should be removed from anovan.

Moreover, the function header also mentions (2nd point under TESTING RESULTS) that it gives identical results to Matlab. This is true except for the factors being swapped (very misleading!). So this must be fixed as well.