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

[Pedro Rafael]

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 em 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 em 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-tp4656799p4656804.html
>>  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>
>>
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121101/1ae7efb9/attachment.html>