
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.
2011/6/30 Walmes Zeviani <walmeszeviani@gmail.com>
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 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.