Dear all,
I have follwing factors in my data set (G) names(G) relative ht, biomass, treatment, variety, soil group relative height biomass treatment (4 levels)- control+lime, control+ash, contaminated soil+lime, contaminated soil+ash variety (4 levels)- Salix swherinii, Salix mursinifolia, Klara, Karin soilgroup (2 level)- control and contaminated if I fit aov function this is ok anova1<-aov(ht~treatment*variety,data=G) summary(anova1) R- output Df Sum Sq Mean Sq F value Pr(>F) treatment 3 14015 4672 65.021 < 2e-16 *** variety 3 5883 1961 27.295 2.95e-12 *** treatment:variety 9 5474 608 8.465 8.93e-09 *** Residuals 80 5748 72 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 above I can see my treatment and variety interaction effect anova3<-aov(ht~treatment*variety*soilgroup,data=G) summary(anova3) R-ouput Df Sum Sq Mean Sq F value Pr(>F) treatment 3 14015 4672 65.021 < 2e-16 *** variety 3 5883 1961 27.295 2.95e-12 *** treatment:variety 9 5474 608 8.465 8.93e-09 *** soilgroup -missing treatment:variety:soilgroup - missing Here, I am missing value of soil group and interation of treatment: variety: soilgroup- Can any one please tell me why? Thanks for help Best regards Salam [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.