[R-br] desdobramento de uma analise fatorial

walmes . walmeszeviani em gmail.com
Segunda Setembro 29 11:52:32 BRT 2014


Segue um CMR.

data(ex5, package="ExpDes")
str(ex5)

xtabs(~trat+genero, ex5)

m0 <- aov(sabor~bloco+trat/genero, data=ex5)

names(coef(m0))

m0$assign
gINt.effects <- names(coef(m0))[m0$assign==3]

generoINtrat <- lapply(paste0("trat", levels(ex5$trat)),
                       grep, gINt.effects)

summary(m0, split=list("trat:genero"=generoINtrat))

m1 <- aov(sabor~bloco+genero/trat, data=ex5)

names(coef(m1))

m1$assign
tINg.effects <- names(coef(m1))[m1$assign==3]

tratINgenero <- lapply(paste0("genero", levels(ex5$genero)),
                       grep, tINg.effects)

summary(m1, split=list("genero:trat"=tratINgenero))

À disposição.
Walmes.

​
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140929/cec58128/attachment.html>


Mais detalhes sobre a lista de discussão R-br