标签：plt python DBETA DFBETA 序列图 sm np import dt

数据来自Regression Analysis by Example第12章的表12.1，一个破产概率的估计的例子。

检查了一下发现statsmodels的Logic中还没有诊断工具，但可以调用GLM，所以就实现了一下。

先加载上需要用到的包

import pandas as pd
import numpy as np
import statsmodels.api as sm
import matplotlib.pyplot as plt

读取数据

dt = pd.read_csv('D://python/regression/P322.csv')
print(dt)

作图

x = dt[['X1', 'X2', 'X3']]
y = dt[['Y']]

X1 = sm.add_constant(x)
m1 = sm.GLM(y, X1, family=sm.families.Binomial(sm.families.links.logit())).fit()
outliers = m1.get_influence(observed=False)
DBETA = outliers.dfbetas
DBETAi = []
for i in range(0, len(y)):
    DBETAi.append(np.maximum(np.sum(DBETA[i]), -np.sum(DBETA[i])))

index = []
for i in range(1, len(y)+1):
    index.append(i)

plt.figure(figsize=(6, 5))
plt.scatter(index, DBETAi)
plt.xlabel('Index', fontsize=20)
plt.ylabel('$DBETA_i$', fontsize=20)
n = np.array([9, 14, 52, 53])
plt.show()

 DFBETA序列图就完美呈现了。

标签：plt,python,DBETA,DFBETA,序列图,sm,np,import,dt
来源： https://blog.csdn.net/weixin_61532136/article/details/121714384

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。