[R-br] Problema com optim() em modelos não lineares
ASANTOS
alexandre.santos em cas.ifmt.edu.br
Sex Mar 1 12:03:55 -03 2019
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)
e-mails:alexandresantosbr em yahoo.com.br
alexandre.santos em cas.ifmt.edu.br
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/
======================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190301/7a0d1c14/attachment.html>
Mais detalhes sobre a lista de discussão R-br