让statsmodels使用heteroskedasticity修正系数t检验中的标准误差

Question

我一直在深入研究API,并且已经找到了如何检索不同类型的异方差校正标准误差(通过属性等等)但是,我无法弄清楚如何获得系数的t-测试使用这些更正的标准错误.有没有办法在API中执行此操作,还是必须手动执行此操作？如果是后者,您是否可以就如何使用statsmodels结果建议任何指导？ statsmodels.regression.linear_model.RegressionResultsHC0_se

Answer 1

所述fit的线性模型,离散模型和GLM的方法,采取了cov_type和cov_kwds参数用于指定鲁棒协方差矩阵.这将附加到结果实例,并用于摘要表中报告的所有推理和统计信息.

不幸的是,文档还没有以适当的方式显示这一点.实际根据选项选择三明治的辅助方法显示选项和必需参数:http: //statsmodels.sourceforge.net/devel/generated/statsmodels.regression.linear_model.OLS.fit.html

例如,HC3可以使用估计OLS模型和使用协方差矩阵

model_ols = OLS(...)
result = model_ols.fit(cov_type='HC3')
result.bse
result.t_test(....)

一些三明治需要额外的参数,例如群集健壮的标准错误,可以通过以下方式选择,假设mygroups是包含组标签的数组:

results = OLS(...).fit(cov_type='cluster', cov_kwds={'groups': mygroups}
results.bse
...

一些稳健的协方差矩阵在不进行检查的情况下对数据做出了额外的假设.例如,异方差性和自相关鲁棒标准误差或Newey-West HAC,标准误差假设顺序时间序列结构.一些面板数据稳健的标准误差也假设个人堆叠时间序列.

有一个单独的选项use_t可用于指定默认情况下是否应使用t和F或普通和chisquare分布进行Wald测试和置信区间.