<div><span style="font-size: large; color: #0000ff;">Não consegui instalar o pacote minipack.lm, diz que não existe na versão R3.4.3. Estou usando o Windows, por acaso é só no Linux que êle existe?</span></div>
<div> </div>
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin: 0 0 0 .8ex; border-left: 1px #ccc solid; padding-left: 1ex;">
<div class="gmail_extra">
<div class="gmail_quote">
<blockquote class="m_4719827550724991136quote" style="margin: 0 0 0 .8ex; border-left: 1px #ccc solid; padding-left: 1ex;">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr"> </div>
<div dir="ltr"> </div>
<div dir="ltr">Boa tarde pessoal,
<div> </div>
<div>Estou querendo fazer a comparação de duas curvas ajustadas utilizando a função "nlsLM" do pacote "minpack.lm". No entanto, não consigo ajustar as duas curvas em um único nlsLM, para depois usar a anova para realizar um teste da razão de verossimilhança para testar a restrição em A, B e C. Alguém pode me ajudar a solucionar esse problema. Desde já agradeço muito a todas as sugestões. Além disso, preciso usar esse pacote por causa do algoritmo de Levenberg-Marquardt<span style="color: #000000; font-family: 'Times New Roman'; font-size: medium;">.</span>  Segue o script reproduzível. </div>
<div> </div>
<div>
<div>if(!require("minpack.lm")){install.packages("minpack.lm",dependencies=T);require("minpack.lm")}</div>
<div>if(!require("lattice")){install.packages("lattice",dependencies=T);require("lattice")}</div>
<div> </div>
<div>d <- expand.grid(parA=c(1,1.99,3,4,5,6,6), sistema = factor(c("TRAT1","TRAT2")))</div>
<div>d$diversidade <- c(55.000, 19.572, 10.231, 7.225, 6.004, 5.4022, 5.0566,</div>
<div>                   62.000, 31.6952, 21.4178, 17.6125, 15.7654, 14.7072, 14.0334)</div>
<div> </div>
<div>#ajusta as duas curvas em um unico nlsLM</div>
<div>n0 <- nlsLM(diversidade ~ A[sistema]/(1+B[sistema]*exp(-C[sistema]*parA)),</div>
<div>            data = d, start = list(A = c(-2.790,11.106), B = c(-1.050,-0.821),</div>
<div>                                       C = c(-0.090, 0.248))) </div>
<div>#Mensagem de erro:</div>
<div># Error in dimnames(x) <- dn : </div>
<div># length of 'dimnames' [2] not equal to array extent</div>
<div># In addition: Warning message:</div>
<div># In matrix(out$hessian, nrow = length(unlist(par))) :</div>
<div># data length [9] is not a sub-multiple or multiple of the number of rows [6]</div>
<div> </div>
<div> </div>
<div>#ajusta o modelo em que A, B e C seja comum aos dois tratamentos</div>
<div>n1 <- nlsLM(diversidade ~ A/(1+B*exp(-C*parA)),</div>
<div>            data = d, start = list(A = 4.158, B = -0.9355, C = 0.079)) </div>
<div> </div>
<div>anova(n0,n1)</div>
</div>
<div> </div>
<div><br />
<div> </div>
-- <br />
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr">
<div>
<div>
<div style="text-align: center;">Diego dos Santos Vieira</div>
</div>
</div>
<div style="text-align: center;">Engenheiro Florestal - UFRA</div>
<div style="text-align: center;">Mestre em Ciência Florestal - UFVJM</div>
<div style="text-align: center;">Doutorando em Ciência Florestal - UFVJM</div>
<div style="text-align: center;"> </div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</blockquote>
</div>
<div> </div>
<div> </div>
<div>-- </div>
<div class="gmail_signature" dir="ltr" data-smartmail="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr">
<div>
<div>
<div style="text-align: center;">Diego dos Santos Vieira</div>
</div>
</div>
<div style="text-align: center;">Engenheiro Florestal - UFRA</div>
<div style="text-align: center;">Mestre em Ciência Florestal - UFVJM</div>
<div style="text-align: center;">Doutorando em Ciência Florestal - UFVJM</div>
<div style="text-align: center;"> </div>
</div>
</div>
</div>
</div>
</div>
<div>_______________________________________________<br />R-br mailing list<br />R-br@listas.c3sl.ufpr.br<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.</div>