原文鏈接：http://tecdat.cn/?p=23934

最近我們被客戶要求撰寫關(guān)于ARIMA-GARCH的研究報(bào)告，包括一些圖形和統(tǒng)計(jì)輸出。

在本文中，我們將嘗試為蘋果公司的日收益率尋找一個(gè)合適的 GARCH 模型

波動(dòng)率建模需要兩個(gè)主要步驟。

• 指定一個(gè)均值方程（例如 ARMA，AR，MA，ARIMA 等）。

• 建立一個(gè)波動(dòng)率方程（例如 GARCH, ARCH，這些方程是由 Robert Engle 首先開發(fā)的）。

要做(1)，你需要利用著名的Box-Jenkins方法，它包括三個(gè)主要步驟。

• 識(shí)別

• 估算

• 診斷檢查

這三個(gè)步驟有時(shí)會(huì)有不同的名稱，這取決于你讀的是誰(shuí)的書。在本文中，我將更多地關(guān)注（2）。

我將使用一個(gè)名為quantmod的軟件包，它代表量化金融建模框架。這允許你在R中直接從各種在線資源中抓取金融數(shù)據(jù)。

#install.packages("quantmod")?-需要先安裝該軟件包getSymbols(Symbols?=?"AAPL", ???????????src="yahoo", #其他來源包括：谷歌、FRED等。

收益通常有一個(gè)非常簡(jiǎn)單的平均數(shù)方程，這導(dǎo)致了簡(jiǎn)單的殘差。

我們首先要測(cè)試序列依賴性，這是條件異方差的一個(gè)指標(biāo)（序列依賴性與序列相關(guān)不同）。這是通過對(duì)原始序列的平方/絕對(duì)值進(jìn)行測(cè)試，并使用Ljung和Box（1978）的Ljung-Box測(cè)試等聯(lián)合假設(shè)進(jìn)行測(cè)試，這是一個(gè)Portmentau檢驗(yàn)，正式檢驗(yàn)連續(xù)自相關(guān)，直到預(yù)定的滯后數(shù)，如下所示。

其中T是總的周期數(shù)，m是你要測(cè)試的序列相關(guān)的滯后期數(shù)，ρ2k是滯后期k的相關(guān)性，Q?(m)～χ2α有m個(gè)自由度。

檢查

下面是AAPL對(duì)數(shù)收益時(shí)間序列及其ACF，這里我們要尋找顯著的滯后期（也可以運(yùn)行pacf）或存在序列自相關(guān)。

通過觀察ACF，水平序列（對(duì)數(shù)收益）并不是真正的自相關(guān)，但現(xiàn)在讓我們看一下平方序列來檢查序列依賴性。

點(diǎn)擊標(biāo)題查閱往期相關(guān)內(nèi)容

R語(yǔ)言用多元ARMA,GARCH ,EWMA, ETS,隨機(jī)波動(dòng)率SV模型對(duì)金融時(shí)間序列數(shù)據(jù)建模

左右滑動(dòng)查看更多

01

02

03

04

我們可以看到，平方序列的ACF顯示出顯著的滯后。這是一個(gè)信號(hào)，說明我們應(yīng)該在某個(gè)時(shí)候測(cè)試ARCH效應(yīng)。

平穩(wěn)性

我們可以看到，AAPL的對(duì)數(shù)回報(bào)在某種程度上是一個(gè)平穩(wěn)的過程，所以我們將使用Augmented Dicky-Fuller檢驗(yàn)（ADF）來正式檢驗(yàn)平穩(wěn)性。ADF是一個(gè)廣泛使用的單位根檢驗(yàn)，即平穩(wěn)性。我們將使用12個(gè)滯后期，因?yàn)楦鶕?jù)文獻(xiàn)的建議，我們有每日數(shù)據(jù)。何：存在單位根（系列是非平穩(wěn)的

##?##?Title:##??Augmented?Dickey-Fuller?Test##?##?Test?Results:##???PARAMETER:##?????Lag?Order:?12##???STATISTIC:##?????Dickey-Fuller:?-14.6203##???P?VALUE:##?????0.01?##?##?Description:##??Mon?May?25?16:45:37?2020?by?user:?Florian

上面的P值為0.01，表明我們應(yīng)該拒絕Ho，因此，該系列是平穩(wěn)的。

結(jié)構(gòu)突變_檢驗(yàn)_

