acho pouco provável que seja a versão do seu R. O meu é o 2.14.1 e rodou perfeitamente o código do exemplo. Dê um dput(dados) e envie-o junto com os seus comandos para que outros possam tentar rodar para você. Olhando assim fica difícil encontrar um erro. Em 20 de abril de 2012 14:21, isabel de sousa Amorim < isabelsousaamorim@yahoo.com.br> escreveu:
Obrigada pela ajuda, Walmes e Fernando.
O meu modelo esse mesmo Fernando. Eu estou utilizando o mesmo modelo do livro do Bates. Porém eu tenho vários atributos, então eu criei uma função para rodar todos os atributos de só vez. (segue a função abaixo da mensagem).
O problema está no comando "VarCorr". Eu consigo fazer "Summary(Model2), mas VarCorr(Model2) não.
Eu tentei rodar o exemplo do help do R e deu a mesma mensagem de erro que dá quando eu rodo meu modelo.
fm1 <- lme(distance ~ age, data = Orthodont, random = ~age) VarCorr(fm1)
Erro em function (classes, fdef, mtable) : unable to find an inherited method for function "VarCorr", for signature "lme"
Será que o problema está na minha versão do R? Eu estou trabalhando com a versão 2.15.0 (2012-03-30). Isso tem alguma coisa a ver?
Obrigada pela ajuda Isabel
Minha função: attrnames<- names(panel)[4:(L+3)] compvar=list() for (j in 1:L){ const<-rep(1, length(panel2$atr1)) panel2$prod<-factor(panel2$prod) panel2$ass<-factor(panel2$ass) form2<- formula(paste(attrnames[j],"~prod")) model2<-lme(form2, random=list(const=pdBlocked(list(~ass-1, ~ass:prod-1),pdClass="pdIdent")),data=panel2) compvar<-c(compvar, list(unique(VarCorr(model2)[,1]))) }
------------------------------ *De:* Walmes Zeviani <walmeszeviani@gmail.com> *Para:* r-br@listas.c3sl.ufpr.br *Enviadas:* Sexta-feira, 20 de Abril de 2012 8:06
*Assunto:* Re: [R-br] VarCorr - nlme
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 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.
_______________________________________________ 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.