calculates p-value based on Bayesian z-score or MCMC sampling
Accessory function to HPD... function series
mcmc.pval(dat, testlim = 0, sided = 2, ptype="z")
dat |
a table of MCMC samples |
testlim |
positive value to test whether the posterior crosses it (for variance components) |
sided |
sidedness of the test |
ptype |
'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.
Mikhail Matz, UT Austin <matz@utexas.edu>
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.