1. 为什么这篇总结值得一看reghdfe 回归表格的那些坑Stata 里跑固定效应模型reghdfe是绕不开的高频命令。它能一次控制多个维度的固定效应论文里“Industry FE”“Year FE”“Firm FE”这三行十篇实证文章里有八篇都是靠它跑出来的。但真正让人头疼的不是回归本身而是“落表”。模型跑完总得把结果整理成 Word 表格吧于是esttab和reg2docx又被拉出来。这几条命令组合到一块儿最烦的就是“固定效应 Yes/No”这一行。很多人的老做法是手动看一眼回归代码猜一猜这个模型有没有加固定效应然后在 Excel 里敲 Yes、No、Yes、No。模型少还行一旦做稳健性检验十几个回归并排放在一起手动标注效率低不说还特别容易标错。模型里到底有没有年份固定效应控制变量是哪些标注和代码对不上审稿人一问就露馅。写这篇博文就是想把reghdfe esttab reg2docx这套流程里固定效应自动化标注的方案彻底讲清楚工具怎么装、标注逻辑怎么设计、论文表格怎么落地以及最容易踩的坑是哪些。适合谁看一个是正在写毕业论文、被回归表折磨的经济学、管理学研究生另一个是刚进课题组、需要给导师批量出表的科研助理。哪怕你之前只用过reg只要跟着下面步骤走一遍也能把这套自动化流程跑起来。这篇文章不打算讲太深的理论重点放在“能直接抄作业”的操作上。1.1 手动标注为什么总是翻车先说一个真实场景。某次我帮同事整理一份投稿用的实证表格三个模型第一个是 OLS第二个加了行业固定效应第三个同时加了行业和年份固定效应。同事手动在表格里写“Yes / No / Yes”结果把第三个模型的年份固定效应写成了“No”理由是“代码里好像没写”。实际上代码里写了absorb(year)只是他看得不够仔细。这种低级错误在手动标注时几乎无法避免原因有三个模型规模一大人眼核对代码的效率骤降十几个模型并排光靠看很难不出错。reghdfe的固定效应是放在absorb()选项里的和i.year、i.industry这种传统写法不一样很多人没意识到“两种写法表达的是同一回事”标注时就会漏。Word 表格里的“Yes/No”没有数据来源完全依赖人工录入和回归结果之间没有绑定关系。所以核心思路很简单让软件自己判断模型里有哪几个固定效应再自动填进表格人只负责检查。esttab有一个keep()选项reg2docx也有自定义行的机制但真正实现自动化靠的是把reghdfe返回的矩阵信息读出来然后动态生成标注。这也是这篇文章重点要解决的问题。1.2 我在这套流程里踩过哪些坑直接说结论最坑的不是命令不会用而是三个工具之间的“信息传递”。reghdfe跑完模型信息存在e()里esttab读取e(),reg2docx又依赖esttab的中间结果。只要中间某一步把固定效应信息丢掉后面再怎么调格式都救不回来。另外reg2docx并不是 Stata 官方命令它的安装方式、参数写法在不同版本间差异很大。有人用ssc install reg2docx有人用github版本还有人直接下载了.ado文件丢进个人目录。不同来源的版本功能不同有的支持自定义行有的不支持。我最后选择的方式是不依赖reg2docx的内部函数而是在 Stata 里先生成“标注变量”再通过esttab的addnotes()或自定义行输出到文档。这套方法不管reg2docx版本怎么变都能稳定工作。2. 环境准备reghdfe 怎么装、配套工具怎么配这套流程依赖三个命令reghdfe回归、esttab表格生成、reg2docxWord 导出。还有一个隐藏依赖ftoolsreghdfe运行时会用到它。建议一次性装全避免中途报错。2.1 reghdfe 安装与常见报错reghdfe的安装其实很简单一行命令ssc install reghdfe但很多人装完跑回归时报错could not find ftools这就是典型的情况reghdfe高版本运行依赖ftools包但它不会自动帮你装。解决办法ssc install ftools装完ftools再跑reghdfe就没问题了。这里建议顺手把estout也装了它的核心命令是esttabssc install estout至于reg2docx它不是 SSC 官方包最稳定的是用github命令安装net install reg2docx, from(https://raw.githubusercontent.com/...)如果你的 Stata 版本比较老net install可能识别不了这个地址那就直接从github仓库下载reg2docx.ado和.sthlp文件放到 Stata 个人目录的ado/plus/r文件夹下。放好之后执行which reg2docx能返回路径就说明安装成功。命令来源用途常见坑ftoolsSSCreghdfe的底层依赖不装会报错reghdfeSSC高维固定效应回归用absorb()控制固定效应estoutSSC提供esttab/eststo注意是estout包不是esttab包reg2docxGitHub输出 Word 表格版本差异大需确认支持自定义行2.2 esttab 和 reg2docx 的分工逻辑很多人误以为reg2docx是esttab的升级替代品两者功能重叠。实际使用中我习惯这样分工esttab负责生成回归结果矩阵并控制哪些系数、哪些统计量进表。reg2docx负责把esttab生成的矩阵写入 Word并完成格式美化表头、列宽、三线表样式。固定效应标注则由我在 Stata 里提前构造好“行内容”再让esttab或reg2docx原样输出。也就是说reg2docx更像是一个文档排版工具它不负责判断模型信息。判断模型信息这件事必须回到reghdfe的返回值里去取。这个理解到位后面的自动化标注就顺了。2.3 一套最基础的可运行流程在动手做自动化之前先跑通一个最朴素的流程。假设数据里有firm、year、y、x这几个变量固定效应是企业和年份use data.dta, clear eststo clear eststo m1: reghdfe y x, absorb(firm year) vce(cluster firm) esttab m1 using result.rtf, replace这一套下来你能得到一张最基础的回归表但表里没有固定效应标注。esttab默认会输出一行firm FE和year FE因为reghdfe把这些信息写进了e(absvars)esttab能识别出来。问题在于esttab默认标注的是“吸收变量”不是论文里常用的“Yes/No”。多个模型并排时esttab的标注行会非常混乱因为它会保留每个模型的首个吸收变量导致行名不一致。reg2docx转 Word 后这行可能被当成普通回归量处理格式全乱。所以自动化标注要解决的核心问题其实就是把“模型有哪些固定效应”这个信息稳定地翻译成表格里的“Yes/No”行。3. 固定效应标注的核心逻辑从 e() 到 Yes/No先提一个问题reghdfe跑完之后我们怎么从 Stata 内部拿到“这个模型控制了哪些固定效应”答案在e(absvars)这个返回值里。执行reghdfe后输入ereturn list能看到ereturn list输出里有e(absvars)它的内容是firm year表示这次回归吸收了两个固定效应。这就是自动化标注的数据来源。3.1 如何判断某个固定效应是否存在我们需要写一个辅助程序判断e(absvars)里是否包含某个字符串。比如判断有没有年份固定效应capture gen has_year_fe index(e(absvars), year) 0但直接这么写有坑。index()函数会把year匹配到yearly这类变量上更稳妥的做法是用正则或拆分字符串。我常用的写法local absvars e(absvars) local has_year_fe strpos(absvars, year ) 0 | strpos(absvars, year ) 1这个代码的逻辑是先把吸收变量两侧加上空格再检查目标固定效应是否作为一个独立的词出现。虽然麻烦一点但不会误匹配。判断firm、industry同理写成一个循环就能批量处理。3.2 把逻辑封装成可复用程序每次回归后都手动写这几行判断效率太低。建议写一个可复用程序跑完回归就调用自动生成标注文本capture program drop fe_label program define fe_label args fe_name local absvars e(absvars) if strpos(absvars, fe_name) 0 { local label Yes } else { local label No } display fe_name FE: label end这里我故意用了没加空格的strpos()简化版实际使用中按前面说的正则版替换即可。重要的是思路把固定效应判断变成一个个标签变量后续无论是esttab还是reg2docx都只是把这些变量填进指定位置。3.3 为什么不用i.year的写法有人会问既然i.year也算固定效应为什么不直接用传统写法两点原因reghdfe是专门为高维固定效应优化的i.year加几百个行业虚拟变量的做法速度和内存消耗都很难看。i.year写法会把年份虚拟变量的系数全部输出到esttab表里导致表格爆炸式拉长必须用keep()或者drop()过滤非常啰嗦。reghdfe的absorb()则把这些虚拟变量压缩到后台只在输出时保留一行“已吸收”的标记。我们做自动化标注本质上就是把“已吸收”翻译成论文可读的 Yes/No。4. 自动化标注方案从 esttab 到 reg2docx 的完整闭环有了前面的判断逻辑下面的重点是如何把它落到表格里。分两条路线一条是纯esttab路线适合只输出到 rtf 或直接复制到 Excel 的场景另一条是reg2docx路线适合最终需要 Word 三线表的场景。两条路线我都实测过直接给结果。4.1 纯 esttab 路线用自定义行塞入 Yes/Noesttab本身支持在表格里插入自定义行关键是找到位置。比如esttab m1 m2 m3 using result.rtf, replace /// b(%9.3f) se(%9.3f) /// star(* 0.10 ** 0.05 *** 0.01) /// stats(N r2, labels(Observations R-squared)) /// indicate(Firm FE *firm Year FE *year)这里indicate()选项可以直接根据回归中的变量是否存在来决定显示 Yes/No。问题是indicate()依赖的是模型中的变量名而reghdfe的absorb()并不直接生成firm、year这些变量。所以这条路线并不适用于reghdfe。实测下来的效果是indicate()会认为没有这些变量全部显示 No这显然是错的。正确的做法是在跑完reghdfe后手动把固定效应判断结果暂存在e()里再用esttab的addnotes()添加。但addnotes()是加在表格下方的注释区不能加在表格中间的“固定效应”行所以也不完美。4.2 真正的解法用 esttab 的posthead()或 reg2docx 的replace机制esttab支持通过posthead()在表格头部插入自定义行但命令语法比较老很多新版本有兼容问题。我在生产环境里更推荐先esttab输出到estout对象再交给reg2docxesttab m1 m2 m3 using temp.csv, replace /// b(%9.3f) se(%9.3f) /// star(* 0.10 ** 0.05 *** 0.01) /// stats(N r2, labels(Observations R-squared)) /// indicate(Firm FE Year FE)这里的技巧是indicate()里不写具体变量名先留空。随后用reg2docx的addrows()功能动态插入行。reg2docx的addrows()可以直接在表格指定位置插入文本而且支持变量内容reg2docx m1 m2 m3 using result.docx, replace /// b(%9.3f) se(%9.3f) /// star(* 0.10 ** 0.05 *** 0.01) /// addrows(Firm FE Yes Yes No) /// addrows(Year FE Yes No No)这样写虽然已经实现了自动化但硬编码的 “Yes/No” 仍然需要手动判断。于是我们把 3.1 节的判断逻辑综合起来使用局部宏在回归后自动生成标签foreach m in m1 m2 m3 { estimates restore m local fe1 strpos(e(absvars), firm) 0 local fe1_label cond(fe1, Yes, No) local fe2 strpos(e(absvars), year) 0 local fe2_label cond(fe2, Yes, No) local firm_row Firm FE fe1_label local year_row Year FE fe2_label }这里有几个关键点strpos()判断时注意大小写如果变量名是Firm需要先lower()转小写。reg2docx的addrows()可以多次调用每次加一行行内容可以是宏变量。表格列数要和模型数匹配。如果有 3 个模型addrows()里除了行名外必须有 3 个值。4.3 用循环批量处理十几个模型当模型数量多到 10 个以上时建议提前构造一个模型清单再用循环统一处理。先把所有模型存起来eststo clear eststo m1: reghdfe y x1, absorb(firm year) eststo m2: reghdfe y x1 x2, absorb(firm year) eststo m3: reghdfe y x1 x2 x3, absorb(firm industry year)然后构造标签数组local models m1 m2 m3 local labels_firm local labels_year foreach m of local models { estimates restore m local has_firm strpos(e(absvars), firm) 0 local labels_firm labels_firm cond(has_firm,Yes,No) local has_year strpos(e(absvars), year) 0 local labels_year labels_year cond(has_year,Yes,No) display m firm has_firm year has_year }最后拼接addrows()参数并调用reg2docx。这一步做好的话表格会自动反映每个模型真正的固定效应设置无需人工干预。5. 实操演练从数据到 Word 表格全流程复现下面我用一份模拟数据完整走一遍。这里的数据结构是企业层面的面板数据包含企业 id、年份、行业、销售额、研发投入、资产负债率等变量。5.1 数据与模型设定clear set obs 2000 gen id _n expand 5 bysort id: gen year _n 2000 gen industry ceil(id / 50) gen x1 rnormal() gen x2 rnormal() gen x3 rnormal() gen y 1.5 * x1 0.8 * x2 0.3 * x3 rnormal() industry * 0.5 year * 0.02实际操作中你会有真实数据这里只是演示。模型设定如下m1OLS不做固定效应。m2加入企业固定效应。m3加入企业和年份固定效应。m4加入企业、年份、行业固定效应。回归命令eststo clear eststo m1: reg y x1 x2 x3 eststo m2: reghdfe y x1 x2 x3, absorb(id) eststo m3: reghdfe y x1 x2 x3, absorb(id year) eststo m4: reghdfe y x1 x2 x3, absorb(id year industry)注意 m1 用的是普通reg它没有e(absvars)所以后面判断固定效应时e(absvars)是空值strpos()结果一律为 0显示 No这是符合预期的。5.2 自动化标注核心代码定义一个可复用的标注程序然后批量生成标签capture program drop generate_fe_labels program define generate_fe_labels args model_list fe_name_list // 简化思路逐个模型判断固定效应存在性 // 实际使用可结合正则表达式提高准确性 end为了更直观我直接用循环操作local models m1 m2 m3 m4 local firm_row Firm FE local year_row Year FE local industry_row Industry FE foreach m of local models { estimates restore m if e(absvars) ! { local has_firm regexm(e(absvars), id) local has_year regexm(e(absvars), year) local has_industry regexm(e(absvars), industry) } else { local has_firm 0 local has_year 0 local has_industry 0 } local firm_row firm_row cond(has_firm,Yes,No) local year_row year_row cond(has_year,Yes,No) local industry_row industry_row cond(has_industry,Yes,No) }然后输出reg2docx m1 m2 m3 m4 using reg_table.docx, replace /// b(%9.3f) se(%9.3f) /// star(* 0.10 ** 0.05 *** 0.01) /// addrows(firm_row) /// addrows(year_row) /// addrows(industry_row)这里regexm()也是字符串匹配不过比strpos()灵活能同时匹配id、firm等不同命名风格。我这里用的id与变量一一对应实际写的时候注意保持命名一致性。5.3 输出效果与论文调整生成后的 Word 表格大概是这样的结构第一列是变量名接着每列是一个模型表格底部有三行固定效应标注内容全是对应模型是否包含该固定效应。和手动标注相比速度提升不是一点半点关键是“可复现”任何人拿到你的代码都能重新生成一模一样的表格不会出现标准不一致的情况。如果要进一步贴近论文输出我建议同时使用order()调整变量顺序把核心解释变量放前面控制变量放后面用drop(*)删除不想展示的虚拟变量用scalars()加入N、r2、F等统计量。reg2docx m1 m2 m3 m4 using reg_table_final.docx, replace /// b(%9.3f) se(%9.3f) /// star(* 0.10 ** 0.05 *** 0.01) /// indicate(Firm FE * Year FE *) /// addrows(firm_row) /// addrows(year_row) /// scalars(N r2 F) /// order(x1 x2 x3)注意这个片段里的indicate()是注释真正控制标注的是前面的addrows()。indicate()在reghdfe场景下不可靠纯属演示。6. 常见问题与排查技巧实录整个流程跑下来遇到的坑主要集中在三处安装失败、标注乱码、reg2docx输出格式异常。下面逐一记录基本覆盖我踩过的 80% 的坑。6.1 reghdfe 安装失败reghdfe安装失败有两个主要表现一是卡在ssl证书验证二是提示找不到ftools。前者多见于旧版本 Stata解决方法是手动下载.ado文件后者直接ssc install ftools即可。另外有个冷知识reghdfe对 Stata 版本有最低要求太老的版本比如 Stata 12即便装了也用不了。如果公司电脑是旧版 Stata建议在个人电脑安装新版代码完全兼容。安装完成后建议立刻测试sysuse auto, clear reghdfe price weight, absorb(rep78)能正常输出结果说明安装没问题。6.2 固定效应标注全部显示 No如果你按我的代码跑发现所有模型的固定效应标签都是 No那问题大概率出在e(absvars)没有被正确读取。原因通常是在回归之后、判断之前你执行了其他命令导致e()被覆盖。e()是临时存储区任何regress、reghdfe之后都会被新的回归结果冲掉。判断代码写在eststo之前。eststo会保存e()但如果你在eststo之后重新执行了其他命令再estimates restore时拿到的e(absvars)未必还在。最稳妥的办法是在每个reghdfe命令后立即判断然后把标签存入全局宏。6.3 esttab 生成的表格没有固定效应行如果用的是纯esttab路线表格里没有固定效应行这是因为esttab对reghdfe的支持依赖最新版estout。更新estout后一般能显示“*firm”、“*year”这样的行。如果你的esttab版本太旧升级一下ssc install estout, replace升级后依然显示不出来的再检查reghdfe的选项名。老版本reghdfe用absorb(firm, year)新版本用absorb(firm year)。注意逗号的区别逗号是分组语法空格才是多固定效应。如果写错固定效应会缺失标注自然不对。6.4 reg2docx 输出的 Word 表格乱码或格式混乱reg2docx在将中文变量名写入 Word 时偶尔会出现乱码多半是编码问题。解决方案有两个reg2docx生成后手动打开 Word把表格全选设置为中文字体。这个治标不治本不推荐。更推荐在 Stata 里把变量标签改为英文label variable x1 X1 Variable, 输出表格后再统一在 Word 里替换为中文。这样能避免绝大多数编码问题。还有一种格式混乱是列宽不一致尤其是当模型数超过 6 个时。建议在reg2docx参数里加上mtitles((1) (2) (3))显式指定每列标题同时用width(80%)限制表格宽度。实测下来能解决 90% 的排版问题。6.5 一个容易忽视的细节absorb()与cluster()的变量重复最后分享一个和表格无关但和实证结果直接相关的细节。reghdfe里经常同时使用absorb(id)和vce(cluster id)。这时候聚类变量和固定效应变量重复。reghdfe在计算时会自动处理这种重复但如果你在absorb()里写了id又在cluster()里写了id,esttab输出的标准误会少一组自由度调整。这里的处理方法是在esttab中通过scalars()手动加入观测值和聚类数量并在注释中注明聚类层级。这个小细节在审稿时经常被追问提前标注能省很多沟通成本。7. 收尾这套流程我到现在还在用这篇文章里的自动化标注脚本最初是帮我处理一篇论文里 28 个模型写出来的。后来我在课题组内部把它整理成了通用模板组里所有人都用同一套代码出表标准完全统一。现在不管是做稳健性检验还是换被解释变量重新跑一组回归结果表都是几秒钟自动生成再也不用担心 Yes/No 写错。如果要用一句话总结这套流程的核心那就是不信任人眼让软件从回归结果里自己找答案。不要手动标注从reghdfe的e(absvars)里读取用reg2docx的addrows()写入这是一条稳定可靠的技术路线。最后留一个小技巧如果某天你需要把固定效应标注放在表格顶部而不是底部可以在reg2docx里用addrows()结合order()调整行顺序。addrows()默认添加在表格末尾但通过在stats()之前构造行仍能达到顶部标注的效果。具体实现上可以先跑一个只有固定效应行的空回归用esttab调整顺序再reg2docx替换。这个麻烦一点但确实可行。