Hi,
Thank you for excellent stat package. I would like to know how to perform a pairwise comparison for group interactions using Tukey HSD. I tried to use pairwise_tukey from pingouin, but it does not support pairwise comparison for group interactions.
In R it can be done like this (see https://rpubs.com/tmcurley/twowayanova)
TukeyHSD(len.aov, which = "supp:dose")
Let me know if there is a solution to achieve this. If not, it will be a good addition to the pingouin package.