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).

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

`node` |
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 |

`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. |

`thin` |
Only use every |

`bins` |
Number of blocks |

`plot` |
Logical, whether to plot the BGR statistics or only return the values.
If |

`mfrow, ask, ann` |
Graphical parameters, see |

`...` |
Further graphical parameters as in |

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

`Iteration` |
end iteration of corresponding bin |

`pooledChain80pct)` |
80pct interval (normalized) of pooled chains |

`withinChain80pct` |
80pct interval (normalized) of mean within chain |

`bgrRatio` |
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.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.