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

[Pedro Rafael]

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/20121103/9761e3c6/attachment.html>