当前位置：首页 > news >正文

python实战：数据分析基础知识

news 2025/8/14 2:48:02

当涉及到数据分析和统计建模时，Python 提供了强大的工具和库，如 pandas、numpy、statsmodels 和 matplotlib。本文将以一个实际的案例为例，介绍如何利用这些工具进行回归分析，并通过可视化工具进行结果展示和解释。

1. 背景介绍

在本文中，我们将探索如何使用 Python 进行回归分析，具体来说，我们将使用的数据集包含了关于不良贷款（y）和一些可能影响不良贷款的因素（x）的信息。我们将通过简单线性回归模型来探索这些因素对不良贷款的影响程度。

2. 数据准备与加载

首先，我们需要加载数据并进行初步的检查，确保数据的完整性和正确性。

import pandas as pd
import numpy as np
import statsmodels.api as sm
import matplotlib.pyplot as plt# 设置中文显示和正常显示负号
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False# 读取数据集
example10_1 = pd.read_csv('exercise10_1.csv', encoding='gbk')# 打印数据集的前几行和列名，确保数据读取正确
print(example10_1.head())
print(example10_1.columns)

3. 数据预处理

在了解数据结构后，我们需要根据实际情况调整列名，并准备用于回归分析的自变量（X）和因变量（y）。

# 根据实际情况调整列名，这里假设实际列名为 '不良贷款','贷款余额','应收贷款','贷款项目个数','固定资产投资'
x = example10_1[['贷款余额', '应收贷款', '贷款项目个数', '固定资产投资']]  # 根据实际列名修改
y = example10_1['不良贷款']  # 根据实际列名修改# 添加常数项
X = sm.add_constant(x)# 拟合线性回归模型
model = sm.OLS(y, X).fit()# 打印回归结果摘要
print(model.summary())

4. 回归结果分析与可视化

完成模型拟合后，我们可以通过图形化方式来评估模型的适配程度和残差的分布情况。

# 创建一个包含两个子图的图像
fig, axes = plt.subplots(1, 2, figsize=(9, 4))# 绘制残差与拟合值的散点图
axes[0].scatter(model.fittedvalues, model.resid)
axes[0].set_xlabel('拟合值')
axes[0].set_ylabel('残差')
axes[0].set_title('(a) 残差值与拟合值图', fontsize=15)
axes[0].axhline(0, ls='--')# 绘制正态 Q-Q 图
sm.qqplot(model.resid, line='r', ax=axes[1])
axes[1].set_xlabel('期望正态值')
axes[1].set_ylabel('标准化的观测值')
axes[1].set_title('正态Q-Q图', fontsize=15)# 调整布局
plt.tight_layout()# 显示图像
plt.show()