<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div># Gamma inversa </div><div><br></div><div>g <- 1/gamma(n,a,b)</div><div> </div><div>        Fábio Mathias Corrêa<br><br></div><div>     Universidade Estadual de Santa Cruz<br>Departamento de Ciências Exatas e da Terra - DCET</div><br><br><div>Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna<br>CEP 45662-900. Ilhéus-Bahia</div><div><br><br></div><div>Tel.: 73-3680-5076<br><div style="font-size: 12pt; font-family: 'times new roman', 'new york', times, serif; "><div style="font-size: 12pt; font-family: 'times new roman', 'new york', times, serif; "><font size="2" face="Arial"><hr size="1"><b><span style="font-weight:bold;">De:</span></b> Vinicius Brito Rocha <viniciusbritor@gmail.com><br><b><span style="font-weight: bold;">Para:</span></b>
 r-br@listas.c3sl.ufpr.br<br><b><span style="font-weight: bold;">Enviadas:</span></b> Segunda-feira, 23 de Maio de 2011 21:09<br><b><span style="font-weight: bold;">Assunto:</span></b> Re: [R-br] Dúvida Inferência Bayesiana<br></font><br><div id="yiv1577682507">Vitor,<br><br>vc pode gerar de uma Gamma comum e inverter, depois pode olhar os quantis.<br><br>Att.<br><br><div class="yiv1577682507gmail_quote">Em 23 de maio de 2011 21:03, Greice Laureano <span dir="ltr"><<a rel="nofollow" ymailto="mailto:greice.laureano@gmail.com" target="_blank" href="mailto:greice.laureano@gmail.com">greice.laureano@gmail.com</a>></span> escreveu:<br>
<blockquote class="yiv1577682507gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">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><div><div class="yiv1577682507h5"><br><div class="yiv1577682507gmail_quote">2011/5/23 Benilton Carvalho <span dir="ltr"><<a rel="nofollow" ymailto="mailto:beniltoncarvalho@gmail.com" target="_blank" href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>></span><br>
<blockquote class="yiv1577682507gmail_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 rel="nofollow" ymailto="mailto:victorduca08@gmail.com" target="_blank" href="mailto:victorduca08@gmail.com">victorduca08@gmail.com</a>>:<br>
<div><div>> 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 rel="nofollow" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> <a rel="nofollow" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><span class="yiv1577682507wrc0" style="padding-right:16px;width:16px;height:16px;"></span><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 rel="nofollow" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a rel="nofollow" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><span class="yiv1577682507wrc0" style="padding-right:16px;width:16px;height:16px;"></span><br>

</font></blockquote></div><br></div></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a rel="nofollow" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a rel="nofollow" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><span class="yiv1577682507wrc0" style="padding-right:16px;width:16px;height:16px;"></span><br>

<br></blockquote></div><br><br clear="all"><br>-- <br>Vinicius Brito Rocha.<br>Estatístico e Atuário,<br>Mestre em Pesquisa Operacional<br><br><a rel="nofollow" target="_blank" href="http://www.aplicademic.blogspot.com">www.aplicademic.blogspot.com</a><span class="yiv1577682507wrc11" style="padding-right:16px;width:16px;height:16px;"></span><br>
<a rel="nofollow" target="_blank" href="http://twitter.com/viniciusbritor">http://twitter.com/viniciusbritor</a><span class="yiv1577682507wrc13" style="padding-right:16px;width:16px;height:16px;"></span><br>
<br>"Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.<br><br><br>
</div><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" 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></div></div></div></div></body></html>