基于 Python 的事件研究法

项目获取途径

• 下载下面压缩包，然后解压缩即可
💡
Event Study Analysis.zip
1.2 GiB

• 也可以点击此处通过百度网盘下载，提取码:ywuh

事件研究法（Event Study）

什么是事件研究法？

短期事件研究法依赖的三个基本假设：

• 根据有效市场假说 (Efficient Markets Hypothesis, EMH) , 金融市场是有效的, 即股票价格反映所有已知的公共信息

• 所研究的事件是市场未预期到的，因此这种异常收益可以度量股价对事件发生或信息披露异常反应的程度

• 在事件发生的窗口期间无其他事件的混合效应

事件研究法的基本步骤：

• 整理事件列表，事件列表中包含了公司名称和事件发生日期。

• 整理事件列表中的公司在事件发生日前后的收益率数据和同一时期的市场收益率数据。

• 根据模型估计事件窗口期的预期收益率并比较公司在事件窗口期的实际收益率与预期收益率的差异。

• 通过回归来检验特定事件对公司股票收益率是否有显著的影响（t检验）。

本代码相关信息

✨ 功能特点

核心功能

• 异常收益率计算：基于管理学中最常见的资本资产定价模型（CAPM）计算个股的异常收益率（AR）和累积异常收益率（CAR）

• 多进程并行计算：利用 Python multiprocessing 模块实现多核心并行处理，显著提升计算速度

• 断点续传：自动跳过已处理的事件，支持中断后继续运行，避免重复计算

• 实时进度显示：使用 tqdm 显示实时进度条，包括处理速度、成功/失败统计等信息

数据处理

• 窗口配置：支持自定义事件窗口（Event Window）和估计窗口（Estimation Window）的长度及间隔. 如果将事件日在时间轴上视为0点，那么默认事件窗口为[-3, 5]（即事件日前3与后5个交易日），默认估计窗口为[-128, -9]（即120个交易日）.
💡
当事件日当天不是股市交易日时（事件日在节假日），自动使用事件日后的第一个交易日作为0点

• 多种指数支持：默认是沪深300，可以自动根据股票代码匹配对应的市场指数（上证综指、深证综指、沪深300），也支持自定义指数

• 数据持久化：使用 SQLite 数据库存储股票日行情数据，提高查询效率

统计分析

• 统计检验：提供详细的统计检验脚本，包括：
• 检验（检验均值是否为0）
• Wilcoxon 符号秩检验（检验中位数是否为0）
• 二项式符号检验（检验正负比例是否为0.5）

• 结果可视化：生成 AR 和 CAR 的趋势图，支持显示均值和中位数

🏗️ 技术框架

核心算法

• 市场模型（Market Model）：使用 OLS 回归估计正常收益率
其中：
• ：股票 在 日的收益率
• ：市场指数在 日的收益率
• ：通过估计窗口的 OLS 回归得到
• ：残差项

• 异常收益率（AR）计算：

• 累计异常收益率（CAR）计算：

技术栈

• 数据处理：pandas, numpy

• 统计分析：statsmodels, scipy

• 数据库：SQLite3

• 并行计算：multiprocessing

• 可视化：matplotlib

• 进度显示：tqdm

📁 项目结构

🔧 环境依赖

Python 版本

• Python 3.7+

必需的 Python 库

内置模块

• sqlite3（Python 标准库）

• multiprocessing（Python 标准库）

• logging（Python 标准库）

📦 安装指南

1. 克隆或下载项目

2. 安装依赖

3. 准备数据文件

• 事件数据：Event study_data/event info.csv
• 必须包含列：company_code（公司代码）、event_date（事件时间）
• 可选列：company_year（公司年份）、industry_code（公司行业代码）、company_province（公司省份）、company_city（公司城市）（会被添加到结果中）

• 股票数据：公司数据集/日个股回报率/TRD_Daily_20150101-20250101.csv
• 从 CSMAR 数据库下载的日个股回报率数据（数据量太大，故使用csv格式加快读取速度）

• 指数数据：公司数据集/指数文件1990-12-19 - 2025-04-24/TRD_Index.xlsx
• 从 CSMAR 数据库下载的市场指数数据

🚀 使用说明

步骤 1：导入数据到 SQLite 数据库

• 读取 公司数据集/日个股回报率/TRD_Daily_20150101-20250101.csv

• 创建 SQLite 数据库：Event study_data/stock_daily_data.db

• 建立索引以优化查询性能

• 显示导入进度

• 如果数据库已存在，脚本会先删除旧数据库再创建新的

• 导入过程可能需要较长时间，取决于数据量大小

步骤 2：运行主程序计算 AR 和 CAR

• 读取事件列表：Event study_data/event info.csv

• 自动检测已处理的事件，跳过重复计算（支持断点续传）

• 使用多进程并行计算（默认使用 CPU 核心数 - 6）

• 实时显示进度条，包括处理速度、成功/失败统计

• 结果实时保存到 CSV 文件

• Event study_results/异常收益率（AR）&累积异常收益率(CAR).csv

• Event study_results/平均累计异常收益率（Mean CAR）_统计检验结果_按公司与事件日分类.csv

步骤 3：运行统计检验

• 读取 AR 和 CAR 数据

• 按相对日期（relative_day）分组计算统计量

• 进行三种统计检验：t 检验、Wilcoxon 符号秩检验、二项式符号检验

• 计算均值、中位数、正负比例等统计量

• Event study_results/异常收益率（AR）&累积异常收益率(CAR)_统计检验结果_按相对日期分类.csv

步骤 4：生成可视化图表（可选）

• 绘制平均异常收益率（AAR）和累积平均异常收益率（CAAR）的趋势图

