1source("../../R/Rcheck.R") 2load.module("glm") 3d <- read.jagsdata("seeds-data.R") 4inits <- read.jagsdata("seeds-init.R") 5m <- jags.model("seedsuni.bug", d, inits, n.chains=2, n.adapt=2500) 6check.data(m, d) 7update(m, 2500) 8x <- coda.samples(m, c("alpha0", "alpha1","alpha2","alpha12","sigma"), 9 n.iter=10000, thin=10) 10source("bench-test4.R") 11check.fun() 12