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