我正在使用statmodel OLS进行迭代离群值消除。我已经使用了模型。
ols_result = sm.OLS(y,X).fit()
然后我可以使用]获得外部和bonferroni的学生化删除残差>
ols_result.outlier_test(method="bonf")
我要删除的样本超过bonferroni p>%10,在该样本中,每次迭代的烹饪距离也是最高的。直到我没有bonf(p)>%10的样本,然后我才得到原始样本的子集。
假设我有400个样本,而在异常删除后,我得到了380个样本。现在我想再次找到学生化的缺失残差和bonferroni,以及针对380个样本拟合的400个样本。查看删除的离群值是否真的是离群值。
这是问题开始的地方。我一直在寻找使用statmodels OLS模型获取拟合值的残差和库克距离的简便方法,而不是自己编写这些函数。但是.outlier_test()
和.get_influence()
似乎可以与OLS Result对象一起使用。
你们有没有任何简单的方法就可以在没有太多代码的情况下完成这些测试?
提前感谢
我正在使用statmodel OLS进行迭代离群值消除。我已经使用了模型。 ols_result = sm.OLS(y,X).fit()然后我可以从外部获得学生化的删除残差,然后...
经过深入研究,得出了这种实用性。我找不到使用statmodels库执行此操作的简便方法。相反,我编写了自己的代码,还从statmodels源代码获得帮助。如果有人需要我在下面共享的代码