萌糍粑 / STATA / Stata結果輸出:outreg2 命令詳解

分享

   

Stata結果輸出:outreg2 命令詳解

2020-01-13  萌糍粑

作者:云鋒 (華南師范大學) 

2020.3.26-29

目錄

  • 1.引言

  • 2. outreg2 命令代碼及實例

    • 2.1 描述性分析

    • 2.2 回歸結果輸出

  • 3.小結

1.引言

一篇實證論文中,最基本也是最重要的部分就是展示 Stata 中得出的統計分析、回歸結果等表格。但自己動手做表格往往非常繁瑣,Word 排版也常常令人抓狂。而 outreg2 命令可以讓 Stata 自動輸出我們想要的表格,為你解決所有結果輸出的煩惱。因此,熟練掌握 outreg2 命令對我們快速導出 Stata 結果,一步到位的完成實證結果展示有莫大幫助。

2. outreg2 命令代碼及實例

2.1 描述性分析

首先我們介紹描述性統計分析表格的輸出,一般實證論文都選擇展示全部變量的一般統計指標,這里我們將各種情況都羅列出來,并利用 Stata 自帶的 1978 年美國汽車交易數據 auto.dta 為大家進行實例演示,方便大家全面了解 outreg2 命令。

2.1.1 全部變量一般統計指標

  • 代碼

outreg2 using xxx.doc,replace sum(log) title(Decriptive statistics)

  • 實例
sysuse auto, clear
outreg2 using xxx.doc, replace sum(log) title(Decriptive statistics)
全部變量一般統計指標輸出

命令解析

  1. xxx.doc 為輸出文件名為 xxx 的 word 文檔命令。類似得,如果需要輸出為 excel 表格則更改代碼為 xxx.xls 即可。
  2. sum(log) 即輸出一般統計指標命令,一般統計指標包括樣本數、中值、標準誤、最大值和最小值。
  3. title(Decriptive statistics) 為自定義輸出表格的名稱命令,可在括號內自行編輯表格名稱。

2.1.2 部分變量一般統計指標輸出

  • 代碼

outreg2 using xxx.doc, replace sum(log) keep(var1 var2 var3) title(Decriptive statistics)

  • 實例
sysuse auto, clear
outreg2 using xxx.doc, replace sum(log) keep(**price** **weight** **foreign**) title(Decriptive statistics)
部分變量一般統計指標輸出結

命令解析

  1. keep(var1 var2 var3) 為保留部分變量輸出命令,其中 var1 var2 var3 是我們所需要保留的變量名稱。
  2. xxx.doc 為輸出文件名為 xxx 的 word 文檔命令。類似得,如果需要輸出為 excel 表格則更改代碼為 xxx.xsl 即可。
  3. title(Decriptive statistics) 為自定義輸出表格的名稱命令,可在括號內自行編輯表格名稱。

2.1.3 全部統計指標輸出

  • 代碼

outreg2 using xxx.doc, replace sum(detail) title(Decriptive statistics)

  • 實例
sysuse auto, clear
outreg2 using xxx.doc, replace sum(detail) title(Decriptive statistics)
全部統計指標輸出

命令解析

  1. sum(detail) 為輸出全部統計指標命令,除去一般統計指標外還包括方差、偏度、峰度、極差等指標。
  2. xxx.doc 為輸出文件名為 xxx 的 word 文檔命令。類似得,如果需要輸出為 excel 表格則更改代碼為 xxx.xsl 即可。
  3. title(Decriptive statistics) 為自定義輸出表格的名稱命令,可在括號內自行編輯表格名稱。

2.1.4 部分變量部分統計指標輸出

  • 代碼

outreg2 using xxx.doc, replace sum(log) keep(var1 var2 var3) eqkeep(N min mix) title(Decriptive statistics)

  • 實例
sysuse auto, clear
outreg2 using xxx.doc, replace sum(log) keep(price weight foreign) eqkeep(N min mix) title(Decriptive statistics)
部分變量部分統計指標輸出結果

