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 <beniltoncarvalho@gmail.com>

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" <pedro.rafael.marinho@gmail.com> 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] <ml-node+s2285057n4656817h4@n4.nabble.com>
Date: 2012/11/1
Subject: Re: [R-br] [Dúvida] Função Gama Incompleta.
To: Pedro Rafael Diniz Marinho <pedro.rafael.marinho@gmail.com>


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
R-br@listas.c3sl.ufpr.br
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.