[R-br] Dúvida Inferência Bayesiana

[Vinicius Brito Rocha]

Vitor,

vc pode gerar de uma Gamma comum e inverter, depois pode olhar os quantis.

Att.

Em 23 de maio de 2011 21:03, Greice Laureano
<greice.laureano em gmail.com>escreveu:

> Oi Victor,
>
> Pode usar se inspirar no exemplo abaixo:
>
> library(TeachingDemos)
> x=c(62,58,61,60,51,61, 66, 63, 71, 59, 68, 50, 62, 62, 58, 54, 56, 58, 49,
> 65)
> p <- 0.95
> s2=var(x)
> gl=length(x)-1
>
> pinvchisq <- function(x, df=gl, scale=s2) {
>     prob <- pchisq(scale/x, df)
>     return(1-prob)
> }
>
>
> qinvchisq <- function(p, df=gl, scale) {
>     qprob <- qchisq(1-p, df)
>     return(scale/qprob)
> }
>
> hpd(qinvchisq,df=gl,scale=s2,conf=.95)
>
> Att,
> Greice.
>
> 2011/5/23 Benilton Carvalho <beniltoncarvalho em gmail.com>
>
>> teoria de probabilidade: se X eh Gamma(a, b), entao Y=1/X eh Gamma(a,
>> 1/b)...
>>
>> 2011/5/24 Victor Eduardo <victorduca08 em gmail.com>:
>> > Galera, qual comando que me da a distribuição gama inversa? Como que
>> calculo
>> > intervalo de credibilidade?
>> >
>> >
>> > Abraços,
>> >
>> > Victor Eduardo
>> >
>> > _______________________________________________
>> > R-br mailing list
>> > R-br em listas.c3sl.ufpr.br
>> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> >
>> >
>>
>>
>>
>> --
>> Successful people ask better questions, and as a result, they get
>> better answers. (Tony Robbins)
>> _______________________________________________
>> R-br mailing list
>> R-br em listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>


-- 
Vinicius Brito Rocha.
Estatístico e Atuário,
Mestre em Pesquisa Operacional

www.aplicademic.blogspot.com
http://twitter.com/viniciusbritor

"Não se preocupe muito com as suas dificuldades em Matemática, posso
assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110523/6fbb68af/attachment.html>