model {
y <- 10
y ~ dt(mu, 1, 4)
mu ~ dt(0, 1, 4)
}
Dbar = post.mean of -2logL; Dhat = -2LogL at post.mean of stochastic nodes
Dbar Dhat pD DIC
y 10.791 12.015 -1.224 9.568
total 10.791 12.015 -1.224 9.568
node mean sd MC error 2.5% median 97.5% start sample
mu 4.914 4.286 0.2148 -1.164 4.724 11.2 1001 10000