ANOVA¶
Analysis of Variance models
Examples¶
In [1]: import statsmodels.api as sm
In [2]: from statsmodels.formula.api import ols
In [3]: moore = sm.datasets.get_rdataset("Moore", "car",
...: cache=True) # load data
...:
In [4]: data = moore.data
In [5]: data = data.rename(columns={"partner.status":
...: "partner_status"}) # make name pythonic
...:
In [6]: moore_lm = ols('conformity ~ C(fcategory, Sum)*C(partner_status, Sum)',
...: data=data).fit()
...:
In [7]: table = sm.stats.anova_lm(moore_lm, typ=2) # Type 2 ANOVA DataFrame
In [8]: print(table)
sum_sq ... PR(>F)
C(fcategory, Sum) 11.614700 ... 0.759564
C(partner_status, Sum) 212.213778 ... 0.002874
C(fcategory, Sum):C(partner_status, Sum) 175.488928 ... 0.022572
Residual 817.763961 ... NaN
[4 rows x 4 columns]
A more detailed example can be found here: