[R] Model simplification using anova()

[ChCh]

Hello all,

I've become confused by the output produced by a call to
anova(model1,model2). First a brief background.  My model used to predict
final tree height is summarised here:

                            Df Sum Sq Mean Sq F value    Pr(>F)    
Treatment              2   748.35   374.17 21.3096   7.123e-06 ***
HeightInitial            1   0.31      0.31     0.0178    0.89519    
DiameterInitial         1   0.52      0.52     0.0298    0.86460    
Frost                     1   38.29    38.29    2.1807    0.15392    
HeightInitial:Frost    1   85.83     85.83   4.8882    0.03774 *  
DiameterInitial:Frost 1   97.90    97.90    5.5754    0.02749 *  
Residuals                22 386.30   17.56                      
---

Based on this, I should not remove either of the interaction terms, so I
turned my attention to the main factors. Based on p-values, I removed
HeightInitial and used a call to anova(model1,model2) to see if this
resulted in a weaker model.  Here is the output:

Model 1: HeightFinal ~ Treatment + HeightInitial + DiameterInitial + Frost +
HeightInitial:Frost + DiameterInitial:Frost

Model 2: AbsoluteDiameterDiff ~ Treatment + DiameterInitial + Frost +
HeightInitial:Frost + DiameterInitial:Frost

  Res.Df   RSS Df Sum of Sq F Pr(>F)
1     22 386.3                      
2     22 386.3  0 -1.08e-12        

This is not the output that I'm used to seeing.  Typically, a pvalue would
be provided to suggest if the models differ significantly from one another.  

Can anyone explain why there's no pvalue in this situation and whether or
not removing HeightInitial was justified.

Thanks,
Jacob Freel
-- 
View this message in context: http://www.nabble.com/Model-simplification-using-anova%28%29-tp17636012p17636012.html
Sent from the R help mailing list archive at Nabble.com.