Stata：毕业论文大礼包 A——实证结果输出命令大比拼

Stata 连享会   主页 || 视频 || 推文

温馨提示： 定期 清理浏览器缓存，可以获得最佳浏览体验。

✌ 课程详情： https://gitee.com/arlionn/Course   |   lianxh.cn

✨ 课程主页： https://gitee.com/arlionn/Course

作者： 王美庭 (中南民族大学经济学院)
Email： 2017110097@mail.scuec.edu.cn

目录

啥也不说，上图……

Stata 结果输出系列推文：

• Stata：毕业论文大礼包 A——实证结果输出命令大比拼
• Stata：毕业论文大礼包 B——神速实证结果输出之搜狗短语

1. 本文目的

目前 Stata 有着众多的实证结果输出命令，连享会对于 asdoc、xxx2docx 系列、outreg2、logout 和 esttab 都有讲解，不熟悉的同学可以先去了解。然而，这么多的命令，各有其特点，往往让我们眼花缭乱。到底哪一种真正适合自己？这也许是很多同学都会提出的问题。

我们的测评结果呈现于文首的表格中，下文展开详细分析。

一般而言，我们将实证结果导入到 Word 的最多，如果要做学术报告或者发 SSCI，那可能还需要用到 LaTeX，所以本文以 Word 和 LaTeX 为基础，目的在于：从我们实证论文常用到的描述性统计、分组 T 均值检验、相关系数矩阵和回归结果四个方面出发，去评判以上几种常见输出命令的优劣，最后再给出一个总结。

整体说明和注意事项：

由于以上所有命令都支持 Word 输出，而只有部分命令支持 LaTeX 输出，且 LaTeX 输出命令中只有少量与 Word 输出的命令不一致，所以我们在 3 - 6 节均以输出到 Word 的角度去阐述和评价，而在第 7 节再统一以 LaTeX 的角度去阐述和评价。

另外，outreg2 只支持描述性统计输出和回归结果输出，而 logout 支持分组 T 均值检验和相关系数矩阵的输出，且两者命令结构相似，且都来自同一个作者（Roy Wada），于是本文将这两个命令合并到一起看作一个命令进行讲解（简称 out_log 命令）。于是本文要评价的输出命令为：asdoc、xxx2docx、out_log 和 esttab 四个命令。

以上命令都可以通过 ssc install 命令名称, replace 安装更新，本文在测试时均采用了行文时的最新版本。

2. 数据的前期处理

sysuse "nlsw88.dta", clear
tabulate race, gen(race_num)
drop race_num1

3. 描述性统计输出

3.1 命令展示与描述

3.1.1 asdoc

local varlist "wage age race married grade collgrad south union occupation"
asdoc sum `varlist', save(Myfile.doc) replace ///
	stat(N mean sd min p50 max)  dec(3) ///
	title(asdoc_Table: Descriptive statistics)

命令特点及存在的问题：

1. 在将结果导入到 Word 的同时，Stata 界面也能看到相应的结果(下同)。
2. 该命令不支持中文。
3. 若有字符串变量，命令会报错。

3.1.2 sum2docx

local varlist "wage age race married grade collgrad south union occupation"
sum2docx `varlist' using Myfile.docx,replace ///
	stats(N mean(%9.2f) sd(%9.3f) min(%9.2f) median(%9.2f) max(%9.2f)) ///
	title(sum2docx_Table: Descriptive statistics)

命令特点及存在的问题：

1. 在将结果导入到 Word 的同时，Stata 界面不能看到相应的结果(下同)。
2. 该命令支持中文。
3. 若有字符串变量，命令会报错。
4. 能分别设置每个统计量的小数点位数。

3.1.3 outreg2

