/* Due to very high autocorrelation for sigma, we set a very long thinning
   interval
*/
model in "oxford.bug"
data in "oxford-data.R"
load glm
compile, nchains(2)
inits in "oxford-init.R"
initialize
samplers to foo-samplers.txt
update 1000
monitor alpha, thin(10)
monitor beta1, thin(10)
monitor beta2, thin(10)
monitor sigma, thin(10)
update 20000
coda *