cd '.'
load dic
model in 'example1model.txt'
data in 'example1data.jags.txt'
compile, nchains(3)
parameters in 'example1inits1.jags.txt', chain(1)
parameters in 'example1inits2.jags.txt', chain(2)
parameters in 'example1inits3.jags.txt', chain(3)
initialize
monitor beta
monitor psi
monitor errvar
monitor pD, type(mean)
monitor deviance, type(mean)
update 5000
coda *, stem('jagsoutputCODA')
coda pD, stem('PD')
coda deviance, stem('DEV')
exit
