全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
8108 14
2021-11-24
求教各位大佬,用rugarch包拟合偏t分布时,结果中的skew和shape参数应该是偏t分布的参数吧?如下图

那应该怎么求该偏t分布(skew=1.143593,shape=8.87965)下的95%分位数或临界值呢?
我看有人推荐gamlss包里的qSST函数,但是函数需要4个参数(mu,sigma,nu,tau),对不上啊?或且应该怎么算出来呢?提前感谢大佬的解答
附件列表
1.jpg

原图尺寸 111.61 KB

1.jpg

二维码

扫码加我 拉你入群

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

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

全部回复
2021-11-25 17:50:11
你要计算ARMA(1,1)-EGARCH(1,1)极大似然估计结果下这个sstd分布的分位数是不是? rugarch包的sstd用的是Fernande and Steel的偏t分布公式,如果你用gamlss包里的qSST函数,需要先选择正确的偏t分布类型,因为gamlss包里的这个函数给了5个不同的偏t分布类型设定,都是不同的学者提出的,不要用错了,那理论和计算就会出现偏差。rugarch包用的是Fernande and Steel的偏t分布。你要计算偏t分布的分位数,是不是需要知道这个偏t分布的mean均值,sd方差,自由度, 偏度?所以这里一定要有4个参数输入,前两个是均值方差,默认分别是0和1,后面两个我猜应该分别是偏度nu,自由度tau,具体你要去看gamlss包的说明文档。rugarch模型估计完只会得到skew偏度还有shape自由度,均值和方差都是预设进建模过程了(0, 1),变量叫什么名字只是个皮,形参名字随便叫,但是输进函数里的实参变量本身是偏度还是峰度一定要正确。我建议你用fGARCH包的sstd函数中的qsstd()函数就可以,这个就是Fernande and Steel的规格,比较方便,这个函数的nu是自由度,xi是偏度,详细请参照fGARCH包说明文档。

btw,你这个计算结果不能用,EGARCH(1,1)里alpha1这个负责sign effect的变量估计值不够显著,普通的hessian inverse标准误和稳健标准误之下的t检验都是这样,建议你换一个不考虑波动率不对称性的GARCH去使用,你的数据可能不适合使用非对称GARCH进行建模分析。
二维码

扫码加我 拉你入群

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

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

2021-11-26 17:51:25
719812133 发表于 2021-11-25 17:50
你要计算ARMA(1,1)-EGARCH(1,1)极大似然估计结果下这个sstd分布的分位数是不是? rugarch包的sstd用的是Fer ...
感谢大神,我去试一下
二维码

扫码加我 拉你入群

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

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

2021-11-26 18:13:14
719812133 发表于 2021-11-25 17:50
你要计算ARMA(1,1)-EGARCH(1,1)极大似然估计结果下这个sstd分布的分位数是不是? rugarch包的sstd用的是Fer ...
大神,我还想问一下,alpha1系数不显著我可不可以直接在结果中将该项系数写成零(或者我看到有文章直接把alpha1系数写在结果中,即使显著性检验未通过)?另一个问题我也不懂,是否用对称GARCH或非对称GARCH不是应该看gamma1系数是否通过显著性检验么?
二维码

扫码加我 拉你入群

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

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

2021-11-26 18:49:04
Aixir 发表于 2021-11-26 18:13
大神,我还想问一下,alpha1系数不显著我可不可以直接在结果中将该项系数写成零(或者我看到有文章直接把 ...
alpha1系数不显著不可以直接在结果里写0,EGARCH(1,1)又不是linear regression model解释变量不显著可以drop掉然后再换其他的,alpha1你要是写0那EGARCH这个模型规格就整个就不存在了。如果这个EGARCH(1,1)的结果不是你最优的模型,即准备拿来分析的模型结果,那你写进文章里当个robustness check还说得过去,若是准备当主要模型结果来用的,那就不可以,关键变量t检验通不过,审稿肯定问你。

在rugarch包的EGARCH(1,1)里,gamma1这个变量负责的是size effect,当然不能只看这一个变量就决定EGARCH用与不用,alpha1也很重要。

是否用对称与非对称GARCH要取决于具体的模型规格的参数估计结果,非对称GARCH有好几个,gamma1具体指哪个模型规格里的哪一个参数要自己去看,不能只认死gamma1这个变量名,具体模型规格变量命名一切以你用的程序包定义为准,你可以去这个链接Introduction_to_the_rugarch_package.pdf - 下载 - 人大经济论坛 (pinggu.org)下载这份文档,里面详细说明了rugarch所用的模型公式,看完你就明白了。
二维码

扫码加我 拉你入群

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

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

2021-11-28 11:22:38
719812133 发表于 2021-11-26 18:49
alpha1系数不显著不可以直接在结果里写0,EGARCH(1,1)又不是linear regression model解释变量不显著可以d ...
哦哦,感谢大佬解惑,我知道了
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群