powered by
Proposal for burn-in and thin rate
diagnostic(results, random)
Bayes.fit class object
one out of 1, 2, c(1,2)