[R-br] VarCorr - nlme

[Fernando Antonio de souza]

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 em 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 em gmail.com>
> *Para:* r-br em 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 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.
>
>
> _______________________________________________
> 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/20120420/8630159a/attachment.html>