<html><head></head><body><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 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 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><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><br>try<-read.csv(”try.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE)</pre><pre><a 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 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<br>original$area<-rep(1,nrow(original))<br>original$area[original$lat<=-15&original$lat>=-45]<-2<br>original$A=as.factor(original$area)<br>summary(original)</pre><pre>original$S<-factor(original$id.group.fleet)<br>original$Q<-factor(original$Quarter)<br>original$Y<-factor(original$Year)</pre><pre>original$logEf<-log(original$effort1)</pre><pre>summary(original)<br></pre><pre># Criar Fatores para o banco try</pre><pre>try$area<-rep(1,nrow(try))<br>try$area[try$lat<=-15&try$lat>=-45]<-2<br>try$A=as.factor(try$area)<br>summary(try)</pre><pre>try$S<-factor(try$id.group.fleet)<br>try$Q<-factor(try$Quarter)<br>try$Y<-factor(try$Year)<br>try$logEf<-try$effort1*0<br>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 data-focusfrompointer="true"> </pre><pre data-focusfrompointer="true">Humberto </pre><pre> </pre><pre> </pre><pre> </pre><pre> </pre></div><div data-signatureblock="true">Enviado do Email do Windows</div><div> </div></div></div></div></body></html>