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