[R] How to calculate the correct SE in a nested or spliplot anova?

From: Ronaldo Reis-Jr. <chrysopa_at_gmail.com>
Date: Thu 02 Jun 2005 - 21:43:25 EST


Hi!

How to calculate the correct SE of mean in a nested or spliplot anova?

Nested example:



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

Call:
aov(formula = Glycogen ~ Treatment + Error(Treatment/Rat/Liver))

Grand Mean: 142.2222

Stratum 1: Treatment

Terms:

                Treatment
Sum of Squares   1557.556
Deg. of Freedom         2

Estimated effects may be unbalanced

Stratum 2: Treatment:Rat

Terms:

                Residuals
Sum of Squares   797.6667
Deg. of Freedom         3

Residual standard error: 16.3061

Stratum 3: Treatment:Rat:Liver

Terms:

                Residuals
Sum of Squares        594
Deg. of Freedom        12

Residual standard error: 7.035624

Stratum 4: Within

Terms:

                Residuals
Sum of Squares        381
Deg. of Freedom        18

Residual standard error: 4.600725
>



I need to make a barplot with error bar for treatment. How to calculate the correct SE?

SplitPLot example:



> m <-

aov(yield~irrigation*density*fertilizer+Error(block/irrigation/density/fertilizer))
> m

Call:
aov(formula = yield ~ irrigation * density * fertilizer + Error(block/irrigation/density/fertilizer))

Grand Mean: 99.72222

Stratum 1: block

Terms:

                Residuals
Sum of Squares   194.4444
Deg. of Freedom         3

Residual standard error: 8.050765

Stratum 2: block:irrigation

Terms:

                irrigation Residuals
Sum of Squares    8277.556  1411.778
Deg. of Freedom          1         3

Residual standard error: 21.69315
8 out of 9 effects not estimable
Estimated effects are balanced

Stratum 3: block:irrigation:density

Terms:

                 density irrigation:density Residuals
Sum of Squares  1758.361           2747.028  2787.944
Deg. of Freedom        2                  2        12

Residual standard error: 15.24233
8 out of 12 effects not estimable
Estimated effects may be unbalanced

Stratum 4: block:irrigation:density:fertilizer

Terms:

                fertilizer irrigation:fertilizer density:fertilizer
Sum of Squares   1977.4444              953.4444           304.8889
Deg. of Freedom          2                     2                  4
                irrigation:density:fertilizer Residuals
Sum of Squares                       234.7222 3108.8333
Deg. of Freedom                             4        36

Residual standard error: 9.292819
Estimated effects may be unbalanced


I need the SE for each combination. How to calculate this?

Thanks for all.
Ronaldo

--
|>   // | \\   [***********************************]
|   ( §   § )  [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 21:53:50 2005

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