Oi Walmes
obrigada por responder.
Então, o modelo é estimado com seguinte código:
model2<-lme(form2, random=list(const=pdBlocked(list(~ass-1, ~ass:prod-1),pdClass="pdIdent")),data=panel2)
onde -1 indica remover o intercepto.
O comando summary(model2) não retorna os componentes de variância.
O erro deve estar na forma como eu especifiquei meu modelo.
Vou tentar corrigir aqui.
Obrigada
Isabel
Isabel,
O seu modelo é estimado com o seu código? A princípio, se o summary() retorna os componentes de variância o VarCorr() deveria fazê-lo também. Só não estou seguro sobre a forma que você está usando no argumento random=, o que significa o "-1" (remover intercepto?), e certifique-se se a lme() entende o ":" pois é usual declarar com "/", ou seja, "A/B" ao invés de "A+A:B". Se não for nada disso, envie nos um CMR.
À 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
==========================================================================
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-brLeia o guia de postagem (
http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.