Calculate the proportion of variables in an mcmc object that have "good" gelman.diag

Calculate the proportion of variables in an mcmc object that have an [coda::gelman.diag()] value below or equal psrf_limit (defalt = 1.10)

Usage

gelman_good_proportion(x, psrf_limit = 1.1)

Arguments

x

an mcmc object

psrf_limit

numeric the limit to classify the diag as good

Value

A numeric vector of length two. The first number is the proportion of gelman.diag() <= `psrf_limit` the second number is the number of variables in `x`.

Examples

require(coda)
#> Loading required package: coda
data(line) # data from package coda
gelman_good_proportion(line)
#>  Proportion N_variables 
#>           1           3