[R] Ancova doesn't return test statistics
    Richard M. Heiberger 
    rmh at temple.edu
       
    Thu Nov 15 18:11:46 CET 2007
    
    
  
Your model is fully saturated.  It specifies terms that use
up all degrees of freedom.  There are no degrees of freedom left
over for a Residual term and therefore there is no denominator for 
the tests.
When you drop one term, then those degrees of freedom are left over,
that is they form the Residual, and are used as the denominator for
the tests.  The usual practice is to suppress the high-order interactions,
in your example by
model <- aov(GR ~ SR*HS*Pop*Popsize*Year - SR:HS:Pop:Popsize:Year)
Please use spaces around the arrow, tilde, and + and - signs for
legibility.
Rich
    
    
More information about the R-help
mailing list