R語言基于ARMA-GARCH-VaR模型擬合和預(yù)測實證研究分析案例
原文鏈接:http://tecdat.cn/?p=3186
?
本文顯示了如何基于潛在的ARMA-GARCH過程(當(dāng)然也涉及更廣泛意義上的QRM)來擬合和預(yù)測風(fēng)險價值(VaR)。
1?從ARMA-GARCH進程模擬(log-return)數(shù)據(jù)
我們考慮使用\(t \)分布的ARMA(1,1)-GARCH(1,1)過程。
模擬一條路徑(用于說明目的)。
?
nu <- 3
fixed.p <- list(mu = 0, # ? mu (intercept)
ar1 = 0.5, # ? phi_1 (AR(1) 參數(shù) of mu_t)
ma1 = 0.3, # ? theta_1 (MA(1) 參數(shù) of mu_t)
omega = 4, # ? alpha_0 (intercept)
alpha1 = 0.4, # ? alpha_1 (GARCH(1) 參數(shù) of sigma_t^2)
beta1 = 0.2, # ? beta_1 (GARCH(1) 參數(shù) of sigma_t^2)
shape = nu) #
armaOrder <- c(1,1) # ARMA 參數(shù)
garchOrder <- c(1,1) # GARCH 參數(shù)
varModel <- list(model = "sGARCH", garchOrder = garchOrder)
spec <- ugarchspec(varModel, mean.model = list(armaOrder = armaOrder),
fixed.pars = fixed.p, distribution.model = "std") # t 標準殘差
?
作為一個完整性檢查,讓我們繪制模擬路徑,條件標準偏差和殘差。
?
plot(X, ? type = "l", xlab = "t", ylab = expression(X[t]))
plot(sig, type = "h", xlab = "t", ylab = expression(sigma[t]))

plot(eps, type = "l", xlab = "t", ylab = expression(epsilon[t]))
?

2將ARMA-GARCH模型擬合到(模擬)數(shù)據(jù)
擬合ARMA-GARCH流程 。
讓我們再考慮一些健全性檢查。
## 擬合 an ARMA(1,1)-GARCH(1,1) model
spec <- ugarchspec(varModel, mean.model = list(armaOrder = armaOrder),
distribution.model = "std") #
fit <- ugarchfit(spec, data = X) # fit
##
mu. <- fitted(fit) # 擬合 hat{mu}_t (= hat{X}_t)
sig. <- sigma(fit) # 擬合 hat{sigma}_t
##
stopifnot(all.equal(as.numeric(mu.), ?fit@fit$fitted.values),
all.equal(as.numeric(sig.), fit@fit$sigma))



3計算VaR時間序列
計算VaR估計值。請注意,我們也可以在這里使用基于GPD的估算器。
4 Backtest VaR估計值
讓我們回顧一下VaR的估計。
## [1] 10
## [1] 12
## [1] "Correct Exceedances"
## [1] "Fail to Reject H0"
## [1] "Correct Exceedances & Independent"
## [1] "Fail to Reject H0"
5基于擬合模型預(yù)測VaR
現(xiàn)在預(yù)測VaR。
6模擬\((X_t)\)的未來軌跡并計算相應(yīng)的VaR
模擬路徑,估計每個模擬路徑的VaR(注意quantile()
這里不能使用,因此我們必須手動構(gòu)建VaR)并計算\(\ mathrm {VaR} _ \ alpha \)的bootstrap置信區(qū)間。
7
最后,讓我們顯示所有結(jié)果。

?
非常感謝您閱讀本文,有任何問題請在下面留言!

最受歡迎的見解
1.R語言基于ARMA-GARCH-VaR模型擬合和預(yù)測實證研究
2.R語言時變參數(shù)VAR隨機模型
3.R語言時變參數(shù)VAR隨機模型
4.R語言基于ARMA-GARCH過程的VAR擬合和預(yù)測
5.GARCH(1,1),MA以及歷史模擬法的VaR比較
6.R語言時變參數(shù)VAR隨機模型
7.R語言實現(xiàn)向量自動回歸VAR模型
8.R語言隨機搜索變量選擇SSVS估計貝葉斯向量自回歸(BVAR)模型
9.R語言VAR模型的不同類型的脈沖響應(yīng)分析