[R] RES: R Mixed Anova

[Rodrigo Aluizio]

Thanks Mark, it was exactly what I was looking for.

-----Mensagem original-----
De: r-help-bounces em r-project.org [mailto:r-help-bounces em r-project.org] Em nome de Mark Difford
Enviada em: quinta-feira, 6 de novembro de 2008 09:22
Para: r-help em r-project.org
Assunto: Re: [R] R Mixed Anova


Hi Rodrigo,

Here are two options; for each type, the second version gives 2nd order
interactions

## aov
T.aovmod <- aov(response ~ Season + Beach + Line + Error(Block/Strata))
T.aovmod <- aov(response ~ (Season + Beach + Line)^2 + Error(Block/Strata))

## lme
library(nlme)
T.lmemod <- lme(response ~ Season + Beach + Line, random = ~ 1 |
Block/Strata)
T.lmemod <- lme(response ~ (Season + Beach + Line)^2, random = ~ 1 |
Block/Strata)
anova(T.lmemod)
?anova.lme

## Use package multcomp for doing post-hoc analysis.

Regards, Mark.


Rodrigo Aluizio wrote:
> 
> Hi list, I was searching how to properly write a command line for a mixed
> ANOVA. Well honestly, there are so many material on the older post of the
> list that just confused me.
> 
> I have five factors.
> 
>  
> 
> Season (fixed)
> 
> Beach (fixed)
> 
> Line (fixed)
> 
> Block (random)
> 
> Strata (random) nested in Block
> 
>  
> 
> And for each of the tree strata per block I got 3 replicates. 
> 
>  
> 
> I saw lots of things about “different” linear models commands… lm, lme… On
> R
> help I found aov and anova…
> 
>  
> 
> Any one knows what should be more adequate and how I write and appropriate
> command line? I need an ANOVA table, a p value, and be able to run post
> hoc
> tests.
> 
>  
> 
> Thanks for your attention and patience.
> 
>  
> 
> ___________________________________
> MSc.  <mailto:r.aluizio em gmail.com> Rodrigo Aluizio
> Centro de Estudos do Mar/UFPR
> Laboratório de Micropaleontologia
> Avenida Beira Mar s/n - CEP 83255-000
> Pontal do Paraná - PR - BRASIL
> Fone: (0**41) 3455-1496 ramal 217
> Fax: (0**41) 3455-1105
> 
> 
> 
> 
> 	[[alternative HTML version deleted]]
> 
> 
> ______________________________________________
> R-help em r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
> 
> 

-- 
View this message in context: http://www.nabble.com/R-Mixed-Anova-tp20358811p20359309.html
Sent from the R help mailing list archive at Nabble.com.

______________________________________________
R-help em r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.