display("log")
#model: Simeoni_PAGE_Estimation_BUGS
check('model_BUGS.txt')
data('data_BUGS.txt')
compile(1)
gen.inits()
#thin.updater(1)
#thin.samples(1)
update(1000)
set(logLAMBDA0_POP)
set(logLAMBDA1_POP)
set(logK1_POP)
set(logK2_POP)
set(logW0_POP)
set(invCV2)
set(Y_res)
set(Y_mean)

update(5000)
coda(*,'output')
quit('yes')



