1source("../../R/Rcheck.R")
2data <- read.jagsdata("birats-data.R")
3inits <- read.jagsdata("birats-inits.R")
4m <- jags.model("birats4.bug", data, inits, n.chains=2)
5check.data(m, data)
6update(m, 1000)
7x <- coda.samples(m, c("mu.beta","sigma.beta","r", "alpha0"), thin=10,
8                  n.iter=10000)
9source("bench-test4.R")
10check.fun()
11
12