命令解析

  1. eqkeep(N min mix) 為保留部分統計指標輸出命令,這里我們以樣本數、最大值和最小值為例。
  2. keep(var1 var2 var3) 為保留部分變量輸出命令,同時利用 keep()eqkeep() 命令進行限定即可進行部分變量和部分統計指標輸出。
  3. xxx.doc 為輸出文件名為 xxx 的 word 文檔命令。類似得,如果需要輸出為 excel 表格則更改代碼為 xxx.xsl 即可。
  4. title(Decriptive statistics) 為自定義輸出表格的名稱命令,可在括號內自行編輯表格名稱。

2.1.5 分組描述

  • 代碼

bysort var1: outreg2 using xxx.doc, replace sum(log) title(Decriptive statistics)

  • 實例
sysuse auto, clear
bysort foreign:outreg2 using xxx.doc, replace sum(log) title(Decriptive statistics)
分組描述輸出結果

命令解析

  1. bysort var1 為按照變量 var1 進行分組命令,在實例中我們以是否為國外品牌進行分組描述。
  2. 若論文中需要描述部分變量或部分統計指標時,同樣可以添加前文所述 keep()eqkeep() 命令,此處不再贅述。
  3. xxx.doc 為輸出文件名為 xxx 的 word 文檔命令。類似得,如果需要輸出為 excel 表格則更改代碼為 xxx.xsl 即可。
  4. title(Decriptive statistics) 為自定義輸出表格的名稱命令,可在括號內自行編輯表格名稱。.

2.1.6 頻數

頻數是指落入不同組的個體數。

  • 代碼

outreg2 var1 using xxx.doc, replace cross title(Decriptive statistics)

  • 實例
sysuse auto, clear
outreg2 foreign using xxx.doc, replace cross title(Decriptive statistics)
頻數輸出結果

命令解析

  1. 代碼中 var1 即所需描述頻數的分組變量。
  2. xxx.doc 為輸出文件名為 xxx 的 word 文檔命令。類似得,如果需要輸出為 excel 表格則更改代碼為 xxx.xsl 即可。
  3. title(Decriptive statistics) 為自定義輸出表格的名稱命令,可在括號內自行編輯表格名稱。

2.2 回歸結果輸出

2.2.1 線性回歸結果輸出

  • 代碼

outreg2 using xxx.doc,replace tstat bdec(3) tdec(2) ctitle(y)

  • 實例
sysuse auto, clear
reg price mpg
outreg2 using xxx.doc,replace tstat bdec(3) tdec(2) ctitle(y)
線性回歸結果輸出

命令解析

  1. ctitle 為自定義表格內標題命令,如果不進行設定則直接輸出為被解釋變量名。
  2. xxx.doc 為輸出文件名為 xxx 的 word 文檔命令。類似得,如果需要輸出為 excel 表格則更改代碼為 xxx.xsl 即可。
  3. 按照 outreg2 命令輸出的表格內相關系數下括號內數字為標準誤,因此我們一般利用命令 tstat 將其更改為 t 值。
  4. outreg2 命令輸出時默認相關系數和 t 值都保留 3 位有效數字,而一般期刊要求相關系數保留 3 位有效數字,t 值保留 2 位有效數字,因此我們利用 bdec(3)tdec(2) 命令限定。

2.2.2 面板固定效應回歸結果輸出

  • 代碼

outreg2 using xxx.doc,replace tstat ctitle(y) bdec(3) tdec(2) keep(var1 var2 var3) addtext(Company FE, YES )

  • 實例
webuse grunfeld,clear
xtset company year
xtreg invest mvalue kstock,fe robust
outreg2 using xxx.doc,replace tstat bdec(3) tdec(2) ctitle(y) keep(invest mvalue kstock) addtext(Company FE, YES )
面板固定效應回歸結果輸出

