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@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@gmail.com To: r-br@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@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@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@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@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.

Por curiosidade eu tentei todos os recursos que conheço para ajustar esses modelos aos dados passados no CMR e não tive sucesso. Você está seguro de que esse modelo é adequado, é capaz de dar a forma que os dados apresentam? Viu aplicações desse modelo em artigos dessa área? O modelo é uma soma de funções exponenciais e a literatura já aponta que são modelos que apresentam problemas de identificabilidade. Ao todo são 5 parâmetros e considerando que você tem 12 observações, talvez um polinômio com menos de 5 termos consiga um bom ajuste ou até um gam. A interpretação dos parâmetros desse modelo, de inicio considerando que não explorei com devido tempo, é meio vaga, então abordagens mais empíricas não devem ser descartadas. Nessa página www.leg.ufpr.br/mrnl2013 tem todo o material do Curso em Modelos de Regressão Não Linear oferecido durante 58º RBRAS. Um zip contém todos os scripts. Veja o script05-mgi.R. O rpanel não funciona dentro do RStudio. Use outro editor de scripts, ex o emacs, gedit, tinn-R, vi, etc. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================