Com respeito a função beta incompleta. Existe um pacote chamado zipfR. A função Igamma (gama imcompleta) bate com a resposta do mathematica assim como a função implementada pelo Benilton. No pacote a definição da gamma incompleta é:<br>
<br>B(x;a,b) = integral_0^x t^(a-1) (1-t)^(b-1) dt<br><br>A definição do mathematica é: <br><br>B(z;a,b) = integral_0^z u^(a-1) (1-u)^(b-1) du<br><br>Ver em:<br><a href="http://mathworld.wolfram.com/IncompleteBetaFunction.html">http://mathworld.wolfram.com/IncompleteBetaFunction.html</a><br>
<br>Concordam que são as mesmas integrais?<br><br>Resposta do R: Ibeta(3,1,2) = 0.5<br>Resposta d Mathematica: N[Beta[3, 1, 2], 5] = -1.5000<br><br>Seria um problema da função ou de documentação ou um problema meu ainda?<br>
<br><br><br><br>
<div class="gmail_extra"><br clear="all">[ ],<br>Pedro Rafael Diniz Marinho.<br>
<br><br><div class="gmail_quote">Em 3 de novembro de 2012 10:44, Pedro Rafael <span dir="ltr"><<a href="mailto:pedro.rafael.marinho@gmail.com" target="_blank">pedro.rafael.marinho@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Valeu Benilton, agora está funcionando. Obrigado.<div class="im"><br clear="all">[ ],<br>Pedro Rafael Diniz Marinho.<br>
<br><br></div><div class="gmail_quote">---------- Mensagem encaminhada ----------<br>De: <b class="gmail_sendername">beniltoncarvalho [via R-br]</b> <span dir="ltr"><<a href="mailto:ml-node%2Bs2285057n4656829h76@n4.nabble.com" target="_blank">ml-node+s2285057n4656829h76@n4.nabble.com</a>></span><br>
Data: 2 de novembro de 2012 14:21<br>Assunto: Re: [R-br] Fwd: [Dúvida] Função Gama Incompleta.<br>Para: Pedro Rafael Diniz Marinho <<a href="mailto:pedro.rafael.marinho@gmail.com" target="_blank">pedro.rafael.marinho@gmail.com</a>><div>
<div class="h5"><br>
<br><br>
o q complicou e' q a cada post vc modificou a pergunta... de qq forma, a recomendacao inicial continua valida (e' so' verificar a definicao da funcao)...<div><br></div><div>para sua conveniencia e exercicio de fim de semana, esta e' a solucao q vc quer:</div>
<div><br></div><div><div>gammainc <- function(x, a, b){</div><div><span style="white-space:pre-wrap"> </span>r1 <- gamma(x)</div><div><span style="white-space:pre-wrap"> </span>if (missing(b)){</div><div><span style="white-space:pre-wrap"> </span>if (isTRUE(all.equal(a, 0)))</div>
<div><span style="white-space:pre-wrap"> </span>return(r1)</div><div><span style="white-space:pre-wrap"> </span>return(r1*pgamma(a, x, lower=FALSE))</div><div><span style="white-space:pre-wrap"> </span>}</div>
<div><span style="white-space:pre-wrap"> </span>gammainc(x, a)-gammainc(x, b)</div><div>}</div></div><div><br></div><div>boa sorte.</div><div class="gmail_extra"><br><br><div class="gmail_quote">2012/11/1 Benilton Carvalho <span dir="ltr"><<a href="http://user/SendEmail.jtp?type=node&node=4656829&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a>></span><div>
<div><br>
<blockquote style="border-left:2px solid #cccccc;padding:0 1em" class="gmail_quote"><p dir="ltr">Bom, se no seu primeiro email vc postou o número errado... Difícil entender... Tento olhar nisso depois. Mas a dica tá aí... b</p>
<div class="gmail_quote"><div><div>On 1 Nov 2012 18:59, "Pedro Rafael" <<a href="http://user/SendEmail.jtp?type=node&node=4656829&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>> wrote:<br type="attribution">
</div></div><blockquote style="border-left:2px solid #cccccc;padding:0 1em" class="gmail_quote"><div><div>
Na verdade foi erro meu, copiei o mesmo valor.<br><br>Mathematica:<br><br>Gamma[1.5, 0, 0.10] = 0.019861<br><br>R:<br><br>> pgamma(1.5, .10, lower=FALSE)*gamma(.10)<br>[1] 0.1079214<br><br>Vendo a formula achei que o certo seria: pgamma(1.5, .10, lower=FALSE)*gamma(1.5) mesmo assim não é a mesma resposta do mathematica.<br>
<br><br clear="all">[ ],<br>Pedro Rafael Diniz Marinho.<br>
<br><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">beniltoncarvalho [via R-br]</b> <span dir="ltr"><<a href="http://user/SendEmail.jtp?type=node&node=4656829&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a>></span><br>
Date: 2012/11/1<br>Subject: Re: [R-br] [Dúvida] Função Gama Incompleta.<br>To: Pedro Rafael Diniz Marinho <<a href="http://user/SendEmail.jtp?type=node&node=4656829&i=3" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
<br>
<br>
Acho q nao deu resultado diferente nao... ou eu estou enxergando coisas estranhas... vc postou no seu email original:<div><div><br></div><div><div><span style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7273px">No Mathematica: N[Gamma[.10, 1.5], 10] = 0.107921</span></div>
</div><div><span style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7273px"><br></span></div></div><div>e eu postei:</div><div>pgamma(1.5, .10, lower=FALSE)*gamma(.10)<br></div><div><br></div><div>que vc confirmou como:</div>
<div>
<div><br></div><div><span style="font-family:arial,sans-serif;font-size:12.7273px">pgamma(1.5, .10, lower=FALSE)*gamma(.10) = 0.1079214</span><br></div><div><span style="font-family:arial,sans-serif;font-size:12.7273px"><br>
</span></div></div><div><span style="font-family:arial,sans-serif;font-size:12.7273px">acho q nao consigo ver a diferenca?</span></div><div><div><br></div>
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4656817&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
<hr color="#cccccc" noshade size="1">
</div><div style="color:rgb(68,68,68);font:12px tahoma,geneva,helvetica,arial,sans-serif"><div>
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
</div><a href="http://r-br.2285057.n4.nabble.com/R-br-Duvida-Funcao-Gama-Incompleta-tp4656799p4656817.html" rel="nofollow" link="external" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Duvida-Funcao-Gama-Incompleta-tp4656799p4656817.html</a>
</div><div><div>
<div style="color:rgb(102,102,102);font:11px/1.5em tahoma,geneva,helvetica,arial,sans-serif;margin-top:0.4em">
To unsubscribe from R-br, <a rel="nofollow" link="external">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" style="font:9px serif" link="external" target="_blank">NAML</a>
</div></div></div></div><br>
<br></div></div><div>_______________________________________________<br>
R-br mailing list<br>
<a href="http://user/SendEmail.jtp?type=node&node=4656829&i=4" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
</div><div><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></div></blockquote></div>
</blockquote></div></div></div><br></div><div><div>
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4656829&i=5" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
<hr color="#cccccc" noshade size="1">
</div></div><div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><div><div>
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
</div></div><a href="http://r-br.2285057.n4.nabble.com/R-br-Duvida-Funcao-Gama-Incompleta-tp4656799p4656829.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Duvida-Funcao-Gama-Incompleta-tp4656799p4656829.html</a>
</div><div><div>
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">
To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=cGVkcm8ucmFmYWVsLm1hcmluaG9AZ21haWwuY29tfDMzNTc5ODJ8NTAyMjI0MDYw" target="_blank">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" style="font:9px serif" target="_blank">NAML</a>
</div></div></div></div></div></div><br>
</blockquote></div><br></div>