Accessory function to HPD... function series
mcmc.pval(dat, testlim = 0, sided = 2, ptype="z")a table of MCMC samples
positive value to test whether the posterior crosses it (for variance components)
sidedness of the test
'z' for z-score based p-values, 'mcmc' for standard sampling based p-values
A vector of calculated Bayesian p-values. For standard 'mcmc' p=values, the minimal possible value is 1/(MCMC sample size). z-score based p-values are useful to approximate very low p-values with limited MCMC sample sizes.