The function runMCMC_btadjust() returns a mcmc.list object which is the output of a Markov Chain Monte Carlo obtained - from either 'JAGS', 'nimble' or 'greta' - after adjusting burn-in and thinning parameters to meet pre-specified criteria in terms of convergence & effective sample size. Used with 'nimble', runMCMC_btadjust() allows extra calculations (e.g. information criteria for model comparison and goodness-of-fit p-values for model diagnosis).