使用FGARCH包。
fit=garchFit(formula=~garch(2,1), data=3,cond.dist="QMLE")
结果如下：
Series Initialization:
ARMA Model:                arma
Formula Mean:              ~ arma(0, 0)
GARCH Model:               garch
Formula Variance:          ~ garch(2, 1)
ARMA Order:                0 0
Max ARMA Order:            0
GARCH Order:               2 1
Max GARCH Order:           2
Maximum Order:             2
Conditional Dist:          QMLE
h.start:                   3
llh.start:                 1
Length of Series:          1
Recursion Init:            mci
Series Scale:              NA

错误于lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :
  0 (non-NA) cases

不知为什么会这样。