statsmodels双因素方差分析

'''
1.正态性
2.方差齐性
3.独立性
'''
#双因素方差分析
import pandas as pd
import numpy as np
from statsmodels.formula.api import ols
from statsmodels.stats.anova import anova_lm
data = pd.DataFrame([[1, 1, 32],
                     [1, 2, 35],
                     [1, 3, 35.5],
                     [1, 4, 38.5],
                     [2, 1, 33.5],
                     [2, 2, 36.5],
                     [2, 3, 38],
                     [2, 4, 39.5],
                     [3, 1, 36],
                     [3, 2, 37.5],
                     [3, 3, 39.5],
                     [3, 4, 43]], 
                    columns=['A', 'B', 'value'])

model = ols('value~C(A) + C(B)', data=data[['A', 'B', 'value']]).fit()
#数据表导入的数据不需要指定数据
#model = ols('value~C(group1) + C(group2)',data).fit()  
anovat = anova_lm(model)
print(model.summary())
print(anovat)

结果：