<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span>Olá Sérgio.</span></div><div><span>Tb aconteceu comigo.</span></div><div><span>Estou ajudando uma aluna da Ensp tb e usando a mesma modelagem com dados espaço-temporais.</span></div><div><span>Até achei que fosse problema no pacote, e estou tentando rodar diretamente no WinBugs.</span></div><div><span>Alguns modelos tive êxito.</span></div><div><span>Abs.</span></div><div><span></span></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> Sérgio Henrique almeida da silva ju <sergio.edfisica@gmail.com><br> <b><span style="font-weight: bold;">Para:</span></b>
r-br@listas.c3sl.ufpr.br <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Sexta-feira, 11 de Maio de 2012 11:55<br> <b><span style="font-weight: bold;">Assunto:</span></b> [R-br] Modelo Bayesiano<br> </font> </div> <br>Olá<br><br>Estou tentando estimar a prevalência de uma doença (True Prevalence),<br>através de um modelo Bayesiano (Joseph L, Gyorkos TW, Coupal L<br>Bayesian estimation of<br>disease prevalence and the parameters of diagnostic tests in the<br>absence of a gold standard Am J Epidemiol 1995,141:263-72.)<br><br>Para isso estou usando a biblioteca (R2WinBUGS) e modelando através de<br>distribuição Binomial Negativa<br><br>##Parameterized by p and r:<br>model {<br> for( i in 1 : N ) {<br> y[i] ~ dnegbin( p , r )<br> }<br> p ~ dbeta(1.001,1.001)<br> r ~ dgamma(0.01,0.01)<br> m <- r*(1-p)/p<br> v <-
r*(1-p)/(p*p)<br>}<br><br><br>Porém quando rodo o modelo:<br><br>N <-nrow(dados1)<br>y <-dados1$esqui<br>data <- list("N", "y")<br><br>inits <- function(){list(p=rbeta(N,1.001,1.001),r=rgamma(1,0.01,0.01))}<br><br>esqui.sim <- bugs(data, inits, model.file = "esqui.bugs",<br>parameters = c("p", "r"), n.chains = 2, n.iter = 1000,<br>bugs.directory = "C:/Users/sergio/Downloads/winbugs14",debug=T)<br><br>Aparece o seguinte erro no log:<br><br>display(log)<br>check(C:/Users/sergio/AppData/Local/Temp/Rtmp4SGBkz/esqui.bugs.txt)<br>model is syntactically correct<br>data(C:/Users/sergio/AppData/Local/Temp/Rtmp4SGBkz/data.txt)<br>expected collection operator c<br>compile(2)<br>inits(1,C:/Users/sergio/AppData/Local/Temp/Rtmp4SGBkz/inits1.txt)<br>command #Bugs:inits cannot be executed (is greyed out)<br>inits(2,C:/Users/sergio/AppData/Local/Temp/Rtmp4SGBkz/inits2.txt)<br>command #Bugs:inits cannot be executed (is greyed
out)<br>gen.inits()<br>command #Bugs:gen.inits cannot be executed (is greyed out)<br>thin.updater(1)<br>update(500)<br>command #Bugs:update cannot be executed (is greyed out)<br>set(p)<br>command #Bugs:set cannot be executed (is greyed out)<br>set(r)<br>command #Bugs:set cannot be executed (is greyed out)<br>set(deviance)<br>command #Bugs:set cannot be executed (is greyed out)<br>dic.set()<br>command #Bugs:dic.set cannot be executed (is greyed out)<br>update(500)<br>command #Bugs:update cannot be executed (is greyed out)<br>coda(*,C:/Users/sergio/AppData/Local/Temp/Rtmp4SGBkz/coda)<br>command #Bugs:coda cannot be executed (is greyed out)<br>stats(*)<br>command #Bugs:stats cannot be executed (is greyed out)<br>dic.stats()<br><br>DIC<br>history(*,C:/Users/sergio/AppData/Local/Temp/Rtmp4SGBkz/history.odc)<br>command #Bugs:history cannot be executed (is greyed
out)<br>save(C:/Users/sergio/AppData/Local/Temp/Rtmp4SGBkz/log.odc)<br>save(C:/Users/sergio/AppData/Local/Temp/Rtmp4SGBkz/log.txt)<br><br>O que estou fazendo de errado?<br><br>Abraços<br><br>-- <br>Sérgio Henrique Almeida da Silva Junior<br>Doutorando em Epidemiologia em Saúde Pública<br>Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ<br><a href="http://lattes.cnpq.br/1611345552843383" target="_blank">http://lattes.cnpq.br/1611345552843383</a><br>Tel: (21) 94429486/78101651 id: 123*20942<br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" 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><br> </div> </div> </div></body></html>