[R-br] matematica com o R

[Benilton Carvalho]

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
<samukajm em yahoo.com.br>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
> R-br em listas.c3sl.ufpr.br
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130912/141d8bd7/attachment.html>