Python：使用statsmod找到logit回归模型的拐点

Question

我在Python中使用statsmodel：

import statsmodels.formula.api as smf

我有一个数据帧，有x1值和d值。d列是我的因变量。我画出d与x1，得到：

现在，我使用statsmodel拟合logit回归模型：

^{pr2}$

结果是：

我如何得到那个logit曲线的拐点的值？在

谢谢！在

Answer 1

import pandas as pds
df1  = pds.DataFrame({'lgt' : logithat,'x':x1}).set_index('x')
grid = np.linspace(df1.index.min(),df1.index.max(),df1.shape[0]+1)
df2 = pds.DataFrame({'grid':grid}).reset_index().set_index('grid')
df2 = pds.concat([df1,df2],join='outer',axis=1)
df2['lgt_'] = df2['lgt'].interpolate()


inflexion       = df2['lgt_'].diff().argmax()

fig = plt.figure()
ax = fig.add_subplot(111)
df2['lgt_'].plot(ax=ax)
ax.vlines(inflexion,df2['lgt_'].min(),df2['lgt_'].ix[inflexion])