BD | R Documentation |
An internal function, for use in sixway
, which calculates the
Brooks-Draper diagnostic, based on an unpublished paper by David Draper. It
estimates the length of a Markov chain required to produce a mean estimate
to k significant figures with a given accuracy (alpha). See Browne (2012)
for further details.
BD(est, var, rho, k = 2, alpha = 0.05)
est |
Numeric scalar for the mean of the distribution |
var |
Numeric scalar for the variance of the distribution |
rho |
The first lag (i.e. after zero) of the auto-correlation function (ACF) diagnostic |
k |
Integer scalar corresponding to the number of significant figures (defaults to |
alpha |
Numeric scalar indicating the desired accuracy (defaults to |
The Brooks-Draper diagnostic statistic is returned.
Zhang, Z., Charlton, C.M.J., Parker, R.M.A., Leckie, G., and Browne, W.J. (2016) Centre for Multilevel Modelling, University of Bristol.
Browne, W.J. (2012) MCMC Estimation in MLwiN, v2.26. Centre for Multilevel Modelling, University of Bristol.
sixway
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.