作者typhoon20025 (typhoon20025)
看板Statistics
標題[程式] Stata的三個問題
時間Mon May 18 00:55:43 2015
------------------------------------------------------------------------
[軟體程式類別]:
Stata
[程式問題]:
資料處理、迴歸、敘述統計
[軟體熟悉度]:
低(1~3個月)
[問題敘述]:
問題一：
需要將資料分產業(ind，8個)及年份(year，12年)跑回歸式，並保留回歸式的殘差項作
為變數，利用產生新變數yearind=year*100+ind，寫迴圈方式if yearind=k，就run
regression
程式如程式範例的問題一，
出現了unrecognized command: estpostr(199)的問題，
請問語法應該要怎麼做修正呢?
問題二：
(1)敘述性統計
(2)相關矩陣
以上兩個程式如範例問題二，同樣出現unrecognized command: estpostr(199)的問題
，請問語法應該要怎麼做修正呢?
問題三：
(1)t檢定及(2)異質修正(heteroscedasticity corrected)程式如範例問題三
tse變數出現ambiguous abbreviationr(111)訊號，tse24 tse 25 tse26 tse27 tse28
tse29 tse30為虛擬變數不知道語法應該要如何做修正呢?
[程式範例]:
問題一：
est clear
estpost tab ind
esttab using regression.rtf, cells("b(label(freq)) pct(fmt(2))
cumpct(fmt(2))") varlabels(, blist(Total)) nonumber noobs replace
esttab using regression.csv, cells("b(label(freq)) pct(fmt(2))
cumpct(fmt(2))") varlabels(, blist(Total)) nonumber noobs replace
est clear
estpost tab year
esttab using regression, cells("b(label(freq)) pct(fmt(2)) cumpct(fmt(2))")
varlabels(, blist(Total)) nonumber noobs replace
esttab using regression, cells("b(label(freq)) pct(fmt(2)) cumpct(fmt(2))")
varlabels(, blist(Total)) nonumber noobs replace
gen dc_a= .
foreach k of local values {
reg ta_lag1a vlag1a drev_lag1_drec_lag1a ppe_lag1a if yearind== `k'
predict residuals, resid, if yearind== `k'
replace dc_a=residuals if resid!=.
drop residuals
問題二：
(1)敘述性統計
est clear
local vars "troi_a treps trroe eso ccq"
estpost sum `vars', detail
esttab using descriptive_statistics.rtf, cells("mean(fmt(3)) sd(fmt(2))
min(fmt(1)) max(fmt(0)) p25(fmt(2)) p50(fmt(2)) p75(fmt(2)) skewness(fmt(2))
kurtosis(fmt(2))") replace
esttab using descriptive_statistics.csv, cells("mean(fmt(3)) sd(fmt(2))
min(fmt(1)) max(fmt(0)) p25(fmt(2)) p50(fmt(2)) p75(fmt(2)) skewness(fmt(2))
kurtosis(fmt(2))") replace
(2)相關矩陣
est clear
pwcorr `vars'
qui estpost correlate `vars', matrix
esttab using correlation_matrix.rtf, unstack not noobs compress star(* 0.10
** 0.05 *** 0.01) b(3) replace
esttab using correlation_matrix.csv, unstack not noobs compress star(* 0.10
** 0.05 *** 0.01) b(3) replace
問題三：
(1)t檢定:
reg troi_a eso_a s_a tse24 tse 25 tse26 tse27 tse28 tse29 tse30
estat hettest
tsset ind year
xtserial troi_a eso_a s_a tse24 tse 25 tse26 tse27 tse28 tse29 tse30
est clear
eststo: reg troi_a eso_a s_a tse24 tse 25 tse26 tse27 tse28 tse29 tse30
esttab , replace starlevels(* .1 ** 0.05 *** 0.01)
*輸出在excel ".csv" / 輸出在word ".rtf"
esttab * using test_H1-1.rtf ,replace
esttab * using test_H1-1.csv ,replace
(2)異質修正:
newey troi_a eso_a s_a tse24 tse 25 tse26 tse27 tse28 tse29 tse30
不好意思是第一次發文 排版有點亂請見諒 再麻煩高手幫忙解答 謝謝~~
-----------------------------------------------------------------------------
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 123.195.202.61
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1431881746.A.E29.html
※ 編輯: typhoon20025 (123.195.202.61), 05/18/2015 00:58:11
1F：→ porfu: 第三個問題看來出在"tse 25"。Stata猜不出tse是哪個變數 05/18 09:40
2F：→ porfu: 前面兩點的話，錯誤訊息說的是有一個estpostr的命令被輸入 05/18 09:42
3F：→ porfu: 檢查一下語法是不是都有拼對。 05/18 09:43
4F：推 evilove: 會不會沒有裝 estpost , 先 help estpost看一下 05/19 07:42
5F：→ evilove: tse 25 中間有空格 05/19 07:43
6F：→ typhoon20025: !!!問題解決了~~~非常謝謝樓上兩位幫忙~~~謝謝你們~ 05/21 01:20