【Stata】Stata论文实证分析的基础代码分享

最新推荐文章于 2024-06-27 14:27:10 发布

旅途中的宽~

最新推荐文章于 2024-06-27 14:27:10 发布

阅读量1.4w

点赞数 27

分类专栏：研一上学期应用回归课程论文文章标签： stata 实证分析

本文链接：https://blog.csdn.net/wzk4869/article/details/128881238

版权

研一上学期应用回归课程论文专栏收录该内容

6 篇文章

订阅专栏

文章提供了使用Stata进行论文实证分析的步骤，包括描述性统计、数据处理（如取对数、缩尾）、相关性分析、单位根检验、协整检验、主回归模型设定、稳健性测试以及异质性检验的详细代码示例。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

一、前言

论文实证分析的一整套基础流程：描述性、相关性、相关检验(单位根、协整检验)、主回归模型、稳健性分析、异质性检验。

二、基础代码

2.1 描述性分析

ssc install asdoc
asdoc sum y x1 x2 x3 x4 x5 x6 x7

2.2 数据处理

取对数：

foreach var of varlist y x1 x2 x3 x4 x5 x6 x7{
gen ln`var'=log(`var')
}

如果有缩尾：

winsor2 lny lnx1 lnx2 lnx3 lnx4 lnx5 lnx6 lnx7, replace cuts(1 99)

2.3 相关性分析

asdoc pwcorr_a lny lnx1 lnx2 lnx3 lnx4 lnx5 lnx6 lnx7, star1(0.01) star5(0.05) star10(0.1)

2.4 主回归模型之前的一些数据检验

单位根检验：

1、ADF检验

基本语法：xtunitroot fisher varname, trend dfuller demean lags(1)

2、LLC检验

基本语法：xtunitroot llc varname, lag(n) 一般使用T比较大

xtunitroot llc varname, trend demean lags(bic 12) --含个体效应和线性时间趋势项

3、IPS检验

基本语法：xtunitroot ips 变量名, lag(n) 其他语法同LLC一致

注：使用xtunitroot，help一下，便会出现所有的检验命令及其用法和实例

协整检验：

1、Kao 检验

基本语法：xtcointtest kao y x1 x2

2、Pedroni 检验

基本语法：.xtcointtest pedroni y x1 x2，trend

3、Westerlund 检验

基本语法： xtcointtest westerlund y x1 x2，trend

2.5 主回归模型

xtset province year
xtabond2 lny L.lny lnx1 lnx2 lnx3 lnx4 lnx5 lnx6 lnx7, gmmstyle(L.lny lnx1 lnx2) ivstyle(lnx3 lnx4)twostep robust small orthogonal
est store m

2.6 稳健性检验

xtreg lny lnx1 lnx2 lnx3 lnx4 lnx5 lnx6 lnx7, fe （fe是固定效应模型）
est store fe （保存回归结果）
esttab fe using reg2.rtf, replace b(%6.3f) se(%6.3f) se ar2(3) star(* 0.1 ** 0.05 *** 0.01) compress nogap mtitles("model1" ) title("Table1")
Ols 
reg lny lnx1 lnx2 lnx3 lnx4 lnx5 lnx6 lnx7,robust （最小二乘法）
est store reg
esttab fe reg using reg2.rtf, replace b(%6.3f) se(%6.3f) se ar2(3) star(* 0.1 ** 0.05 *** 0.01) compress nogap mtitles("model1""model2" ) title("Table1")

2.7 异质性检验

sort x7 
gen group_x7 = group(4) 

forvalues i = 1/4{
qui reg lny lnx1 lnx2 lnx3 lnx4 lnx5 lnx6 lnx7 if group_x7 == `i', r 

qui outreg2 using "groupsave.xls", append bdec(3) tdec(2) ctitle(`y')
est store groupx7_`i'
}
esttab groupx7_*, replace nogap compress b(%6.3f) s(N r2_a) se star(* 0.1 ** 0.05 *** 0.01) addnotes("*** 1% ** 5% * 10%")