model clear
data clear
model in "model.txt"
data in "data.txt"
compile
inits in "initial_1.txt"
initialize
update 10000
monitor mu, thin(2)
monitor psi, thin(2)
monitor gam, thin(2)
update 100000
coda *