Simultaneous credible band computation
scrBesag(samples, level=0.95)
matrix with ‘lower’ and ‘upper’ rows
Thomas Kneib
m by n matrix where m is the number of parameters,
n is the number of samples and hence each (multivariate) sample is a column
in the matrix samples
the credible level (default: 0.95)
Calculates a series of simultaneous credible bounds for one parameter type, following section 6.3 in the seminal paper "Bayesian computation and stochastic systems". The corresponding algorithm was invented by Peter Green.
J. Besag, P. Green, D. Higdon, K. Mengersen (1995): Bayesian computation and stochastic systems, Statistical Science 10/1, 3--41