R语言用GARCH模型波动率建模和预测、回测风险价值 (VaR)分析股市收益率时间序列|附代码数据

原文链接：http://tecdat.cn/?p=26897

最近我们被客户要求撰写关于GARCH的研究报告，包括一些图形和统计输出。

风险价值 (VaR) 是金融风险管理中使用最广泛的市场风险度量，也被投资组合经理等从业者用来解释未来市场风险

风险价值 (VaR)

• VaR 可以定义为资产在给定时间段内以概率 θ 超过的市场价值损失。对于收益率 rt 的时间序列，VaRt将是这样的

其中 It-1表示时间 t-1 的信息集。

• 尽管 VaR 在提供资产组合下行风险的简单总结时具有吸引人的简单性，但没有单一的计算方法。

1% 风险价值

• 将价格转换为收益

library(ggplot2)
# 计算收益率的正态密度

# 价格与收益的关系

bp2 = Close
# 转换收益率

bret = dailyReturn
# 改变列名
colnames(data_rd) = c("x", "y")
# 正态分位数
vr1 = quantile
 ggplot(data, aes(x = x, y = y))

图 ：1% VaR

---

点击标题查阅往期内容

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

左右滑动查看更多

01

02

03

04

• 在分布术语中，对于分布 F，VaR 可以定义为它的第 p 个分位数，由下式给出

其中 F−1是分布函数的倒数，也称为分位数函数。因此，一旦可以定义收益序列的分布，VaR 就很容易计算。

使用 GARCH 进行波动率建模和预测

• 广义自回归条件异方差 (GARCH) 模型 ，用于预测条件波动率的最流行的时间序列模型。

• 这些模型是条件异方差的，因为它们考虑了时间序列中的条件方差。GARCH 模型是在金融风险建模和管理中用于预测 VaR 和条件 VaR 等金融风险度量的最广泛使用的模型之一。

• GARCH 模型是 ARCH 模型的广义版本。具有旨在捕获波动率聚类的 p 滞后项的标准 ARCH(p) 过程可以编写如下

  

  其中，第 t 天的收益为 Yt=σtZt和 Zt∼iid(0,1)，即收益的创新是由随机冲击驱动的

• GARCH(p,q) 模型在 ARCH(p) 模型中包含滞后波动率，以纳入历史收益的影响

  

• GARCH(1,1) 每个阶数只使用一个滞后，是实证研究和分析中最常用的版本。

   GARCH(1,1) 预测 VaR

• 其中最通用和最有能力的一种是 rugarch 包。在这里，我们使用数据集来演示使用 rugarch 包中可用的函数和方法对 GARCH 进行建模。

• 具有恒定均值方程的 GARCH(1,1) 模型 可以指定如下：

  ugarchspec(variance.model = list(model = "sGARCH", garchOrder = c(1,
      1)), mean.model = list(armaOrder = c(0, 0)))

• 上面存储的规范 garch_spec 现在可用于将 GARCH(1,1) 模型拟合到我们的数据。以下代码使用该函数将 GARCH(1,1) 模型拟合到 BHP 对数收益并显示结果。

  

  

• 使用对象类可用的各种方法获得选定的拟合统计量

  par1 = par() #保存图形参数

 标准化残差

plot(figarch, which = 10)

 2. 条件SD 

plot(fiarch, which = 3)

