[R-br] Fwd: [Dúvida] Função Gama Incompleta.
Pedro Rafael
pedro.rafael.marinho em gmail.com
Domingo Novembro 4 17:50:07 BRST 2012
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 é:
B(x;a,b) = integral_0^x t^(a-1) (1-t)^(b-1) dt
A definição do mathematica é:
B(z;a,b) = integral_0^z u^(a-1) (1-u)^(b-1) du
Ver em:
http://mathworld.wolfram.com/IncompleteBetaFunction.html
Concordam que são as mesmas integrais?
Resposta do R: Ibeta(3,1,2) = 0.5
Resposta d Mathematica: N[Beta[3, 1, 2], 5] = -1.5000
Seria um problema da função ou de documentação ou um problema meu ainda?
[ ],
Pedro Rafael Diniz Marinho.
Em 3 de novembro de 2012 10:44, Pedro Rafael <pedro.rafael.marinho em gmail.com
> escreveu:
> Valeu Benilton, agora está funcionando. Obrigado.
>
> [ ],
> Pedro Rafael Diniz Marinho.
>
>
> ---------- Mensagem encaminhada ----------
> De: beniltoncarvalho [via R-br] <ml-node+s2285057n4656829h76 em n4.nabble.com
> >
> Data: 2 de novembro de 2012 14:21
> Assunto: Re: [R-br] Fwd: [Dúvida] Função Gama Incompleta.
> Para: Pedro Rafael Diniz Marinho <pedro.rafael.marinho em gmail.com>
>
>
>
> 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)...
>
> para sua conveniencia e exercicio de fim de semana, esta e' a solucao q vc
> quer:
>
> gammainc <- function(x, a, b){
> r1 <- gamma(x)
> if (missing(b)){
> if (isTRUE(all.equal(a, 0)))
> return(r1)
> return(r1*pgamma(a, x, lower=FALSE))
> }
> gammainc(x, a)-gammainc(x, b)
> }
>
> boa sorte.
>
>
> 2012/11/1 Benilton Carvalho <[hidden email]<http://user/SendEmail.jtp?type=node&node=4656829&i=0>
> >
>
> 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
>> On 1 Nov 2012 18:59, "Pedro Rafael" <[hidden email]<http://user/SendEmail.jtp?type=node&node=4656829&i=1>>
>> wrote:
>>
>>> Na verdade foi erro meu, copiei o mesmo valor.
>>>
>>> Mathematica:
>>>
>>> Gamma[1.5, 0, 0.10] = 0.019861
>>>
>>> R:
>>>
>>> > pgamma(1.5, .10, lower=FALSE)*gamma(.10)
>>> [1] 0.1079214
>>>
>>> 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.
>>>
>>>
>>> [ ],
>>> Pedro Rafael Diniz Marinho.
>>>
>>>
>>> ---------- Forwarded message ----------
>>> From: beniltoncarvalho [via R-br] <[hidden email]<http://user/SendEmail.jtp?type=node&node=4656829&i=2>
>>> >
>>> Date: 2012/11/1
>>> Subject: Re: [R-br] [Dúvida] Função Gama Incompleta.
>>> To: Pedro Rafael Diniz Marinho <[hidden email]<http://user/SendEmail.jtp?type=node&node=4656829&i=3>
>>> >
>>>
>>>
>>> Acho q nao deu resultado diferente nao... ou eu estou enxergando coisas
>>> estranhas... vc postou no seu email original:
>>>
>>> No Mathematica: N[Gamma[.10, 1.5], 10] = 0.107921
>>>
>>> e eu postei:
>>> pgamma(1.5, .10, lower=FALSE)*gamma(.10)
>>>
>>> que vc confirmou como:
>>>
>>> pgamma(1.5, .10, lower=FALSE)*gamma(.10) = 0.1079214
>>>
>>> acho q nao consigo ver a diferenca?
>>>
>>>
>>> _______________________________________________
>>> R-br mailing list
>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4656817&i=0>
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>>> código mínimo reproduzível.
>>>
>>> ------------------------------
>>> If you reply to this email, your message will be added to the
>>> discussion below:
>>>
>>> http://r-br.2285057.n4.nabble.com/R-br-Duvida-Funcao-Gama-Incompleta-tp4656799p4656817.html
>>> To unsubscribe from R-br, click here.
>>> NAML<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>
>>>
>>>
>>> _______________________________________________
>>> R-br mailing list
>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4656829&i=4>
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>>> código mínimo reproduzível.
>>>
>>
>
> _______________________________________________
> R-br mailing list
> [hidden email] <http://user/SendEmail.jtp?type=node&node=4656829&i=5>
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
> ------------------------------
> If you reply to this email, your message will be added to the discussion
> below:
>
> http://r-br.2285057.n4.nabble.com/R-br-Duvida-Funcao-Gama-Incompleta-tp4656799p4656829.html
> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=cGVkcm8ucmFmYWVsLm1hcmluaG9AZ21haWwuY29tfDMzNTc5ODJ8NTAyMjI0MDYw>
> .
> NAML<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>
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121104/5eb547cf/attachment.html>
Mais detalhes sobre a lista de discussão R-br