Davi,<div>Agora sim, acho que pode ser resolvido assim:</div><div><div>MU <- 0.3</div><div>f <- function(p){exp(-p) * (1 - p) - exp( 1 - (1/MU))}</div><div><br></div><div>require(rootSolve)</div><div>#?rootSolve</div>
<div>p <- uniroot.all(f,lower=0,upper=1)</div><div>p</div><div>curve(f,0,1)</div><div>abline(h=0,lty=2)</div><div>abline(v=p,lty=2)</div><div><br></div><div>Att</div><br><div class="gmail_quote">Em 1 de julho de 2012 17:07, Davi <span dir="ltr"><<a href="mailto:davi.butturi@gmail.com" target="_blank">davi.butturi@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Na verdade estava querendo reparametrizar a distribuição logarítmica pela<br>
sua média.<br>
<br>
Então se X é uma v.a. com função de probabilidade logarítimica então E(X) =<br>
- ( 1 / log (1-p) ) * ( p / (1-p) )<br>
<br>
A função rlog(n,p) {VGAM} gera n realizações dessa distribuição, mas<br>
parametrizada pela probabilidade p.<br>
Assim, se E(X) é um valor MU qualquer conhecido, eu quero obter essa<br>
probabilidade.<br>
<br>
Resolvendo para p:<br>
<br>
exp( -p ) * (1 - p) = exp( 1 - (1/MU) )<br>
<br>
onde f(p) = exp( -p ) * (1 - p) é uma função implícita de p<br>
e g(MU)  = exp( 1 - (1/MU) ) uma função que pode ser avaliada pq conheço MU.<br>
<br>
Agora preciso tornar essa igualdade verdadeira fazendo f( p ) = g ( MU ), e<br>
obter p.<br>
<br>
Ficou mais claro?<br>
<br>
Obrigado! Abraços,<br>
<br>
Davi.<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
-----Mensagem original-----<br>
De: <a href="mailto:r-br-bounces@listas.c3sl.ufpr.br">r-br-bounces@listas.c3sl.ufpr.br</a><br>
[mailto:<a href="mailto:r-br-bounces@listas.c3sl.ufpr.br">r-br-bounces@listas.c3sl.ufpr.br</a>] Em nome de FHRB Toledo<br>
Enviada em: domingo, 1 de julho de 2012 11:59<br>
Para: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
Assunto: Re: [R-br] Função implícita<br>
<div class="HOEnZb"><div class="h5"><br>
Lagrange?<br>
<br>
2012/7/1 Eder David Borges da Silva <<a href="mailto:eder@leg.ufpr.br">eder@leg.ufpr.br</a>>:<br>
> Davi,<br>
> Sem um CMR fica difícil.....Mas segue algumas funções que podem de dar uma<br>
> ajuda:<br>
> ?optim<br>
> ?rootSolve<br>
> Essa são genéricas  e numéricas , a sua função não tem solução analítica?<br>
> Att<br>
><br>
><br>
><br>
> Em 1 de julho de 2012 00:43, Davi <<a href="mailto:davi.butturi@gmail.com">davi.butturi@gmail.com</a>> escreveu:<br>
>><br>
>> Bom dia pessoal,<br>
>><br>
>><br>
>><br>
>> O problema que estou é o seguinte: preciso encontrar uma (ou mais, não<br>
>> sei) soluções para uma equação do tipo<br>
>><br>
>><br>
>><br>
>> f(x) - c =0<br>
>><br>
>><br>
>><br>
>> onde f(x) é uma função implícita e c uma constante qualquer. Tem algum<br>
>> método?<br>
>><br>
>><br>
>><br>
>> Muito obrigado,<br>
>><br>
>><br>
>><br>
>> Davi<br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>> ___________________________________________________<br>
>><br>
>> Davi Butturi-Gomes<br>
>><br>
>> Doutorando em Estatística e Experimentação Agronômica - ESALQ /USP -<br>
>> Piracicaba<br>
>><br>
>> Mestre em Biometria - UNESP - Botucatu<br>
>><br>
>> Ecólogo - UNESP - Rio Claro<br>
>><br>
>> MSN: <a href="mailto:davi_bgomes@hotmail.com">davi_bgomes@hotmail.com</a><br>
>><br>
>> Skype: davibg87<br>
>><br>
>> <a href="tel:%2819%299809.5753" value="+551998095753">(19)9809.5753</a><br>
>><br>
>><br>
>><br>
>><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<br>
>> código mínimo reproduzível.<br>
><br>
><br>
><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<br>
código<br>
> mínimo reproduzível.<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<br>
mínimo reproduzível.<br>
<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>
</div></div></blockquote></div><br></div>