Model GARCH di R
Kris Boudt
Professor of finance and econometrics
Berdasarkan model GARCH terestimasi, kita memiliki:

Implementasi
e <- residuals(tgarchfit)
mean(e ^ 2)
Model GARCH menghasilkan:

Implementasi
e <- residuals(tgarchfit)
d <- e ^ 2 - sigma(tgarchfit) ^ 2
mean(d ^ 2)
tgarchspec <- ugarchspec(mean.model = list(armaOrder = c(0, 0)),
variance.model = list(model = "sGARCH", variance.targeting = TRUE),
distribution.model = "std")
tgarchfit <- ugarchfit(data = EURUSDret, spec = tgarchspec)
# Hitung galat kuadrat rata-rata untuk mean
e <- residuals(tgarchfit) ^ 2
mean(e ^ 2) # 3.836205e-05
# Hitung galat kuadrat rata-rata untuk varians
d <- e ^ 2 - sigma(tgarchfit) ^ 2
mean(d ^ 2) # 5.662366e-09
likelihood(tgarchfit) # returns 18528.58
Analisis dengan membandingkan model lain:
# Model kompleks dengan banyak parameter
flexgarchspec <- ugarchspec(mean.model = list(armaOrder = c(1, 0)),
variance.model = list(model = "gjrGARCH"),
distribution.model = "sstd")
flexgarchfit <- ugarchfit(data = EURUSDret, spec = flexgarchspec)
likelihood(flexgarchfit) # returns 18530.49
Perhatian: Kita menggunakan evaluasi in-sample, di mana sampel estimasi dan evaluasi bertepatan.
Risiko overfitting:
Overfitting adalah memilih model terlalu kompleks yang cocok untuk return pada sampel estimasi, tetapi tidak untuk return di luar sampel.
kriteria informasi = - likelihood + penalti(jumlah parameter)
Aturan praktis:
Pilih model dengan kriteria informasi terendah.
Metode infocriteria() menampilkan kriteria informasi untuk berbagai penalti
infocriteria(tgarchfit)
out
Akaike -7.468081
Bayes -7.462833
Shibata -7.468083
Hannan-Quinn -7.466241
`
Interpretasi perlu membandingkan dengan kriteria informasi dari model lain.
tgarchspec <- ugarchspec(mean.model = list(armaOrder = c(0, 0)),
variance.model = list(model = "sGARCH", variance.targeting = TRUE),
distribution.model = "std")
tgarchfit <- ugarchfit(data = EURUSDret, spec = tgarchspec)
length(coef(tgarchfit)) # only 5 parameters
likelihood(tgarchfit) # equals 18528.58
flexgarchspec <- ugarchspec(mean.model = list(armaOrder = c(1, 0)),
variance.model = list(model = "gjrGARCH"), distribution.model = "sstd")
flexgarchfit <- ugarchfit(data = EURUSDret, spec = flexgarchspec)
length(coef(flexgarchfit)) # we now have 8 parameters
likelihood(flexgarchfit) # 18530.49: likelihood increased
Likelihood lebih tinggi lebih baik. Kriteria informasi lebih rendah lebih baik.
infocriteria(tgarchfit) # Model sederhana
Akaike -7.468435
Bayes -7.464499
infocriteria(flexgarchfit) # Model kompleks
Akaike -7.467239
Bayes -7.456742
Model sederhana memiliki kriteria informasi terendah dan sebaiknya dipilih di sini.
tgarchfit <- ugarchfit(data = msftret, spec = tgarchspec)
flexgarchfit <- ugarchfit(data = msftret, spec = flexgarchspec)
infocriteria(tgarchfit)
Akaike -5.481895
Bayes -5.477833
infocriteria(flexgarchfit)
Akaike -5.489087
Bayes -5.478255
Model kompleks memiliki kriteria informasi terendah dan sebaiknya dipilih di sini.
Model GARCH di R