A model-adequacy diagnostic that simulates replicated datasets from the posterior predictive distribution and compares summary statistics (mean, variance, quantiles, test quantities) of the replications against the observed data. Systematic discrepancies indicate model misspecification. The standard goodness-of-fit procedure in the Bayesian workflow.