[R-br] Modelo não linear

[Fernando Antonio de souza]

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 <
tiagosouzamarcal em hotmail.com> 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: ribas.aca em gmail.com
> To: r-br em listas.c3sl.ufpr.br
> 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 <
> tiagosouzamarcal em hotmail.com> 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
> 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.
>
>
>
>
> --
> 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
> 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/20130731/78752c12/attachment.html>