[R-br] Parametrização com lme4
walmes .
walmeszeviani em gmail.com
Quarta Novembro 20 11:05:27 BRST 2013
Fernando,
Acho que é assim:
## fatores cruzados
xtabs(~local+trat, dados)
## fatores aninhados
xtabs(~trat+eras, dados)
## local
## +bloco dentro de local
## +tratamento
## +interação local tratamento
mf1 <- lm(terms(resp~local/rep+trat+local:trat, keep=TRUE),
data=dados)
anova(mf1)
## local
## +bloco dentro de local
## +eras
## +tratamento dentro de eras
## +interação local eras
## +interação local tratamento dentro de eras
mf2 <- lm(terms(resp~local/rep+eras/trat+local:eras/trat, keep=TRUE),
data=dados)
anova(mf2)
require(lme4)
mm1 <- lmer(resp~(1|local/rep)+(1|trat)+(1|local:trat), data=dados)
summary(mm1)
mm2 <- lmer(resp~(1|local/rep)+(1|eras/trat)+
(1|local:eras)+(1|local:eras:trat), data=dados)
summary(mm2)
À disposição.
Walmes.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131120/7c3011d6/attachment.html>
Mais detalhes sobre a lista de discussão R-br