|
你好,不知道你这个问题解决了没有。
要得到残差的话就follow 下面的例子就可以了
garch(p=1,q=1,mv=bek,dist=ged,iters=5000, robust,hmatrices=hd,rvectors=rd) / x1 x2
set resid1 %regstart() %regend() = rd(t)(1)
set resid2 %regstart() %regend() = rd(t)(2)
以上是得到残差,如果想得到estimate出来的方差就是如下命令:
set variance1 %regstart() %regend() = hd(t)(1,1)
set variance2 %regstart() %regend() =hd(t)(2,2)
如果想得到标准化了的残差,就是如下命令:
set z1 = rd(t)(1)/sqrt(hd(t)(1,1))
set z2 = rd(t)(2)/sqrt(hd(t)(2,2))
命令的原理就是存有你残差和方差-协方差的矩阵里调出你要的残差和方差。
所以如果你想要协方差的话,那命令就是:
set cov12 %regstart() %regend() = hd(t)(1,2)
同理,相关系数就是
set corr12 %regstart() %regend() = hd(t)(1,2)/sqrt(hd(t)(1,1)*hd(t)(2,2))
希望对你帮助
另外要作图的话,就在rats里面直接graph就好了,命令栏里面可以找到
|