stata结果导出

stata结果导出

资料来源微信公众号：爬⾍俱乐部、 stata连享会

结果输出##

clear all //清空内存

set more off // 关闭more选项。

webuse nlswork,clear //导⼊⽹上数据

xtset idcode year //设置⾯板标识

tab year,gen(yd) //⽣成年份虚拟变量

*回归1：未控制年份的OLS，对标准误进⾏聚类修正

*回归2：控制年份的OLS，对标准误进⾏聚类修正

*回归3：控制年份的⾯板固定效应模型，对标准误进⾏聚类修正

注1：star选项可以对显著性进⾏调整，此处

表⽰10%，表⽰5%，*表⽰1%。

表⽰5%，*表⽰1%。

*注2：b()，t()选项分别为调整输出系数和t值的格式，%6.3f表⽰，长度为6，⼩数3位。

*注3：compress和nogap选项可以压缩⾏距和列距。

*注4：由于year是控制变量，我们也不想观察其系数，故在输出时使⽤drop()选项删除。

*注5：stats()选项可以在输出结果中增加⾏注释或者回归存储值；这⾥我们输出了观测值个数N和调整后的R2，并为其设置显⽰格式。

     reg ln_w age ttl_exp tenure not_smsa south, vce(cluster idcode)

     est store m1

     reg ln_w age ttl_exp tenure not_smsa south yd*, vce(cluster idcode)

     est store m2

      xtreg ln_w age ttl_exp tenure not_smsa south yd*, fe vce(cluster idcode)

      est store m3

   local s "using $Out\Table11.csv"     // 

指定存储结果的

 Excel  

⽂档名称

   local m " m1 m2 m3 "               //

指定

"m1 m2 m3 " 

为

m

   esttab  `m' `s' , star(* 0.1 ** 0.05 *** 0.01) b(%6.3f) t(%6.3f) 

   compress nogap ///

   drop(yd*) stats(N r2_a, fmt(%12.0f  %9.3f)) varwidth(20) ///

    title("Table1 Wage") mtitle("OLS" "OLS" "FE")