Em 23 de maio de 2011 21:03, Greice Laureano <greice.laureano@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@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@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@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@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br