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.

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 at 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 at 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.

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.

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.