2015-09-15 49 views
1

我想从一个特定的值:statsmodels.api收敛失败

result = logit.fit() 
print result.summary() 

      Logit Regression Results 
============================================================================== 
Dep. Variable:      y No. Observations:     8039 
Model:       Logit Df Residuals:      8033 
Method:       MLE Df Model:       5 
Date:    Tue, 15 Sep 2015 Pseudo R-squ.:     0.01873 
Time:      10:54:35 Log-Likelihood:    -1851.0 
converged:      True LL-Null:      -1886.4 
             LLR p-value:     7.422e-14 

我想从融合(对于以上是True) 我无法得到的值找到一种方法来从总结中获得价值。 有没有方法可以判断模型是否收敛?

回答

2

优化程序提供的附加信息存储在mle_retvals中。这应该适合你:

from statsmodels.discrete.discrete_model import Logit 
import numpy as np 
np.random.seed(42) 

n, d = 20, 10 
y = np.random.randint(2, size=n) 
X = np.random.rand(n, d) 

res = Logit(y, X).fit() 

did_converge = res.mle_retvals["converged"]