在这篇文章中，我们将学习一种在价格序列中建立波动性模型的标准方法，即广义自回归条件异方差（GARCH）模型。
价格波动的 GARCH 模型的思想是利用误差结构的近期实现来预测误差结构的未来实现。更简单地说，我们经常看到在高波动性或低波动性时期的聚类，因此我们可以利用近期的波动性来预测近期未来的波动性。
我们将使用SPY价格来说明波动率的模型。下面的图显示了SPY收益率。
colnames(SPYRet) <- c('SPY')
SPYRet <- tidy(SPYRet)
ggplot(SPYRet) +
geom_line() +
labs(title = "SPY收益率"
这篇文章的首要目标是对价格收益序列的变化的平均值和方差进行建模。
可下载资源
作者
接下来，我们将通过两种常用的方式来显示收益率的变化方差。这些是绘制价格收益的绝对值。
或价格收益的平方
这两种情况都是有意义的，因为方差总是一个正数，并且受到偏离平均值的影响。当然这是在我们知道收益序列的平均值为0的情况下。
在价格数据中，百分比回报率的平均值几乎都非常接近于0。 如果平均回报率非零，那么我们可以直接绘制
add_column( value^2, abs(value))
ggplot(SPYRet) +
geom_line() +
labs(title = "SPY绝对收益值"
ggplot(Ret, aes(y = SquaredReturns) +
geom_line() +
labs(title = "SPY平方收益率"
波动率的GARCH模型
普通的（GARCH模型有很多变体）GARCH模型如下。
第一行是为均值建模的方程。这里没有ARMA效应，但如果你发现它们很重要，可以很容易地把它们放进去。只有一个截距和一个误差项。接下来的三行将更多的结构放在误差项上，
第二行为什么我们要把两个项相乘来得到ϵt？
要看到这一点，重要的是要牢记这里的目标。我们正在寻找一个能给我们带来rSPYt变化方差的模型
随时关注您喜欢的主题
因此，如果基本收益模型是
接下来的步骤依赖于随机变量方差的属性。具体来说，如果
因此，如果我们得到一个εt的模型
考虑GARCH模型中的第二行。
请注意 σtσ2ϵ2
最后一行是由于
估计GARCH模型
下面的代码使用rugarch R包来估计GARCH(p = 1, q = 1)模型。请注意，p和q表示σ2tϵ2t的滞后数。
第一条命令要求它用model = “sGARCH “指定一个普通的GARCH。它要求它使用ARMA(1, 1)作为收益模型，即armaOrder = c(1, 1), include.mean = TRUE。我们要求它使用N(0,1)，distribution.model=”norm”，模型输出显示在 “最优参数 “下。mu、ar1和ma1系数来自均值模型（ARMA（1，1）），omega、alpha1和beta1是来自
garchspec(variance.m = list(model = "sGARCH", mean.m= list(armaOrder = c(1, 1)
dist = "norm")
garchfit
现在让我们用这个估计的模型来产生滚动预测，即
可以看到，这个模型很好的预测了波动率峰值将保持多长时间，或者说是模拟了波动率峰值回落到长期平均水平的路径。
forecast(spec, n.ahead = 1, n.roll = 2499, out = 2500)
plot(forecast)
由于所有的计量经济学模型都是用过去的数值来预测当前的数值，所以它无法预见波动率最初上升的情况。
可下载资源
关于作者
Kaizong Ye是拓端研究室（TRL）的研究员。Kaizong Ye是拓端研究室（TRL）的研究员。在此对他对本文所作的贡献表示诚挚感谢，他在上海财经大学完成了统计学专业的硕士学位，专注人工智能领域。擅长Python.Matlab仿真、视觉处理、神经网络、数据分析。
本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。
非常感谢您阅读本文，如需帮助请联系我们！
【视频讲解】线性时间序列原理及混合ARIMA-LSTM神经网络模型预测股票收盘价研究实例
R语言股价跳跃点识别：隐马尔可夫hmm和GARCH-Jump对sp500金融时间序列分析
R 语言GJR-GARCH、GARCH-t、GARCH-ged分析金融时间序列数据波动性预测、检验、可视化
