保姆级教程：用Stata处理2000-2021年A股上市公司控制变量（附完整代码与数据）

Stata实战：A股上市公司控制变量构建全流程解析

第一次接触实证研究时，最让我头疼的不是模型设定，而是数据清洗。记得研一那年，导师扔给我一份从CSMAR导出的原始数据，要求两周内完成控制变量构建。面对密密麻麻的Excel表格和缺失值警告，我在图书馆熬了三个通宵才勉强交差——如果当时有这样一份手把手教程，至少能省下50%的时间。

这份教程将用最直白的语言，带你走完从原始数据到回归-ready控制变量的完整流程。不同于单纯的数据说明文档，我们会聚焦实际操作的每个技术细节：怎么处理行业分类的特殊规则？为什么我的缩尾处理结果和文献不一致？这些在论文方法部分从不交代的"黑箱操作"，正是新手最容易踩坑的地方。

1. 数据准备与环境配置

在打开Stata之前，有几个前期准备工作直接影响后续效率。我习惯在D盘根目录创建/research/raw_data和/research/clean_data两个文件夹，分别存放原始数据和清洗后的文件。这种简单的目录管理能避免后期路径混乱——相信我，当你的do文件里出现十几种路径时，绝对会感谢这个决定。

1.1 原始数据检查

从CSMAR/Wind导出的Excel通常包含这些关键表：

• 基础财务表（BalanceSheet）
• 公司治理表（CorporateGovernance）
• 股票交易表（StockTrading）
• 行业分类表（IndustryClassification）

用这个命令快速查看数据结构：

import excel using "raw_data/balance_sheet.xlsx", firstrow clear describe list in 1/5

常见问题排查表：

1.2 Stata环境预设

这些初始化设置能避免90%的常见报错：

set more off // 取消分页暂停 set varabbrev off // 禁用变量缩写 set excelxlsxlargefile on // 处理大Excel文件

内存优化技巧：

• 对于2000-2021年A股数据（约50万条观测值），建议设置：

clear all set maxvar 30000 // 最大变量数 set matsize 11000 // 矩阵维度

2. 数据清洗核心步骤

2.1 缺失值与异常值处理

金融数据的缺失往往有规律可循：

// 标记财务数据缺失 foreach var of varlist 资产总计 负债合计 营业收入 { gen miss_`var' = missing(`var') } // 按年度统计缺失比例 table year, c(mean miss_资产总计 mean miss_负债合计)

行业处理需要特别注意：

// 剔除金融行业(J开头的代码)和ST公司 drop if substr(行业代码,1,1)=="J" drop if strmatch(证券名称,"*ST*") | strmatch(证券名称,"ST*")

提示：证监会2012版行业分类中，制造业使用两位代码（如C27），其他行业用一位字母代码。合并时需特别注意这个差异。

2.2 关键变量生成

以最常用的Size变量为例：

// 总资产自然对数 gen Size = ln(资产总计) // 分年度行业调整版本 bys year 行业代码: egen ind_avg_size = mean(Size) gen adj_Size = Size - ind_avg_size

资产负债率需要处理分母为零的情况：

gen Lev = 负债合计/资产总计 replace Lev = . if 资产总计<=0 // 处理异常分母

2.3 缩尾处理的艺术

很多论文对缩尾处理语焉不详，实际操作中要注意：

// 按年度-行业分组缩尾（1%水平） foreach var of varlist Size Lev ROA { bys year 行业代码: egen p99 = pctile(`var'), p(99) bys year 行业代码: egen p1 = pctile(`var'), p(1) replace `var' = p99 if `var' > p99 & !missing(`var') replace `var' = p1 if `var' < p1 & !missing(`var') }

常见误区对比：

• 全局缩尾 vs 分组缩尾
• 缩尾比例选择（1% vs 5%）
• 是否对负值变量缩尾

3. 面板数据架构搭建

3.1 长面板转换

从宽格式（每个公司一行）转为长格式（公司-年度观测值）：

reshape long 资产总计 负债合计, i(stkcd) j(year)

检查平衡面板：

xtset stkcd year xtdes // 显示面板结构

3.2 行业-年度效应生成

创建虚拟变量时，这个技巧能节省大量时间：

// 制造业细分到二位代码，其他行业用一位 gen ind_code = substr(行业代码,1,1) replace ind_code = substr(行业代码,1,2) if substr(行业代码,1,1)=="C" // 生成虚拟变量 tab ind_code, gen(IND_)

4. 最终数据校验与输出

4.1 描述性统计检查

这些指标必须人工验证：

tabstat Size Lev ROA, stats(mean sd p50 min max N) by(year)

与文献对比参考值：

4.2 数据保存规范

使用这种命名规则便于版本管理：

save "clean_data/A股控制变量_$DATE.dta", replace

同时输出代码本：

label data "A股上市公司2000-2021控制变量面板" notes: 最后更新日期 $DATE

5. 效率提升技巧

5.1 自动化do文件结构

一个标准的处理流程应该包含：

/*=========================================== 项目：A股控制变量构建 作者：YourName 日期：$DATE ===========================================*/ // 第一节：环境设置 version 17 set more off ... // 第二节：数据导入 import excel "raw_data/balance.xlsx", firstrow clear ... // 第三节：变量生成 gen Size = ln(资产总计) ... // 第四节：结果输出 save "clean_data/final.dta", replace

5.2 常见报错解决方案

遇到这些问题时不要慌：

• "variable not found"：检查describe确认变量名拼写
• "no observations"：用browse查看数据过滤条件是否过严
• "matrix too big"：调整set matsize数值

最后分享一个血泪教训：永远在do文件开头注明数据来源和处理日期。三个月后当审稿人问起某个变量的计算方式时，这行注释能救你的命。