[R-br] Modelo misto

[Greice Laureano]

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 em gmail.com>

> Vou fazer pela lmer() mesmo então.
> Obridado novamente,
> Greice.
>
>  2011/6/30 Walmes Zeviani <walmeszeviani em 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 em ufpr.br
>> twitter: @walmeszeviani
>> homepage: http://www.leg.ufpr.br/~walmes
>> linux user number: 531218
>> ==========================================================================
>>
>>
>> _______________________________________________
>> R-br mailing list
>> R-br em 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.
>>
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110630/82cd74b4/attachment.html>