• 支持显示均值和中位数（可在脚本中配置）

• 使用科研风格的图表样式

• Event study_results/异常收益率（AR）&累积异常收益率(CAR)_graph.png

⚙️ 参数配置

EventStudy 类参数

• event_window_before：事件窗口开始日相对于事件日的天数（负数，如 -3 表示事件日前 3 天）

• event_window_after：事件窗口结束日相对于事件日的天数（正数，如 5 表示事件日后 5 天）

• estimation_window_length：估计窗口的长度（交易日数）

• estimation_window_gap：估计窗口结束日与事件窗口开始日之间的间隔天数

多进程配置

• 根据 CPU 核心数调整保留的核心数

• 如果系统资源充足，可以减少保留的核心数

• 如果系统资源紧张，可以增加保留的核心数

市场指数配置

• '000300'：沪深300指数（默认）

• '000001'：上证综合指数

• '399106'：深证综合指数

• None：自动根据股票代码匹配对应的市场指数

📊 输出结果

1. AR 和 CAR 数据文件

• company_code：公司代码

• event_date：事件日期

• relative_day：相对日期（事件日为 0，事件日前为负数，事件日后为正数）

• abnormal_return：异常收益率（AR）

• cumulative_abnormal_return：累积异常收益率（CAR）

• 公司年份、公司行业代码、公司省份、公司城市：从事件数据中继承的额外信息

2. 统计检验结果（按相对日期分类）

• AR_Mean / CAR_Mean：均值

• AR_Median / CAR_Median：中位数

• AR_Min / CAR_Min：最小值

• AR_Max / CAR_Max：最大值

• AR_t_stat / CAR_t_stat： 检验统计量

• AR_p_value / CAR_p_value： 检验 值

• AR_w_stat / CAR_w_stat：Wilcoxon 检验统计量

• AR_w_p_value / CAR_w_p_value：Wilcoxon 检验 值

• AR_Positive_Ratio / CAR_Positive_Ratio：正收益率比例

• AR_Negative_Ratio / CAR_Negative_Ratio：负收益率比例

• AR_Binom_p_value_Positive / CAR_Binom_p_value_Positive：二项式检验 值（正收益率）

• AR_Binom_p_value_Negative / CAR_Binom_p_value_Negative：二项式检验 值（负收益率）

3. 统计检验结果（按公司与事件日分类）

• company_code：公司代码

• event_date：事件日期

• mean_CAR：平均累积异常收益率

• t_stat：t 检验统计量

• p_value：t 检验 p 值

• significant：是否显著（基于 值 < 0.05）

4. 可视化图表

• 上图：平均异常收益率（AAR）趋势

• 下图：累积平均异常收益率（CAAR）趋势

⚠️ 注意事项

数据要求

1. 事件数据格式：
• CSV 文件必须包含 company_code（公司代码）与event_date（事件时间）列
• 事件时间 格式应为 YYYY-MM-DD 或 YYYY/MM/DD的datetime格式

2. 股票数据要求：
• 数据应包含股票代码、交易日期、收益率等必要字段
• 确保数据时间范围覆盖所有事件日期

3. 指数数据要求：
• 确保指数数据的时间范围覆盖所有事件日期
• 指数代码应与股票代码匹配（沪市股票对应上证指数，深市股票对应深证指数）

运行环境

1. Windows 系统：
• 使用 multiprocessing.freeze_support() 确保多进程正常工作
• 确保在 if __name__ == '__main__': 块中运行主程序

2. 内存管理：
• 大规模数据处理时，注意内存使用情况
• 程序采用实时写入 CSV 的方式，避免内存溢出

3. 日志文件：
• 日志文件保存在 Event study_code/Event study program/event_study.log
• 如果程序异常，可以查看日志文件获取详细信息

性能优化

1. 多进程设置：
• 根据 CPU 核心数合理设置并行进程数
• 过多进程可能导致系统资源竞争，反而降低效率

2. 数据库索引：
• SQLite 数据库已建立索引，查询速度较快
• 如果数据量特别大，可以考虑使用更专业的数据库（如 PostgreSQL）

3. 断点续传：
• 程序支持断点续传，可以安全地中断和重启
• 已处理的事件会自动跳过，不会重复计算

❓ 常见问题

Q1: 程序运行很慢怎么办？

• 检查是否使用了多进程（查看日志中的进程数）

• 减少保留给系统的 CPU 核心数

• 确保数据已导入 SQLite 数据库（而不是直接读取 CSV）

Q2: 某些事件计算失败，显示"数据获取失败"？

• 检查股票代码是否正确

• 确认股票数据的时间范围是否覆盖事件日期

• 检查指数数据是否完整

• 查看日志文件获取详细错误信息

Q3: 如何修改事件窗口和估计窗口？

• 在 process_event_wrapper 函数中修改 EventStudy 类的初始化参数

• 注意：修改后需要重新运行程序，已处理的事件不会自动更新

Q4: 如何更换市场指数？

• 在 process_event_wrapper 函数中修改 custom_index 参数

• 或者设置为 None 以自动匹配市场指数

Q5: 程序中断后如何继续？

• 直接重新运行程序即可

• 程序会自动检测已处理的事件并跳过

• 已处理的结果会保留在 CSV 文件中

Q6: 如何查看详细的处理日志？

• 查看 Event study_code/Event study program/event_study.log 文件

• 日志级别设置为 INFO，包含详细的处理信息

Q7: 统计检验的 p 值如何解释？

• 值 < 0.05：在 5% 的显著性水平下拒绝原假设（认为异常收益率显著不为 0）

• 值 < 0.01：在 1% 的显著性水平下拒绝原假设（高度显著）

• 值 >= 0.05：不能拒绝原假设（异常收益率不显著）