Hi, I tried to do a model simplification in a repeated anova with 3 factors (population,treatment,site). I tried to get rid of some interactions. Obviously it's not possible to use the function "anova" straight away to compare the different models. I got the suggestion to use the function proj (see below). In my data, however, I have in "block" not the third order interaction (n:p:k), but two factors and their interaction (population, site, population:site).
I somehow didn't manage to get the right aov after the projection. > From: "Richard M. Heiberger" <[EMAIL PROTECTED]> > > > npk.aovE <- aov(yield ~ N*P*K + Error(block), npk) > > np.k.aovE <- aov(yield ~ N*P+K + Error(block), npk) > > > npk.proj <- proj(npk.aovE) > npk.x <- npk.proj$block[,"N:P:K"] > > npk.aov <- aov(terms(yield ~ npk.x + block + N+P+K+N:P+N:K+P:K, > keep.order=TRUE), npk) > > np.k.aov <- aov(terms(yield ~ block + N+P+K+N:P, keep.order=TRUE), > npk) > > summary(npk.aov) > summary(np.k.aov) > > anova(npk.aov, np.k.aov) > ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