同理,若要同時控制時間固定效應代碼如下:

  • 代碼

outreg2 using xxx.doc,replace tstat bdec(3) tdec(2) ctitle(y) keep(var1 var2 var3) addtext(Company FE, YES,Year FE, YES)

  • 實例
webuse grunfeld,clear
xtset company year
xtreg invest mvalue kstock i.year,fe robust
outreg2 using xxx.doc,replace tstat bdec(3) tdec(2) ctitle(y) keep(invest mvalue kstock) addtext(Company FE, YES,Year FE, YES)
同時控制時間固定效應結果輸出

命令解析

  1. addtext 為在表中增加信息命令,由于 Stata 進行固定效應回歸后單純利用 outreg2 命令輸出不會展示是否控制固定效應,因此我們需要利用 addtext 命令追加。
  2. keep() 命令的作用是保持輸出表格簡潔,括號中是所有需要展示的變量。
  3. ctitle 為自定義表格內標題命令,如果不進行設定則直接輸出為被解釋變量名。
  4. xxx.doc 為輸出文件名為 xxx 的 word 文檔命令。類似得,如果需要輸出為 excel 表格則更改代碼為 xxx.xsl 即可。
  5. 按照 outreg2 命令輸出的表格內相關系數下括號內數字為標準誤,因此我們一般利用命令 tstat 將其更改為 t 值。
  6. outreg2 命令輸出時默認相關系數和 t 值都保留 3 位有效數字,而一般期刊要求相關系數保留 3 位有效數字,t 值保留 2 位有效數字,因此我們利用 bdec(3)tdec(2) 命令限定。

2.2.3 逐步增加變量回歸合并輸出

  • 代碼
outreg using xxx.doc,replace tstat bdec(3) tdec(2) ctitle(y)
outreg using xxx.doc,append tstat bdec(3) tdec(2) ctitle(y)
  • 實例
reg price mpg
outreg2 using xxx.doc,replace tstat bdec(3) tdec(2) ctitle(y)
reg price mpg rep78
outreg2 using xxx.doc,append tstat bdec(3) tdec(2) ctitle(y)
逐步增加變量回歸結果合并輸出

命令解析

  1. append 為合并命令,在實例中我們先進行了 pricempg 兩個變量回歸,之后增加變量 rep78 回歸后再利用 outreg2 + append 命令進行合并輸出。
  2. ctitle 為自定義表格內標題命令,如果不進行設定則直接輸出為被解釋變量名。
  3. xxx.doc 為輸出文件名為 xxx 的 word 文檔命令。類似得,如果需要輸出為 excel 表格則更改代碼為 xxx.xsl 即可。
  4. 按照 outreg2 命令輸出的表格內相關系數下括號內數字為標準誤,因此我們一般利用命令 tstat 將其更改為 t 值。
  5. outreg2 命令輸出時默認相關系數和 t 值都保留 3 位有效數字,而一般期刊要求相關系數保留 3 位有效數字,t 值保留 2 位有效數字,因此我們利用 bdec(3)tdec(2) 命令限定。

2.2.4 不同模型回歸結果合并輸出

  • 代碼
outreg2 using xxx.doc,replace tstat bdec(3) tdec(2) ctitle(model1)
outreg2 using xxx.doc,append tstat bdec(3) tdec(2) ctitle(model2)
outreg2 using xxx.doc,append tstat bdec(3) tdec(2) ctitle(model3)
  • 實例
webuse grunfeld,clear
xtset company year
reg invest mvalue kstock,robust
outreg2 using xxx.doc,replace tstat bdec(3) tdec(2) ctitle(OLS)
xtreg invest mvalue kstock,fe robust
outreg2 using xxx.doc,append tstat bdec(3) tdec(2) ctitle(FE) addtext(Company FE, YES)
xtreg invest mvalue kstock i.year,fe robust
outreg2 using xxx.doc,append tstat bdec(3) tdec(2) ctitle(FE) keep(invest mvalue kstock) addtext(Company FE, YES,Year FE, YES)
不同模型回歸結果合并輸出

