用R语言做GARCH-MIDAS

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

想用mfGARCH包做GARCH-MIDAS，变量是收益率，协变量是月度已实现方差，想据此识别出收益率序列的长期波动和短期波动。代码是fit_mfgarch(data = mylist2, y = "return", x = "rv",low.freq = "month", K = 24)。做出来结果tau（长期波动率）和g（短期波动率）前面都是NA，请问有哪位老师知道原因吗？焦急的菜鸟博士请求高人指导

扫码加我 拉你入群

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

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

分享0 收藏3 回帖

关键词：garch-m GARCH MIDAS ARCH ARC R语言 GARCH-MIDAS

因为需要用滞后的低频因子去计算长期波动率，为了保证高频收益率能够对应一定滞后阶数的低频因子，前若干期的高频收益率样本需要被丢弃掉。如果低频因子滞后12个月，那么最开始那一年的高频收益率样本就没有被利用到，对应的长期波动率与短期波动率自然就为NA。

因此改变滞后阶数就会改变高频数据有效样本，这样一来，不同滞后阶数下的AIC与BIC就不具有可比性了。这个程序包无法根据信息准则去确定最优滞后阶数，这是它的缺陷，也是经常会被误用的地方。

1258225671 发表于 2022-4-21 11:50
因为需要用滞后的低频因子去计算长期波动率，为了保证高频收益率能够对应一定滞后阶数的低频因子，前若干期 ...

明白了，特别感谢您的回复。看论文中有一个疑惑，设定短期波动率服从GJR-GARCH(1,1)，长期波动率是正负已实现半方差的加权和tau=tau1+θ1∑RS(+)+θ2∑RS(-)。把波动率分解成长期（和短期）好的波动和坏的波动了。我的理解是正、负半方差相当于两个协变量。但是fit_mfgarch输出只有一个tau和一个g，不知道怎么分成好的波动和坏的波动的。这块仍然很疑惑。请问这方面您有所研究吗？希望得到您的指教。

于果1992 发表于 2022-4-21 14:48
明白了，特别感谢您的回复。看论文中有一个疑惑，设定短期波动率服从GJR-GARCH(1,1)，长期波动率是正负已 ...

在mfGARCH包中，长期波动方程不能使用tau=tau1+θ1∑RS(+)+θ2∑RS(-)，只能使用tau=exp(tau1+θ1∑RS(+)+θ2∑RS(-))，这样tau被分解成三部分：tau = exp(tau1)*exp(θ1∑RS(+))*exp(θ2∑RS(-))。在fit_mfgarch的返回值中，有权重函数值，可以基于此编写代码对两个低频因子分别做midas滤波并取exp。

1258225671 发表于 2022-4-21 15:11
在mfGARCH包中，长期波动方程不能使用tau=tau1+θ1∑RS(+)+θ2∑RS(-)，只能使用tau=exp(tau1+θ1∑RS(+) ...

原来如此，豁然开朗。真是太感谢您了！还想问一下，我看r package‘mfGARCH’里只有哪些函数和参数代表什么意思。您说的这个包里只能用指数形式计算长期波动率是看package里Description里的原文知晓的吗？还有就是tau=tau1+θ1∑RS(+)+θ2∑RS(-)这种形式可以在别的包里实现吗，还是只能自己编写代码？

于果1992 发表于 2022-4-21 15:33
原来如此，豁然开朗。真是太感谢您了！还想问一下，我看r package‘mfGARCH’里只有哪些函数和参数代表什 ...

我看过里面每一行源代码。需要你事先构造出这两个因子，然后用mfGARCH包估计双因子模型。其他免费的包暂时不能实现，因为它们只适用于单因子模型。

1258225671 发表于 2022-4-21 16:56
我看过里面每一行源代码。需要你事先构造出这两个因子，然后用mfGARCH包估计双因子模型。其他免费的包暂时 ...

明白了。感谢您耐心的回答，帮我厘清了思路，算是基本上了解这个模型了。我真是太幸运了，祝您工作顺利，万事如意。

通过修改mfGARCH包的源代码，可以很容易将exp去掉。取exp主要是为了防止tau变成零或负数。

1258225671 发表于 2022-4-21 17:06
通过修改mfGARCH包的源代码，可以很容易将exp去掉。取exp主要是为了防止tau变成零或负数。

知道啦，再次感谢您[em23][em23]。我再看看怎么用mfGARCH估计双因子模型