在 R 中构建 GARCH 模型
Kris Boudt
Professor of finance and econometrics
Robert Engle

Tim Bollerslev






关于均值的 AR(MA) 模型,参见 Datacamp 的时间序列分析课程。



为使 GARCH 过程更贴近现实,需要:
$\omega$、$\alpha$ 和 $\beta$ 均 $>0$:确保任意时点 $\sigma^2_t >0$。
$\alpha + \beta < 1$:确保预测方差 $\sigma^2_t$ 回归长期方差:
$$ \sigma^{2}_{t} = \omega + \alpha e^{2}_{t-1} \beta \sigma^{2}_{t-1} $$
# Set parameter values
alpha <- 0.1
beta <- 0.8
omega <- var(sp500ret) * (1 - alpha - beta)
# Then: var(sp500ret) = omega / (1 - alpha - beta)
# Set series of prediction error
e <- sp500ret - mean(sp500ret) # Constant mean
e2 <- e ^ 2
# We predict for each observation its variance.
nobs <- length(sp500ret)
predvar <- rep(NA, nobs)
# Initialize the process at the sample variance
predvar[1] <- var(sp500ret)
# Loop starting at 2 because of the lagged predictor
for (t in 2:nobs){
predvar[t] <- omega + alpha * e2[t - 1] + beta * predvar[t-1]
}
# Volatility is sqrt of predicted variance
predvol <- sqrt(predvar)
predvol <- xts(predvol, order.by = time(sp500ret))
# We compare with the unconditional volatility
uncvol <- sqrt(omega / (1 - alpha-beta))
uncvol <- xts(rep(uncvol, nobs), order.by = time(sp500ret))
# Plot
plot(predvol)
lines(uncvol, col = "red", lwd = 2)

在 R 中构建 GARCH 模型