best_models

Assesses the precision in estimating the ranking of the \(k\) best-performing models.

Estimates the precision of transdimensional Markov chain Monte Carlo
(MCMC) output, which is often used for Bayesian analysis of models with different
dimensionality (e.g., model selection). Transdimensional MCMC (e.g., reversible
jump MCMC) relies on sampling a discrete model-indicator variable to estimate
the posterior model probabilities. If only few switches occur between the models,
precision may be low and assessment based on the assumption of independent
samples misleading. Based on the observed transition matrix of the indicator
variable, the method of Heck, Overstall, Gronau, & Wagenmakers (2019,
Statistics & Computing, 29, 631-643) <doi:10.1007/s11222-018-9828-0> draws
posterior samples of the stationary distribution to (a) assess the uncertainty
in the estimated posterior model probabilities and (b) estimate the effective
sample size of the MCMC output.

Daniel W. Heck

MCMCprecision

Precision of Discrete Parameters in Transdimensional MCMC

best_models function

<dl><dt>samples</dt>
<dd>a matrix with posterior samples (one per row) for the model posterior probabilities (one model per column). Can be estimated using <code>stationary</code> with the argument <code>summary = FALSE</code>.</dd>
<dt>k</dt>
<dd>number of best-performing models to be considered</dd>
<dt>ties.method</dt>
<dd>a character string specifying how ties are treated, see <code><a href="/link/rank?package=MCMCprecision&version=0.4.2" data-mini-rdoc="MCMCprecision::rank">rank</a></code></dd></dl>

Arguments

Precision for the k Best-Performing Models — best_models

<dl>

<dt>samples</dt>
<dd>a matrix with posterior samples (one per row) for the model posterior probabilities (one model per column). Can be estimated using <code>stationary</code> with the argument <code>summary = FALSE</code>.</dd>


<dt>k</dt>
<dd>number of best-performing models to be considered</dd>


<dt>ties.method</dt>
<dd>a character string specifying how ties are treated, see <code><a href='https://rdrr.io/r/base/rank.html'>rank</a></code></dd>

</dl>

best_models: Precision for the k Best-Performing Models

Description

Usage

Arguments

Examples