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

[Ronaldo Reis-Jr.]

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