[R-br] Modelo não linear

[Tiago Souza Marçal]

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/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. 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130731/29a86b4e/attachment.html>