[R] anova and lm results differ

Iasonas Lamprianou lamprianou at yahoo.com
Tue Aug 31 18:11:02 CEST 2010


Dear all
   
 I have found that the two "equivalent" commands do not produce the same results.
    1. (I wrote this command by hand, this is what I would do usually)
    >summary(aov(eduyrs ~ cntry * edf, data=ESS1))
                  Df Sum Sq Mean Sq  F value    Pr(>F)   
    cntry          1    257  256.65  21.2251 4.243e-06 ***
    edf            4  11010 2752.42 227.6296 < 2.2e-16 ***
    cntry:edf      4     24    6.03   0.4987    0.7367   
    Residuals   3205  38754   12.09                       
    ---
    Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
   
 2. (I run this from Rcmdr, this is what my students would normally do)
    >AnovaModel.2 <- (lm(eduyrs ~ cntry*edf, data=ESS1))
    > Anova(AnovaModel.2)
    Anova Table (Type II tests)
    Response: eduyrs
              Sum Sq   Df  F value  Pr(>F)   
    cntry         38    1   3.1158 0.07763 . 
    edf        11010    4 227.6296 < 2e-16 ***
    cntry:edf     24    4   0.4987 0.73672   
    Residuals  38754 3205                     
    ---
    Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
   
    Unfortunately, the results are different. Could anyone, please, explain why?
    

Dr. Iasonas Lamprianou


Assistant Professor (Educational Research and Evaluation)
Department of Education Sciences
European University-Cyprus
P.O. Box 22006
1516 Nicosia
Cyprus 
Tel.: +357-22-713178
Fax: +357-22-590539


Honorary Research Fellow
Department of Education
The University of Manchester
Oxford Road, Manchester M13 9PL, UK
Tel. 0044  161 275 3485
iasonas.lamprianou at manchester.ac.uk






More information about the R-help mailing list