1source("../../R/Rcheck.R")
2d <- read.jagsdata("inhaler-data.R")
3inits <- read.jagsdata("inhaler-inits.R")
4load.module("glm")
5m <- jags.model("inhaler1.bug", d, inits, n.chains=2)
6names(list.samplers(m))
7update(m, 1000)
8x <- coda.samples(m, c("a","beta","kappa","log.sigma","pi","sigma"),
9                  n.iter=20000, thin=10)
10source("bench-test1.R")
11check.fun()
12