[R] Error in JAGS, cannot monitor z
alon.benari
alon.benari at gmail.com
Fri Jan 27 15:51:23 CET 2012
Hello.
I am running a state space model (the example from Marc kery's Bayesian
population analysis using WinBUGS )
The code runs fine for the examle as shown in the book.
However I played around and changed the nodes monitored to the latent
variable z I get the following error- even though it is monitored in the
code:
Initializing model
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100%
Error in jags.samples(model, variable.names, n.iter, thin, type = "trace",
:
Failed to set trace monitor for z
Node not found
The code is:
model {
# Priors and constraints
for (i in 1:nind){
for (t in f[i]:(n.occasions-1)){
phi[i,t] <- alpha[t]
p[i,t] <- beta[t]
} #t
} #i
for (t in 1:n.occasions-1)
{
alpha[t] ~ dbeta(1, 1) # Prior for mean survival
beta[t] ~ dbeta(1, 1) # Prior for mean recapture
}
# Prior for mean recapture
# Likelihood
for (i in 1:nind){
# Define latent state at first capture
z[i,f[i]] <- 1
for (t in (f[i]+1):n.occasions){
# State process
z[i,t] ~ dbern(mu1[i,t])
mu1[i,t] <- phi[i,t-1] * z[i,t-1]
# Observation process
y[i,t] ~ dbern(mu2[i,t])
mu2[i,t] <- p[i,t-1] * z[i,t]
} #t
} #i
}
ni <- 10000
nt <- 6
nb <- 5000
nc <- 3
cjs.ran <- jags.model(file="cjs22time.jags",data=bugs.data, n.chains =
nc,n.adapt=2000)
coda.cjs<-coda.samples(cjs.ran,c("z"),n.iter=5000,thin=5)
Any idea why would not trace z node?
Thank you.
--
View this message in context: http://r.789695.n4.nabble.com/Error-in-JAGS-cannot-monitor-z-tp4333685p4333685.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list