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