Modelagem utilizando a função nls() 

# vamos fazer um exemplo simples com os dados simulados abaixo

 

x = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10)

y = c(5, 11, 13, 22, 35, 40, 65, 80, 123, 175)

 

# modelo: y=a*eb*x

 

 

# utilizando a função nls() ajustamos o modelo com dois starts a=1 e b=1

 

modelo=nls(y~a*exp(b*x) ,start=c(a=1,b=1))

 

# para ver o modelo ajustado

 

modelo

 

# significância dos coeficientes

 

summary(modelo)

 

 

# criando uma função a partir dos parâmetros estimados

 

f=function(x){5.12*exp(0.35*x)}

 

 

# criando um gráfico com os pontos e inserindo a equação (função f criada acima)

 

plot(x,y)

curve(f, 1, 10, add=T,col=2, lty=3)

#1=limite inferior, 10=limite superior

 

 

# veja abaixo como ficou o gráfico