<span class=""><span name="beniltoncarval." class=""></span> </span><div class=""><img class="" id=":1m1" tabindex="0" src="https://mail.google.com/mail/images/cleardot.gif" alt=""></div><div id=":1mf" class=""><div id=":1mg">
Mesmo assim deu valores diferentes. <br>hy<br><br><br>Mathematica:<br>Gamma[1.5, 0, 0.10] = 0.019861<br><br>OBS: O argumento zero é para dizer que estou calculando a gamma inferior, integral de 0 à x;<br><br>R:<br><br>pgamma(1.5, .10, lower=FALSE)*gamma(.10) = 0.1079214<br>
<div class="gmail_extra"><br><br>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 é:<br><br>fda_bllg <- function(par,x){<br>
a <- par[1]<br> b <- par[2]<br> alpha <- par[3]<br> Beta <- par[4]<br> (((x^(-Beta))^(-a))*((1/alpha)^(a*Beta))*Re(hypergeo(a, a+b, 1+a,<br> -(x/alpha)^Beta, tol = 0, maxiter=20000)))/(a*beta(a,b))<br>
}<br>
<br>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. <br>
<br>As EMV da distribuição KLL são:<br><br>est = c(6.474587e+02,2.043228e+02, 5.867625e-03,1.210073e-01)<br><br>Seja x = 10. Temos que:<br><br>fda_bllg(est,10)<br><br>A mensagem de erro é:<br>[1] Inf<br>Mensagens de aviso perdidas:<br>
1: In i15.3.8(A, B, C) : valor fora de limites em 'gammafn'</div></div></div><div class="gmail_extra"><br clear="all">[ ],<br>Pedro Rafael Diniz Marinho.<br>
<br><br><div class="gmail_quote">Em 1 de novembro de 2012 14:41, Pedro Rafael <span dir="ltr"><<a href="mailto:pedro.rafael.marinho@gmail.com" target="_blank">pedro.rafael.marinho@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Mesmo assim deu valores diferentes. <br>hy<br><br><br>Mathematica:<br>Gamma[1.5, 0, 0.10] = 0.019861<br><br>OBS: O argumento zero é para dizer que estou calculando a gamma inferior, integral de 0 à x;<br>
<br>R:<br><br>pgamma(1.5, .10, lower=FALSE)*gamma(.10) = 0.1079214<br>
<div class="gmail_extra"><br><br>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 é:<br><br>fda_bllg <- function(par,x){<br>
a <- par[1]<br> b <- par[2]<br> alpha <- par[3]<br> Beta <- par[4]<br> (((x^(-Beta))^(-a))*((1/alpha)^(a*Beta))*Re(hypergeo(a, a+b, 1+a,<br> -(x/alpha)^Beta, tol = 0, maxiter=20000)))/(a*beta(a,b))<br>
}<br>
<br>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. <br>
<br>As EMV da distribuição KLL são:<br><br>est = c(6.474587e+02,2.043228e+02, 5.867625e-03,1.210073e-01)<br><br>Seja x = 10. Temos que:<br><br>fda_bllg(est,10)<br><br>A mensagem de erro é:<br>[1] Inf<br>Mensagens de aviso perdidas:<br>
1: In i15.3.8(A, B, C) : valor fora de limites em 'gammafn'<div class="im"><br clear="all">[ ],<br>Pedro Rafael Diniz Marinho.<br>
<br><br></div><div class="gmail_quote">Em 1 de novembro de 2012 06:44, beniltoncarvalho [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4656804h1@n4.nabble.com" target="_blank">ml-node+s2285057n4656804h1@n4.nabble.com</a>></span> escreveu:<div>
<div class="h5"><br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
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...<div><br></div><div>o que vc quer e':</div>
<div><br></div><div>pgamma(1.5, .10, lower=FALSE)*gamma(.10)<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">2012/11/1 Pedro Rafael <span dir="ltr"><<a href="http://user/SendEmail.jtp?type=node&node=4656804&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a>></span><br>
<blockquote style="border-left:2px solid rgb(204,204,204);padding:0px 1em" class="gmail_quote"><div><div>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. <br>
<br><br>Por exemplo: <br><br>No R: pgamma(.10,1.5) = 0.0224107<br>
No Mathematica: N[Gamma[.10, 1.5], 10] = 0.107921<br><br>Qual seria o problema que não estou enxergando?<br><br clear="all">[ ],<br>Pedro Rafael Diniz Marinho.<br>
<br></div></div>_______________________________________________<br>
R-br mailing list<br>
<a href="http://user/SendEmail.jtp?type=node&node=4656804&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4656804&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
<hr color="#cccccc" noshade size="1">
<div style="color:rgb(68,68,68);font:12px tahoma,geneva,helvetica,arial,sans-serif">
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
<a href="http://r-br.2285057.n4.nabble.com/R-br-Duvida-Funcao-Gama-Incompleta-tp4656799p4656804.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Duvida-Funcao-Gama-Incompleta-tp4656799p4656804.html</a>
</div>
<div style="color:rgb(102,102,102);font:11px/1.5em tahoma,geneva,helvetica,arial,sans-serif;margin-top:0.4em">
To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=cGVkcm8ucmFmYWVsLm1hcmluaG9AZ21haWwuY29tfDMzNTc5ODJ8NTAyMjI0MDYw" target="_blank">click here</a>.<br>
<a href="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" rel="nofollow" style="font:9px serif" target="_blank">NAML</a>
</div></blockquote></div></div></div><br></div>
</blockquote></div><br></div>