arma模型谱估计matlab_R语言基于ARMA-GARCH过程的VaR拟合和预测

原文链接

R语言基于ARMA-GARCH过程的VaR拟合和预测​tecdat.cn
6e9abd1650ac79b4a1c1631cdbc77f17.png

本文展示了如何基于基础ARMA-GARCH过程(当然这也涉及广义上的QRM)来拟合和预测风险价值(Value-at-Risk,VaR)。

library(qrmtools)# for qq_plot() library(rugarch)

模拟数据

我们考虑具有t的ARMA(1,1)-GARCH(1,1)过程

95754920b9366142ad6ac01789fc8243.png

e2ce092fd2cc2b50ad3ff51cc60e36c2.png

将ARMA-GARCH模型拟合到(模拟的)数据

拟合一个ARMA-GARCH过程。

961b276e677f4242ef66056e79aefd7d.png

计算VaR时间序列

计算风险价值估计值。请注意,我们也可以在这里使用基于GPD的估计器。

通过随机性检查进行后测

我们来回溯一下VaR估计值。

## Backtest VaR_0.99btest <- VaRTest(alpha,actual =X,VaR =VaR,conf.level =0.95)btest$expected.exceed# 0.99 * n ## [1] 990 btest$actual.exceed ## [1] 988 btest$uc.Decision# unconditional test decision (note: cc.Decision is NA here) ## [1] "Fail to Reject H0"

基于拟合模型预测VaR

现在预测风险价值。

模拟(X)的未来轨迹并计算相应的VaR

模拟路径,估算每个模拟路径的VaR(注意,quantile()这里不能使用,所以我们必须手动构建VaR)。

deb58610d9e1e818fd8e00b0a54082b9.png

相关文章

R语言中的风险价值模型度量指标TVaR与VaR

href="https://blog.csdn.net/qq_19600291/article/details/103491160?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522158652610019724811809250%2522%252C%2522scm%2522%253A%252220140713.130056874..%2522%257D&request_id=158652610019724811809250&biz_id=0&utm_source=distribute.pc_search_result.none-task-blog-blog_SOOPENSEARCH-2">R语言VAR模型的不同类型的脉冲响应分析

R语言用向量自回归(VAR)进行经济数据脉冲响应研究...

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

matlab使用Copula仿真优化市场风险数据VaR分析

GARCH(1,1),MA以及历史模拟法的VaR比较

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

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、付费专栏及课程。

余额充值