load glm
model in "stagnant.bug"
data in "stagnant-data.R"
compile, nchains(2)
parameters in "stagnant-inits1.R"
initialize
update 1000
monitor alpha, thin(10)
monitor beta, thin(10)
monitor k, thin(10)
update 10000 
coda *, stem("stag1")