<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span><div style="font-size: 12pt; font-family: times, serif; ">Obrigada pela ajuda, Walmes e Fernando.</div><div style="font-size: 12pt; font-family: times, serif; "><br></div><div style="font-size: 12pt; font-family: times, serif; ">O meu modelo esse mesmo Fernando. Eu estou utilizando o mesmo modelo d<span style="font-size: 12pt; ">o livro do Bates.</span></div><div style="font-size: 12pt; font-family: times, serif; "><span>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).</span></div><div style="font-size: 12pt; font-family: times, serif; "><br></div><div style="font-size: 12pt; font-family: times, serif; ">O problema está no comando "VarCorr". Eu consigo fazer "Summary(Model2), mas VarCorr(Model2) não.</div><div
style="font-size: 12pt; font-family: times, serif; "><br></div><div style="font-size: 12pt; font-family: times, serif; ">Eu tentei rodar o exemplo do help do R e deu a mesma mensagem de erro que dá quando eu rodo meu modelo.</div><div><pre style="font-family: times, serif; font-size: 12pt; ">fm1 <- lme(distance ~ age, data = Orthodont, random = ~age)
VarCorr(fm1)</pre><pre style="font-family: 'Courier New'; ">Erro em function (classes, fdef, mtable) :
unable to find an inherited method for function "VarCorr", for signature "lme"<font face="'times new roman', 'new york', times, serif" size="4">
</font></pre><div style="font-family: times, serif; font-size: 12pt; ">Será que o problema está na minha versão do R? </div><div style="font-family: times, serif; font-size: 12pt; ">Eu estou trabalhando com a versão 2.15.0 (2012-03-30).</div><div style="font-family: times, serif; font-size: 12pt; ">Isso tem alguma coisa a ver?</div><div style="font-family: times, serif; font-size: 12pt; "><br></div><div style="font-family: times, serif; font-size: 12pt; ">Obrigada pela ajuda</div><div style="font-family: times, serif; font-size: 12pt; ">Isabel</div></div><div style="font-size: 12pt; font-family: times, serif; "><br></div><div style="font-size: 12pt; font-family: times, serif; ">Minha função:</div><div><div>attrnames<- names(panel)[4:(L+3)]</div><div>compvar=list()</div><div>for (j in 1:L){</div><div>const<-rep(1,
length(panel2$atr1))</div><div>panel2$prod<-factor(panel2$prod)</div><div>panel2$ass<-factor(panel2$ass)</div><div>form2<- formula(paste(attrnames[j],"~prod"))</div><div>model2<-lme(form2, random=list(const=pdBlocked(list(~ass-1, ~ass:prod-1),pdClass="pdIdent")),data=panel2) </div><div>compvar<-c(compvar, list(unique(VarCorr(model2)[,1])))</div><div>} </div></div></span></div><div><br></div> <div style="font-size: 12pt; font-family: 'times new roman', 'new york', times, serif; "> <div style="font-size: 12pt; font-family: 'times new roman', 'new york', times, serif; "> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1"> <b><span style="font-weight:bold;">De:</span></b> Walmes Zeviani <walmeszeviani@gmail.com><br> <b><span style="font-weight: bold;">Para:</span></b> r-br@listas.c3sl.ufpr.br <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Sexta-feira, 20 de Abril de 2012 8:06<br>
<b><span style="font-weight: bold;">Assunto:</span></b> Re: [R-br] VarCorr - nlme<br> </font> </div> <br><div id="yiv824103454"><div><font face="trebuchet ms,sans-serif">Isabel,</font></div><div><font face="trebuchet ms,sans-serif"><br></font></div><div><font face="trebuchet ms,sans-serif">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.</font></div>
<div><font face="trebuchet ms,sans-serif"><br></font></div><div><font face="trebuchet ms,sans-serif">À disposição.</font></div><div><font face="trebuchet ms,sans-serif">Walmes.</font></div><font face="trebuchet ms,sans-serif"><br clear="all">
</font><span style="font-family: 'trebuchet ms', sans-serif; ">==========================================================================</span><br style="font-family: 'trebuchet ms', sans-serif; "><span style="font-family: 'trebuchet ms', sans-serif; ">Walmes Marques Zeviani</span><br style="font-family: 'trebuchet ms', sans-serif; ">
<span style="font-family: 'trebuchet ms', sans-serif; ">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family: 'trebuchet ms', sans-serif; "><span style="font-family: 'trebuchet ms', sans-serif; ">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family: 'trebuchet ms', sans-serif; ">
<span style="font-family: 'trebuchet ms', sans-serif; ">fone: (+55) 41 3361 3573</span><br style="font-family: 'trebuchet ms', sans-serif; "><span style="font-family: 'trebuchet ms', sans-serif; ">VoIP: (3361 3600) 1053 1173</span><br style="font-family: 'trebuchet ms', sans-serif; ">
<span style="font-family: 'trebuchet ms', sans-serif; ">e-mail: <a rel="nofollow" ymailto="mailto:walmes@ufpr.br" target="_blank" href="mailto:walmes@ufpr.br">walmes@ufpr.br</a></span><br style="font-family: 'trebuchet ms', sans-serif; "><span style="font-family: 'trebuchet ms', sans-serif; ">twitter: @walmeszeviani</span><br style="font-family: 'trebuchet ms', sans-serif; ">
<span style="font-family: 'trebuchet ms', sans-serif; ">homepage: <a rel="nofollow" target="_blank" href="http://www.leg.ufpr.br/%7Ewalmes">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family: 'trebuchet ms', sans-serif; "><span style="font-family: 'trebuchet ms', sans-serif; ">linux user number: 531218</span><br style="font-family: 'trebuchet ms', sans-serif; ">
<span style="font-family: 'trebuchet ms', sans-serif; ">==========================================================================</span><br>
</div><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br><br> </div> </div> </div></body></html>