作者alexyerray (如是我聞)
看板Statistics
標題[程式] stata loop 問題
時間Wed Mar 20 22:54:01 2019
[軟體程式類別]:
stata
[程式問題]:
做bsample 跑1000次回歸取平均beta和pvalue
[軟體熟悉度]:
新手
[問題敘述]:
1. 想對100萬筆資料抽1000次，每次一萬筆資料出來跑回歸並同時記錄
1000次回歸的beta 和 p-value，前500多次都很正常，但到後面突然independent
variable 的beta都會變成0(原本想這樣算正常，因為我的independent variable是dummy
也許剛好都抽到0，都之後到1000次都是0，實在很怪
[程式碼]:
ge coef=0
ge z=0
ge pvalue=0
ge b1=0
ge p=0
forvalues i=1(1)1000 {
bsample 10000
di `i'
reg sbuy_turn_w lagevent2 i.month_num i.zip i.dow , cluster(broker_id_num)
replace b1=_b[lagevent2]
replace z = _b[lagevent2] /_se[lagevent2]
replace p = 2*normal(-abs(z))
replace coef=coef + b1
replace pvalue = pvalue + p
}
replace coef=coef/1000
replace pvalue = pvalue/1000
實在想不通為什麼，還請stata的高手幫忙，非常感謝
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.112.87.155
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1553093643.A.842.html
※ 編輯: alexyerray (140.112.87.155), 03/20/2019 23:04:11