The superdiag function takes MCMC samples as input. It provides a comprehensive test suite for Markov chain nonconvergence, which integrates four standard empirical MCMC convergence diagnostics: Gelman-Rubin, Geweke, Heidelberger-Welch, and Raftery-Lewis.