source("../../R/Rcheck.R")
data <- read.jagsdata("ice-data.R")
load.module("glm")
m <- jags.model("icear.bug", data, n.chains=2)
check.data(m, data)
update(m, 10000)
x <- coda.samples(m, c("sigma", "logRR"), n.iter=50000, thin=50)
source("bench-test1.R")
check.fun()