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
