cd '.'
load dic
model in 'example6model.txt'
data in 'example6data.jags.txt'
compile, nchains(3)
initialize
monitor lam
monitor errorvar
monitor latcov
monitor latcor
monitor gamma
monitor a
monitor vareta
monitor pD, type(mean)
monitor deviance, type(mean)
update 5000
coda *, stem('jagsoutputCODA')
coda pD, stem('PD')
coda deviance, stem('DEV')
exit