連享會計量方法專題……[1],https://gitee.com/arlionn/Course

命令解析

  1. 與逐步增加變量合并輸出相類似,不同模型回歸結果輸出也是利用 append 命令。
  2. 利用 ctitle 可將標題命名為模型名稱。
  3. 固定效應回歸后需利用 addtext 增加是否控制固定效應信息。
  4. xxx.doc 為輸出文件名為 xxx 的 word 文檔命令。類似得,如果需要輸出為 excel 表格則更改代碼為 xxx.xsl 即可。
  5. 按照 outreg2 命令輸出的表格內相關系數下括號內數字為標準誤,因此我們一般利用命令 tstat 將其更改為 t 值。
  6. outreg2 命令輸出時默認相關系數和 t 值都保留 3 位有效數字,而一般期刊要求相關系數保留 3 位有效數字,t 值保留 2 位有效數字,因此我們利用 bdec(3)tdec(2) 命令限定。

2.2.5 包含兩階段回歸的結果輸出——以 2SLS 為例

  • 第一階段回歸結果輸出代碼
ivregress2 2sls y var3 (var1=var2), first
est restore first
outreg2 using xxx.doc, cttop(first) tstat bdec(3) tdec(2) replace
  • 第二階段回歸結果輸出代碼
ivregress2 2sls y var3 (var1=var2), first
outreg2 using xxx.doc, cttop(two) tstat bdec(3) tdec(2) replace
  • 實例
sysuse auto
ivregress2 2sls mpg weight (length=displacement),first
est restore first
outreg2 using xxx.doc,cttop(first) tstat bdec(3) tdec(2) replace

ivregress2 2sls mpg weight (length=displacement), first
outreg2 using xxx.doc,cttop(two) tstat bdec(3) tdec(2)
2SLS 兩階段回歸結果輸出
  • 命令解析
    1. 兩階段回歸輸出時,若只展示第二階段回歸,則直接利用第二階段回歸輸出代碼即可。如需展示第一階段回歸結果,則可以利用實例中所示方法,利用命令先將第一階段回歸結果保存輸出,再輸出第二階段回歸結果。實例中去掉了第二階段輸出時的 replace 這樣就可以將兩階段結果輸出在一張表上。
    2.  按照 outreg2 命令輸出的表格內相關系數下括號內數字為標準誤,因此我們一般利用命令 tstat 將其更改為 t 值。
    3.  outreg2 命令輸出時默認相關系數和 t 值都保留 3 位有效數字,而一般期刊要求相關系數保留 3 位有效數字,t 值保留 2 位有效數字,因此我們利用 bdec(3) 和 tdec(2) 命令限定。

連享會計量方法專題……[2],https://gitee.com/arlionn/Course

2.2.6 自定義統計量的輸出

前文所介紹的輸出命令下 Stata 默認輸出相關系數、標準誤、樣本數和 R 平方,并保留 3 位有效數字。但根據實際情況,有些論文需要報告偽 R 平方、F 值,或者更改相關系數和 t 值有效數字位數,這里我們進一步為大家介紹自定義統計量輸出命令。

  • 更改有效數字位數實例
sysuse auto
reg price mpg rep78
outreg2 using xxx.doc,replace tstat ctitle(y) bdec(4) tdec(5)
更改有效數字位數

命令解析 

bdec(4) 是指相關系數保留 4 位有效數字,tdec(5) 指 t 值保留 5 位有效數字。

  • 在 logit 估計之后報告偽 R 平方實例
