ANOVA

線形 OLSModel による ANOVA 分析のための anova_lm を含む分散分析モデルと、バランスの取れたデータに対する ANOVA の within のための AnovaRM。

例

In [1]: import statsmodels.api as sm

In [2]: from statsmodels.formula.api import ols

In [3]: moore = sm.datasets.get_rdataset("Moore", "carData",
   ...:                                  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    df          F    PR(>F)
C(fcategory, Sum)                          11.614700   2.0   0.276958  0.759564
C(partner_status, Sum)                    212.213778   1.0  10.120692  0.002874
C(fcategory, Sum):C(partner_status, Sum)  175.488928   2.0   4.184623  0.022572
Residual                                  817.763961  39.0        NaN       NaN

anova_lm に適した詳細な例は こちら でご覧いただけます。

モジュールリファレンス

anova_lm(*args, **kwargs)

1 つ以上の適合済み線形モデルのための Anova 表。

AnovaRM(data, depvar, subject[, within, ...])

最小二乗回帰を使用した Repeated measures Anova
最終更新: Oct 03, 2024