local varlist "wage age race married grade collgrad south union occupation"
outreg2 using Myfile, sum(detail) replace word eqkeep(N mean sd min p50 max) ///
	fmt(f) keep(`varlist') sortvar(wage age grade) ///
	title(outreg2_Table: Descriptive statistics)

命令特点及存在的问题：

1. 在将结果导入到 Word 的同时，Stata 界面能看到相应的结果(下同)。
2. 该命令不支持中文。
3. 若有字符串变量，命令会在窗口说明什么变量是字符型，并在报告列表中自动剔除该变量。
4. 支持变量排序。

3.1.4 esttab

local varlist "wage age race married grade collgrad south union occupation"
estpost summarize `varlist', detail
esttab using Myfile.rtf, ///
	cells("count mean(fmt(2)) sd(fmt(2)) min(fmt(2)) p50(fmt(2)) max(fmt(2))") ///
	noobs compress replace title(esttab_Table: Descriptive statistics)

命令特点及存在的问题：

1. 在将结果导入到 Word 的同时，Stata 界面能看到相应的结果(下同)。
2. 该命令不支持中文(有些教程说可以对中文进行 gb18030 转码后导入，但本文试验不成功)。
3. 若有字符串变量，命令仍会直接运行，在输出列表中字符型变量名会写上去，但后面的统计量为空白。
4. 能分别设置每个统计量的小数点位数。

3.2 小结

• 在将结果导入到 Word 的同时，命令asdoc、outreg2、esttab能在 Stata 界面能看到相应的结果，而sum2docx不行(下同)。
• 仅sum2docx支持中文，其余命令不支持。
• 若变量里有字符串变量，outreg2命令的处理最智能化。
• sum2docx和esttab能分别设置每个统计量的小数点位数，其他的命令不行。
• 仅outreg2支持变量排序。

4. 分组 T 均值检验输出

4.1 命令展示与描述

4.1.1 asdoc

local common_exp "save(Myfile.doc) by(south) stat(obs mean p)"
asdoc ttest wage, `common_exp' replace title(asdoc_Table: T_test by group)
asdoc ttest age, `common_exp' rowappend
asdoc ttest race, `common_exp' rowappend
asdoc ttest married, `common_exp' rowappend
asdoc ttest grade, `common_exp' rowappend
asdoc ttest collgrad, `common_exp' rowappend
asdoc ttest union, `common_exp' rowappend

命令特点及存在的问题：

1. 每次只能检验一个变量, 当然它可以通过rowappend选项进行叠加，但这样麻烦不少。
2. 内容不完整，本身没有也无法添加 MeanDiff 一列。

4.1.2 t2docx

local varlist "wage age race married grade collgrad union"
t2docx `varlist' using Myfile.docx,replace ///
	not by(south) title(t2docx_Table: T_test by group)

命令特点及存在的问题：

1. 表格报告的内容剔除了所输入变量中包含缺漏值的观测值，且该操作无法被禁止。
2. 支持中文。
3. 有 MeanDiff 一列，并默认以* p < 0.1 ** p < 0.05 *** p < 0.01的方式标注星号。

4.1.3 logout

local varlist "wage age race married grade collgrad union"
logout, save(Myfile) word replace: ttable2 `varlist', by(south)

命令特点及存在的问题：

1. 能一次性分组 T 均值检验所有变量
2. 有 MeanDiff 一列，并默认以* p < 0.1 ** p < 0.05 *** p < 0.01的方式标注星号
3. 该命令不能书写表格标题。

4.1.4 esttab

local varlist "wage age race married grade collgrad union"
estpost ttest `varlist', by(south)
esttab using Myfile.rtf, ///
	cells("N_1 mu_1(fmt(3)) N_2 mu_2(fmt(3)) b(star fmt(3))") starlevels(* 0.10 ** 0.05 *** 0.01) ///
	noobs compress replace title(esttab_Table: T_test by group)

命令特点及存在的问题：

1. 能一次性分组 T 均值检验所有变量。
2. 有 MeanDiff 一列，可设置以* p < 0.1 ** p < 0.05 *** p < 0.01 的方式标注星号
3. 能分别设置每一列的小数点位数。

4.2 小结

• t2docx 表格报告的内容剔除了所输入变量中包含缺漏值的观测值，且该操作无法被禁止，研究者在查看输出结果时需要注意。其他的命令无此问题。
• 在将结果导入到 Word 的同时，asdoc、logout、esttab命令能在 Stata 界面能看到相应的结果，而t2docx不行。
• logout、esttab 和 t2docx 能一次性针对所有变量进行分组 T 均值检验，而asdoc需要多行命令，不太方便。
• asdoc 输出的内容不完整，没有 MeanDiff 一列，更无法设置相应的星号标注，而其他命令无此问题。
• logout 无法书写表格标题，而其他命令无此问题。
• esttab 能设置每一列的小数点位数，而其他命令不行。
• t2docx 支持中文，而其他的命令不行。

5. 相关系数矩阵输出

5.1 命令展示与描述

5.1.1 asdoc

local varlist "wage age race married grade collgrad"
asdoc cor `varlist', save(Myfile.doc) replace nonum dec(3) ///
	title(asdoc_Table: correlation coefficient matrix)

命令特点及存在的问题：

1. 无法使用 asdoc pwcorr 命令，而简易的 asdoc cor 命令不能报告 p 值(显著性)，不能标星号。

5.1.2 corr2docx

local varlist "wage age race married grade collgrad"
corr2docx `varlist' using Myfile.docx, replace spearman(ignore) pearson(pw) ///
	star title(corr2docx_Table: correlation coefficient matrix)

命令特点及存在的问题：

1. 支持中文
2. 不能报告 p 值，但能以* p < 0.1 ** p < 0.05 *** p < 0.01的方式标注星号

5.1.3 logout

local varlist "wage age race married grade collgrad"
logout, save(Myfile) word replace : pwcorr_a `varlist', ///
	star1(0.01) star5(0.05) star10(0.1)

命令特点及存在的问题：

1. 可设置以* p < 0.1 ** p < 0.05 *** p < 0.01的方式标注星号。
2. 表格会有串行的问题，需要后期自己手动调整。
3. 不能设置表格标题。

5.1.4 esttab

local varlist "wage age race married grade collgrad"
estpost correlate `varlist', matrix
esttab using Myfile.rtf, ///
	unstack not noobs compress nogaps replace star(* 0.1 ** 0.05 *** 0.01) ///
	b(%8.3f) p(%8.3f) title(esttab_Table: correlation coefficient matrix)

命令特点及存在的问题：

1. 能报告 p 值，可设置以* p < 0.1 ** p < 0.05 *** p < 0.01的方式标注星号。

5.2 小结

• corr2docx支持中文，而其他命令不行。
• esttab既能报告 p 值，也能以一般文献的要求报告星号，logout也可以，但有串行的问题；corr2docx只能标注星号；asdoc都不能做。
• logout不能设置表格标题，而其他的命令可以。

连享会现场课程报名中……

6. 回归结果输出

6.1 命令展示与描述

6.1.1 asdoc

asdoc reg wage age married occupation, save(Myfile.doc) nest replace ///
	cnames(OLS-1) rep(se) add(race, no)
asdoc reg wage age married collgrad occupation, save(Myfile.doc) nest append ///
	cnames(OLS-2) add(race, no)
asdoc reg wage age married collgrad occupation race_num*, save(Myfile.doc) nest append ///
	add(race, yes) cnames(OLS-3) dec(3) drop(occupation race_num*) ///
	stat(r2_a, F, rmse, rss) title(asdoc_Table: regression result)

命令特点及存在的问题：

1. 可添加不同列的列名。
2. 没有诸如indicate("race=race_num*")的选项，需用add(race, yes)代替。
3. 不能对变量进行排序。
4. 表格布局与一般的文献不同，且不能用命令调整。

6.1.2 reg2docx

reg wage age married occupation
est store m1
reg wage age married collgrad occupation
est store m2
reg wage age married collgrad occupation race_num*
est store m3
reg2docx m1 m2 m3 using Myfile.docx, replace indicate("race=race_num*") ///
	b(%9.2f) se(%7.2f) scalars(r2(%9.3f) r2_a(%9.2f) N) ///
	drop(occupation) order(married) ///
	title(reg2docx_Table: regression result) ///
	mtitles("OLS-1" "OLS-2" "OLS-3")

命令特点及存在的问题：

1. 可添加不同列的列名
2. 有诸如indicate("race=race_num*")的选项。
3. 命令运行之后会将所有已经储存的回归结果清除，导致后面无法再继续使用前面的回归结果。
4. 支持中文。

6.1.3 outreg2

local subexp "nor2 noobs e(r2 r2_a F N) fmt(f) bdec(2) sdec(2) drop(occupation race_num*) sortvar(married)"
reg wage age married occupation
outreg2 using Myfile, word replace title(outreg2_Table: regression result) ///
	ctitle(OLS-1) `subexp' addtext(race, no)

reg wage age married collgrad occupation
outreg2 using Myfile, word append ctitle(OLS-2) `subexp' addtext(race, no)

reg wage age married collgrad occupation race_num*
outreg2 using Myfile, word append ctitle(OLS-3) `subexp' addtext(race, yes)

命令特点及存在的问题：

1. 可添加不同列的列名。
2. 没有诸如indicate("race=race_num*")的选项，需用addtext(race, yes)代替。

6.1.4 esttab

reg wage age married occupation
est store m1
reg wage age married collgrad occupation
est store m2
reg wage age married collgrad occupation race_num*
est store m3
esttab m1 m2 m3 using Myfile.rtf, ///
	replace star( * 0.10 ** 0.05 *** 0.01 ) nogaps compress ///
	order(married) drop(occupation) b(%20.3f) se(%7.2f)  ///
	r2(%9.3f) ar2 aic bic obslast scalars(F)  ///
	indicate("race=race_num*") mtitles("OLS-1" "OLS-2" "OLS-3") ///
	title(esttab_Table: regression result)

命令特点及存在的问题：

1. 可添加不同列的列名。
2. 有诸如indicate("race=race_num*")的选项。
3. 表格形式最符合一般文献报告表格的格式。

6.2 小结

• 都可添加不同列的列名。
• reg2docx和esttab都有诸如indicate("race=race_num*")的选项，而其他命令都需要手动添加。
• reg2docx运行之后会将所有已经储存的回归结果清除，导致后面无法再继续使用前面的回归结果。而esttab无此问题。
• reg2docx支持中文，而其他命令不能。
• 相对于其他命令，esttab输出到 Word 的表格形式最符合一般文献报告表格的格式。

7. LaTeX 代码输出

以上四个命令中，只有out_log和esttab支持 LaTeX 代码输出，所以本部分从这两个命令入手，阐述它们的特点和优劣。

由于 LaTeX 代码编译之后，表格与 Word 的一致，所以本部分不再展示导出后的表格。

7.1 out_log

7.1.1 代码部分

*描述性统计输出
local varlist "wage age race married grade collgrad south union occupation"
outreg2 using Myfile, sum(detail) replace tex eqkeep(N mean sd min p50 max) ///
	fmt(f) keep(`varlist') sortvar(wage age grade) ///
	title(outreg2_Table: Descriptive statistics)
/*
默认：输出一个完整的可以编译的但格式较为简易 .tex 文件，表格默认居左对齐。
outreg2 的 tex() 的专有选项：
1. fragment：只输出表格本身。
2. pretty：在默认的基础上增加了更为丰富的格式，并将表格居中对齐。
3. landscape：横置页面。
*/


*分组 T 均值检验输出
local varlist "wage age race married grade collgrad union"
logout, save(Myfile) tex replace: ttable2 `varlist', by(south)
/*
1. LaTeX 输出不能增加表格标题。
2. logout 的 tex 没有选项，只有默认形式。
*/


*相关系数矩阵输出
local varlist "wage age race married grade collgrad"
logout, save(Myfile) tex replace : pwcorr `varlist', ///
	star(0.05) listwise
/*
1. 表格会有串行的问题，需要后期自己手动调整。
2. LaTeX 输出不能增加表格标题。
3. logout 的 tex 没有选项，只有默认形式。
*/


*回归结果输出
local subexp "tex nor2 noobs e(r2 r2_a F N) fmt(f) bdec(2) sdec(2) drop(occupation race_num*) sortvar(married)"
reg wage age married occupation
outreg2 using Myfile, replace title(outreg2_Table: regression result) ///
	ctitle(OLS-1) `subexp' addtext(race, no)
reg wage age married collgrad occupation
outreg2 using Myfile, append ctitle(OLS-2) `subexp' addtext(race, no)
reg wage age married collgrad occupation race_num*
outreg2 using Myfile, append ctitle(OLS-3) `subexp' addtext(race, yes)

7.1.2 小结

• out_log针对每个输出部分都能导出结果，但每个部分前后不能append在一起。本文测试过，如果用命令强行append，表格内容会乱套。所以，如果要使用该命令导出 LaTeX 代码，建议自己先写一个总的 LaTeX 框架，然后针对每个导出的部分采用\input{}命令即可。
• 表格本身也存在着一些小问题，如采用logout命令的部分无法书写表格标题，表格可能会出现串行。

7.2 esttab

7.2.1 代码部分

*描述性统计输出
local varlist "wage age race married grade collgrad south union occupation"
estpost summarize `varlist', detail
esttab using Myfile.tex, ///
	cells("count mean(fmt(2)) sd(fmt(2)) min(fmt(2)) p50(fmt(2)) max(fmt(2))") ///
	noobs compress replace title(esttab_Table: Descriptive statistics) ///
	booktabs  page(array, makecell) alignment(cccccc) width(\hsize)
/*
esttab 的 LaTeX 输出的专有选项：
1. booktabs: 用 booktabs 宏包输出表格(三线表格)。
2. page[(packages)]: 创建完成的 LaTeX 文档以及添加括号里的宏包
3. 如果写了 booktabs 选项，则 page[(packages)] 将自动添加\usepackage{booktabs}。
4. alignment(cccccc)：定义从第二列开始的列对齐方式(默认居中)。
5. width(\hsize)：可以使得表格宽度为延伸至页面宽度
6. fragment：不输出表头表尾，只输出表格本身内容，其不能与 page[(packages)] 选项共存。
*/


*分组 T 均值检验输出
local varlist "wage age race married grade collgrad union"
estpost ttest `varlist', by(south)
esttab using Myfile.tex, ///
	cells("N_1 mu_1(fmt(3)) N_2 mu_2(fmt(3)) b(star fmt(3))") starlevels(* 0.10 ** 0.05 *** 0.01) ///
	noobs compress append title(esttab_Table: T_test by group) ///
	booktabs  page width(\hsize)
/*
1. 这里的 page 选项可以将本部分附着在先前的 document 环境中，从而形成一个整体。
*/


*相关系数矩阵输出
local varlist "wage age race married grade collgrad"
estpost correlate `varlist', matrix listwise
esttab using Myfile.tex, ///
	unstack not noobs compress nogaps append star(* 0.1 ** 0.05 *** 0.01) ///
	b(%8.3f) p(%8.3f) title(esttab_Table: correlation coefficient matrix) ///
	booktabs  page width(\hsize)


*回归结果输出
reg wage age married occupation
est store m1
reg wage age married collgrad occupation
est store m2
reg wage age married collgrad occupation race_num*
est store m3
esttab m1 m2 m3 using Myfile.tex, ///
	append star( * 0.10 ** 0.05 *** 0.01 ) nogaps compress ///
	order(married) drop(occupation) b(%20.3f) se(%7.2f)  ///
	r2(%9.3f) ar2 aic bic obslast scalars(F)  ///
	indicate("race=race_num*") mtitles("OLS-1" "OLS-2" "OLS-3") ///
	title(esttab_Table: regression result) booktabs  page width(\hsize)

7.2.2 小结

• 这里所有的表格都可以生成标题，也没有串行的问题。
• 与out_log不一样，esttab能形成一个整体，所有的表格都能完美的append在一起，并且可以设置生成导言区和正文区的document环境。所以以上命令生成的 .tex 文件可以直接拿去编译，并生成所有表格，一步到位。之后我们再统一做一些微调就可以放在 LaTeX 正文中了。

8. 总结

上面从各个角度对四个实证结果输出命令，这里做一个总结。

• 从描述性统计输出的角度，四个命令的效果差不多，只是sum2docx和esttab能分别设置每个统计量的小数点位数而占据了上风。
• 从分组 T 均值检验输出 的角度，asdoc输出内容不完整且操作最麻烦；logout和t2docx输出表格虽然完整但各有优缺点（logout不能设置表格标题；t2docx表格报告的内容剔除了所输入变量中包含缺漏值的观测值，且该操作无法被禁止）；esttab功能最完整。
• 从相关系数矩阵输出的角度，esttab既能报告 p 值，也能以一般文献的要求报告星号，logout也可以，但有串行的问题；corr2docx只能标注星号；asdoc都不能做。
• 从回归结果输出的角度，asdoc的效果最差，不能对变量进行排序，表格布局与一般的文献不同，且不能用命令调整；其他三个命令功能都是完整的，而对于outreg2没有 indicate()选项使得操作更麻烦，reg2docx运行之后会将所有已经储存的回归结果清除，导致后面无法再继续使用前面的回归结果（这应该是该命令近期更新后的一个 bug）。
• 从 Word 整体性的角度（即将所有表格都append到一个 Word 里的程度），asdoc、xxx2docx和esttab都能完美的自动实现一个 Word 文件装上所有表格，而out_log无法实现这一点，主要是由于其中的logout命令没有append选项。
• 从 LaTeX 代码输出的角度，只有out_log和esttab支持 LaTeX 输出。其中esttab能自动的形成包含所有表格的单个 .tex 文件，且该 .tex 能直接进行编译。这使得esttab在这方面极大的优于out_log。
• 从综合评价的角度，功能越完整，自动化程度越高，那么这样的命令往往也越会获得我们的青睐。通过前面的论述，本文认为，asdoc的 Word 整体性虽然很好，但由于其他的缺陷实在太多，于是将其归为最差的位置；out_log既支持 Word 输出也支持 LaTeX 输出，在复杂的回归结果输出中有着强悍的力量，但是，在一般的应用中，其整体效果还是要差于xxx2docx；xxx2docx是唯一能够完美支持中文和将表格导出到 .docx 文件的命令，命令表达式也通俗易懂，但是其不支持 LaTeX 代码的输出和近期更新后出现的一些 bug，使得该命令的综合表现要差于esttab；esttab在某些方面虽然语句复杂一些，但也造就了其功能的完整性，使得其在每一方面都能有很好的表现，因此本文对于该命令的综合评价位于四大输出命令之首。

为了能够更形象的表示出本文对于以上四个实证结果输出命令在各个方面的评价，本文将以推荐指数的形式，列于下方的表格中。

9. 相关附件

本文使用的所有 Dofiles：「全体 Stata 实证结果输出命令.do」

相关课程

连享会-直播课 上线了！
http://lianxh.duanshu.com

免费公开课：

• 直击面板数据模型 - 连玉君，时长：1小时40分钟
• Stata 33 讲 - 连玉君, 每讲 15 分钟.
• 部分直播课 课程资料下载 (PPT，dofiles等)

课程一览

支持回看，所有课程可以随时购买观看。

Note: 部分课程的资料，PPT 等可以前往 连享会-直播课 主页查看，下载。

关于我们

• Stata连享会 由中山大学连玉君老师团队创办，定期分享实证分析经验。直播间 有很多视频课程，可以随时观看。
• 连享会-主页 和 知乎专栏，300+ 推文，实证分析不再抓狂。
• 公众号推文分类： 计量专题 | 分类推文 | 资源工具。推文分成 内生性 | 空间计量 | 时序面板 | 结果输出 | 交乘调节 五类，主流方法介绍一目了然：DID, RDD, IV, GMM, FE, Probit 等。
• 公众号关键词搜索/回复 功能已经上线。大家可以在公众号左下角点击键盘图标，输入简要关键词，以便快速呈现历史推文，获取工具软件和数据下载。常见关键词：课程, 直播, 视频, 客服, 模型设定, 研究设计, stata, plus, 绘图, 编程, 面板, 论文重现, 可视化, RDD, DID, PSM, 合成控制法 等

连享会小程序：扫一扫，看推文，看视频……

扫码加入连享会微信群，提问交流更方便

✏ 连享会学习群-常见问题解答汇总：
✨ https://gitee.com/arlionn/WD