So if I ran an experiment with multiple treatment groups and a control, I would analyse the results using Statsmodel ols to see if any of the treatment group were statistically different from the control group: y ~ C(treatment_group, Treatment(‘Control’) I would then run results.t_test_pairwise() to find out if the coefficients of each treatment group were equal. I.e. to know