<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Prezado Wagner,<div>O teste de Vuong está implementado no pacote "pscl".</div><div>Para comparar dois modelos m1 e m2, utiliza-se o código vuong (m1,m2).</div><div><br></div><div>Nos links abaixo você pode consultar alguns exemplos de aplicação do teste.</div><div><br></div><div>http://www.ats.ucla.edu/stat/r/dae/zinbreg.htm</div><div>http://www.ats.ucla.edu/stat/r/dae/zipoisson.htm</div><div><br></div><div>Att,</div><div>José Vilton Costa<br><br><div><div id="SkyDrivePlaceholder"></div><hr id="stopSpelling">Date: Fri, 21 Dec 2012 05:54:04 -0800<br>From: wtassinari@yahoo.com.br<br>To: r-br@listas.c3sl.ufpr.br<br>Subject: [R-br] Teste de Vuong para comparação entre os modelos Poisson, Binomial Negativo e ZIP<br><br><div style="color:#000;background-color:#fff;font-family:times new roman, new york, times, serif;font-size:12pt"><div><span>Saudações,</span></div><div style="color:rgb(0, 0, 0);font-size:16px;font-family:times new roman,new york,times,serif;background-color:transparent;font-style:normal"><br><span></span></div><div style="color:rgb(0, 0, 0);font-size:16px;font-family:times new roman,new york,times,serif;background-color:transparent;font-style:normal"><span>Aproveitando essa discussão dos modelos ZINB e ZIP, gostaria de saber se já existe implementado no R o teste de Vuong, para comparar o ajuste entre os modelos zero-inflated (ZIP), Poisson e os Binomiais Negativos. A referência é a </span>do próprio Vuong Q.H. 1989. "Likelihood ratio tests for model
selection and non-nested hypotheses." <em>Econometrica</em>. 57:307-333. <br></div><div style="color:rgb(0, 0, 0);font-size:16px;font-family:times new roman,new york,times,serif;background-color:transparent;font-style:normal"><br></div><div style="color:rgb(0, 0, 0);font-size:16px;font-family:times new roman,new york,times,serif;background-color:transparent;font-style:normal">Abs,</div><div style="color:rgb(0, 0, 0);font-size:16px;font-family:times new roman,new york,times,serif;background-color:transparent;font-style:normal"><br></div><div style="color:rgb(0, 0, 0);font-size:16px;font-family:times new roman,new york,times,serif;background-color:transparent;font-style:normal">Wagner Tassinari<br></div><div><br></div> <div style="font-family:times new roman, new york, times, serif;font-size:12pt"> <div style="font-family:times new roman, new york, times, serif;font-size:12pt"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1"> <b><span style="font-weight:bold">De:</span></b> Humberto Hazin <hghazin@hotmail.com><br> <b><span style="font-weight:bold">Para:</span></b> R_br <r-br@listas.c3sl.ufpr.br>; R_br <r-br@listas.c3sl.ufpr.br> <br> <b><span style="font-weight:bold">Enviadas:</span></b> Quinta-feira, 20 de Dezembro de 2012 9:30<br> <b><span style="font-weight:bold">Assunto:</span></b> Re: [R-br] Cálculo do percentage prediction error em ZINB<br> </font> </div> <br><div id="ecxyiv776560573"><div><div style="font-family:Calibri, "Segoe UI", Meiryo, "Microsoft YaHei UI", "Microsoft JhengHei UI", "Malgun Gothic", "Khmer UI", "Nirmala UI", Tunga, "Lao UI", Ebrima, sans-serif;font-size:16px"><div>Opa Leonard,</div><div> </div><div>Muito obrigado vou dar uma olhada! Sim estou usando o pacote pscl para rodar os modelos
ZINB! Se eu conseguir resolver o problema coloco na lista. </div><div>Cara se não me engano tive esse problema também a um tempo atrás, veja se código abaixo para ver se resolve o teu problema. Só uma observação, com esse código não precisas rodar o pacote pscl. </div><div> </div><div><a rel="nofollow" href="https://www.dropbox.com/s/8l480on83dvm3m7/zeroinflCorreto.R" target="_blank">https://www.dropbox.com/s/8l480on83dvm3m7/zeroinflCorreto.R</a></div><div> </div><div>Um abraço,</div><div> </div><div>Humberto</div><div> </div><div>Enviado do Email do Windows</div><div> </div> <div style="border-top-color:rgb(229, 229, 229);border-top-width:2px;border-top-style:solid"> <b>De:</b> Leonard de Assis<br> <b>Enviado:</b> 19 de dezembro de 2012 22:04<br> <b>Para:</b> r-br@listas.c3sl.ufpr.br<br> <b>Assunto:</b> Re: [R-br] Cálculo do percentage
prediction error em ZINB<br> </div> <div> </div>
<div class="ecxyiv776560573moz-cite-prefix">Humberto<br>
<br>
vc está usando o pscl?<br>
<br>
Eu estou no momento com um problema envolvendo ZINB, só que o meu
é um pouco mais chato, pois meu modelo é multivariado (E pra
complicar, está dando que uma das dimensões é singular, hehehe)<br>
<br>
Sugiro ler os tutoriais, vc acha um monte via rseek.<br>
<br>
Esse aqui é um deles (De antemão, não tem exatamente o que vc
precisa, mas tem muita coisa)<br>
<br>
<a rel="nofollow" class="ecxyiv776560573moz-txt-link-freetext" href="http://www.ats.ucla.edu/stat/r/dae/zinbreg.htm" target="_blank">http://www.ats.ucla.edu/stat/r/dae/zinbreg.htm</a><br>
<pre class="ecxyiv776560573moz-signature">[]s
Leonard de Assis
<a rel="nofollow" class="ecxyiv776560573moz-txt-link-freetext" href="http://about.me/ldeassis" target="_blank">http://about.me/ldeassis</a></pre>
Em 19/12/2012 17:44, Humberto Hazin escreveu:<br>
</div>
<blockquote>
<div>
<div>
<div>
<div>
<pre>Prezados,
Estou com uma dúvida que pode parecer besteira, porém gostaria da ajuda
de vocês! A dúvida é a seguinte
como calcular o percentual de "Bias" (percentage prediction error) por
linha de um modelo ZINB (package- pscl)? O meu objetivo é fazer um boxplot com esse percentual por ano </pre>
<pre>Eu iniciei da seguinte forma:</pre>
<pre>dois bancos de dados um de try.csv (10% retirados do banco original) e outro original.csv</pre>
<pre>try<-read.csv(”try.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE)</pre>
<pre><a rel="nofollow" href="https://www.dropbox.com/s/v7vruztus11nnbi/try.csv" target="_blank">https://www.dropbox.com/s/v7vruztus11nnbi/try.csv</a></pre>
<pre> </pre>
<pre>original<-read.csv(”original.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE)</pre>
<pre> <a rel="nofollow" href="https://www.dropbox.com/s/7htqc2bqi45rkku/original.csv" target="_blank"><u><font color="#0066cc">https://www.dropbox.com/s/7htqc2bqi45rkku/original.csv</font></u></a></pre>
<pre># Criar Fatores para o banco original
original$area<-rep(1,nrow(original))
original$area[original$lat<=-15&original$lat>=-45]<-2
original$A=as.factor(original$area)
summary(original)</pre>
<pre>original$S<-factor(original$id.group.fleet)
original$Q<-factor(original$Quarter)
original$Y<-factor(original$Year)</pre>
<pre>original$logEf<-log(original$effort1)</pre>
<pre>summary(original)
</pre>
<pre># Criar Fatores para o banco try</pre>
<pre>try$area<-rep(1,nrow(try))
try$area[try$lat<=-15&try$lat>=-45]<-2
try$A=as.factor(try$area)
summary(try)</pre>
<pre>try$S<-factor(try$id.group.fleet)
try$Q<-factor(try$Quarter)
try$Y<-factor(try$Year)
try$logEf<-try$effort1*0
summary(try)</pre>
<pre> </pre>
<pre>zeroinfl.fit_original<-zeroinfl(BSH~Y+S+Q+A+offset(logEf)|Y+S+Q+A,dist="negbin",data=original)</pre>
<pre>pred<-predict(zeroinfl.fit_original,newdata=try,type="response")</pre>
<pre> </pre>
<pre>Daqui em diante não tenho mais ideia de como gerar esse percentage prediction error</pre>
<pre> </pre>
<pre>Alguem poderia me ajudar!</pre>
<pre>Desde já Agradeço</pre>
<pre> </pre>
<pre>Humberto </pre>
<pre> </pre>
<pre> </pre>
<pre> </pre>
<pre> </pre>
</div>
<div>Enviado do Email do Windows</div>
<div> </div>
</div>
</div>
</div>
<br>
<fieldset class="ecxyiv776560573mimeAttachmentHeader"></fieldset>
<br>
<pre>_______________________________________________
R-br mailing list
<a rel="nofollow" class="ecxyiv776560573moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a rel="nofollow" class="ecxyiv776560573moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a rel="nofollow" class="ecxyiv776560573moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
</div></div></div><br>_______________________________________________<br>R-br mailing list<br><a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br><br> </div> </div> </div><br>_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
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.</div></div> </div></body>
</html>