arma模型谱估计matlab_【大数据部落】R语言多元Copula GARCH 模型时间序列预测

原文链接

R语言多元Copula GARCH 模型时间序列预测​tecdat.cn
19984283817a34fa3028fce75ce1485b.png

和宏观经济数据不同,金融市场上多为高频数据,比如股票收益率序列直观的来说 ,后者要比前者“抖动”多了有漂移且随机波动的序列,在一元或多元的情况下,构建Copula函数模型和GARCH模型是最好的选择。

多元GARCH家族中,种类非常多,需要自己多推导理解,选择最优模型。本文使用R软件对3家上市公司近十年的每周回报率为例建立模型。

9bf2e446ff70c76f8e04367ac645decb.png

首先我们可以绘制这三个时间序列。

57a4458a1ba4094240c8fe14cd604617.png

在这里使用多变量的ARMA-GARCH模型。

本文考虑了两种模型
1 ARMA模型残差的多变量GARCH过程

2 ARMA-GARCH过程残差的多变量模型(基于Copula)

1 ARMA-GARCH模型

> fit1 = garchFit(formula = ~arma(2,1)+ garch(1,1),data = dat [,1],cond.dist =“std”)

可视化波动

cf50b75be5f6789ad26c94aeb9510c14.png

隐含的相关性

> emwa_series_cor = function(i = 1,j = 2){+ if((min(i,j)== 1)&(max(i,j)== 2)){+ a = 1; B = 5; AB = 2} +}

c34a6e8cf1844702a7074ac0c7fff37c.png

2 BEKK(1,1)模型:

BEKK11(dat_arma)

99257ef1f73ac8ffc2d1a79535da76b1.png

隐含的相关性

42e5cff13c0b6f10c4915791fbb82737.png

对单变量GARCH模型残差建模

第一步可能是考虑残差的静态(联合)分布。单变量边际分布是

b36bb549d665d5cc5f0116869beaba9a.png

而联合密度为

d915b00356a402069b7cd36ebf575522.png

可视化 密度

c2b130f43b3a4b46a80db79cdb0505e7.png

d944bef8eb3ad20628cbd28558da53a4.png

查看相关性是否随着时间的推移而稳定。

46c7eb111fe2656fcd920a97f149da6e.png

斯皮尔曼相关性

fe5ecbe3a38203b6155df308b4ea974d.png

肯德尔相关性

b1f1b9699fb70254e3852b284804f8fd.png

对相关性建模,考虑DCC模型

6f5a70c7b96e0780c8d5969da783f363.png

对数据进行预测

> fcst = dccforecast(dcc.fit,n.ahead = 200)

c4e017b8e1f35a564af8a16386962b76.png

我们已经完全掌握了多元GARCH模型的使用,接下来就可以放手去用[R处理时间序列了!

相关推荐:

R语言基于ARMA-GARCH-VaR模型拟合和预测实证研究分析

ef="https://blog.csdn.net/qq_19600291/article/details/79690330?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522158652518719726867831337%2522%252C%2522scm%2522%253A%252220140713.130056874..%2522%257D&request_id=158652518719726867831337&biz_id=0&utm_source=distribute.pc_search_result.none-task-blog-blog_SOOPENSEARCH-2">GARCH(1,1),MA以及历史模拟法的VaR比较

ef="https://blog.csdn.net/qq_19600291/article/details/90519433?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522158652518719726867831337%2522%252C%2522scm%2522%253A%252220140713.130056874..%2522%257D&request_id=158652518719726867831337&biz_id=0&utm_source=distribute.pc_search_result.none-task-blog-blog_SOOPENSEARCH-3">R语言: GARCH模型股票交易量的研究道琼斯股票市场

R语言基于ARMA-GARCH过程的VaR拟合和预测

matlab实现MCMC的马尔可夫切换ARMA - GARCH模型估计

href="https://blog.csdn.net/qq_19600291/article/details/101350137?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522158652518719726867831337%2522%252C%2522scm%2522%253A%252220140713.130056874..%2522%257D&request_id=158652518719726867831337&biz_id=0&utm_source=distribute.pc_search_result.none-task-blog-blog_SOOPENSEARCH-8">R语言GARCH-DCC模型和DCC(MVT)建模估计

matlab预测ARMA-GARCH 条件均值和方差模型 https://blog.csdn.net/qq_19600291/article/details/81225151?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522158652518719726867831337%2522%252C%2522scm%2522%253A%252220140713.130056874..%2522%257D&request_id=158652518719726867831337&biz_id=0&utm_source=distribute.pc_search_result.none-task-blog-blog_SOOPENSEARCH-9

r语言预测波动率的实现:ARCH模型与HAR-RV模型

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Matlab功率谱估计的详尽分析——绝对原创 功率谱估计是信息学科中的研究热点,在过去的30多年里取得了飞速的发展。现代谱估计主要是针对经典谱估计(周期图和自相关法)的分辨率低和方差性能不好的问题而提出的。其内容极其丰富,涉及的学科和领域也相当广泛,按是否有参数大致可分为参数模型估计和非参数模型估计,前者有AR模型、MA模型ARMA模型、PRONY指数模型等;后者有最小方差方法、多分量的MUSIC方法等。 ARMA谱估计叫做自回归移动平均谱估计,它是一种模型化方法。由于具有广泛的代表性和实用性,ARMA谱估计在近十几年是现代谱估计中最活跃和最重要的研究方向之一。 二: AR参数估计及其SVD—TLS算法。 谱分析方法要求ARMA模型的阶数和参数以及噪声的方差已知.然而这类要求在实际中是不可能提供的,即除了一组样本值x(1),x(2),…,x(T)以供利用(有时会有一定的先验知识)外,再没有其它可用的数据.因此必须估计有关的阶数和参数,以便获得谱密度的估计.在ARMA定阶和参数之估计中,近年来提出了一些新算法,如本文介绍的SVD—TLS算法便是其中之一。 三:实验结果分析和展望 1,样本数多少对估计误差的影响。(A=[1,0.8,-0.68,-0.46]) 图1 上部分为N=1000;下部分为取相同数据的前N=50个数据产生的结果。 图1 N数不同:子图一N=1000,子图二N=200,子图三 N=50 由图可知,样本数在的多少,在对功率谱估计的效果上有巨大的作用,特别在功率谱密度函数变化剧烈的地方,必须有足够多的数据才能完整的还原原始功率谱密度函数。 2,阶数大小对估计误差的影响。 A=[1,-0.9,0.76] A=[1,-0.9,0.76,-0.776] 图二 阶数为二阶和三阶功率密度函数图 A=[1,-0.9,0.86,-0.96,0.7] A=[1,-0.9,0.86,-0.96,0.7,-0.74] 图三 阶数为三阶和四阶功率密度函数图 如图所示,阶数相差不是很大时,并不能对结果产生较大的影响。但是阶数太低,如图二中二阶反而不能很好的估计出原始值。 3,样本点分布对估计误差 对于相同的A=[1,-0.9,0.86,-0.96,0.7];样本的不同,在估计时的误差是不可避免的。因此,我们在取得样本时,应该尽可能的减少不必要的误差。 图四:不同的样本得到不同的估计值 4,奇异值的阈值判定范围不同对结果的影响。 上图是取奇异值的阈值大于等于0.02,而下图是取阈值大于等于0.06,显然在同种数据下,阈值的选取和最终结果有密切关系。由于系数矩阵和其真实值的逼近的精确度取决于被置零的那些奇异值的平方和。所以选取太小,导致阶数增大,选取太大会淘汰掉真实的系数。根据经验值,一般取0.05左右为最佳。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值