Dear Bugs Users,
I am running a spatial random effects model in WinBUGS 1.4.1 with GeoBUGS (code below). The model compiled, I loaded most of the inits and generated the remaining ones, which all worked as usual. However, when I started to update, I got a Rejection1 error in the bottom left-hand corner of the screen. I did some poking around in the manual and on the list, but there were no answers there. If I continue to hit the update buttom, WinBUGS will continue to update, but each iteration I continue to get this message. Can anyone tell me what this means?
Thanks in advance,
Dave.
model{
for(i in 1:152){
for(j in 1:24){
ycs[j,i] ~ dpois(p.ycs[j,i])
log(p.ycs[j,i]) <- exp.kill[j,i] + alpha[1,j] + S[j,i]
ed[j,i] ~ dpois(p.ed[j,i])
log(p.ed[j,i]) <- exp.kill[j,i] + alpha[2,j] + S[j,i]
ar[j,i] ~ dpois(p.ar[j,i])
log(p.ar[j,i]) <- exp.kill[j,i] + alpha[3,j] + S[j,i]
hrw[j,i] ~ dpois(p.hrw[j,i])
log(p.hrw[j,i]) <- exp.kill[j,i] + alpha[4,j] + S[j,i]
mean[j,i] <- (1/4)*(p.ycs[j,i] + p.ed[j,i] + p.hrw[j,i] + p.ar[j,i])
}
}
for(k in 1:152){
mu[k] <- 0
}
for(j in 1:24){
S[j,1:152] ~ spatial.exp(mu[], x[], y[], tau, phi, 1)
alpha[1,j] ~ dnorm(0, 0.001)
alpha[2,j] ~ dnorm(0, 0.001)
alpha[3,j] ~ dnorm(0,
0.001)
alpha[4,j] ~ dnorm(0, 0.001)
}
tau ~ dgamma(0.01, 0.01)
phi ~ dunif(0.33, 6)
}
--
Dave Armstrong
University of Oxford
Department of Politics and International Relations
Center for Research Methods in the Social Sciences
Nuffield College
e-mail: [log in to unmask]
web: www.davearmstrong-ps.com
Facts are meaningless. You can use facts to prove anything that's even remotely true. - Homer Simpson
-------------------------------------------------------------------
This list is for discussion of modelling issues and the BUGS software.
For help with crashes and error messages, first mail [log in to unmask]
To mail the BUGS list, mail to [log in to unmask]
Before mailing, please check the archive at www.jiscmail.ac.uk/lists/bugs.html
Please do not mail attachments to the list.
To leave the BUGS list, send LEAVE BUGS to [log in to unmask]
If this fails, mail [log in to unmask], NOT the whole list
