Rで学ぶGARCHモデル
Kris Boudt
Professor of finance and econometrics
names(coef(flexgarchfit))"mu" "ar1" "omega" "alpha1" "beta1" "gamma1" "skew" "shape"
ar1パラメータが0であれば、定数平均モデルを使用できます。gamma1パラメータが0であれば、GARCHインミーン効果はなく、GJRの代わりに標準GARCHモデルを使用できます。ar1とgamma1のパラメータは0ですか?round(coef(flexgarchfit), 6)
mu ar1 omega alpha1 beta1 gamma1 skew shape
-0.000021 0.000150 0.000000 0.034281 0.968688 -0.010093 1.013487 9.139252
方法:推定パラメータを標準誤差と比較します。
標準誤差はパラメータ推定値の標準偏差です。
t統計量 = 推定パラメータ / 標準誤差
モデルの特定と推定
flexgarchspec <- ugarchspec(mean.model = list(armaOrder = c(1, 0)),
variance.model = list(model = "gjrGARCH"),
distribution.model = "sstd")
flexgarchfit <- ugarchfit(data = msftret, spec = flexgarchspec)
パラメータ推定値、標準誤差、t統計量の表を以下で出力します:
round(flexgarchfit@fit$matcoef, 6)
round(flexgarchfit@fit$matcoef, 6)
Estimate Std. Error t value Pr(>|t|)
mu 0.000610 0.000189 3.220843 0.001278
ar1 -0.037799 0.013718 -2.755532 0.005860
omega 0.000002 0.000001 2.617696 0.008853
alpha1 0.034574 0.003395 10.182558 0.000000
beta1 0.935927 0.007163 130.667531 0.000000
gamma1 0.055483 0.009772 5.677857 0.000000
skew 1.059959 0.020676 51.264435 0.000000
shape 4.392327 0.256700 17.110745 0.000000
注: 3.220843 = 0.000610 / 0.000189、-2.755532 = -0.037799 / 0.013718、...
t統計量 = 推定パラメータ / 標準誤差
判断の目安:
もし |t統計量| > 2
ならば 推定パラメータは統計的に有意
したがって 真のパラメータ $\neq$ 0 と結論します。
すべてのt統計量が2より大きい:すべてのパラメータ推定値は統計的に有意です。
Estimate Std. Error t value Pr(>|t|)
mu 0.000610 0.000189 3.220843 0.001278
ar1 -0.037799 0.013718 -2.755532 0.005860
omega 0.000002 0.000001 2.617696 0.008853
alpha1 0.034574 0.003395 10.182558 0.000000
beta1 0.935927 0.007163 130.667531 0.000000
gamma1 0.055483 0.009772 5.677857 0.000000
skew 1.059959 0.020676 51.264435 0.000000
shape 4.392327 0.256700 17.110745 0.0000000
最終列のp値を用いても同じ結論が得られます。
判断の目安:
もし p値 < 5%
ならば 推定パラメータは統計的に有意
したがって 真のパラメータ $\neq$ 0 と結論します。
Rで学ぶGARCHモデル