Extract MSE from ANOVA (within-group variance), compute df = N − k to find the appropriate q critical value for the chosen alpha and k, calculate each q = |mean_i − mean_j| / sqrt(MSE / n) (for equal n), compare each q to qcrit, and report which pairs are significant along with their q statistics.