第15周 方差分析实践

引入数据

import pandas as pd
data_value = { ‘无促销’:[23,19,17,26,28,23,24,30],
‘被动促销’:[26,22,20,30,36,28,30,32],
‘主动促销’:[30,23,25,32,48,40,41,46]}# 因变量
da = pd.DataFrame( data_value ).stack()
da.columns = [‘水平’,‘观测值’]

ols模块进行分析

from statsmodels.formula.api import ols
from statsmodels.stats.anova import anova_lm

formula = ‘{} ~ {}’.format(da.columns[1], da.columns[0])
model = ols( formula, da ).fit()
anovat = anova_lm(model)
print(anovat)