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

[Greice Laureano]

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
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110523/30362ced/attachment.html>