GARCH模型的估计

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

bysort id: gen t = _n
order id t
qui tab id
local N = r(r) (这一句的作用是什么？)
gen ret = firm_r - hush300
     //这样处理后，ARCH模型中就不用加常数项了，否则很难收敛

*local N = 20  // test  （为什么要这样做？）
mat R = J(`N', 2, .)  // 记录系数的矩阵
forvalues i=1/`N' {
preserve
  qui keep if id==`i'
  qui tsset t
  cap arch ret, arch(1) garch(1)  
        // 有些公司的估计可能无法收敛,所以我附加了cap命令
  if _rc == 0{   （这一句的作用是什么？）
    mat b = e(b)
    mat R[`i', 1] = b[1,2..3]
  }
restore
}
  mat colnames R = arch  garch
  mat list R

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：GARCH模型 ARCH模型 GARCH ARCH ARC local 模型 记录

我如果想计算波动率的话，该怎么改？

bysort id: gen t = _n
order id t
qui tab id
local N = r(r) (这一句的作用是什么？)
答复：记录公司数目。
gen ret = firm_r - hush300
     //这样处理后，ARCH模型中就不用加常数项了，否则很难收敛

*local N = 20  // test  （为什么要这样做？）
答复：这一句是注释，建议你看看Stata简介中的内容。
mat R = J(`N', 2, .)  // 记录系数的矩阵
forvalues i=1/`N' {
preserve
  qui keep if id==`i'
  qui tsset t
  cap arch ret, arch(1) garch(1)  
        // 有些公司的估计可能无法收敛,所以我附加了cap命令
  if _rc == 0{   （这一句的作用是什么？）
答复：上一句的注释中说明了这个处理的目的，至于 _rc，请输入 help _rc 查看说明。
    mat b = e(b)
    mat R[`i', 1] = b[1,2..3]
  }
restore
}
  mat colnames R = arch  garch
  mat list R

  if _rc == 0{   （这一句的作用是什么？）

这个呢？

help _rc

另外，想删除那些连续 T （例如 T=6）年没有资料的公司，请参考：
https://bbs.pinggu.org/thread-1171693-1-1.html
只需把其中的
gen p6 = strpos(pp, "111111")
修改为：
gen p6 = strpos(pp, "......")

进而把
drop if p6==0               // 仅保留连续六年有资料的公司
修改为：
drop if p6>0                 // 仅保留连续六年有资料的公司

谢谢，我再试试。