[R] help in anova

Mir Salam mir.salam at uef.fi
Wed Feb 11 15:47:18 CET 2015


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]]



More information about the R-help mailing list