ARCH模型类

ARCH模型类

1 ARCH模型

异方差问题一般出现在带有截面的观测数据中，但对于某些时间序列也存在异方差问题。例如股票收益率波动较大的观测一般集聚在一起，而波动较小的观测集聚在一起，这表现出一种“波动集聚”。

1.1 ARCH模型

考虑模型
$$y_t={\mathbf{x}}_t^{\prime }\mathbf{\beta }+{\epsilon }_t\text{\hspace{1em}(1)}$$
记${\epsilon }_t$条件方差为${\sigma }_t^2\equiv \mathrm{Var}\left({\epsilon }_t\mid {\epsilon }_{t-1},\cdots \right)$，这表明扰动项不是同方差，它随着时间变动而变动。模型(1)中，${\epsilon }_t$是引起$y_t$变化的系统扰动的部分，$y_t$的波动程度取决于${\epsilon }_t$的分布情况，将扰动项的方差视为扰动项的函数
$${\sigma }_t^2={\alpha }_0+{\alpha }_1{\epsilon }_{t-1}^2\text{\hspace{1em}(2)}$$
式表明，上一期的扰动项的会对本期扰动项的方差产生影响。称模型(1)为均值方程，模型(2)为方差方程。均值方程的自变量如果视为被解释变量的滞后值，则是标准的ARCH模型，记作
$$\{\begin{array}{l}{x}_t={\beta }_0+{\beta }_1{x}_{t-1}+{\beta }_2{x}_{t-2}+\dots +{\beta }_p{x}_{t-p}+u_t\\ {\sigma }_t^2=\mathrm{E}\left(u_t^2\right)={\alpha }_0+{\alpha }_1{u}_{t-1}^2+{\alpha }_2{u}_{t-2}^2+\dots +{\alpha }_q{u}_{t-q}^2\end{array}\text{\hspace{1em}(3)}$$
式中第一个方程为均值方程，即服从自回归过程$AR(p)$，第二个方程为条件方差方程，${\mu }_t\sim ARCH(q)$ 。由于自回归方程为平稳序列建模方法，因此随机过程 { x t } \{x_t\} {xt​}必须是平稳的，且特征方程满足
$$1-{\beta }_{1}L-{\beta }_2{L}^2-\dots -{\beta }_p{L}^p=0$$
的特征根全部在单位圆外。 { x t } \{x_t\} {xt​}条件期望为
$$\mathbf{E}\left(x_t\mid x_{t-1},\dots ,x_{t-p}\right)={\beta }_0+{\beta }_1{x}_{t-1}+{\beta }_2{x}_{t-2}+\dots +{\beta }_p{x}_{t-p}$$
当时间$t\to T$时，$x_t$的无条件期望为
$$\mathbf{E}\left({\mathbf{x}}_t\right)=\frac{{\beta }_0}{1-{\beta }_1-\cdots -{\beta }_p}$$

---

1.2 ARCH模型性质

• 非负性

ARCH方程${\sigma }_t^2=\mathrm{E}\left(u_t^2\right)={\alpha }_0+{\alpha }_1{u}_{t-1}^2+{\alpha }_2{u}_{t-2}^2+\dots +{\alpha }_q{u}_{t-q}^2$中所有扰动项平方项均非负，因此各个参数${\alpha }_i\ge 0(i=1,2\dots q)$,且${\alpha }_0>0$。当${\alpha }_i=0(i=1,2\dots q)$，则${\sigma }_t^2={\alpha }_0>0$

• 平稳性

为了保证方差序列也是平稳的，需要满足ARCH的特征方程
$$1-{\alpha }_{1}L-{\alpha }_2{L}^2-\dots -{\alpha }_q{L}^q=0$$
的所有特征根全部在单位圆外。

• 归一性

ARCH方程的所有参数${\sum }_i^q{\alpha }_i\in [0,1)(i=1,2\dots q)$,

证明：ARCH方程两端求期望
$$\begin{array}{rl}{\sigma }_t^2& ={\alpha }_0+{\alpha }_1\mathbf{E}\left(u_{t-1}^2\right)+{\alpha }_2\mathbf{E}\left(u_{t-2}^2\right)+\dots +{\alpha }_q\mathbf{E}\left(u_{t-q}^2\right)\\ & ={\alpha }_0+{\alpha }_1{\sigma }_{t-1}^2+{\alpha }_2{\sigma }_{t-2}^2+\dots +{\alpha }_q{\sigma }_{t-q}^2\end{array}$$
当$t\to \infty$时，
$${\sigma }^2={\alpha }_0+{\alpha }_1{\sigma }^2+{\alpha }_2{\sigma }^2+\dots +{\alpha }_q{\sigma }^2$$
解得
$${\sigma }^2=\frac{1}{1-{\sum }_{i=1}^q{\alpha }_i}{\alpha }_0$$
由于${\alpha }_0\ge 0$,故${\sum }_i^q{\alpha }_i\in [0,1)(i=1,2\dots q)$，才能保证${\sigma }_t^2$非负。

---

1.3 ARCH模型估计

ARCH模型用极大似然方法估计。设均值方程为
$$y_t={\mathbf{x}}_t^{\prime }\mathbf{\beta }+{\epsilon }_t$$
其中参数向量$\beta ={\left({\beta }_0{\beta }_1,\dots ,{\beta }_{k-1}\right)}^{\prime }$，数据向量$x_t={\left(1x_1,\dots ,x_{k-1}\right)}^{\prime }$。注意$x_t$即包括$y$的滞后项(AR过程)也可以包括其他外生变量(ADL过程)。扰动项${\epsilon }_t\sim ARCH(q)$。由于ARCH过程存在滞后$q$项，为了保证样本容量为$T$，假设总的样本容量对应的时间观测点为$1,2,3\dots T,\dots T+q$。设方程(1)的扰动项
$${\mathbf{\epsilon }}_t=\sqrt{h_t}{\mathbf{v}}_t$$
其中$v_t$为服从$iid(0,1)$,$Cov(x_t,v_t)=0$,$h_t$为扰动项${\epsilon }_t$的方差,即
$$h_t={\alpha }_0+{\alpha }_1{u}_{t-1}^2+{\alpha }_2{u}_{t-2}^2+\dots +{\alpha }_q{u}_{t-q}^2$$
显然${\sigma }_t^2=\mathbf{E}\left(u_t^2\right)=h_t,\mathbf{E}\left(u_t\right)=0$。如果${\epsilon }_t$服从正态分布，则$y_t$的条件概率密度函数为
$$f\left(y_t\mid x_t,{\alpha }_i,\beta \right)=\frac{1}{\sqrt{2\pi h_t}}\exp \left(-\frac{{\left(y_t-x_t^{\prime }\beta \right)}^2}{2h_t}\right)$$
另外
$$h_t={\alpha }_0+{\alpha }_1{\left(y_{t-1}-x_{t-1}^{\prime }\beta \right)}^2+{\alpha }_2{\left(y_{t-2}-x_{t-2}^{\prime }\beta \right)}^2+\dots +{\alpha }_q{\left(y_{t-q}-x_{t-q}^{\prime }\beta \right)}^2$$
将参数$\beta$有$\alpha ={\left(\begin{array}{lllll}{\alpha }_0& {\alpha }_1& {\alpha }_2& \dots & {\alpha }_q\end{array}\right)}^{\prime }$记作$\gamma =(\beta \alpha {)}^{\prime }$，考虑
$$\{\begin{array}{l}{y}_t=x_t^{\prime }\beta +u_t\\ u_t=\sqrt{h_t}{v}_t\\ h_t={\alpha }_0+{\alpha }_1{u}_{t-1}^2+{\alpha }_2{u}_{t-2}^2+\dots +{\alpha }_q{u}_{t-q}^2\end{array}$$
的极大似然函数为
$$\begin{array}{rl}\log L(\gamma )& =\sum _{t=1}^T\log f\left(y_t\mid x_t,\gamma \right)\\ & =-\frac{T}{2}\log (2\pi )-\frac{1}{2}\sum _{t=1}^T\log \left(h_t\right)-\frac{1}{2}\sum _{t=1}^T\frac{{\left(y_t-x_t^{\prime }\beta \right)}^2}{h_t}\end{array}$$
接下来寻找$\gamma =\hat{\gamma }$时，使得$logL(\gamma )$值最大，

下面是求解步骤：
$$\begin{array}{rl}\frac{\partial \log L(\gamma )}{\partial \gamma }& =-\frac{1}{2}\sum _{t=1}^T\frac{\partial \log h_t}{\partial \gamma }-\frac{1}{2}\sum _{t=1}^T\left[\frac{1}{h_t}\frac{\partial {\left(y_t-{\mathbf{x}}_t^{\prime }\mathbf{\beta }\right)}^2}{\partial \gamma }-\frac{{\left(y_t-{\mathbf{x}}_t^{\prime }\mathbf{\beta }\right)}^2}{h_t^2}\frac{\partial h_t}{\partial \mathbf{\gamma }}\right]\\ & =\frac{1}{2}\sum _{t=1}^T\left(-\frac{1}{h_t}\frac{\partial h_t}{\partial \mathbf{\gamma }}-\frac{1}{h_t}\frac{\partial {\left(y_t-{\mathbf{x}}_t^{\prime }\mathbf{\beta }\right)}^2}{\partial \gamma }+\frac{u_t^2}{h_t^2}\frac{\partial h_t}{\partial \gamma }\right)\\ & =\frac{1}{2}\sum _{t=1}^T(\frac{u_t^2}{h_t^2}\frac{\partial h_t}{\partial \gamma }-\frac{1}{h_t}\frac{\partial h_t}{\partial \gamma }\partial \mathbf{\gamma }-\frac{1}{h_t}\frac{\partial {\left(y_t-{\mathbf{x}}_t^{\prime }\mathbf{\beta }\right)}^2}{\partial \gamma })\\ & =\frac{1}{2}\sum _{t=1}^T\left(\frac{u_t^2-h_t}{h_t^2}\frac{\partial h_t}{\partial \gamma }-\frac{1}{h_t}\frac{\partial {\left(y_t-{\mathbf{x}}_t^{\prime }\mathbf{\beta }\right)}^2}{\partial \gamma }\right)\end{array}$$
整理得
$$\frac{\partial \log L(\gamma )}{\partial \gamma }=\frac{1}{2}\sum _{t=1}^T\left\{\frac{u_t^2-h_t}{h_t^2}\left[\begin{array}{c}-2{\sum }_{j=1}^q{\alpha }_j{\mathbf{x}}_{t-j}{u}_{t-j}\\ 1\\ u_{t-1}^2\\ ⋮\\ u_{t-q}^2\end{array}\right]-\frac{1}{h_t}\left[\begin{array}{c}-2x_t{u}_t\\ 0\end{array}\right]\right\}$$
在上式为0的条件下求出$\hat{\gamma }$即为极大似然估计量，具有一致性。

---

1.4 ARCH模型检验

• 先检验均值方程的变量是否平稳（单位根检验系列）

• ARCH过程存在性检验（LM检验,F检验，LR检验,Q检验）

LM检验

建立原假设：
$$\begin{array}{rl}& H_0:{\alpha }_i=0(i=1,2\dots q)\\ & H_1:{\alpha }_i不全为0\end{array}$$
估计均值方程$y_t={\mathbf{x}}_t^{\prime }\mathbf{\beta }+{\epsilon }_t$，得到估计值${\hat{\epsilon }}_t$,并计算${\hat{\epsilon }}_t^2$，构造辅助回归
$${\hat{\epsilon }}_t^2={\alpha }_0+{\alpha }_1{\hat{\epsilon }}_{t-1}^2+{\alpha }_2{\hat{\epsilon }}_{t-2}^2+\dots +{\alpha }_q{\hat{\epsilon }}_{t-q}^2+v_t$$
得到可决系数$R^2$,构造LM统计量$LM=T{R}^2\sim {\chi }^2(q)$。若$LM>{\chi }_{\alpha }^2(q)$，则拒绝原假设，存在ARCH效应。

F检验

建立原假设：
$$\begin{array}{rl}& H_0:{\alpha }_i=0(i=1,2\dots q)\\ & H_1:{\alpha }_i不全为0\end{array}$$
估计均值方程$y_t={\mathbf{x}}_t^{\prime }\mathbf{\beta }+{\epsilon }_t$，得到估计值${\hat{\epsilon }}_t$,并计算${\hat{\epsilon }}_t^2$，构造两个辅助回归
$$\begin{array}{l}{\hat{\epsilon }}_t^2={\alpha }_0+v_t\\ {\hat{\epsilon }}_t^2={\alpha }_0+{\alpha }_1{\hat{\epsilon }}_{t-1}^2+{\alpha }_2{\hat{\epsilon }}_{t-2}^2+\dots +{\alpha }_q{\hat{\epsilon }}_{t-q}^2+v_t\end{array}$$
第一个辅助回归为约束回归，即$H_0$,第二个为无约束回归即$H_1$,利用有约束的检验方法构建F统计量
$$\mathbf{F}=\frac{\left({\mathrm{SSE}}_r-SS{E}_u\right)/q}{{\mathrm{SSE}}_u/(T-q-1)}\sim \mathbf{F}(\mathbf{q},\mathbf{T}-\mathbf{q}-1)$$
若 $F$<$F_{\alpha }(q,T-q-1)$则不拒绝原假设，反之拒绝。

LR检验

与F检验类似，先构建有约束与无约束的辅助回归，分别利用极大似然估计法得到有约束与无约束的极大似然函数最大值$Log{L}_r$与$Log{L}_u$，构建LR统计量
$$LR=-2\left(\log L_r-\log L_u\right)\sim {\chi }^2(m)$$
若$LR<{\chi }^2{}_{\alpha }(m)$不拒绝原假设，反之拒绝。

Q检验：对均值方程的残差平方（即扰动项方差）进行自相关检验，若存在自相关意味着存在ARCH。

---

2 GARCH模型

在ARCH模型中由于ARCH方程是关于扰动项平方的分布滞后模型，为避免扰动项平方滞后过多（样本容量损失），可以引入扰动项方差的滞后项。例如
$${\sigma }_t^2={\alpha }_0+{\alpha }_1{\epsilon }_{t-1}^2+{\lambda }_1{\sigma }_{t-1}^2\text{\hspace{1em}(4)}$$
此模型称为广义自回归条件异方差模型，记作GARCH(1,1)。其中${\epsilon }_{t-1}$为ARCH项，${\sigma }_{t-1}^2$为GARCH项。${\sigma }_t^2$为${\epsilon }_t$的条件方差。显然
$${\alpha }_0>0,{\alpha }_1\ge 0,{\lambda }_1\ge 0$$
若${\lambda }_1\in [0,1)$则方程(4)可以改写为
$$\begin{array}{l}\left(1-{\lambda }_{1}L\right){\sigma }_t^2={\alpha }_0+{\alpha }_1{u}_{t-1}^2\\ \\ {\sigma }_t^2=\frac{{\alpha }_0}{1-{\lambda }_1}+\frac{{\alpha }_1}{1-{\lambda }_{1}L}{u}_{t-1}^2\\ \\ =\frac{{\alpha }_0}{1-{\lambda }_1}+\left({\alpha }_1+{\alpha }_1{\lambda }_{1}L+{\alpha }_1{\lambda }_1^2{L}^2+{\alpha }_1{\lambda }_1^3{L}^3+\dots \right)u_{t-1}^2\end{array}$$
其中$L$为滞后算子。由此看出，GARCH模型是ARCH模型的无限阶滞后。

---

2.1 GARCH定义

GARCH模型包含$q$个ARCH项与$p$个GARCH项，即GARCH(p,q),
$${\sigma }_t^2={\alpha }_0+{\lambda }_1{\sigma }_{t-1}^2+\dots +{\lambda }_p{\sigma }_{t-p}^2+{\alpha }_1{u}_{t-1}^2+\dots +{\alpha }_q{u}_{t-q}^2$$
其中参数
$$\{\begin{array}{l}{\alpha }_0>0\\ {\alpha }_i\ge 0,i=1,2,\dots q\\ {\lambda }_i\ge 0,i=1,2,\dots p\\ 0\le \left({\sum }_{i=1}^q{\alpha }_i+{\sum }_{i=1}^p{\lambda }_i\right)<1\end{array}$$

---

2.2 GARCH性质

GARCH模型的均值方程的被解释变量的条件期望与方差
{ E { y t ∣ x t } = x t β Var ⁡ { y t ∣ x t } = σ t 2 \left\{\begin{array}{l} \mathbf{E}\left\{y_{t} \mid x_{t}\right\}=x_{t} \beta \\ \operatorname{Var}\left\{y_{t} \mid x_{t}\right\}=\sigma_{t}^{2} \end{array}\right. {E{yt​∣xt​}=xt​βVar{yt​∣xt​}=σt2​​
当然，也可以对GARCH方程求期望，并令$t\to \infty$得到
$${\sigma }^2=\frac{{\alpha }_0}{1-{\sum }_{i=1}^q{\alpha }_i-{\sum }_{i=1}^p{\lambda }_i}$$
要保证${\sigma }^2$非负，需要保证$0\le \left({\sum }_{i=1}^q{\alpha }_i+{\sum }_{i=1}^p{\lambda }_i\right)<1$。在现实建模过程中。GARCH(1,1)与GARCH(2,1)足可以对自回归异方差进行描述。

---

3 I-GARCH模型

在GARCH模型中，
$${\sigma }_t^2={\alpha }_0+{\alpha }_1{\epsilon }_{t-1}^2+\cdots +{\alpha }_q{\epsilon }_{t-q}^2+{\beta }_1{\sigma }_{t-1}^2+\cdots +{\beta }_p{\sigma }_{t-p}^2$$
可能在估计过程中并不能保证
$$\sum _{i=1}^q{\alpha }_i+\sum _{i=1}^p{\beta }_i<1$$
而是出现
$$\sum _{i=1}^q{\alpha }_i+\sum _{i=1}^p{\beta }_i>1$$
的情形。此时扰动项产生的冲击将是永久的。因此I-GARCH模型在GARCH模型的基础上加入限制条件
$$\sum _{i=1}^q{\alpha }_i+\sum _{i=1}^p{\lambda }_i=1$$
以保证平稳性，此时将GARCH模型记作I-GARCH(p,q),
$${\sigma }_t^2={\alpha }_0+\sum _{i=1}^p{\alpha }_i{ϵ}_{t-i}^2+\sum _{i=1}^q{\beta }_i{\sigma }_{t-i}^2$$
限制条件${\sum }_{i=1}^p{\alpha }_i+{\sum }_{i=1}^q{\beta }_i=1$

---

4 M-ARCH模型

定义：将ARCH模型（包括扩展类）的均值方程引入方差作为解释变量，即
$$\{\begin{array}{l}{\mathbf{y}}_t=x_t^{\prime }\beta +\varphi \sqrt{{\sigma }_t^2}+{\epsilon }_t\\ {\sigma }_t^2={\alpha }_0+{\alpha }_1{\epsilon }_{t-1}^2+\cdots +{\alpha }_p{\epsilon }_{t-p}^2\end{array}$$
称为M-ARCH§。当然，也可将GARCH模型均值方程引入方差，得到M-GARCH(p,q)，将IGARCH模型方差方程引入方差，得到M-IGARCH(p,q)模型

---

5 T-GARCH模型

T-GARCH模型又称门限GARCH模型，Glosten, Jagannathan and Runkle (1993)提出了非对称(asymmetric)的“门限GARCH”模型(Threshold GARCH，简记TARCH)。门限GARCH模型根据扰动项的正负在方差方程加入门槛变量，即
$${\sigma }_t^2={\alpha }_0+{\alpha }_1{\epsilon }_{t-1}^2+{\lambda }_1\underset{\text{tarch }}{\underbrace{{\epsilon }_{t-1}^2\cdot 1\left({\epsilon }_{t-1}>0\right)}}+{\beta }_1{\sigma }_{t-1}^2$$
其中$\mathbf{l}({\epsilon }_{t-1}>0)$为示性函数，满足括号里的记作1，反之为0。${\epsilon }_{t-1}>0$表示利好消息，${\epsilon }_{t-1}<0$为利差消息。显然利好与离差的对条件方差的影响是不对称的。

• 利好消息时，扰动项的平方对条件方差的影响为${\lambda }_1+{\alpha }_1$
• 利差消息时，扰动项的平方对条件方差的影响为${\alpha }_1$
• 当${\lambda }_1=0$时，即不存在门槛效应，此时利好与利差的影响是对称的
• 当${\lambda }_1\ne 0$时，条件方差对冲击的反应是非对称的，称这种现象为“杠杆作用”

更一般的TGARCH模型可以记作
$${\sigma }_t^2={\alpha }_0+\sum _{i=1}^q{\alpha }_i{\epsilon }_{t-i}^2+{\lambda }_1{\epsilon }_{t-1}^2\cdot 1\left({\epsilon }_{t-1}>0\right)+\sum _{j=1}^p{\beta }_j{\sigma }_{t-j}^2$$

---

6 E-GARCH模型

标准的GARCH模型存在参数限制，因此可以考虑对数形式的条件方差方程，即E-GARCH模型：
$$\ln {\sigma }_t^2={\alpha }_0+{\alpha }_1\underset{\text{earch }}{\underbrace{\left({\epsilon }_{t-1}/{\sigma }_{t-1}\right)}}+{\lambda }_1\underset{{\text{earch }}_{-}a}{\underbrace{|{\epsilon }_{t-1}/{\sigma }_{t-1}|}}+{\beta }_1\underset{\text{egarch }}{\underbrace{\ln {\sigma }_{t-1}^2}}$$
其中$\left({\epsilon }_{t-1}/{\sigma }_{t-1}\right)$表示${\epsilon }_{t-1}$的标准化，表示非对称效应。$|{\epsilon }_{t-1}/{\sigma }_{t-1}|$表示对称效应。$ln{\sigma }_{t-1}^2$为EGARCH项。无论$ln{\sigma }_{t-1}^2$取何值，都有${\sigma }_t^2=\exp \left(\ln {\sigma }_t^2\right)>0$，故称作“指数GARCH”模型

---

7 ABS-GARCH模型

为了保证方差为正，还可以引入绝对值，即
$${\sigma }_t^2={\alpha }_0+\sum _{i=1}^q{\alpha }_i\left|{\mathbf{\epsilon }}_{t-i}\right|+\sum _{j=1}^p{\lambda }_j{\sigma }^2{}_{t-j}$$
与GARCH模型比较，上述模型用$|{\epsilon }_{t-i}|$代替了${\epsilon }_{t-i}^2$，采用绝对值形式缩小了${\epsilon }_t^2$的幅度。当然也可将ABS-GARCH模型写作
$${\sigma }_t={\alpha }_0+\sum _{i=1}^q{\alpha }_i\left|{\epsilon }_{t-i}\right|+\sum _{j=1}^p{\lambda }_j{\sigma }_{t-j}$$

---

8 P-GARCH模型

考虑冲击的非对称性，还可以用P-garch模型来刻画，即
$${\sigma }_t^k={\alpha }_0+\sum _{j=1}^q{\beta }_j{\sigma }_{t-j}^k+\sum _{i=1}^p{\alpha }_i{\left(\left|{\epsilon }_{t-i}\right|-{\gamma }_i{\epsilon }_{t-i}\right)}^k$$
其中$k>0$,$|{\gamma }_i|\le 1(i=1,2,\dots r)$

---

9 ARMA-GARCH模型

设线性模型为
$$y_t={\mathbf{x}}_t^{\prime }\mathbf{\beta }+u_t$$
其中扰动项服从ARMA(p,q)过程，即
$$u_t=\sum _{i=1}^p{\rho }_i{u}_{t-i}+\sum _{j=1}^q{\theta }_j{\epsilon }_{t-j}+{\epsilon }_t$$
其中$\epsilon$为ARCH或GARCH模型的扰动项，代入方程
$$y_t={\mathbf{x}}_t^{\prime }\mathbf{\beta }+\sum _{i=1}^p{\rho }_i\left(y_{t-i}-{\mathbf{x}}_{t-i}^{\prime }\mathbf{\beta }\right)+\sum _{j=1}^q{\theta }_j{\epsilon }_{t-j}+{\epsilon }_t$$

---

10 其他GARCH模型

ARCH模型的扩展形式非常丰富(还有LM-GARCH,FI-GARCH,FIE-GARCH等)，以均值方程的设定形式又分为线性与非线性类型。总之，ARCH模型及其扩展的选择应从研究目的来做出选择。

---

-END-

参考文献

陈强. 高级计量经济学及Stata应用[M].高等教育出版社，2014