作者drewlin (小卓..)
看板NCCU08_MAT
標題[STATA] regression model、敘述統計、相關係數 輸出製表
時間Fri Apr 23 19:50:11 2010
第一次使用時請先更新estout套件
#ssc install estout,replace
開啟範例檔案
#sysuse auto,clear
輸出迴歸結果:
每次執行之前都先清除之前的estimation
#est clear
跑三個迴歸 然後存下來並命名m1 m2 m3
#reg price rep78 trunk
#est store m1
#reg price rep78 trunk displacement
#est store m2
#reg price rep78 trunk displacement gear_ratio
#est store m3
輸出結果
#esttab * using reg_result.rtf,label r2 ar2 scalars(F) nogaps p(%4.3f) replace
為了輸出結果可以直接使用,不在需要改變數名稱
可以先將變數的標籤(Label)進行設定
設定方法為舉例如下 將ln_mv的Label設為Ln(MV)
label variable ln_mv "Ln(MV)"
上面指令逗號後的label即代表輸出時用標籤命名而不是變數名稱
若要改變預設的星號level,
可以在上行後面加上starlevels(+ 0.2 ++ 0.1 * 0.05 ** 0.01 *** 0.001)
p<0.2會是一個+號 p<0.1則是兩個+號 以次類推 符號可以自訂
另外請利用word的取代功能,將標號更改顏色與格式(ex:紅色 上標)
如果不會的請洽立婷或青雲
輸出敘述統計:
一樣先清estimation
#est clear
輸出sum,price.....,gear_ratio的變數之敘述統計
#eststo: estpost sum price rep78 trunk displacement gear_ratio,detail
#esttab * using descriptive_statistic.rtf,label
cells("mean(fmt(3)) sd(fmt(3)) p25(fmt(3))
p50(fmt(3)) p75(fmt(3))") replace
最後三行為同一行請接在一起
代表輸出平均、標準差、25 50 75百分位數 並取三位小數
如果需要其他敘述統計量 可以參考第一個指令跑完的結果加進去即可
輸出相關係數矩陣:
還是繼續先清estimation
#est clear
輸出price ,rep78....,gear_ratio的相關係數
#estpost correlate price rep78 trunk displacement gear_ratio, matrix listwise
#est store c1
#esttab * using correlation.rtf,label p(%4.3f) replace b(%4.3f)
nogaps unstack not noobs compress nostar
最後兩行為同一行 p(%4.3f)為同時輸出p value 格式為X.XXX 四位數字中有三位小數
b(%4.3f)只是指定係數的格式
若不需要p value只需要星號的話可以改成下面這行
#esttab * using correlation.rtf,label replace b(%4.3f)
nogaps unstack not noobs compress star
有時候會因為變數太多導致出現錯誤: varlist not allowed
解決方法為先輸出成.csv用excel開再複製到word
#esttab * using correlation.csv,label p(%4.3f) replace b(%4.3f)
nogaps unstack not noobs compress nostar
若您在使用上遇到問題請回覆這篇文章
以方便其他人幫忙解決問題
感恩不盡
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.119.202.234
推 erica1215:大師好厲害~~ 04/23 20:05
推 chopper625: 好厲害~~ 感謝大師!! 真的是辛苦你了.... 04/24 00:49
推 sjyshower:高手! 04/24 17:59
推 newfan:真的超感謝大師~~~可以在參考文獻寫Drewlin(2010)嗎? 04/25 16:21
推 stustu:stata魔人 我也好想用stata跑 04/25 17:10
推 happywave:大師~我想輸出spearman相關係數矩陣+p值...試不出來@@ 04/28 00:54
→ drewlin:esttab + spearman 之前用起來有點怪怪的 04/28 01:19
→ drewlin:所以就sperman xx xx,stats(rho p)再到excel整理吧 04/28 01:20
推 happywave:了解! 謝啦~ ^^ 04/28 01:30