如何运行python post hoc测试？

import pandas as pd
from statsmodels.stats.multicomp import pairwise_tukeyhsd

data = pd.read_csv('data.csv')  # 读取数据文件，可以是CSV、Excel等格式

# 假设数据中有一个因变量（dependent variable）'y'和一个自变量（independent variable）'group'
from statsmodels.formula.api import ols
model = ols('y ~ group', data=data).fit()
anova_table = sm.stats.anova_lm(model, typ=2)

posthoc = pairwise_tukeyhsd(data['y'], data['group'])
print(posthoc)