[STATA] regression model、敘述統計、相關係數輸出製表

看板NCCU08_MAT作者drewlin (小卓..)時間15年前 (2010/04/23 19:50)推噓7(7推 0噓 2→)

留言9則, 7人參與討論串1/1

第一次使用時請先更新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, , 1^F

04/23 20:05, 1^F

推

chopper625

04/24 00:49, , 2^F

04/24 00:49, 2^F

推

sjyshower

04/24 17:59, , 3^F

04/24 17:59, 3^F