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]>
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
To unsubscribe from R-br, click here.
NAML