Mesmo assim deu valores diferentes. hy Mathematica: Gamma[1.5, 0, 0.10] = 0.019861 OBS: O argumento zero é para dizer que estou calculando a gamma inferior, integral de 0 à x; R: pgamma(1.5, .10, lower=FALSE)*gamma(.10) = 0.1079214 Uma outra dúvida que tenho é sobre a função hipergeométrica 2F1. Estou utilizando o pacote hypergeo. A minha função que envolva a função hypergeo é: fda_bllg <- function(par,x){ a <- par[1] b <- par[2] alpha <- par[3] Beta <- par[4] (((x^(-Beta))^(-a))*((1/alpha)^(a*Beta))*Re(hypergeo(a, a+b, 1+a, -(x/alpha)^Beta, tol = 0, maxiter=20000)))/(a*beta(a,b)) } fda_bllg trata-se da função de distribuição acumulada da distribuição Beta Log-Logistica KLL, uma distribuição com 4 parâmetros utilizada em análise de sobrevivência e suporte x>=0. Quando substituo os estimadores de máxima verossimilhança ocorre erro. As EMV da distribuição KLL são: est = c(6.474587e+02,2.043228e+02, 5.867625e-03,1.210073e-01) Seja x = 10. Temos que: fda_bllg(est,10) A mensagem de erro é: [1] Inf Mensagens de aviso perdidas: 1: In i15.3.8(A, B, C) : valor fora de limites em 'gammafn' [ ], Pedro Rafael Diniz Marinho. Em 1 de novembro de 2012 14:41, Pedro Rafael <pedro.rafael.marinho@gmail.com
escreveu:
Mesmo assim deu valores diferentes. hy
Mathematica: Gamma[1.5, 0, 0.10] = 0.019861
OBS: O argumento zero é para dizer que estou calculando a gamma inferior, integral de 0 à x;
R:
pgamma(1.5, .10, lower=FALSE)*gamma(.10) = 0.1079214
Uma outra dúvida que tenho é sobre a função hipergeométrica 2F1. Estou utilizando o pacote hypergeo. A minha função que envolva a função hypergeo é:
fda_bllg <- function(par,x){ a <- par[1] b <- par[2] alpha <- par[3] Beta <- par[4] (((x^(-Beta))^(-a))*((1/alpha)^(a*Beta))*Re(hypergeo(a, a+b, 1+a, -(x/alpha)^Beta, tol = 0, maxiter=20000)))/(a*beta(a,b)) }
fda_bllg trata-se da função de distribuição acumulada da distribuição Beta Log-Logistica KLL, uma distribuição com 4 parâmetros utilizada em análise de sobrevivência e suporte x>=0. Quando substituo os estimadores de máxima verossimilhança ocorre erro.
As EMV da distribuição KLL são:
est = c(6.474587e+02,2.043228e+02, 5.867625e-03,1.210073e-01)
Seja x = 10. Temos que:
fda_bllg(est,10)
A mensagem de erro é: [1] Inf Mensagens de aviso perdidas: 1: In i15.3.8(A, B, C) : valor fora de limites em 'gammafn'
[ ], Pedro Rafael Diniz Marinho.
Em 1 de novembro de 2012 06:44, beniltoncarvalho [via R-br] < ml-node+s2285057n4656804h1@n4.nabble.com> escreveu:
o problema e' q vc vai precisar revisitar as definicoes da funcao gama
incompleta... ler a documentacao detalhada de pgamma() e comparar com o que Gamma[] no Mathematica faz...
o que vc quer e':
pgamma(1.5, .10, lower=FALSE)*gamma(.10)
2012/11/1 Pedro Rafael <[hidden email]<http://user/SendEmail.jtp?type=node&node=4656804&i=0>
Pessoal qual o comando certo ou pacote que trabalha com a função gamma incompleta? Estou obtendo valores muito diferentes da gamma incompleta no R e no mathematica.
Por exemplo:
No R: pgamma(.10,1.5) = 0.0224107 No Mathematica: N[Gamma[.10, 1.5], 10] = 0.107921
Qual seria o problema que não estou enxergando?
[ ], Pedro Rafael Diniz Marinho.
_______________________________________________ R-br mailing list [hidden email] <http://user/SendEmail.jtp?type=node&node=4656804&i=1> 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=4656804&i=2> 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>