[R] How to automatic obtain the p value from a nested aov?

[Ronaldo Reis-Jr.]

Hi,

I try this nested anova:

> m <- aov(Glycogen~Treatment+Error(Treatment/Rat/Liver))
> summary(m)

Error: Treatment
          Df  Sum Sq Mean Sq
Treatment  2 1557.56  778.78

Error: Treatment:Rat
          Df Sum Sq Mean Sq F value Pr(>F)
Residuals  3 797.67  265.89               

Error: Treatment:Rat:Liver
          Df Sum Sq Mean Sq F value Pr(>F)
Residuals 12  594.0    49.5               

Error: Within
          Df Sum Sq Mean Sq F value Pr(>F)
Residuals 18 381.00   21.17               
> 

My objective is the effect of treatment.

I know how to calculate the F for this:

778.78/265.89

But, exist any package that make this operation automatic? Is for a class 
work.

Thanks
Ronaldo
-- 
Vou rezar 1/3, para encontrar 1/2 de te levar para 1/4
--
|>   // | \\   [***********************************]
|   ( õ   õ )  [Ronaldo Reis Júnior                ]
|>      V      [UFV/DBA-Entomologia                ]
|    /     \   [36570-000 Viçosa - MG              ]
|>  /(.''`.)\  [Fone: 31-3899-4007                 ]
|  /(: :'  :)\ [chrysopa at insecta.ufv.br            ]
|>/ (`. `'` ) \[ICQ#: 5692561 | LinuxUser#: 205366 ]
|    ( `-  )   [***********************************]
|>>  _/   \_Powered by GNU/Debian Woody/Sarge