Sample difficult target distributions with Metropolis–Hastings, Gibbs sampling, Hamiltonian Monte Carlo, and convergence diagnostics.