Ronaldo Reis-Jr.
2005-Jun-02 11:43 UTC
[R] How to calculate the correct SE in a nested or spliplot anova?
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))> mCall: 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))> mCall: 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