
# These commands fit a univariate t-distribution with d degrees of freedom
# to the data in the file "tdata".  The prior for the location parameter 
# is Normal(0,10^2).  The prior for the log of the scale parameter is 
# Normal(0,1).  Metropolis updates are done, with a stepsize of 1 for 
# w and 2 for u.

dist-spec tlog12 d=2 "u ~ Normal(0,10^2) + w ~ Normal(0,1)" \
                     "w + [(d+1)/2] Log { 1 + [(t-u)/Exp(w)]^2/d }"
data-spec tlog12 0 1 / tdata .
dist-stepsizes tlog12 w=1 u=2
mc-spec tlog12 repeat 10 metropolis 1
dist-mc tlog12 5000