sysuse auto, clear
logit foreign price trunk weight
outreg2 using xxx.doc,replace tstat bdec(3) tdec(2) addstat(Pseudo R-squared, `e(r2_p)')
增加報告偽 R 方回歸結果輸出

命令解析

  1. logit 模型回歸是沒有 R 方的,因此在報告時我們需要輸出偽 R 方的值,在正常輸出回歸表格代碼后增加 addstat() 命令可以增加報告偽 R 方。
  2. outreg2 命令輸出時默認相關系數和 t 值都保留 3 位有效數字,而一般期刊要求相關系數保留 3 位有效數字,t 值保留 2 位有效數字,因此我們利用 bdec(3)tdec(2) 命令限定。
  • 回歸后表格中添加 F 值、調整后 R 方 和 F 檢驗 P 值輸出實例
sysuse auto,clear
reg price length rep78 weight
outreg2 using xxx.doc,replace tstat bdec(3) tdec(2) e(r2_a,F) addstat(F test,e(p))
增加報告調整后 R 方 F 值等結果輸出

命令解析

  1. 與增加偽 R 方類似,增加 F 檢驗同樣是利用 addstat 命令,括號內代碼更換為 F test,e(p) 即可。
  2. outreg2 命令輸出時默認相關系數和 t 值都保留 3 位有效數字,而一般期刊要求相關系數保留 3 位有效數字,t 值保留 2 位有效數字,因此我們利用 bdec(3)tdec(2) 命令限定。

3.小結

通過對 outreg2 命令的學習,我們不難發現它可以輸出除相關系數矩陣外,一般實證論文所需幾乎全部的表格,實用性非常強,可謂一招鮮吃遍天。最后我們整理一下利用 outreg2 輸出主要表格的代碼,方便大家理清 outreg2 命令語法。

*描述性統計表格輸出
outreg2 using xxx.doc, replace sum(detail) keep(var1 var2 var3) eqkeep(N min mix) title(Decriptive statistics) \\部分變量部分統計指標輸出

*-回歸結果輸出
outreg2 using xxx.doc,replace tstat bdec(3) tdec(2) ctitle(y) \\普通 OLS 回歸結果輸出

outreg2 using xxx.doc,replace tstat bdec(3) tdec(2) ctitle(y) keep(var1 var2 var3) addtext(Company FE, YES,Year FE, YES) \\面板固定效應回歸結果輸出(個體固定效應和時間固定效應)


*-逐步增加變量回歸結果合并輸出
outreg2 using xxx.doc,replace tstat ctitle(y)
outreg2 using xxx.doc,append tstat ctitle(y)

*-不同模型回歸結果合并輸出
outreg2 using xxx.doc,replace tstat ctitle(model1)
outreg2 using xxx.doc,append tstat ctitle(model2)
outreg2 using xxx.doc,append tstat ctitle(model3)

*-2SLS 兩階段模型回歸結果輸出
ivregress2 2sls y var3 (var1=var2), first
est restore first
outreg2 using xxx.doc, cttop(first) replace

ivregress2 2sls y var3 (var1=var2), first
outreg2 using xxx.doc, cttop(two)

連享會直播回放:我的特斯拉—實證研究設計(連玉君主講)
課程主頁:https://gitee.com/arlionn/Live

    本站是提供個人知識管理的網絡存儲空間,所有內容均由用戶發布,不代表本站觀點。如發現有害或侵權內容,請點擊一鍵舉報

    0條評論

    發表

    請遵守用戶 評論公約

    類似文章 更多
    喜歡該文的人也喜歡 更多

    ×
    ×

    .00

    微信或支付寶掃碼支付:

    開通即同意《個圖VIP服務協議》

    全部>>
    华东15选5幸运之门彩票网 安徽快三在哪里买 易发游戏每天赠6元 山西扣点点官网 一定牛陕西十一选五 深圳风采2011090 888棋牌官方下载最新版 官网武汉麻将赖子山庄 天天捕鱼赢话费破解版无限金币 广东十一选五计划app 快3青海基本走势图 黑金快乐8登录网址导航 北京麻将app 湖北赖子麻将怎么打 快船vs76人 亿客隆 吉林快三昨天开奖结果