display("log")
#model: Simeoni_PAGE_Estimation_POP_BUGS
check('model_BUGS.txt')
data('data_BUGS.txt')
compile(1)
gen.inits()
#thin.updater(1)
#thin.samples(1)
update(1)
set(logLAMBDA0_POP)
set(logLAMBDA1_POP)
set(logK1_POP)
set(logK2_POP)
set(logW0_POP)
set(TAU_LAMBDA0)
set(TAU_LAMBDA1)
set(invCV2)
set(LAMBDA0_POP)
set(LAMBDA1_POP)
set(K1_POP)
set(K2_POP)
set(W0_POP)
set(LAMBDA0)
set(OMEGA_LAMBDA0)
set(OMEGA_LAMBDA1)
set(CV)
set(Y_res)
set(Y_mean)

update(50)
coda(*,'output')
#quit('yes')



