[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