a sacada e' notar que optim() minimiza (ou maximiza, dependendo dos argumentos) a funcao objet... nao existem garantias de que sua fx assuma apenas valores nao negativos, entao vc tem q transforma-la numa funcao nao-negativa (c.f. sugestao).
Vc pode, entao, minimizar "fx^2" (pois essa funcao e' nao-negativa (o que nao pode ser garantido para fx), sob a pena de erros maiores... ou, mais "acuradamente" minimizar abs(fx). Seguindo o seu exemplo: gx = function(...) abs(fx(...)) optim(c(1), gx, a=3.3,b=-1.2,method="BFGS") b Em 12 de setembro de 2013 13:17, Samuel Carvalho <[email protected]>escreveu: > Caros, boa tarde! > Alguém poderia por gentileza auxiliar com a duvida relativa ao codigo > abaixo? > > # RCMR > y <- expression(exp(a+b/x)) > D(y,'x') #primeira derivada > dy.dx <- expression(-(exp(a + b/x) * (b/x^2))) #resultado da primeira > derivada > D(dy.dx, 'x') #segunda derivada > fx <- function(a,b,x) exp(a + b/x) * (b * (2 * x)/(x^2)^2) + exp(a + b/x) > * (b/x^2) * (b/x^2) #resultado da segunda derivada > optim(c(1), fx, a=3.3,b=-1.2,method="BFGS") > > # A ideia é encontrar o ponto de inflexao para o modelo representado pelo > objeto y. Por definição o ponto de inflexao é o valor de x que faz a > segunda derivada da > # função igual a zero. Montei este rcmr porém não estou certo se a > otimização de fx está igualando a zero, até porque os resultados não estão > condizentes. > # Gostaria se possível também de alguma dica de algum livro que tenha > dicas de cálculos matemáticos com o R pois gostaria de adaptar algumas > rotinas do maple > # para o R, até porque recentemente precisei isolar um termo de um modelo > não linear que manualmente não é tão simples e gostaria de fazer estes > calculos com o > # R e nao mais com o maple. Segue um exemplo simples de um calculo que > gostaria de reproduzir com o R > > y = b0 + b1x > b0=? > b0 = y - b1x # esta é a saída que gostaria de encontrar com o R. Ate > cheguei a ver um pouco sobre o pacote Ryacas mas sem muito sucesso. > > Ja deixo aqui meus agradecimentos > > *====================================* > *Samuel P. C. Carvalho > *Engº Florestal [UFLA] > Mestre em Ciências Florestais [UFLA] > Doutor em Recursos Florestais [ESALQ/USP] > ============================================= > > _______________________________________________ > R-br mailing list > [email protected] > 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. >
_______________________________________________ R-br mailing list [email protected] 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.
