[R-br] Fwd: [Dúvida] Função Gama Incompleta.

[Benilton Carvalho]

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 em 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 em 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 em 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 em 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] <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<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>
>>
>>
>> _______________________________________________
>> R-br mailing list
>> R-br em 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.
>>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121102/82bf57d3/attachment.html>