全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
28495 44
2011-06-21
悬赏 188 个论坛币 已解决

这个是用先从朋友的哪里找到的程序,然后再修改。朋友程序的原出处:
https://bbs.pinggu.org/viewthread.php?tid=1109721&page=1#pid9381852
这个是他在论坛的帖子。
然后,尝试规范地写成 ucsd_garch 里 dcc_mvgarch 的估计函数形式。不过有一定的错误。状态转移概率甚至会超过1,达到1.4多。
附上程序压缩文件,里面还有论文一篇,期待懂的朋友修改下,如果结果理想。奖励100论坛币。

很期待 epoh 的回复,我期待朋友们帮忙解决这个问题。

注意:这个是根据ucsd_garch 工具箱的形式修改,会调用里面的函数,使用前先确定已经装上这个工具箱,要不估计的时候,tgach等估计函数会没用,会导致错误。所以在使用前,请务必先装上Ucsd_garch toolbox。

————————————————————————————————————————————————————————————
最近忙这个程序用了将近一个月,都差点忘记了自己9月得考一个重要的从业资格证书。暂且只能到这里了。我再把最新修改的程序上传,让后来者可以继续。Epoh说得对,第一步的估计就缺少了 Markov regime swithcing 效应。 第一步调用时 stdresid=fattailed_garch(para,data);这里没有加入状态效应。等后面主优化函数,
[dccpara,]=fmincon('mrs_dcc_mvgarch_likelihood',dccstarting,……opion,stdresid,,,);等到这里调用时, Qt(t)=(1-a-b)Qbar+a*(stdresid(state)(t-1)'*stdresid(state)(t-1)+b*Qt(t-1);,此时函数里的 stdresid(state)(t-1),就缺少状态效应。所以这个程序的修改很复杂。
同时的,因为这样书写。一个是缺少一般性,强行定格在“两个状态”,程序代码臃肿,不美观,zhangtao老师也说这个程序不好阅读。这个是需要后来继续修改的人值得考虑的地方。另一个推荐是看 朱钧钧 提到的 3-state_garch 模型的程序,这个用的是 Ggriddy Gibbs sampling。这个算法用的是 Markov chain monte carlo 算法,避免了MLE 算法陷入局部最优解的可能。

————————————————————————————————————————————————————————————

MarkovDccSource 是2011年6月27日上传,也把有 swgarchlike() 这个函数的文件夹包含进去,让后来的人可以参考修改 fattailed_garch()函数,让它可以有状态转移效应。

MarkovDccSource.rar

大小:932.82 KB

 马上下载

最佳答案

epoh 查看完整内容

两个时间序列, 若先以MRS-GARCH-N估计, 除可各估计出10 parameters外 亦可由swgarchlik()估计出 LLF,likelihoods,e,p1,p1t,p1sm(smoothprob),h1,h2,h 有了这些数据, 就可依Caporin and Billio (2005) page 8/22 进行下列步骤: given the filtered probabilities evaluate the regime dependent likelihood ............. [vi] update the correlation matrix using the following approximation [vii] iterate to ...
二维码

扫码加我 拉你入群

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

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

全部回复
2011-6-21 20:19:14
两个时间序列,
若先以MRS-GARCH-N估计,
除可各估计出10 parameters外
亦可由swgarchlik()估计出
LLF,likelihoods,e,p1,p1t,p1sm(smoothprob),h1,h2,h

有了这些数据,
就可依Caporin and Billio (2005) page 8/22
进行下列步骤:
[i] given the filtered probabilities
[ii] evaluate the regime dependent likelihood
.............
[vi] update the correlation matrix using the following approximation
[vii] iterate [i] to [vi] until the end of the sample

这个部份,应该就是你着墨最深,
也是最深入的部分!
因为你的ms_dcc_mvgarch_likelihood.m
ms_dcc_mvgarch_full_likelihood.m
似乎就是依此编程的.
二维码

扫码加我 拉你入群

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

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

2011-6-23 13:57:45
会不会存在路径依赖~
二维码

扫码加我 拉你入群

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

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

2011-6-23 16:16:37
2# yucongy
我对这个也有很多疑惑。本身DCC GARCH 的估计就已经很难,算是比较那一理解的。再加上 Markov regime switching 效应,对模型的理解就都更难了。
这个方面的研究,在非多元GARCH模型里,比如一般的模型 ,GARCH/ ARCH 等,也有做 MRS的。比较著名的,一个是 汉密尔顿在1994年吧,用卡尔曼滤波技术做出了MRS GARCH 模型。这个是开创的。然后有KIM。KIM专门出了一本书,就叫 《State-space regime model》,同时也有他的论文。不过对于动态相关多元GARCH,Kim的滤波技术必须得修正,这个在资料里的论文里提到。
问题是用论文里的方程写出的代码,最终估计出的结果不理想。所以发帖请大家看看,该怎么修改。
这个运行的时候甚至是复数的。包括 logL 对数似然函数的值,不是一个实数,而是复数。
后来我又修改下,可以正常的优化,没有复数了。不过结果依然不理想,所以想请各位看看,看一下该怎么修改。

谢谢你的积极回复。
二维码

扫码加我 拉你入群

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

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

2011-6-23 20:28:38
老大,真佩服你,也只有你有这种能耐!
整个结构都变了,相当费时!


A Bivariate Markov Regime Switching GARCH Approach to
Estimate Time Varying Minimum Variance Hedge Ratios.pdf




1.State transition probabilities are assumed to
  follow a logistic distribution
  page 11/35,formula(20),formula(21)
2.conditional variances and residuals
  formula(23),formula(24)
二维码

扫码加我 拉你入群

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

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

2011-6-23 21:46:20
??? At compilation "CleanUp" was determined to be variable.  As a variable,
it is uninitialized, but function "CleanUp" is now in context.
In MATLAB 7, you cannot use the same name in a single function as
both a function and a variable.

Error in ==> Demo_MS_DCC at 14
CleanUp

??? Error: File: C:\MarkovSwitchingDCC\ms_dcc_mvgarch.m Line: 100 Column: 17
Incomplete or misformed expression or statement.

??? Error: File: C:\MarkovSwitchingDCC\ms_dcc_mvgarch_full_likelihood.m Line: 21 Column: 7
Incomplete or misformed expression or statement.

Error in ==> ms_dcc_mvgarch_likelihood at 10
[t,k]=size(stdresid);
??? Input argument "x" is undefined.

Error in ==> tvp_markov_filter at 53
[n k] = size(x);

运行你提供的程序,matlab7提示以上错误,不知是什么原因?
另外,因为我主要做金融计量,所以多次做过GARCH模型,
我看你的这个问题懂的人不多,估计还得自己多操心解决了。
当然论坛上的朋友也会尽力的。
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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