![图片](https://p3-juejin.byteimg.com/tos-cn-i-k3u1fbpfcp/c3a8a2e609ae453e8d3f7d0d6a8442e1~tplv-k3u1fbpfcp-zoom-1.image)

图 ：GARCH(1,1) 的两个信息图

# 使用样本外的 VaR 预测 

  


-   让我们使用 Student-t 分布，因为收益并不总是遵循正态分布

 学生-T分布的spec2

spc2 = ugarchspec

  


-   rugarch 包对于估计移动窗口模型和预测 VaR 具有非常有用的功能。

garchroll(spec2, data = bpret

  


-   我们可以使用以下例程绘制 1% 和 5% VaR 预测与实际收益的对比。

 注意绘图方法提供了四张图，其中VaR为选项-4

 预测1%的学生-t GARCH风险值

plot(v.t, which = 4, VRaha = 0.01)

 5%学生-t GARCH风险值

plot(var.t, which = 4, Vaalha = 0.05)

![图片](https://p3-juejin.byteimg.com/tos-cn-i-k3u1fbpfcp/10430383887348dba79364f6d7fd22ca~tplv-k3u1fbpfcp-zoom-1.image)

图：实际收益率与 1% VaR 预测

  


-   最后获得回测

 VaR预测的回测

report(va., VaRha = 0.05)  #α的默认值是0.01

![图片](https://p3-juejin.byteimg.com/tos-cn-i-k3u1fbpfcp/f6fee11190d74763a5d5a4c478070bbe~tplv-k3u1fbpfcp-zoom-1.image)

  


* * *

  


![图片](https://p3-juejin.byteimg.com/tos-cn-i-k3u1fbpfcp/6ca38b6a1d5040c7a26f8f6a2d7e7dd6~tplv-k3u1fbpfcp-zoom-1.image)

点击文末 **“阅读原文”**

获取全文完整资料。

  


本文选自《**R语言用GARCH模型波动率建模和预测、回测风险价值 (VaR)分析股市收益率时间序列**》。

  


**点击标题查阅往期内容**

[R语言使用多元AR-GARCH模型衡量市场风险](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247515716&idx=5&sn=a2d70d76364238dd3231e69674d73906&chksm=fd92804fcae5095992611e228b0470d478a25172f06fee0d76a536246605b076d2b764d1b718&scene=21#wechat_redirect)  
[R语言GARCH模型对股市sp500收益率bootstrap、滚动估计预测VaR、拟合诊断和蒙特卡罗模拟可视化](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247513800&idx=1&sn=2881eac305059dfb6d077bad53718f7b&chksm=fd9288c3cae501d5c878ea0e350e5ac16838e60ab863f8acdb7f87a5bbd0a44689966b0e5bc5&scene=21#wechat_redirect)[R语言单变量和多变量（多元）动态条件相关系数DCC-GARCH模型分析股票收益率金融时间序列数据波动率](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247512552&idx=1&sn=cd0f2de9c840592f17bb735714eccbd1&chksm=fd928fe3cae506f56a4cc4a7825d92f43465e00013d2d042a863d3c3b278d095711ec0d261d5&scene=21#wechat_redirect)  
[R语言中的时间序列分析模型：ARIMA-ARCH / GARCH模型分析股票价格](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247512993&idx=3&sn=88c6792a7eb065dc82d1e0eb00f5f72e&chksm=fd928daacae504bc604e9ad42409b5535a7ca12ceb47636bf1a91684b2d9423a502b6b24a0d3&scene=21#wechat_redirect)  
[GARCH-DCC模型和DCC（MVT）建模估计](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247513518&idx=3&sn=a7584ee8a5940a63d5f6852ff052c474&chksm=fd928ba5cae502b33147fc58d0419817ec3d21b9ef7dab57c646f4e7160e7774765c7e6a3282&scene=21#wechat_redirect)  
[R语言预测期货波动率的实现：ARCH与HAR-RV与GARCH，ARFIMA模型比较](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247510231&idx=2&sn=77e9e264e811bd27312703c14f8e1654&chksm=fd9296dccae51fcaef5050b793a056da2b2f93c051639003543fc11cc714c4167e88b2a975d1&scene=21#wechat_redirect)  
[ARIMA、GARCH 和 VAR模型估计、预测ts 和 xts格式时间序列](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247510062&idx=1&sn=943d6d3b043942222bfa2c3527d5626a&chksm=fd929625cae51f33d54e7c8ea01702e366ca5a56561f8bea1e79dfb39637435671ade4713787&scene=21#wechat_redirect)  
[PYTHON用GARCH、离散随机波动率模型DSV模拟估计股票收益时间序列与蒙特卡洛可视化](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247509977&idx=1&sn=68450b5bc3e8f903c83f96cd9f00d1aa&chksm=fd9299d2cae510c46304710b275fe24e4403052cf010f643b12fd84d58c322721dc8bf37d3e4&scene=21#wechat_redirect)  
[极值理论 EVT、POT超阈值、GARCH 模型分析股票指数VaR、条件CVaR：多元化投资组合预测风险测度分析](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247509289&idx=3&sn=e40a8e7a8ea7c48d56f5ab32b6cc9cbd&chksm=fd929b22cae51234553e6288d34cc071945f5fe56d988bab259bc7ffea3199c57eacd85dee75&scene=21#wechat_redirect)  
[Garch波动率预测的区制转移交易策略](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247508741&idx=3&sn=c49f5f8caab501e2980ee51d509c3276&chksm=fd929d0ecae51418e3791dbb9ff0cb27e60ec40e8e0f3a1d156fd9abb54b826c5f8479e706a2&scene=21#wechat_redirect)  
[金融时间序列模型ARIMA 和GARCH 在股票市场预测应用](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247508260&idx=2&sn=8206765e1e5a18a16909c2084ef5b8b7&chksm=fd929f2fcae51639e69c5e1c9e6f764611238bd347e58f5f85fba093e40c88bf75690f7f26f0&scene=21#wechat_redirect)  
[时间序列分析模型：ARIMA-ARCH / GARCH模型分析股票价格](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247508136&idx=2&sn=efcb350e6e73a6b2ce1a5e4ed5816158&chksm=fd929ea3cae517b581e174ce39c37f3319644e2f0a24d4c713f04713d112f2416fdaa767f127&scene=21#wechat_redirect)  
[R语言风险价值：ARIMA，GARCH，Delta-normal法滚动估计VaR（Value at Risk）和回测分析股票数据](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247506706&idx=1&sn=270827c906243e4bd2fe71e2ba1bde10&chksm=fd92e519cae56c0f662e32c906ffef042386ea30deb7b2794b29dc66ce007159e35b6c3eeca6&scene=21#wechat_redirect)  
[R语言GARCH建模常用软件包比较、拟合标准普尔SP 500指数波动率时间序列和预测可视化](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247506372&idx=1&sn=ddc017011491d0f692b09aa6fda26083&chksm=fd92e7cfcae56ed97bc16105e9e88021a5e7a1c4f58d9f60fb54f7f176eb2f3df4d517b0869f&scene=21#wechat_redirect)  
[Python金融时间序列模型ARIMA 和GARCH 在股票市场预测应用](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247506099&idx=1&sn=4fb892131945deaca373fce7f30fae1b&chksm=fd92e6b8cae56fae222abd07001eeddbe2ac3b1d61ea9aaf17186c35852ddbb0ab8426678698&scene=21#wechat_redirect)  
[MATLAB用GARCH模型对股票市场收益率时间序列波动的拟合与预测](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247505674&idx=1&sn=84aa62f91e4b098c7f8d7dcd93a4b1cd&chksm=fd92e901cae56017ad05da51d9fb043db036cbde047231220faea57f608caba7d0ec5a9418df&scene=21#wechat_redirect)  
[R语言极值理论 EVT、POT超阈值、GARCH 模型分析股票指数VaR、条件CVaR：多元化投资组合预测风险测度分析](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247505427&idx=1&sn=04ce1c89c632fed1a52a4daeb783b3c3&chksm=fd92e818cae5610edce52b308a6d60c6c5fbd83f139fb7dc0a2e0cf33b06fe701962f46369b7&scene=21#wechat_redirect)  
[Python 用ARIMA、GARCH模型预测分析股票市场收益率时间序列](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247504283&idx=1&sn=e7737d2c3b6e1014493a56d9001c1c40&chksm=fd92ef90cae56686d6da7104d2d0f1aa62c0927094b9c94acafd40b09002f40c9371c0acc623&scene=21#wechat_redirect)  
[R语言中的时间序列分析模型：ARIMA-ARCH / GARCH模型分析股票价格](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247503984&idx=2&sn=8ad8573f1b609571b4a1a45a2752f794&chksm=fd92ee7bcae5676dce10814f29939ddadd742e2faeee974043bd403f2fc8383b7424f1f2a2db&scene=21#wechat_redirect)  
[R语言ARIMA-GARCH波动率模型预测股票市场苹果公司日收益率时间序列](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247503337&idx=1&sn=79370847577f9fafb3dfe56b40d953c8&chksm=fd92f3e2cae57af40d5a336ae63e1f7f0be36d6d24ebff81de173514afe0c5f0f04471dfd4bb&scene=21#wechat_redirect)  
[Python使用GARCH，EGARCH，GJR-GARCH模型和蒙特卡洛模拟进行股价预测](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247497021&idx=2&sn=ed0ab6559b75da51504e32772f186eb3&chksm=fd92cb36cae54220ef77066e1d0f7470352864d9f14f061041d69b1fe03d42eaf5e356d6a521&scene=21#wechat_redirect)  
[R语言时间序列GARCH模型分析股市波动率](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247496697&idx=1&sn=26545eb6e70e15e99033ee28fc7dc20c&chksm=fd92cdf2cae544e4a97140229506c5b1a2420fb0f156ac161c81e3d94cba8bf44cd545a90412&scene=21#wechat_redirect)  
[R语言ARMA-EGARCH模型、集成预测算法对SPX实际波动率进行预测](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247496554&idx=2&sn=0edf01886135034bcc152660b0c20c54&chksm=fd92cd61cae54477cfe495d7a64c71f198efe3b1c56625193c39c68fc0e048818313dcea790e&scene=21#wechat_redirect)  
[matlab实现MCMC的马尔可夫转换ARMA - GARCH模型估计](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247495866&idx=2&sn=360f932fa15a7bbdc31c9efecc15a861&chksm=fd92ceb1cae547a7036e2137d5fe28db51deeaf15c4c007e402e0d18d1ba602b14360632ebcd&scene=21#wechat_redirect)  
[Python使用GARCH，EGARCH，GJR-GARCH模型和蒙特卡洛模拟进行股价预测](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247494839&idx=1&sn=72686526dda3fb8b85337d30295ae2bc&chksm=fd92d2bccae55baa5f9c9d3a40ce0653c917bde6b0e060de2c662d8820a0045faf9c819f8304&scene=21#wechat_redirect)  
[使用R语言对S＆P500股票指数进行ARIMA + GARCH交易策略](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247494760&idx=2&sn=e12cca3b8bfd7864b6231228a451ff4b&chksm=fd92d263cae55b75148ac82aebf72fe8962d46f26ef98c3ddc7871b45250735e98c558227d97&scene=21#wechat_redirect)  
[R语言用多元ARMA,GARCH ,EWMA, ETS,随机波动率SV模型对金融时间序列数据建模](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247494103&idx=1&sn=c5dcaaa27ad884f83e8e40edb41ebbfd&chksm=fd92d7dccae55eca8a2f2a333544ee9a35f9744ec66f54edeaa82488bc56465ff6d631fe30ff&scene=21#wechat_redirect)  
[R语言股票市场指数：ARMA-GARCH模型和对数收益率数据探索性分析](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247493463&idx=1&sn=b5d3f36fe60e19c36c7ebd81e50847e3&chksm=fd92d95ccae5504a4a87b208ff7f1311716333a8ee48ef48d526868bd0dcf7aaaabf103dcfbb&scene=21#wechat_redirect)  
[R语言多元Copula GARCH 模型时间序列预测](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247493270&idx=3&sn=65f333f31222b17f35717bab7d145c2c&chksm=fd92d89dcae5518b7ab0adeefc873ff40b8d683085dfffddba321576b5f81c895af5f3b93e65&scene=21#wechat_redirect)  
[R语言使用多元AR-GARCH模型衡量市场风险](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247493156&idx=1&sn=e7e2b5dfca9c78e870176d3b5e4a796d&chksm=fd92d82fcae5513930bb452b1cbb784b551b976ad11cafd61bfefcb92c6aeb9a7527a37ae17b&scene=21#wechat_redirect)  
[R语言中的时间序列分析模型：ARIMA-ARCH / GARCH模型分析股票价格](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247492792&idx=1&sn=e729d58a15a86b04d54df8116c0c428a&chksm=fd92dab3cae553a57aa2af43544a36b56ce29e591752facba174d3efc246eb3957fa2dcbd01c&scene=21#wechat_redirect)  
[R语言用Garch模型和回归模型对股票价格分析](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247492215&idx=1&sn=f59aa15b8af6ca1edfc551c0d568b241&chksm=fd92dc7ccae5556aaa3324947b91c0ede914847d9f69498d1f2d06ee9caa063a085679eec396&scene=21#wechat_redirect)  
[GARCH（1,1），MA以及历史模拟法的VaR比较](http://mp.weixin.qq.com/s?__biz=MzU4NTA1MDk4MA==&mid=2247491996&idx=2&sn=8b84ed0248f273f89953c4030cb872a8&chksm=fd92df97cae55681f49ee24ccd88a6beca4b393466384c24c32bd83be62ca03a98a5f1441cc6&scene=21#wechat_redirect)