[R-br] Problema com optim() em modelos não lineares

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

[R-br] Problema com optim() em modelos não lineares

R-br mailing list

Prezados Membros,

    Estou tentando realizar a estimação dos coeficientes de um modelo não linear com a função optim(), mas só recebo a mensagem de erro:

Error in optim(c(0.97, 0.32), function(a) f(a, x, yexp)) : 
  função alvo em optim retorna um objeto de comprimento 11 ao invés de 1

   Segundo o meu CRM:

#Crio as variáveis

x <- c(1 ,10,  20,  30,  40,  50,  60,  70,  80,  90, 100)
y <- c(0.033823,  0.014779,  0.004698,  0.001584, -0.002017, -0.003436,
-0.000006, -0.004626, -0.004626, -0.004626, -0.004626)

#Especifico a função, que no meu caso é y ~ a4*exp(-x/a5)

f = function(a,x,yexp) {yexp= a[1]*exp(-x/a[2]) }

# Utilizando a otimização

mod = optim(c(0.97, 0.32), function(a) f(a, x, yexp))


Alguém poderia me ajudar.

Obrigado


-- 
======================================================================
Alexandre dos Santos
Proteção Florestal 
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial 
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO)
[hidden email] 
        [hidden email] 
Lattes: http://lattes.cnpq.br/1360403201088680 
OrcID: orcid.org/0000-0001-8232-6722   -   ResearcherID: A-5790-2016
Researchgate: www.researchgate.net/profile/Alexandre_Santos10                       
LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635
Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/
======================================================================

_______________________________________________
R-br mailing list
[hidden email]
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.