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]>
_______________________________________________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]
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]
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-tp4656799p4656804.html