O  "-1" é aceito na função lme. Ele significa retirar o intercepto do modelo aleatório. significa dizer que os efeitos aleatórios não estão associados com o intercepto (intercepto fixo) e somente com as inclinações (pelo menos foi isto que entendi do livro do bates). Acredito que o erro está na forma como esta inserido o modelo aleatório. Entretanto ela utilizou duas vezes o -1 e aí sim não tem sentido uma vez que ela já o removeu. talvez isto funcione

model2<-lme(atr1~prod, random=list(const=pdBlocked(list(~-1,~ass, ~ass:prod),pdClass="pdIdent")),data=panel2)

Em 20 de abril de 2012 12:43, Walmes Zeviani <walmeszeviani@gmail.com> escreveu:
Isabel, a menos que eu esteja cometendo um grande engano, não existe isso de remover o intercepto na formula do efeito aleatório. Isso é restrição necessária para efeitos fixos. Os efeitos aleatórios não tem isso, mesmo porque pela própria definição eles tem esperança zero. Eles são preditos individualmente, não tem restrição paramétrica para formação das matrizes (pelo menos o usuário não pode controlar isso). Sendo assim, acho que deveria tentar sem o "-1". E veja sobre o ":" também.

À 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.