【時系列】GARCHモデル（R）

今回はRで実装します。 今回はこちらのサイトを写経します。

ドル円のボラティリティをGARCHで推定 – Momentum

データ

まず、データはFRED先生からドル円レートを落としてきて使います。

Japan / U.S. Foreign Exchange Rate | FRED | St. Louis Fed

Rで実装

library(tseries)
fx.ret <- fx/lag(fx)-1
fx.garch <- garch(fx.ret["1973-04::"]) #デフォルトはGARCH(1,1)
summary(fx.garch)

結果

 
Call:
garch(x = fx.ret)

Model:
GARCH(1,1)

Residuals:
     Min       1Q   Median       3Q      Max 
-5.20882 -0.61636 -0.04313  0.55319  5.88943 

Coefficient(s):
    Estimate  Std. Error  t value Pr(>|t|)    
a0 2.120e-07   8.670e-08    2.445   0.0145 *  
a1 5.313e-02   8.385e-03    6.337 2.34e-10 ***
b1 9.421e-01   8.660e-03  108.779  < 2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Diagnostic Tests:
    Jarque Bera Test

data:  Residuals
X-squared = 335.28, df = 2, p-value < 2.2e-16


    Box-Ljung test

data:  Squared.Residuals
X-squared = 4.8338, df = 1, p-value = 0.02791
'''