【R语言实战】——带有高斯新息的金融时序的GARCH模型拟合预测及VAR/ES风险度量

简介: 【R语言实战】——带有高斯新息的金融时序的GARCH模型拟合预测及VAR/ES风险度量

该篇文章主要展示了应用一个带有高斯新息的GARCH(1,1)模型,对数据进行拟合并且预测风险损失,同时进行了风险价值VaR和局部均值ES的度量,附完整代码及分析。

1 数据读取及预处理

  运行程序:

da=read.table("F:\\ch7data\\d-ibm-0110.txt",header=T)
xt=-log(da$return+1)   # calculate negative log returns.library(fGarch)

2 GARCH模型拟合

  此处为作演示,拟合GARCH(1,1)模型。

  运行程序:

library(fGarch)
m1=garchFit(~garch(1,1),data=xt,trace=F)
m1

  运行结果:

## 
## Title:
##  GARCH Modelling 
## 
## Call:
##  garchFit(formula = ~garch(1, 1), data = xt, trace = F) 
## 
## Mean and Variance Equation:
##  data ~ garch(1, 1)
## <environment: 0x00000000177235d8>
##  [data = xt]
## 
## Conditional Distribution:
##  norm 
## 
## Coefficient(s):
##          mu        omega       alpha1        beta1  
## -6.0097e-04   4.3781e-06   1.0113e-01   8.8412e-01  
## 
## Std. Errors:
##  based on Hessian 
## 
## Error Analysis:
##          Estimate  Std. Error  t value Pr(>|t|)    
## mu     -6.010e-04   2.393e-04   -2.511 0.012044 *  
## omega   4.378e-06   1.160e-06    3.774 0.000161 ***
## alpha1  1.011e-01   1.851e-02    5.463 4.67e-08 ***
## beta1   8.841e-01   1.991e-02   44.413  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Log Likelihood:
##  7114.066    normalized:  2.828654

3 模型预测

  此处预测未来三期情况。

  运行程序:

predict(m1,3)

  运行结果:

##    meanForecast   meanError standardDeviation
## 1 -0.0006009667 0.007824302       0.007824302
## 2 -0.0006009667 0.008043298       0.008043298
## 3 -0.0006009667 0.008253382       0.008253382

4 VAR、ES风险度量

  运行程序:

source("F:\\ch7data\\RMeasure.R")
m11=RMeasure(-.000601,.0078243)

  运行结果:

## 
##  Risk Measures for selected probabilities: 
##        prob        VaR         ES
## [1,] 0.9500 0.01226883 0.01553828
## [2,] 0.9900 0.01760104 0.02025244
## [3,] 0.9990 0.02357790 0.02574412
## [4,] 0.9999 0.02849770 0.03037133

  根据结果得出拟合的模型为:

image.png

  在5%的水平上,所有系数估计是统计上显著的。在时刻T=2515,均值和波动率的超前一步预测分别为6.01×1047.82×103 。相应地,我们有:

image.png

  然后可以应用上述结果计算该金融头寸的风险度量。接下来的一个交易日有:

image.png

相关文章
|
30天前
【R语言实战】——带有新息为标准学生t分布的金融时序的GARCH模型拟合预测
【R语言实战】——带有新息为标准学生t分布的金融时序的GARCH模型拟合预测
|
30天前
|
移动开发 数据可视化
广义线性模型beta二项分布的淋巴结疾病风险预测可视化R语言2实例合集|附数据代码
广义线性模型beta二项分布的淋巴结疾病风险预测可视化R语言2实例合集|附数据代码
|
30天前
|
数据可视化
R语言汇率、股价指数与GARCH模型分析:格兰杰因果检验、脉冲响应与预测可视化
R语言汇率、股价指数与GARCH模型分析:格兰杰因果检验、脉冲响应与预测可视化
|
30天前
|
数据可视化 数据挖掘 API
【R语言实战】聚类分析及可视化
【R语言实战】聚类分析及可视化
|
30天前
|
机器学习/深度学习 数据可视化
R语言逻辑回归logistic模型ROC曲线可视化分析2例:麻醉剂用量影响、汽车购买行为2
R语言逻辑回归logistic模型ROC曲线可视化分析2例:麻醉剂用量影响、汽车购买行为
|
30天前
|
Web App开发 数据可视化 数据挖掘
利用R语言进行聚类分析实战(数据+代码+可视化+详细分析)
利用R语言进行聚类分析实战(数据+代码+可视化+详细分析)
|
30天前
|
数据采集 数据可视化
利用R语言进行因子分析实战(数据+代码+可视化+详细分析)
利用R语言进行因子分析实战(数据+代码+可视化+详细分析)
|
30天前
利用R语言进行典型相关分析实战
利用R语言进行典型相关分析实战
|
30天前
|
机器学习/深度学习 算法
R语言分类回归分析考研热现象分析与考研意愿价值变现
R语言分类回归分析考研热现象分析与考研意愿价值变现
|
30天前
|
数据可视化 定位技术
R语言贝叶斯INLA空间自相关、混合效应、季节空间模型、SPDE、时空分析野生动物数据可视化
R语言贝叶斯INLA空间自相关、混合效应、季节空间模型、SPDE、时空分析野生动物数据可视化