Multiple comparison for n-way ANOVA

Hi all,

I know I can do a multiple comparison test by passing the 'stats' output from the anova1 function to the multcompare function.

Right now, if I use anovan (n-way ANOVA), and I want to do a multcompare for the interaction of a pair of factors (X1 * X2), how do I do that?

If I just pass the stats struct from anovan to multcompare, I would just get a multiple comparison for the linear effect of the first factor X1, not X1*X2.

E.g.:

[p, ~, stats] = anovan(data, 'model', 'interaction');
[c,m] = multcompare(stats); % <-- this will do a multiple comp for X1 only, not X1*X2.