楼主: hante865
36759 42

[R] DCC-GARCH模型R程序实现   [推广有奖]

41
卓德啊涵 在职认证  学生认证  发表于 2020-12-23 17:07:51
您好,图我做出来了,但是横坐标不是时间,而是从1到n的数据量,请问怎么把最后的DCC图横坐标设置成时间呀,自己找了两天了还是弄不好。

图如下: 1608713868(1).jpg






  1. #对收益率序列进行DCC-GARCH模型的拟合,选择0-4绘制相关图示。0表示退出,其中4是最重要的,是两者之间的动态相关系数时序图。
  2. #本次数据有2个样本序列,两两组合,共1种
  3. #install.packages('rmgarch')
  4. library(rmgarch)

  5. meanSpec=list(armaOrder=c(1,1),include.mean=FALSE,archpow=1,arfima = FALSE, external.regressors = NULL, archex = FALSE)
  6. distSpec=c("mvnorm")
  7. varSpec=list(model="sGARCH",garchOrder=c(1,1))
  8. spec1=ugarchspec(mean.model=meanSpec,variance.model=varSpec)

  9. mySpec=multispec(replicate(2,spec1))

  10. mspec=dccspec(mySpec,VAR=F,robust = F,lag=1,lag.max=NULL,lag.criterion = c("AIC"),external.regressors = NULL,robust.control = list(gamma=0.25,delta=0.01,nc=10,ns=500),dccOrder = c(1,1),distribution = distSpec,start.pars = list(),fixed.pars = list(),model = "DCC")
  11. ##上述是方程参数估计准备过程
  12. fdcc12=dccfit(data=p1ld[,c(1,2)],mspec,out.sample = 0,solver = "solnp",solver.control = list(),fit.control = list(eval.se=TRUE,stationary=TRUE,scale=FALSE),parallel=TRUE,parallel.control=list(pkg=c("multicore"),cores=2),fit=NULL,VAR.fit=NULL,cluster = NULL)
  13. show(fdcc12) #dcca1 和dccb1 显著,表明存在DCC效应。(行业联动性、时变相关性)
  14. plot(fdcc12)
复制代码







1608713780(1).jpg (209.41 KB)

1608713780(1).jpg

1608713780(1).jpg (209.41 KB)

1608713780(1).jpg

42
13880480178 学生认证  发表于 2022-4-7 11:09:19
卓德啊涵 发表于 2020-12-23 17:07
您好,图我做出来了,但是横坐标不是时间,而是从1到n的数据量,请问怎么把最后的DCC图横坐标设置成时间呀, ...
我也遇到这个问题,请问您找到解决方法了吗,谢谢

43
卓德啊涵 在职认证  学生认证  发表于 2022-4-9 09:43:38
13880480178 发表于 2022-4-7 11:09
我也遇到这个问题,请问您找到解决方法了吗,谢谢
没有,如果你摸索出来记得告我一声

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群
GMT+8, 2025-12-29 00:20