Complex PERMANOVA mixed model nested design in PRIMER

Hi All,

I am trying to use a complex ANOVA design with the Permanova function in PRIMER. I have been working on this for a while with no significant progress, so any ideas or help would be greatly appreciated.

The SAS code for my design is:

proc mixed data=treatment;
class Gap Wood Collar Year Month Plot;
model Biomass= Gap|Wood|Collar|Year|Month / outp=Biomass2;
random plot(Gap*Wood) plot(Gap*Wood*Collar) plot(Gap*Wood*Collar*Year);
run;

Resulting in denominator DF of 16 for Gap, Wood, Gap*Wood; 32 for collar and collar* combinations of gap and wood; 48 for year and combinations of gap and wood and collar; and 93 for everything else.


The problem I am having in PRIMER is that I can only nest plot in one combination of these terms, resulting in denominator DF that are too large. I think that the answer to this question has something to do with pooling terms. Any help or direction that you can give would be greatly appreciated. Thanks.


Sincerely,
Tera