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

在这篇文章中，我们将学习一种在价格序列中建立波动性模型的标准方法，即广义自回归条件异方差（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)