作者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/cn.aspx?n=bbs/Statistics/M.1553093643.A.842.html
※ 编辑: alexyerray (140.112.87.155), 03/20/2019 23:04:11