[R-br] RES: Função implícita

[Eder David Borges da Silva]

Davi,
Agora sim, acho que pode ser resolvido assim:
MU <- 0.3
f <- function(p){exp(-p) * (1 - p) - exp( 1 - (1/MU))}

require(rootSolve)
#?rootSolve
p <- uniroot.all(f,lower=0,upper=1)
p
curve(f,0,1)
abline(h=0,lty=2)
abline(v=p,lty=2)

Att

Em 1 de julho de 2012 17:07, Davi <davi.butturi em gmail.com> escreveu:

> Na verdade estava querendo reparametrizar a distribuição logarítmica pela
> sua média.
>
> Então se X é uma v.a. com função de probabilidade logarítimica então E(X) =
> - ( 1 / log (1-p) ) * ( p / (1-p) )
>
> A função rlog(n,p) {VGAM} gera n realizações dessa distribuição, mas
> parametrizada pela probabilidade p.
> Assim, se E(X) é um valor MU qualquer conhecido, eu quero obter essa
> probabilidade.
>
> Resolvendo para p:
>
> exp( -p ) * (1 - p) = exp( 1 - (1/MU) )
>
> onde f(p) = exp( -p ) * (1 - p) é uma função implícita de p
> e g(MU)  = exp( 1 - (1/MU) ) uma função que pode ser avaliada pq conheço
> MU.
>
> Agora preciso tornar essa igualdade verdadeira fazendo f( p ) = g ( MU ), e
> obter p.
>
> Ficou mais claro?
>
> Obrigado! Abraços,
>
> Davi.
>
>
>
>
>
>
>
>
> -----Mensagem original-----
> De: r-br-bounces em listas.c3sl.ufpr.br
> [mailto:r-br-bounces em listas.c3sl.ufpr.br] Em nome de FHRB Toledo
> Enviada em: domingo, 1 de julho de 2012 11:59
> Para: r-br em listas.c3sl.ufpr.br
> Assunto: Re: [R-br] Função implícita
>
> Lagrange?
>
> 2012/7/1 Eder David Borges da Silva <eder em leg.ufpr.br>:
> > Davi,
> > Sem um CMR fica difícil.....Mas segue algumas funções que podem de dar
> uma
> > ajuda:
> > ?optim
> > ?rootSolve
> > Essa são genéricas  e numéricas , a sua função não tem solução analítica?
> > Att
> >
> >
> >
> > Em 1 de julho de 2012 00:43, Davi <davi.butturi em gmail.com> escreveu:
> >>
> >> Bom dia pessoal,
> >>
> >>
> >>
> >> O problema que estou é o seguinte: preciso encontrar uma (ou mais, não
> >> sei) soluções para uma equação do tipo
> >>
> >>
> >>
> >> f(x) - c =0
> >>
> >>
> >>
> >> onde f(x) é uma função implícita e c uma constante qualquer. Tem algum
> >> método?
> >>
> >>
> >>
> >> Muito obrigado,
> >>
> >>
> >>
> >> Davi
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >> ___________________________________________________
> >>
> >> Davi Butturi-Gomes
> >>
> >> Doutorando em Estatística e Experimentação Agronômica - ESALQ /USP -
> >> Piracicaba
> >>
> >> Mestre em Biometria - UNESP - Botucatu
> >>
> >> Ecólogo - UNESP - Rio Claro
> >>
> >> MSN: davi_bgomes em hotmail.com
> >>
> >> Skype: davibg87
> >>
> >> (19)9809.5753
> >>
> >>
> >>
> >>
> >> _______________________________________________
> >> 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.
> >
> >
> >
> > _______________________________________________
> > 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.
> _______________________________________________
> 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.
>
> _______________________________________________
> 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/20120701/89221332/attachment.html>