Practitioners of Bayesian statistics often use Markov chain Monte Carlo (MCMC) samplers to sample from a posterior distribution. This package determines whether the MCMC sample is large enough   to yield reliable estimates of the target distribution. In particular, this calculates a Gelman-Rubin convergence diagnostic using stable and consistent estimators of Monte Carlo variance. Additionally, this uses the connection between an MCMC sample's effective sample size and the Gelman-Rubin diagnostic to produce a threshold for terminating MCMC simulation. Finally, this informs the user whether enough samples have been collected  and (if necessary) estimates the number of samples needed for a desired level of accuracy. The theory underlying these methods can be found in "Revisiting the Gelman-Rubin Diagnostic" by Vats and  Knudson (2018) <arXiv:1812:09384>.

Graphical and computational methods that can be used to assess the 
  stability of results from supervised statistical learning.

An exploratory and heuristic approach for specification search in
    Structural Equation Modeling. The basic idea is to subsample the original data
    and then search for optimal models on each subset. Optimality is defined through
    two objectives: model fit and parsimony. As these objectives are conflicting,
    we apply a multi-objective optimization methods, specifically NSGA-II, to obtain
    optimal models for the whole range of model complexities. From these optimal
    models, we consider only the relevant model specifications (structures), i.e.,
    those that are both stable (occur frequently) and parsimonious and use those to
    infer a causal model.

An implementation of many measures for the
    assessment of the stability of feature selection. Both simple measures
    and measures which take into account the similarities between features
    are available, see Bommert (2020) <doi:10.17877/DE290R-21906>.

Resampling procedures to assess the stability of selected variables
    with additional finite sample error control for high-dimensional variable
    selection procedures such as Lasso or boosting. Both, standard stability
    selection (Meinshausen & Buhlmann, 2010, <doi:10.1111/j.1467-9868.2010.00740.x>) 
    and complementary pairs stability selection with improved error bounds 
    (Shah & Samworth, 2013, <doi:10.1111/j.1467-9868.2011.01034.x>) are
    implemented. The package can be combined with arbitrary user specified
    variable selection approaches.