This function calculates and plots the Gelman-Rubin convergence statistic for a scalar variable, as modified by Brooks and Gelman (1998).
plotBgr(node, plot = TRUE, main = NULL, xlab = "iteration",
ylab = "bgr", col = c("red", "blue", "green"), bins = 50,
...)Character, name of a scalar variable in the model.
Logical, whether to plot the BGR statistics or only return the values.
If TRUE, values are returned invisibly.
annotation, see plot.default
Colours, see Details Section in samplesBgr.
Number of blocks
Further graphical parameters as in par.
Data frame with elements
end iteration of corresponding bin
80pct interval (normalized) of pooled chains
80pct interval (normalized) of mean within chain
BGR ratio
Acts on a scalar variable. See the wrapper function samplesBgr for more details.