## NOT RUN
# library(predictmeans)
# Oats$nitro <- factor(Oats$nitro)
# fm <- lmer(yield ~ nitro*Variety+(1|Block/Variety), data=Oats)
#
# # Permutation Test for model terms
# permanova.lmer(fm)
# permanova.lmer(fm, type=2)
#
# # Compare to F test
# fm0 <- lme(yield ~ nitro*Variety, random=~1|Block/Variety, data=Oats)
# anova(fm0)
Run the code above in your browser using DataLab