OLS Breusch языческий тест в Python


Я использовал пакет statsmodels для оценки моей регрессии OLS. Теперь я хочу Breusch Pagan test. Я использовал пакет pysal для этого теста, но эта функция возвращает ошибку:

import statsmodels.api as sm
import pysal

model = sm.OLS(Y,X,missing = 'drop')
rs = model.fit()
pysal.spreg.diagnostics.breusch_pagan(rs)

Возвращенная Ошибка:

AttributeError: объект' OLSResults 'не имеет атрибута' u '

Что мне делать?

1 2

1 ответ:

Проблема заключается в том, что экземпляр результатов регрессии statsmodels не совместим с экземпляром в pysal.

Можно использовать breushpagan из statsmodels, который принимает остатки OLS и кандидатов в объясняющие переменные для гетероскедастичности и поэтому не зависит от конкретной модели или реализации модели.

Документация: http://statsmodels.sourceforge.net/devel/generated/statsmodels.stats.diagnostic.het_breushpagan.html

С примерами здесь http://statsmodels.sourceforge.net/devel/examples/notebooks/generated/regression_diagnostics.html

Я не знаю, есть ли какие-либо существенные различия в реализации теста Бреуша-Пагана.

Похоже, что название в statsmodels написано с ошибкой.