
# These commands sample from the posterior distribution for the random
# effects model using Metropolis updates.

dist-spec elog.met \
 "u~Normal(0,10^2) + v1~Normal(0,2^2) + v2~Normal(0,2^2)" \
 "t0~Normal(u,Exp(v1)+Exp(v2)/i) + \
  (1-Delta(i-1)) * ExpGamma2(Log(t1),i-1,Exp(v2))"
data-spec elog.met 1 2 / edata .
mc-spec elog.met repeat 25 metropolis 0.5
dist-mc elog.met 10000
