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

From: Ronaldo Reis-Jr. <chrysopa_at_gmail.com>
Date: Thu 02 Jun 2005 - 00:18:09 EST


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@insecta.ufv.br            ]
|>/ (`. `'` ) \[ICQ#: 5692561 | LinuxUser#: 205366 ]
|    ( `-  )   [***********************************]
|>>  _/   \_Powered by GNU/Debian Woody/Sarge

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Thu Jun 02 00:23:35 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:32:20 EST