全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
10674 13
2021-04-23
estimation <- mvBEKK.est(bekk,order = c(1,1),params = NULL,method='Nelder-Mead')
请问用这个代码做完GARCH-BEKK后如何进行wald检验????

二维码

扫码加我 拉你入群

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

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

全部回复
2021-4-23 22:52:32
楼主您好,你用的这个函数mvBEKK.est所属的程序包现在是不是叫mgarchBEKK?如果是的话,这个包没有给参数估计后的对象提供wald检验计算的,你得自己写wald检验的代码去构造统计量然后实现检验。所以你需要提取计算对象的参数估计值coefficient estimate,还需要参数协方差矩阵coefficient covariance matrix。第一个参数估计值直接提取即可,输入estimation$est.params,第二个参数协方差矩阵输入solve(estimation\$estimation\$hessian)就能得到(打代码时把\去掉,我回答里写\是避免变成letex语法),剩下的就需要根据自己的需要构建原假设,然后去套wald检验公式写代码,构建统计量实现检验。
二维码

扫码加我 拉你入群

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

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

2021-4-25 22:00:32
719812133 发表于 2021-4-23 22:52
楼主您好,你用的这个函数mvBEKK.est所属的程序包现在是不是叫mgarchBEKK?如果是的话,这个包没有给参数估 ...
谢谢您,我在aod包中找到一个wald.test函数要求输入估计参数vector,bekk估计的参数是matrix,用这个也不行吧?只能按您提示的自己构造了
二维码

扫码加我 拉你入群

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

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

2021-4-26 14:48:40
卖笑有 发表于 2021-4-25 22:00
谢谢您,我在aod包中找到一个wald.test函数要求输入估计参数vector,bekk估计的参数是matrix,用这个也不 ...
您好,这个aod包的wald.test函数可以使用,应该能满足你的要求。你自己先把bekk估计出来的系数矩阵按照coefficient covariance matrix里元素的顺序(里面都是系数对应的协方差)给排成一个vector,可以拿coefficient covariance matrix主对角线的元素开方后的值去和estimation$asy.se.coef里的标准误做配对比较,这样就知道coefficient covariance matrix里的参数是按什么顺序排列的了,你的系数矩阵按着排成vector即可。然后用wald.test函数里的term变量去声明是哪一个参数要和哪一个参数做wald检验。wald.test函数还支持构造linear combinations of the coefficients的原假设。有些时候呢没能找到现成的函数直接计算,但是我们可以改很多中间过程来满足函数输入的要求,这样问题就可以解决了。
二维码

扫码加我 拉你入群

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

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

2021-4-26 19:25:42
719812133 发表于 2021-4-26 14:48
您好,这个aod包的wald.test函数可以使用,应该能满足你的要求。你自己先把bekk估计出来的系数矩阵按照co ...
万分感谢您赐教。根据solve(estimation\$estimation\$hessian)求出的coefficient covariance matrix的对角线元素是负的,没法进行开方处理,请问是什么原因呢?把估计的系数矩阵排序成vector后再假设股票1到股票2的波动溢出a12和b12的系数为零进行检验就可以吧?
二维码

扫码加我 拉你入群

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

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

2021-4-26 19:36:57
sqrt(abs(solve(estimation\$estimation\$hessian))),代码这样写就可以,源代码里作者是用了绝对值的,就按他这个操作吧。得到的结果就是参数估计值们对应的SE。
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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