<div dir="ltr">Obrigado Simone pela ajuda. Era exatamente o que estava procurando.</div><br><div class="gmail_quote"><div dir="ltr">Em seg, 1 de out de 2018 às 19:43, <<a href="mailto:sisartorio@yahoo.com.br">sisartorio@yahoo.com.br</a>> escreveu:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto">Diego, boa noite!<div dir="auto"><br></div><div dir="auto">Veja minha tese...la tem como vc fazer isso...</div><div dir="auto"><br></div><div dir="auto">Banco de teses da usp/ simone daniela sartorio</div><div dir="auto"><br></div><div dir="auto">Att.</div><div dir="auto">Simone</div></div><div class="gmail_extra"><br><div class="gmail_quote">Em 29 de set de 2018 13:49, Diego Vieira via R-br <<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>> escreveu:<br type="attribution"><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">Boa tarde pessoal,<div><br></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<font color="#000000" face="Times New Roman" size="3">.</font> Segue o script reproduzível. </div><div><br></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><br></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><br></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><br></div><div><br></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><br></div><div>anova(n0,n1)</div></div><div><br></div><div><br clear="all"><div><br></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<br></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"><br></div></div></div></div></div></div></div></div></div></div>
</blockquote></div><br></div></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" 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<br></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"><br></div></div></div></div></div></div>