<div dir="ltr">Olá Alessandro e demais colegas.<div><br></div><div>Você pode usar o optim para encontrar a solução aproximada. Acredito que haja uma solução mais simplificada, mas essa resolve o problema e pode ser generalizada para outras situações:</div><div><br></div><div><br></div><div><div># Define os valores dos parametros que voce tem como conhecidos</div><div>b0=1</div><div>b1=1</div><div>b2=1</div><div>b3=1</div><div>b4=1</div><div>b5=1</div><div>s=1</div><div>id1=1</div><div>ab1=1</div><div><br></div><div># Define expressao como funcao de x</div><div>der=function(x) {</div><div> exp(b0 + b1/s + b2/x + b3 * (id1/x) * log(ab1) + b4 * (1 - id1/x) +b5 * (1 - id1/x) * s) * (b4 * (id1/x^2) - (b3 * (id1/x^2) *log(ab1) + b2/x^2) + b5 * (id1/x^2) * s)/x - exp(b0 + b1/s +b2/x + b3 * (id1/x) * log(ab1) + b4 * (1 - id1/x) + b5 *(1 - id1/x) * s)/x^2</div><div>}</div><div><br></div><div># Caso queira fazer o gráfico</div><div>curve(der(x))</div><div><br></div><div># Função para otimizar. 1 é o valor inicial, e o fnscale=-1 serve para ele maximizar</div><div>optim(1, der, method="BFGS", control=list(fnscale=-1))</div></div><div><br></div><div><br></div><div><div>$par</div><div>[1] <span style="background-color:rgb(255,255,0)">0.2928928</span></div><div><br></div><div>$value</div><div>[1] 50.55462</div><div><br></div><div>$counts</div><div>function gradient </div><div> 29 8 </div><div><br></div><div>$convergence</div><div>[1] 0</div><div><br></div><div>$message</div><div>NULL</div></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><b><br></b></div><div dir="ltr"><b>Paulo Dick</b><div>Estatístico</div><div>Mestrando em Epidemiologia em Saúde Pública</div><div>Tel.: (55 21) 99591-2716</div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">Em 28 de abril de 2015 23:35, <span dir="ltr"><<a href="mailto:aalmeidaalessandro@gmail.com" target="_blank">aalmeidaalessandro@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div dir="ltr" style="font-family:'Calibri','Segoe UI','Meiryo','Microsoft YaHei UI','Microsoft JhengHei UI','Malgun Gothic','sans-serif';font-size:12pt"><div>Como posso resolver essa equação sendo que sei os valores dos betas, id1,ab1 e s, quero achar o valor de x. (Não sei se é função solve ou outra função).</div><div><div><div>der= expression(exp(b0
+ b1/s + b2/x + b3 * (id1/x) * log(ab1) + b4 * (1 - id1/x) +b5 * (1 -
id1/x) * s) * (b4 * (id1/x^2) - (b3 * (id1/x^2) *log(ab1) + b2/x^2) + b5
* (id1/x^2) * s)/x - exp(b0 + b1/s +b2/x + b3 * (id1/x) * log(ab1) + b4
* (1 - id1/x) + b5 *(1 - id1/x) * s)/x^2)<br></div><div><div><br></div><div>Desde já agradeço!</div></div><br></div></div></div>
</div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>