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@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]]
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 commandslm, 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. 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.