/* SEEDS example with uniform prior on the standard deviation of the
random effects. The results should be identical to test4, although
the parameterization of the model is different.
*/
model in "seedspar.bug"
data in "seeds-data.R"
load glm
compile, nchains(2)
inits in "seeds-init.R"
initialize
update 5000
monitor alpha0 , thin(10)
monitor alpha1 , thin(10)
monitor alpha2 , thin(10)
monitor alpha12 , thin(10)
monitor sigma , thin(10)
update 10000
coda *