Các mô hình GARCH trong R
Kris Boudt
Professor of finance and econometrics
Đổi đối số distribution.model của ugarchspec() từ "norm" sang "sstd":
garchspec <- ugarchspec( mean.model = list(armaOrder = c(0, 0)),
variance.model = list(model = "sGARCH"),
distribution.model = "norm")
$$ \downarrow$$
garchspec <- ugarchspec( mean.model = list(armaOrder = c(0, 0)),
variance.model = list(model = "sGARCH"),
distribution.model = "sstd")
Với các giả định mô hình
$$ R_{t} = \mu_{t} + e_{t} $$ $$ e_{t} \sim N(0, \sigma^{2}_{t}) $$
suy ra rằng
$$ \frac{R_{t} - \mu_{t}}{\sigma_{t}} \sim N(0, 1) $$
$$ Z_{t} = \frac{R_{t} - \hat{\mu_{t}}}{ \hat{\sigma_{t}}} $$
Tính trong R
# Obtain standardized returns
stdret <- residuals(garchfit, standardize = TRUE)
library(PerformanceAnalytics)
chart.Histogram(sp500ret, methods = c("add.normal", "add.density"),
colorset = c("gray", "red", "blue"))


Vì vậy, một phân phối thực tế cần bao quát:
Trong rugarch, có thể dùng phân phối skewed student t:
garchspec <- ugarchspec(distribution.model = "sstd")
So với phân phối chuẩn, phân phối skewed student t có hai tham số bổ sung:
shape): $\nu$ càng thấp thì đuôi càng dày.skew): khi $\xi=1$: đối xứng. Khi $\xi<1$: lệch âm. Khi $\xi>1$: lệch dương.Trường hợp đặc biệt:




Đặt đối số distribution.model thành "sstd"
garchspec <- ugarchspec(mean.model = list(armaOrder = c(0,0)),
variance.model = list(model = "sGARCH"),
distribution.model = "sstd")
Ước lượng mô hình
garchfit <- ugarchfit(data = sp500ret, spec = garchspec)
coef(garchfit)
mu omega alpha1 beta1 skew shape
5.669200e-04 6.281258e-07 7.462984e-02 9.223701e-01 9.436331e-01 6.318621e+00
Các mô hình GARCH trong R