samples.bgr: Plot the Gelman-Rubin convergence statistic

Description Usage Arguments Details Value Note References See Also


This function calculates and plots the Gelman-Rubin convergence statistic, as modified by Brooks and Gelman (1998).


samplesBgr(node, beg = samplesGetBeg(), end = samplesGetEnd(), 
    firstChain = samplesGetFirstChain(), 
    lastChain = samplesGetLastChain(), thin = samplesGetThin(),
    bins = 50, plot = TRUE, mfrow = c(3, 2), ask = NULL, 
    ann = TRUE, ...)



Character vector of length 1, name of a variable in the model.

beg, end

Arguments to select a slice of monitored values corresponding to iterations beg:end.

firstChain, lastChain

Arguments to select a sub group of chains to calculate the Gelman-Rubin convergence statistics for. Number of chains must be larger than one.


Only use every thin-th value of the stored sample for statistics.


Number of blocks


Logical, whether to plot the BGR statistics or only return the values. If TRUE, values are returned invisibly.

mfrow, ask, ann

Graphical parameters, see par for details. ask defaults to TRUE unless it is plotting into an already opened non-interactive device. The ann parameter is not available in S-PLUS, and will be ignored if it is set.


Further graphical parameters as in par may also be passed as arguments to plotBgr.


The width of the central 80% interval of the pooled runs is green, the average width of the 80% intervals within the individual runs is blue, and their ratio R (= pooled / within) is red. For plotting purposes the pooled and within interval widths are normalised to have an overall maximum of one. The statistics are calculated in bins of length 50: R would generally be expected to be greater than 1 if the starting values are suitably over-dispersed. Brooks and Gelman (1998) emphasise that one should be concerned both with convergence of R to 1, and with convergence of both the pooled and within interval widths to stability.

If the variable of interest is an array, slices of the array can be selected using the notation variable[lower0:upper0, lower1:upper1, ...]. A star ‘*’ can be entered as shorthand for all the stored samples.

If the arguments are left at their defaults the whole sample for all chains will be used for calculation.


A list containing data frames - one for each scalar variable contained in argument node. Each data frames contains elements


end iteration of corresponding bin


80pct interval (normalized) of pooled chains


80pct interval (normalized) of mean within chain


BGR ratio


If the MCMC simulation has an adaptive phase it will not be possible to make inference using values sampled before the end of this phase.


Brooks, S.P. and Gelman A. (1998): Alternative Methods for Monitoring Convergence of Iterative Simulations. Journal of Computational and Graphical Statistics, 7, 434-455.

See Also

plotBgr, BRugs, help.WinBUGS

BRugs documentation built on Sept. 24, 2021, 5:06 p.m.

Related to samples.bgr in BRugs...