Hello All, The PR gh-13002 <https://github.com/scipy/scipy/pull/13002> implements Tukey’s Honestly Significant Difference Test, a multiple sample test used to find pairs of means that are significantly different from each other. We’d appreciate your thoughts! The test reports the pairwise differences, p-values, and confidence intervals in separate attribute arrays of shape (n, n). This allows the client to access their desired results programmatically, which is not nearly as easy in other programs. There are a few areas with which extra input would be helpful. Specifically: - output interface - __str__ / __repr__ methods - confidence interval method More details are in the PR. Many thanks, Sam