請(qǐng)注意，我從2008年底開始研究APPL序列。以避免08年大衰退，通常會(huì)在數(shù)據(jù)中產(chǎn)生結(jié)構(gòu)性突變（即趨勢(shì)的嚴(yán)重下降/跳躍）。我們將對(duì)結(jié)構(gòu)性突變/變化進(jìn)行Chow測(cè)試。AAPL的日收益率沒有結(jié)構(gòu)性突變

該圖顯示，用于估計(jì)斷點(diǎn)（BP）數(shù)量的BIC（黑線）是BIC線的最小值，所以我們可以確認(rèn)沒有結(jié)構(gòu)性斷點(diǎn)，因?yàn)樽钚≈凳橇?，即零斷點(diǎn)。在預(yù)測(cè)時(shí)間序列時(shí)，斷點(diǎn)非常重要。

估計(jì)

在這一節(jié)中，我們?cè)噲D用auto.arima命令來擬合最佳a(bǔ)rima模型，允許一個(gè)季節(jié)性差異和一個(gè)水平差異。

正如我們所知，{Yt}的一般ARIMA(p,d,q)。

根據(jù)auto.arima，最佳模型是ARIMA(3,0,2)，平均數(shù)為非零，AIC為-14781.55。我們的平均方程如下（括號(hào)內(nèi)為SE）。

Auto.arima函數(shù)挑選出具有最低AIC的ARIMA(p,d,q)，其中。

其中Λθ是觀察到的數(shù)據(jù)在參數(shù)的mle的概率。因此，如果Auto.arima函數(shù)運(yùn)行N模型，其決策規(guī)則為AIC?=min{AICi}Ni=1

診斷檢查

我們可以看到，我們的ARIMA(3,0,2)的殘差是良好的表現(xiàn)。它們似乎也有一定的正態(tài)分布

##?##??Ljung-Box?test##?##?data:??Residuals?from?ARIMA(3,0,2)?with?non-zero?mean##?Q*?=?6.7928,?df?=?4,?p-value?=?0.1473##?##?Model?df:?6.???Total?lags?used:?10

現(xiàn)在我們將通過對(duì)我們的ARIMA(3,0,2)模型的平方殘差應(yīng)用Ljung-Box測(cè)試來檢驗(yàn)ARCH效應(yīng)。

##?##??Box-Ljung?test##?##?data:??resid^2##?X-squared?=?126.6,?df?=?12,?p-value?<?2.2e-16

我們可以看到，殘差平方的 ACF 顯示出許多顯著的滯后期，因此我們得出結(jié)論，確實(shí)存在 ARCH 效應(yīng)，我們應(yīng)該對(duì)波動(dòng)率進(jìn)行建模。

使用 GARCH 建立波動(dòng)率模型

上面將我們的平均數(shù)方程中的殘差進(jìn)行了平方，看看大的沖擊是否緊隨在其他大的沖擊之后（無(wú)論哪個(gè)方向，即負(fù)的或正的），如果是這樣，那么我們就有條件異方差，意味著我們有需要建模的非恒定方差。下面是一個(gè)GARCH(m,s)的樣子。

其中{?2t}mt=1是我們通常的特異性沖擊，iid隨機(jī)變量，即?2t～WN(0,σ2?)。我們可以更緊湊地寫成：

其中B是標(biāo)準(zhǔn)的后移算子Bi?2t=?2t-i，Biσ2t=σ2t-i。對(duì)于任何整數(shù)ii，以及α和β分別是度數(shù)為m和s的多項(xiàng)式

請(qǐng)注意，一個(gè)特殊情況是當(dāng)s=0時(shí)，GARCH(m,0)被稱為ARCH(m)。

當(dāng)我說GARCH家族時(shí)，它表明模型有變化。

• SGARCH。普通GARCH

• EGARCH。指數(shù)GARCH，允許波動(dòng)率不為負(fù)值（這迫使模型只輸出正方差

• FGARCH。這是為長(zhǎng)記憶模型準(zhǔn)備的。它使用了被稱為 ARFIMA 的 Fractionaly integrated ARIMA（即非整數(shù)整合）。

• GARCH-M：這是GARCH的均值，適合你的均值方程中有波動(dòng)率例如CAPM的方程中有σ。

• GJR-GARCH。假設(shè)負(fù)面沖擊和正面沖擊之間存在不對(duì)稱性（金融數(shù)據(jù)幾乎都是這樣）。

為收益率序列建立波動(dòng)率模型包括四個(gè)步驟：

1. 通過測(cè)試數(shù)據(jù)中的序列依賴性來指定一個(gè)均值方程，如果有必要，為收益序列建立一個(gè) 計(jì)量經(jīng)濟(jì)學(xué)模型（例如，ARIMA 模型）來消除任何線性依賴。

