Skip to contents

Gelman-Rubin R-hat for a set of MCMC chains

Usage

morie_bayes_rhat(chains)

Arguments

chains

A list of iteration-by-parameter matrices.

Value

Named numeric R-hat per parameter.

Examples

set.seed(1)
morie_bayes_rhat(list(matrix(rnorm(200), 100), matrix(rnorm(200), 100)))
#> [1] 0.9966664 0.9970985