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@gmail.com> escreveu:
Valeu Benilton, agora está funcionando. Obrigado.

[   ],
Pedro Rafael Diniz Marinho.


---------- Mensagem encaminhada ----------
De: beniltoncarvalho [via R-br] <ml-node+s2285057n4656829h76@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@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]>

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]> 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]>
Date: 2012/11/1
Subject: Re: [R-br] [Dúvida] Função Gama Incompleta.
To: Pedro Rafael Diniz Marinho <[hidden email]>


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]
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


_______________________________________________
R-br mailing list
[hidden email]
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]
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.
NAML