<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">Ótimo, Humberto,<br>
vou dar uma lida no código e ver se aproveito alguma coisa<br>
<pre class="moz-signature" cols="72">[]s
Leonard de Assis
<a class="moz-txt-link-freetext" href="http://about.me/ldeassis">http://about.me/ldeassis</a></pre>
Em 20/12/2012 09:30, Humberto Hazin escreveu:<br>
</div>
<blockquote
cite="mid:COL401-EAS26495B24E691B982DCD2ACFB4370@phx.gbl"
type="cite">
<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 moz-do-not-send="true"
href="https://www.dropbox.com/s/8l480on83dvm3m7/zeroinflCorreto.R">https://www.dropbox.com/s/8l480on83dvm3m7/zeroinflCorreto.R</a></div>
<div> </div>
<div data-signatureblock="true">Um abraço,</div>
<div data-signatureblock="true"> </div>
<div data-signatureblock="true">Humberto</div>
<div data-signatureblock="true"> </div>
<div data-signatureblock="true">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;"> <strong>De:</strong> Leonard
de Assis<br>
<strong>Enviado:</strong> 19 de dezembro de 2012 22:04<br>
<strong>Para:</strong> <a class="moz-txt-link-abbreviated" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
<strong>Assunto:</strong> Re: [R-br] Cálculo do percentage
prediction error em ZINB<br>
</div>
<div> </div>
<div class="moz-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 moz-do-not-send="true" class="moz-txt-link-freetext"
href="http://www.ats.ucla.edu/stat/r/dae/zinbreg.htm">http://www.ats.ucla.edu/stat/r/dae/zinbreg.htm</a><br>
<pre class="moz-signature">[]s
Leonard de Assis
<a moz-do-not-send="true" class="moz-txt-link-freetext" href="http://about.me/ldeassis">http://about.me/ldeassis</a></pre>
Em 19/12/2012 17:44, Humberto Hazin escreveu:<br>
</div>
<blockquote
cite="mid:COL401-EAS450C254C2E479AF3AA76B4CB4300@phx.gbl">
<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 moz-do-not-send="true" href="https://www.dropbox.com/s/v7vruztus11nnbi/try.csv">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 moz-do-not-send="true" href="https://www.dropbox.com/s/7htqc2bqi45rkku/original.csv"><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="mimeAttachmentHeader"></fieldset>
<br>
<pre>_______________________________________________
R-br mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a moz-do-not-send="true" class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a moz-do-not-send="true" class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
</body>
</html>