- 阅读权限
- 255
- 威望
- 3 级
- 论坛币
- 34526 个
- 通用积分
- 55890.0311
- 学术水平
- 3538 点
- 热心指数
- 3698 点
- 信用等级
- 2921 点
- 经验
- 308755 点
- 帖子
- 14503
- 精华
- 12
- 在线时间
- 7393 小时
- 注册时间
- 2016-7-9
- 最后登录
- 2024-4-27
大师
还不是VIP/贵宾
- 威望
- 3 级
- 论坛币
- 34526 个
- 通用积分
- 55890.0311
- 学术水平
- 3538 点
- 热心指数
- 3698 点
- 信用等级
- 2921 点
- 经验
- 308755 点
- 帖子
- 14503
- 精华
- 12
- 在线时间
- 7393 小时
- 注册时间
- 2016-7-9
- 最后登录
- 2024-4-27
|
经管之家送您一份
应届毕业生专属福利!
求职就业群
感谢您参与论坛问题回答
经管之家送您两个论坛币!
+2 论坛币
根据 xtqreg 之原作者之建议,可用 - // https://jmcss.som.surrey.ac.uk/MM-QR-JK.do
- webuse nlswork, clear
- xtset idcode
- by idcode: egen c=count(idcode)
- set seed 123
- gen s=runiform()<0.5
- *** Simple MM-QR estimation with fixed effects
- xtqreg ln_w age ttl_exp tenure not_smsa south if c>=10, ls q(0.10(0.10)0.90)
- *** MM-QR results with JK correction and clustered standard errors
- bs if c>=10, cluster(idcode) rep(100): xtqreg ln_w age ttl_exp tenure not_smsa south
- mat V=e(V)
- mat b=e(b_location)
- mat g=e(b_scale)
- mat Q50=e(q)
- scalar N=e(N)
- qui xtqreg ln_w age ttl_exp tenure not_smsa south if c>=10&s==0
- scalar N0=e(N)
- mat Q050=e(q)
- mat g0=e(b_scale)
- qui xtqreg ln_w age ttl_exp tenure not_smsa south if c>=10&s==1
- scalar N1=e(N)
- mat Q150=e(q)
- mat g1=e(b_scale)
- mat bjk=b+(2*g-g1*(N1/N)-g0*(N0/N))*(2*Q50-Q150*(N1/N)-Q050*(N0/N))
- ereturn post bjk V
- ereturn display
复制代码
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
|
|
-
总评分: 学术水平 + 3
热心指数 + 3
信用等级 + 3
查看全部评分
|