<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>Samuel,<BR> <BR>não utilizo a função optim portanto, não sei o procedimento que ela usa. Uma maneira seria as tentativas baseadas em aproximações gráficas embora ao se igualar a segunda derivada a zero a solução seria analítica (x = -b^2/2*b).<BR> <BR>Outra sugestão é você dar uma olhada na função while() e organizar um código que procure o valor de x que zere a segunda derivada em um intervalo de -0.01 > y'' > 0.01 <BR> <BR>Ex:<BR> <BR>x=seq(0.1,2,0.01)<BR> <BR>fx=function(x){exp(3.3-1.2/x)}<BR> <BR>fx2 <- function(x) exp(3.3 - 1.2/x) * (-1.2 * (2 * x)/(x^2)^2) + exp(3.3 - 1.2/x) * (-1.2/x^2) * (-1.2/x^2)<br><BR>y=fx(x)<BR> <BR>y2=fx2(3.3,-1.2,x)<BR> <BR>plot(y~x)<BR> <BR>par (new=T)<BR> <BR> plot(y2~x,axes=F)<BR> <BR>axis(4)<BR> <BR>abline(h=0)<BR> <BR>abline(v=0.6)<BR> <BR>Att.<BR> <BR>Tiago.<BR><br><br> <BR><div>#################################################################</div><div> </div><div>Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES</div><div> </div><div>Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas</div><div> </div><div>################################################################# </div><br> <BR><div><hr id="stopSpelling">Date: Thu, 12 Sep 2013 09:17:39 -0700<br>From: samukajm@yahoo.com.br<br>To: r-br@listas.c3sl.ufpr.br<br>Subject: [R-br] matematica com o R<br><br><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 12pt; background-color: rgb(255, 255, 255);"><div style='font-family: "times new roman", "new york", times, serif; font-size: 12pt;'><span>Caros, boa tarde!</span></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span>Alguém poderia por gentileza auxiliar com a duvida relativa ao codigo abaixo?</span></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span><br></span></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span># RCMR</span></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span>y <- expression(exp(a+b/x))</span></div><div style="background-color: transparent;">D(y,'x') #primeira derivada</div><div style="background-color: transparent;">dy.dx <- expression(-(exp(a + b/x) * (b/x^2))) #resultado da primeira derivada<br></div><div style="background-color: transparent;">D(dy.dx, 'x') #segunda derivada</div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span style="background-color: transparent;">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</span></div><div style="background-color: transparent;"><span style="background-color: transparent;">optim(c(1), fx, a=3.3,b=-1.2,method="BFGS")<br></span></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span style="background-color: transparent;"><br></span></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span style="background-color: transparent;"># 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 </span></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span style="background-color: transparent;"># 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.</span></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span style="background-color: transparent;"># 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 </span></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span style="background-color: transparent;"># 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 </span></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span style="background-color: transparent;"># R e nao </span><span style="background-color: transparent;">mais com o maple. Segue um exemplo simples de um calculo que gostaria de reproduzir com o R</span></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span style="background-color: transparent;"><br></span></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'>y = b0 + b1x</div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'>b0=?</div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'>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.</div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><br></div><div style='color: rgb(0, 0, 0); font-family: "times new roman", "new york", times, serif; font-size: 16px; font-style: normal; background-color: transparent;'><span style="font-size: 12pt;">Ja deixo aqui meus agradecimentos</span><br></div><div style='font-family: "times new roman", "new york", times, serif; font-size: 12pt;'><br></div><div style="font-family: courier, monaco, monospace, sans-serif; font-size: 12pt;"><font size="3"><b>====================================</b></font></div><div style="font-family: courier, monaco, monospace, sans-serif; font-size: 12pt;"><font size="3"><b>Samuel P. C. Carvalho<br></b><font color="#000000" size="2">Engº Florestal [UFLA]</font></font></div><div style="color: rgb(0, 0, 0); font-family: courier, monaco, monospace, sans-serif; font-size: 16px; font-style: normal; background-color: transparent;"><font size="3"><font color="#000000" size="2">Mestre em Ciências Florestais [UFLA]</font></font></div><div style="font-family: courier, monaco, monospace, sans-serif; font-size: 12pt;"><font color="#000000" size="2">Doutor em Recursos Florestais [ESALQ/USP]</font></div><div style="font-family: courier, monaco, monospace, sans-serif; font-size: 12pt;"><font size="2">=============================================</font></div></div><br>_______________________________________________
R-br mailing list
R-br@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.</div> </div></body>
</html>