2. 使用平均值方程的殘差來測(cè)試ARCH效應(yīng)。

3. 如果ARCH效應(yīng)在統(tǒng)計(jì)上是顯著的，就指定一個(gè)波動(dòng)率模型，并對(duì)均值和波動(dòng)率方程進(jìn)行聯(lián)合估計(jì)。

4. 仔細(xì)檢查擬合的模型，必要時(shí)對(duì)其進(jìn)行改進(jìn)。

一個(gè)簡(jiǎn)單的 GARCH 模型有以下成分。

均值:?

波動(dòng)率方程:??

誤差假設(shè):?

#以下命令將計(jì)算GARCH(m,s)。請(qǐng)記住，對(duì)于某些m和s的組合，它可能不會(huì)收斂。garchlist(model="sGARCH", #其他選項(xiàng)有egarch, fgarch等。?????????????????????????????????????????????????????garchOrder=c(1,2)),?#你可以在這里修改GARCH(m,s)的階數(shù)?????????????????????????????? mean.model ?, #指定你的ARMA模型，暗示你的模型應(yīng)該是平穩(wěn)的。???????????????????????????????distribution.model??????????#其他分布是?"std?"代表t分布，"ged?"代表一般誤差分布

我們的波動(dòng)率方程由GARCH(1,2)給出，AIC:-5.5277（注意GARCH可能無(wú)法收斂）。

下面是使用我們的波動(dòng)率模型對(duì)波動(dòng)率進(jìn)行的預(yù)測(cè)。這看起來是一個(gè)合理的波動(dòng)率預(yù)測(cè)，但是你想改進(jìn)你的模型。

現(xiàn)在讓我們使用rugarch的標(biāo)準(zhǔn)功能，使用估計(jì)的GARCH(1,2)模型來產(chǎn)生σt的滾動(dòng)預(yù)測(cè)，并將它們與|rt|作對(duì)比。

最后，我們可以手動(dòng)編寫代碼來查看隨時(shí)間變化的波動(dòng)率和對(duì)數(shù)收益率rt，如下圖。

#?這將有助于在對(duì)數(shù)收益率上繪制sigma隨時(shí)間變化的圖。sigma.t?#這是你的波動(dòng)率序列ggplot() ??geom_line(aes(x=as.numeric( ?????????????? ??theme_bw()+

結(jié)論

事實(shí)證明，GARCH系列是所謂確定性波動(dòng)率模型的一部分。還有一個(gè)家族叫做隨機(jī)波動(dòng)率模型，它允許模型中存在隨機(jī)性，而GARCH假設(shè)我們對(duì)波動(dòng)率進(jìn)行了完美的建模（如果你對(duì)你所分析的序列非常熟悉，這可能是一個(gè)好的假設(shè)，但實(shí)際情況并不總是這樣）。隨機(jī)波動(dòng)率模型通常是用馬爾科夫鏈蒙特卡洛（MCMC）和準(zhǔn)蒙特卡洛方法來估計(jì)的，如果你學(xué)過隨機(jī)過程的相關(guān)內(nèi)容，你會(huì)知道這是什么。

參考文獻(xiàn)

• Tsay, R. (2010).?Analysis of Financial Time Series. (3rd ed., Wiley Series in Probability and Statistics).

• Brockwell, P., & Davis, Richard A. (2016).?Introduction to time series and forecasting?(3rd ed., Springer texts in statistics). New York: Springer.

• Racine, Jeffrey S. (2019)?Reproducible Econometrics Using R?(Oxford)

點(diǎn)擊文末?“閱讀原文”

獲取全文完整資料。

本文選自《R語(yǔ)言ARIMA-GARCH波動(dòng)率模型預(yù)測(cè)股票市場(chǎng)蘋果公司日收益率時(shí)間序列》。

點(diǎn)擊標(biāo)題查閱往期內(nèi)容

