Найти значение F-статистики для первой линейной модели ANOVA

Question

Я просматриваю статистику по темам Python. Я поражен одной рукой.

Постановка задачи:

Выполнить ANOVA для первой линейной модели, полученной при работе с набором данных mtcars. Отображение значения F-статистики .

Что я сделал для постановки задачи:

import statsmodels.api as sm    
import statsmodels.formula.api as smf    
from statsmodels.stats import anova

mtcars_data = sm.datasets.get_rdataset("mtcars").data    
print(mtcars_data.columns)

mt_model1 = smf.ols('mpg ~ cyl', mtcars_data ).fit()    
print(anova.anova_lm(mt_model1))

Как отобразить статистику F для вышеуказанной проблемы?

Answer 1

Поскольку вы уже установили модель с необходимыми переменными в mt_model1, вы можете напрямую вызвать F-статистику с помощью,

print(mt_model1.fvalue)

Это также можно использовать, когда в вашей модели есть несколько предикторов.

Answer 2

После многих испытаний это получилось.

    import statsmodels.api as sm
    from statsmodels.formula.api import ols

    mtcars = sm.datasets.get_rdataset('mtcars').data

    lm = ols('mpg ~ wt', mtcars).fit()
    av = sm.stats.anova_lm(lm,type=2)
    print(av.F.wt)

Answer 3

Так что, если вы хотите получить значение F Statistic из таблицы anova для атрибута cyl, что-то вроде этого

print(anova.anova_lm(mt_model1).F["cyl"])