Stata—基本统计量输出、模型估计和结果输出
1. 基本统计量的输出
tabstat var, stats(mean sd min max p50 max) c(s) f(%6.2f)
输入上述命令之后,会在stata的结果窗口呈现所需的描述性统计,那么如何直接将结果转到word中?方法如下:
logout, save(mytable) word replace: tabstat var, stats(mean sd min max p50 max) c(s) f(%6.2f)
上述命令中,mytable为文件名称;replace代表每次生成的word结果都会覆盖原有的mytable中的结果;var为变量名称;word可以替换为excel,从而将结果输出至excel中。
2. 相关系数矩阵
logout, save(mytable) word replace: pwcorr var1 var2
输入上述命令可见var1和var2的相关系数矩阵会储存到名为mytable的word文件中。
logout, save(mytable) word replace: pwcorr_a var1 var2
此命令中的pwcorr_a是连玉君博士编写的代码,相比于pwcorr,此命令则可以显示相关系数的显著性(带星号)。
3. 估计结果输出
①est—估计模型
reg sas gender:以sas作为因变量,gender作为自变量行回归分析。
est store a1:将估计结果储存在a1。
reg sas gender dp:以sas作为因变量,gender,dp作为自变量行回归分析。
est store a2:将估计结果储存在a2。
esttab a1 a2:在结果显示框中显示a1和a2结果,但是比较简单,未有呈现出回归中的重要信息。
esttab a1 a2, beta:呈现标准化回归系数。
esttab a1 a2, ar2 compress nogap star(*0.1 **0.05 ***0.01):ar2命令用来显示R方;compress使得结果呈现更加紧凑;nogap命令使得两个自变量之间的空行删除;star是将默认的显著性水平进行更改,也就是说将默认的“*0.05 **0.01 ***0.001”改为“*0.1 **0.05 ***0.01”。
esttab a1 a2, ar2 compress nogap star(*0.1 **0.05 ***0.01) b(%6.3f) brackets p:b(%6.3f) 代表估计的回归系数,并以3位小数呈现b估计值;brackets p,将p值放在中括号中并进行显示。
label var gender “性别”:将变量gender加标签为性别。
esttab a1 a2, label:也就是说估计结果中,呈现的是变量的标签名。
esttab a1 a2, margin constant:指的是回归中含有常数项。
esttab a1 a2 using resluts.csv, replace:该命令将统计结果输出至名为resluts的csv文件中。
②logout—储存结果
logout, save(mytable) word dec(3) replace: reg sas gender:该命令将回归分析的结果输出至名为mytable的word文件中,dec(3)表示转换格式时的敏感程度,连博士认为取值3较为合适。同样地,也可以将word替换为excel。
案例:
第一步:首先进行模型估计,并储存模型结果
reg sas gender
est store a1
reg sas gender dp
est store a2
第二步:将结果输出至word文档中
logout, save(resluts) word replace fix(3): esttab a1 a2, mtitle(model1 model2) b(%6.3f) se(%6.2f) star(*0.1 **0.05 ***0.01) scalar(r1 r2_a N P) compress nogap:mtitle(model1 model2)命令将2个模型分别命名为model1和model2;b(%6.3f) 为系数格式,se(%6.2f)为标准误的格式。
③outreg—储存结果
reg sas gender
est store a1
reg sas gender dp
est store a2
outreg2 [a1 a2] using resluts1, seeout replace:将结果呈现在Stata的数据窗口。
outreg2 [a1 a2] using resluts1, word replace:将结果呈现在名为results1的word文档中。同样地,word可以换为excel。
outreg2 [a1 a2] using resluts1, word replace label:呈现变量标签名。
outreg2 [a1 a2] using resluts1, word replace tstat tdec(2):呈现t值并保留2位小数。
outreg2 [a1 a2] using resluts1, word replace pvalue pdec(3):呈现p值并保留3位小数。