R語(yǔ)言風(fēng)險(xiǎn)價(jià)值：ARIMA，GARCH，Delta-normal法滾動(dòng)估計(jì)VaR（Value at Risk）和回測(cè)分析股票數(shù)據(jù)
R語(yǔ)言GARCH模型對(duì)股市sp500收益率bootstrap、滾動(dòng)估計(jì)預(yù)測(cè)VaR、擬合診斷和蒙特卡羅模擬可視化
R語(yǔ)言單變量和多變量（多元）動(dòng)態(tài)條件相關(guān)系數(shù)DCC-GARCH模型分析股票收益率金融時(shí)間序列數(shù)據(jù)波動(dòng)率
R語(yǔ)言中的時(shí)間序列分析模型：ARIMA-ARCH / GARCH模型分析股票價(jià)格
GARCH-DCC模型和DCC（MVT）建模估計(jì)
R語(yǔ)言預(yù)測(cè)期貨波動(dòng)率的實(shí)現(xiàn)：ARCH與HAR-RV與GARCH，ARFIMA模型比較
ARIMA、GARCH 和 VAR模型估計(jì)、預(yù)測(cè)ts 和 xts格式時(shí)間序列
PYTHON用GARCH、離散隨機(jī)波動(dòng)率模型DSV模擬估計(jì)股票收益時(shí)間序列與蒙特卡洛可視化
極值理論 EVT、POT超閾值、GARCH 模型分析股票指數(shù)VaR、條件CVaR：多元化投資組合預(yù)測(cè)風(fēng)險(xiǎn)測(cè)度分析
Garch波動(dòng)率預(yù)測(cè)的區(qū)制轉(zhuǎn)移交易策略
金融時(shí)間序列模型ARIMA 和GARCH 在股票市場(chǎng)預(yù)測(cè)應(yīng)用
時(shí)間序列分析模型：ARIMA-ARCH / GARCH模型分析股票價(jià)格
R語(yǔ)言風(fēng)險(xiǎn)價(jià)值：ARIMA，GARCH，Delta-normal法滾動(dòng)估計(jì)VaR（Value at Risk）和回測(cè)分析股票數(shù)據(jù)
R語(yǔ)言GARCH建模常用軟件包比較、擬合標(biāo)準(zhǔn)普爾SP 500指數(shù)波動(dòng)率時(shí)間序列和預(yù)測(cè)可視化
Python金融時(shí)間序列模型ARIMA 和GARCH 在股票市場(chǎng)預(yù)測(cè)應(yīng)用
MATLAB用GARCH模型對(duì)股票市場(chǎng)收益率時(shí)間序列波動(dòng)的擬合與預(yù)測(cè)
R語(yǔ)言極值理論 EVT、POT超閾值、GARCH 模型分析股票指數(shù)VaR、條件CVaR：多元化投資組合預(yù)測(cè)風(fēng)險(xiǎn)測(cè)度分析
Python 用ARIMA、GARCH模型預(yù)測(cè)分析股票市場(chǎng)收益率時(shí)間序列
R語(yǔ)言中的時(shí)間序列分析模型：ARIMA-ARCH / GARCH模型分析股票價(jià)格
R語(yǔ)言ARIMA-GARCH波動(dòng)率模型預(yù)測(cè)股票市場(chǎng)蘋果公司日收益率時(shí)間序列
Python使用GARCH，EGARCH，GJR-GARCH模型和蒙特卡洛模擬進(jìn)行股價(jià)預(yù)測(cè)
R語(yǔ)言時(shí)間序列GARCH模型分析股市波動(dòng)率
R語(yǔ)言ARMA-EGARCH模型、集成預(yù)測(cè)算法對(duì)SPX實(shí)際波動(dòng)率進(jìn)行預(yù)測(cè)
matlab實(shí)現(xiàn)MCMC的馬爾可夫轉(zhuǎn)換ARMA - GARCH模型估計(jì)
Python使用GARCH，EGARCH，GJR-GARCH模型和蒙特卡洛模擬進(jìn)行股價(jià)預(yù)測(cè)
使用R語(yǔ)言對(duì)S＆P500股票指數(shù)進(jìn)行ARIMA + GARCH交易策略
R語(yǔ)言用多元ARMA,GARCH ,EWMA, ETS,隨機(jī)波動(dòng)率SV模型對(duì)金融時(shí)間序列數(shù)據(jù)建模
R語(yǔ)言股票市場(chǎng)指數(shù)：ARMA-GARCH模型和對(duì)數(shù)收益率數(shù)據(jù)探索性分析
R語(yǔ)言多元Copula GARCH 模型時(shí)間序列預(yù)測(cè)
R語(yǔ)言使用多元AR-GARCH模型衡量市場(chǎng)風(fēng)險(xiǎn)
R語(yǔ)言中的時(shí)間序列分析模型：ARIMA-ARCH / GARCH模型分析股票價(jià)格
R語(yǔ)言用Garch模型和回歸模型對(duì)股票價(jià)格分析
GARCH（1,1），MA以及歷史模擬法的VaR比較
matlab估計(jì)arma garch 條件均值和方差模型