拓端tecdat|R語言時變波動率和ARCH,GARCH,GARCH-in-mean模型分析股市收益率時間序列
?原文鏈接:http://tecdat.cn/?p=23225?
原文出處:拓端數(shù)據(jù)部落公眾號
自回歸條件異方差(ARCH)模型涉及具有時變異方差的時間序列,其中方差是以特定時間點(diǎn)的現(xiàn)有信息為條件的。
ARCH模型
ARCH模型假設(shè)時間序列模型中誤差項的條件均值是常數(shù)(零),與我們迄今為止討論的非平穩(wěn)序列不同),但其條件方差不是。這樣一個模型可以用公式1、2和3來描述。

方程4和5給出了測試模型和假設(shè),以測試時間序列中的ARCH效應(yīng),其中殘差e^t來自于將變量yt回歸一個常數(shù),如1,或回歸一個常數(shù)加上其他回歸因子;方程4中的測試可能包括幾個滯后項,在這種情況下,無效假設(shè)(方程5)是所有這些項都不顯著。

無效假設(shè)是不存在ARCH效應(yīng)。檢驗統(tǒng)計量為
?下面的例子使用了數(shù)據(jù)集,它包含了500個股票收益率的生成觀測值。圖顯示了數(shù)據(jù)的時間序列圖和柱狀圖。
plot.ts(r)
hist(r)
圖: 變量 的水平和柱狀圖
讓我們首先對數(shù)據(jù)集中的變量r一步一步地進(jìn)行公式4和5中描述的ARCH檢驗。
summary(yd)
ehsq <- ts(resid(mean)^2)
summary(ARCH)

Rsq <- glance(ARCH)[[1]]
LM <- (T-q)*Rsq
Chicr <- qchisq(1-alpha, q)

結(jié)果是LM統(tǒng)計量,等于62.16,與α=0.05和q=1自由度的臨界卡方值進(jìn)行比較;這個值是χ2(0.95,1)=3.84;這表明拒絕了無效假設(shè),結(jié)論是該序列具有ARCH效應(yīng)。
如果我們不使用一步步的程序,而是使用R的ARCH檢驗功能之一,也可以得出同樣的結(jié)論。
ArchTest

函數(shù)garch(),當(dāng)使用order=參數(shù)等于c(0,1)時,成為一個ARCH模型。這個函數(shù)可以用來估計和繪制方程3中定義的方差ht,如以下代碼和圖所示。?
garch(r,c(0,1))

summary(arch)

ts(2*fitted.values^2)
plot.ts(hhat)

圖 對數(shù)據(jù)集的ARCH(1)方差的估計
GARCH模型
# 使用軟件包`garch`來建立GARCH模型
fit(spec=garch, data=r)
coef(Fit)

fitted.values
fit$sigma^2)
plot.ts(hhat)

圖: 使用數(shù)據(jù)集的標(biāo)準(zhǔn)GARCH模型(sGARCH)。
# tGARCH
garchfit(spec, data=r, submodel="TGARCH")
coef(garchfit)

fitted.values
fit$sigma^2)
plot.ts(hhat)

圖: 數(shù)據(jù)集的tGARCH模型
# GARCH-IN-MEAN模型
fit( data=r,
distribution="std",variance=list(model="fGARCH")
coef(garchFit)

fit$fitted.values
fit$sigma^2)
plot.ts(hhat)

圖:使用數(shù)據(jù)集的GARCH-in-mean模型的一個版本
圖顯示了GARCH模型的幾個版本。預(yù)測結(jié)果可以通過ugarchboot()來獲得。

最受歡迎的見解
1.在python中使用lstm和pytorch進(jìn)行時間序列預(yù)測
2.python中利用長短期記憶模型lstm進(jìn)行時間序列預(yù)測分析
3.使用r語言進(jìn)行時間序列(arima,指數(shù)平滑)分析
4.r語言多元copula-garch-模型時間序列預(yù)測
5.r語言copulas和金融時間序列案例
6.使用r語言隨機(jī)波動模型sv處理時間序列中的隨機(jī)波動
7.r語言時間序列tar閾值自回歸模型
8.r語言k-shape時間序列聚類方法對股票價格時間序列聚類
9.python3用arima模型進(jìn)行時間序列預(yù)測