#!/usr/bin/env python
# coding: utf-8
import pandas as pd
from statsmodels.formula.api import ols
from statsmodels.stats.anova import anova_lm
from statsmodels.stats.multicomp import pairwise_tukeyhsd
 
file = r'D:\\价格.xlsx'
data = pd.read_excel(file)
#print(data)
data['年份'] = data['年份'].astype(dtype='str',copy=True  )
data['年份']
data

formula = '影子价格~ 年份 +  地区  年份 * 地区'
anova_results = anova_lm(ols(formula,data).fit())
anova_results