Calculates an R-hat (potential scale reduction factor) value.
# S3 method for mcmcrs
rhat(x, by = "all", as_df = FALSE, na_rm = FALSE, bound = FALSE, ...)An object.
A string indicating whether to determine by "term", "parameter" or "all".
A flag indicating whether to return the results as a data frame versus a named list.
A flag specifying whether to ignore missing values.
flag specifying whether to bind mcmcrs objects by their chains before calculating rhat.
Other arguments passed to methods.
A number >= 1 indicating the rhat value.
By default the uncorrected, unfolded, univariate, split R-hat value.
Gelman, A., and Rubin, D.B. 1992. Inference from Iterative Simulation Using Multiple Sequences. Statistical Science 7(4): 457–472.
Other convergence:
converged_pars(),
converged_terms(),
converged(),
esr_pars(),
esr_terms(),
esr(),
rhat_pars(),
rhat_terms()
rhat(mcmcrs(mcmcr_example, mcmcr_example))
#> $mcmcr1
#> [1] 2.002
#>
#> $mcmcr2
#> [1] 2.002
#>
rhat(mcmcrs(mcmcr_example, mcmcr_example), bound = TRUE)
#> Warning: `rhat.mcmcrs(x, bound = TRUE)` returns scalar was deprecated in mcmcr
#> 0.6.1.9001.
#> ℹ Please use `rhat(x, bound = TRUE)$bound` for previous behaviour instead.
#> ℹ `rhat.mcmcrs(bound = TRUE)` now returns named list of scalars.
#> $mcmcr1
#> [1] 2.002
#>
#> $mcmcr2
#> [1] 2.002
#>
#> $bound
#> [1] 1.891
#>