Tive mais um problema...
Na hora de estimar os componentes de variância por RELM a função convergiu para zero.
fit1 <- lmer(dados$Nota ~ (1|Experimento/Degustador) + (1|Experimento:Leite) +
(1|Experimento:Achocolatado) + (1|Experimento:Café) + (1|Experimento:Leite:Achocolatado:Café) + Café + Leite + Achocolatado +Leite*Achocolatado*Café)
summary(fit1)
AIC BIC logLik deviance REMLdev
1823 1885 -896.7 1789 1793
Random effects:
Groups Name Variance Std.Dev.
Degustador:Experimento (Intercept) 0.945490 0.97236
Experimento:Leite:Achocolatado:Café (Intercept) 0.059272 0.24346
Experimento:Café (Intercept) 0.000000 0.00000
Experimento:Achocolatado (Intercept) 0.148591 0.38548
Experimento:Leite (Intercept) 0.034867 0.18673
Experimento (Intercept) 0.000000 0.00000
Residual 2.847312 1.68740
Isso pode interferir nos meus resíduos marginais e condicionais?
Tem como corrigir?
Obrigado,
Greice.
2011/6/30 Greice Laureano
<greice.laureano@gmail.com>
Vou fazer pela lmer() mesmo então.
Obridado novamente,
Greice.
Greice,
Eu suspeito que não seja possível declarar a sua estrutura de efeitos aleatórios na lme().
À disposição.
Walmes.
==========================================================================Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573VoIP: (3361 3600) 1053 1173
e-mail: walmes@ufpr.brtwitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmeslinux user number: 531218
==========================================================================