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

1 Nov 2012

      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>

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

Pedro Rafael

Benilton Carvalho

Pedro Rafael

Benilton Carvalho

Benilton Carvalho

Pedro Rafael

Pedro Rafael

etiquetas

participantes (2)