GARCH-BEKK
10 次查看(过去 30 天)
显示 更早的评论
[EDIT: 20110617 09:06 CDT - reformat - WDR]
Hi,
I want to evaluate the volatility spill over between bonds, cds and equity using company data.
However, I have a problem with my GARCH BEKK model. I used UCSD toolbox, and followed the following steps for the estimation of the model. Built a ARMA model and obtained the residuals, then demeaned the residuals and run the GARCH BEKK model. Everything is fine so far, but the problem is that I get insignificant results for the coefficients that reflect the volatility spillover.
This is my code:
%ARMA fiting
[parameters, errors, LLF , SEregression, stderrors, robustSE, scores, likelihoods]=armaxfilter(rstock,1,1,1)
%Ljung Box test for serial correlation
[H,pValue,Qstat,CriticalValue] = lbqtest(errors,[5 ,10],0.01)
resstock=errors
resstock=resstock-mean(resstock)
%ARMA fiting
[parameters, errors, LLF , SEregression, stderrors, robustSE, scores, likelihoods]=armaxfilter(rbond,1,1,1)
%Ljung Box test for serial correlation
[H,pValue,Qstat,CriticalValue] = lbqtest(errors,[5 ,10],0.01)
resbond=errors
resbond=resbond-mean(resbond)
%ARMA fiting
[parameters, errors, LLF , SEregression, stderrors, robustSE, scores, likelihoods]=armaxfilter(rcds,1,1,1)
%Ljung Box test for serial correlation
[H,pValue,Qstat,CriticalValue] = lbqtest(errors,[5 ,10],0.01)
rescds=errors
rescds=rescds-mean(rescds)
x=[resbond rescds resstock]
%GARCH BEKK estimation
[parameters, loglikelihood,likelihoods, stdresid, stderrors, A, B, scores] = full_bekk_mvgarch(x,1,1)
%%%%%Test for adequacy
stdresid2=stdresid.^2 %square of the std resid
%Ljunk Box Test
[h,pValue,stat,cValue] = lbqtest(stdresid2,20)
qqplot(stdresid2)
%%%%Extracting std errors from A
st=diag(A,0)
st=sqrt(st)
Could anybody help me?
0 个评论
回答(3 个)
Javier
2012-12-3
编辑:Javier
2012-12-10
Hello Kaloyan
I just start checking UCSD GARCH toolbox. Once you get the H value for the lbqtest of the square residuals equal to 0, it means that the model is ok (UNIVARIATE). To get the Significance of the parameters I use this formula:
parameters/sqrt(diag(A)).
Now you have to test significance in the multivariate sense. There is no function in Matlab or MFE toolbox for that purpose. In this book you can find a multivariate test (New introduction to multiple time series analysis (Luetkepohl 2005) Hope this helps
0 个评论
Javier
2012-12-4
Hello Kaloyan
I just talk with professor Sheppard. There is a new toolbox MFE with the same functionalities. Erase UCSD Toolbox and prove with the new one.
Best regards
Mauricio
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Conditional Variance Models 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!