display('window')

##################### Model Checking, compiling

check('D:/alewin/Bugs-spline/BayesDE.txt')
data('D:/alewin/Bugs-spline/quadsim_100.txt')
data('D:/alewin/Bugs-spline/constants_data.txt')
data('D:/alewin/Bugs-spline/dummy_data_100.txt')

#compile(1)
compile(2)

inits(1,'D:/alewin/Bugs-spline/init1_100.txt')
inits(2,'D:/alewin/Bugs-spline/init2_100.txt')
gen.inits()

refresh(1)

###################### Start updates

### adapting for 4050 iterations, can't set monitors yet

update(4050)

set(alpha[1:3])
set(dstar[1:3])
set(c1[1:3,])
set(c2[1:3,])
set(sig21[1:3])
set(sig22[1:3])
set(mm1)
set(mm2)
set(tt1)
set(tt2)

update(950)
update(5000)

####################### Check convergencee

history(*)
autoC(*)
gr(*)

###################### Updates to save

set.summary(alpha)
set.summary(dstar)
set.summary(c1)
set.summary(c2)
set.summary(sig21)
set.summary(sig22)
set.summary(mm1)
set.summary(mm2)
set.summary(tt1)
set.summary(tt2)
set.summary(pval1)
set.summary(pval2)

set(dstar)

dic.set()

thin.updater(10)
update(1000)

####################### Results

##### results displayed in windows
dic.stats()
stats.summary(*)
stats(*)

###### results saved to files
#display('log')
#dic.stats()
#stats.summary(*)
#stats()
#save('D:/alewin/Bugs-spline/bugs_testsave2.txt')

###### output whole posterior of dstar
#coda(dstar,'D:/alewin/Bugs-spline/coda_dstarprime_spline6_c3_exch_simdata42_1')