Oi Victor,<div><br></div><div>Pode usar se inspirar no exemplo abaixo:</div><div><br></div><div>library(TeachingDemos)</div><div><div>x=c(62,58,61,60,51,61, 66, 63, 71, 59, 68, 50, 62, 62, 58, 54, 56, 58, 49, 65)</div><div>
p <- 0.95</div><div>s2=var(x)</div><div>gl=length(x)-1</div><div><br></div><div>pinvchisq <- function(x, df=gl, scale=s2) {</div><div> prob <- pchisq(scale/x, df)</div><div> return(1-prob)</div><div>}</div>
<div><br></div><div><br></div><div>qinvchisq <- function(p, df=gl, scale) {</div><div> qprob <- qchisq(1-p, df)</div><div> return(scale/qprob)</div><div>}</div><div><br></div><div>hpd(qinvchisq,df=gl,scale=s2,conf=.95)</div>
<div><br></div><div>Att,</div><div>Greice.</div><br><div class="gmail_quote">2011/5/23 Benilton Carvalho <span dir="ltr"><<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
teoria de probabilidade: se X eh Gamma(a, b), entao Y=1/X eh Gamma(a, 1/b)...<br>
<br>
2011/5/24 Victor Eduardo <<a href="mailto:victorduca08@gmail.com">victorduca08@gmail.com</a>>:<br>
<div><div></div><div class="h5">> Galera, qual comando que me da a distribuição gama inversa? Como que calculo<br>
> intervalo de credibilidade?<br>
><br>
><br>
> Abraços,<br>
><br>
> Victor Eduardo<br>
><br>
</div></div>> _______________________________________________<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>
><br>
><br>
<font color="#888888"><br>
<br>
<br>
--<br>
Successful people ask better questions, and as a result, they get<br>
better answers. (Tony Robbins)<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>
</font></blockquote></div><br></div>