Olá Thiago Porque você quer outra função para gera chutes iniciais. O manipulate não lhe atende não? No blog do Walmes (Ridículas) na página inicial ele fornece um código para utilizar a função rpanel . De uma olhada lá e veja se é isso que você quer. http://ridiculas.wordpress.com/
Em 31 de julho de 2013 06:31, Tiago Souza Marçal < [email protected]> escreveu: > Obrigado pela dica Augusto, > > eu conheço a nls(), entretanto, o problema são os chutes iniciais em > start=list(), pois, um ajuste pode existir mas dependendo dos chutes > iniciais não se consegue a convergência. Com outros modelos que tenho > trabalhado a linearização tem dado certo para estimar os chutes iniciais, > mas não consegui linearizar este modelo. > > > Att. > > Tiago. > > ################################################################# > > Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES > > Bolsista de Iniciação Científica da área de Genética e Melhoramento de > Plantas > > ################################################################# > > > ------------------------------ > Date: Wed, 31 Jul 2013 00:23:26 -0400 > From: [email protected] > To: [email protected] > Subject: Re: [R-br] Modelo não linear > > > Opa, eu não entendo muito bem de regressão não linear, mas talvez a função > nls possa fazer o que você tem interesse. > > Tem um livro muito bom da serie use R explicando como usar, vale a pena a > consulta caso a função nls se encaixe nas suas necessidades. > http://www.amazon.com/Nonlinear-Regression-R-Use/dp/0387096159 > > Fora isso tem muitos manuais pela internet também. > Mas mais ou menos é isso que você quer, ajustar esse parâmetros aos dados? > > > ######################################################################################## > #exemplo > tempo=1:12 > > y=c(8.81790512674818,8.83724295007473,8.84652971800299,8.85625726150244,8.86597721721444,8.86647798384983, > > 8.87625601672093,8.89503521901844,8.91478268591035,8.93378059512321,8.97371264813241,9.00355384268932) > > plot(y~tempo) > > modelo<-nls(y~A*(1-B*exp(-C*tempo)) + exp(D + > E*(tempo-tempo[1])),start=list(A=1,B=1,C=1,D=1,E=1)) > summary(modelo) > > curve(coef(modelo)[1]*(1-coef(modelo)[2]*exp(-coef(modelo)[3]*x)) + > exp(coef(modelo)[4] + coef(modelo)[5]*(x-1)),1,12, > add=T,lty=2,col="red") > > > ########################################################################################### > Mas é um bocado de parâmetros para 12 pontos. > Bem, espero ter ajudado. > > > > > Em 30 de julho de 2013 19:39, Tiago Souza Marçal < > [email protected]> escreveu: > > Boa noite pessoal, > > estou querendo modelar uma curva de embebição de sementes através do > seguinte modelo não linear: > > y=A*(1-B*exp(-C*tempo)) + exp(D + E*(tempo-t0)) > > t0 = 6 > > Gostaria de saber se existe um procedimento parecido com o do pacote > manipulate do Rstudio no R para ajustar este modelo. > > CMR > > tempo=1:12 > y=c(8.81790512674818,8.83724295007473,8.84652971800299,8.85625726150244, > > 8.86597721721444,8.86647798384983,8.87625601672093,8.89503521901844,8.91478268591035,8.93378059512321,8.97371264813241,9.00355384268932) > > Att. > > Tiago. > > > > > > > > > > ################################################################# > > Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES > > Bolsista de Iniciação Científica da área de Genética e Melhoramento de > Plantas > > ################################################################# > > _______________________________________________ > R-br mailing list > [email protected] > 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. > > > > > -- > Grato > Augusto C. A. Ribas > > Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> > Github: https://github.com/Squiercg > Lattes: http://lattes.cnpq.br/7355685961127056 > > _______________________________________________ R-br mailing list > [email protected] > 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 > [email protected] > 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 [email protected] 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.
