自回归移动平均模型解析

第二章自回归移动平均模型

一些金融时间序列的变动往往呈现出一定的平稳特征，由 模型就是借助时间序列的随机性来描述平稳序列的相关性信息, 行建模和预测。

第一节ARMA 模型的基本原理

ARMA 模型由三种基本的模型构成：自回归模型( AR, Auto-regressive Model )，移动平 均模型(MA ，Moving Average Model )以及自回归移动平均模型 (ARMA ，Auto-regressive Moving Average Model )。

2.1.1自回归模型的基本原理 1. AR 模型的基本形式

AR 模型的一般形式如下：

办乂「1办」? 2办上 ....... \%申；t

其中，c 为常数项，'1, 2^ \模型的系数，；t 为白噪声序列。我们称上述方程为 p

阶自回归模型，记为 AR(p )。 2. AR 模型的平稳性

此处的平稳性是指宽平稳，即时间序列的均值，方差和自协方差均与时刻无关。 即若时

间序列｛%｝是平稳的，即 E(y t)=^, Var (y t

2

, Cov(y t , y —) =

。

为了描述的方便，对式(2.1 )的滞后项引入滞后算子。若

y t

二x t

j ，定义算子“ L ”，

r .

k

使得y t =Lx t =为4 L 称为滞后算子。由此可知，

L 人=X t±。

对于式子(2.1),可利用滞后算子改写为：

y t =c 丄％

2

L 2%

p

L P %

t

移项整理，可得：

(1- 丄- 2L 2

- - p L p

)y t 二c ；t

Box 和 Jenkins 创立的 ARMA 并由此对时间序列的变化进

AR(p)的平稳性条件为方程1 - 1L - qL2 -…-pL^0的解均位于单位圆外。

3.AR模型的统计性质

(1)AR模型的均值。

假设AR(p )模型是平稳的，对AR(p )模型两边取期望可得：

根据平稳序列的定义知，E(yJ = ? I ,由于随即干扰项为白噪声序列，所以E(；J = 0, 因此上式可化简为：

(1 - 1 - 2 -…- ')亠八0

1 一％-％―…_?p

(2)AR模型的方差。

直接计算AR( p)模型的方差较困难，这里引入Green函数。

AR(p )模型可以改写成如下形式：

y t

：(L)

设■ j…■ p为平稳AR(p )模型的反特征根，则

p

:：J(L) =1一丄一2L2-川一p L p"【(1- 丄)。

i =1

进一步，

p k p p 『：：

y t 亠.t k i(' i L V t k i 订；t「j = ' G j ；t_j

i=11—，-i L imj=o j=0im j =0

p

其中，k i为常数，G j八.k i■ i j，称为Green函数，因为'' p均在单位圆内，所

i =1

以Green函数是呈负指数下降的。

对上式两边取方差，可得：

C3O

2

var(y t) =、G j var( ；t_j)

j=0

由于随机干扰项为白噪声序列，所以var(；t_j)=匚2。因为Green函数是呈负指数下降,

-- 2

所以G j ：：：：:，这说明平稳时间序列方差有界，且等于常数、G j务2。

j =0 j=P

(3)自协方差函数。

假设将原序列已经中心化，贝V E(y t) =0，则对AR( p)模型等号两边同时乘以

y t上(-k _1)，两边取期望得：

E(ytytjJ = 土卜少―) jEWt^ytjJ …? pE(yt_pytjJ E(心…)

因为当期的随机干扰项与过去的时间序列值无关，所以：E(；t y t上)=0。因此，上式可以化为：

G二】「丄2—'……

其中r k，表示k阶自协方差。

2.1.2 移动平均模型的基本原理

1.MA模型的基本形式

MA模型的一般形式如下：

办=U J 」二2 ；2…乜2

其中，u为常数项，RC2…Jp为模型的系数，；t为白噪声序列。我们称上述方程为q 阶移动平均模型，记为MA(q)。

2、MA模型的可逆性

对于一个MA( q)模型：

y t = U ? ；t 7 2 宀2 2 …T q ；tT

将其写成滞后算子的形式：

y t-u=(1,丄F2，q L q)；t

若方程1 ?? 1L2? Tq L q=0的根全部落在单位圆外，则称MA模型是可逆的。可逆性可以保证MA模型可以改写成：

'■■(L)(y t -u^ ’

即MA模型可以转化为AR模型，同时可以保证参数估计的唯一性。

3、MA模型的数字特征

（1）均值

当q ：：:：:时，对于一般的MA （q）模型：

y t = U " ^-t J ^2 -t _2 L 入討_q

两边取期望，可得：

E（yJ = E（u ? “；2 I ；t, * L Vq ；t^）二u

即一般的MA（q）模型的期望值即为模型中的常数项。

（2）方差

对MA（q）模型，两边取方差：

Var （y t） =Var （u ? t ：：勺；t」：：*2 t _2 ?……为戈卫）=（1 W -……V ）；『

（3）协方差函数

r k 二Elyy上）=E[（u ?上T t」二t 2 ?…F ；t』）（u y T 2丄 V '… 為化简可得：匚2（1 V2；；?川v：）,k =0 k二F卞1 ?川启山）,0你"

0,k q

2.1.3自回归移动平均模型的基本原理

1、ARMA模型的基本形式

ARMA模型的一般形式如下:

y t 二C「2丫2 一? ；t 二1 ；2 二2 2 一二p ；tT

显然ARMA（p,q）模型可看成是AR（p）模型和MA（q）模型相结合的混合形式。

2、ARMA模型的平稳性和可逆性

对于一个ARMA（p，q）模型，

y t =C「1丫2「2人,…；t 二1 ；t J 二2 2 f p

将其写为滞后算子的形式：

（1 一丄一\L2-1_一p L p）y t 二c （1 JL ^L2L jq L q）t

两边同时除以（1 - ]L - ；L2-L-\L p）

y t 二亠冋（L）；t

其中:

叫)二

由此可以看出，ARMA模型的平稳性完全取决于AR(p)模型的参数，与MA (q)模型的参数无关。

类似地，ARMA模型的可逆性完全取决于MA( q)模型的参数，与AR ( p)模型的参数无关。

3、ARMA模型的数字特征

(1)期望

对于一个一般的ARMA(p,q)模型两边同时取期望，化简得:

E（yJ =

(2)自协方差函数

厲=E（yt y 十）=E[（艺G J）（无Gj fl t*J]

<=0 j=0

=E[.一G i ?_ G j ;t _i ；t k _j ] i ￡ j ￡

□0

=二2 ''G i G i k

i ￡

第二节时间序列的相关性分析与平稳性

2.2.1时间序列的自相关系数

2.2.1.1自相关函数(ACF

1、AR(p)的自相关函数

在上一节中已经介绍了AR( p)模型的协方差函数满足下式：

r

k =人―2「心'.……「P—

由于自相关系数坯二土，因此:

r o

人二1凡丄-2

；k 2……

「P 凡』

该式表示自相关系数满足 p 阶差分方程。根据差分方程解的性质， 上差分方程的通解可 以写为:

p

珥k ）八G

i=1

其中，C i 为任意不全为 0的常数， 是滞后多项式的反特征根。根据平稳性的性质, 闊灯：【。从自相关系数的一般形式可看出， 闻［:始终不为0，但是随着滞后阶数的增加，自 相关系数慢慢逼近 0，在图形上表现出一定的拖尾性。

2、MA 模型的自相关函数

根据上一节推导的 MA 模型的自协方差函数的表达式，

MA 模型的自相关函数表示为：

1, k =0

忙=丄=1 Fl ?川启/q 4

截尾的。

3、ARMA

模型的自相关函数

根据ARMA 模型的自协方差函数，不难得到 ARMA 模型的自相关函数:

QO

■./ 为 G j G i*

4 =」=

-----

k :.

、G i 2

i =0

由此可以看出，ARMA 模型的自相关函数不具有截尾性。事实上， ARMA 模型若满足可

逆性，其形式相当于一个无穷阶的 AR 模型，因此自相关函数与 AR 模型一样具有拖尾性。

2.2.1.2 偏自相关函数(PACF

1、偏自相关函数的定义

自相关函数卜；不能纯粹地表示 霜与之间的相关性，两者的相关性还会受到：辽-，亠..

;加护……■-的间接影响，为了单纯地表示 .■与 : 之间的相关性，这里引入偏自相关函数。 偏自相关函数表示在固定

……沁I 的情况下豁与聘出之间的相关性。下面介绍

,0 : k _q

1 V ￡ ?川応 0,k q

因此，当k>q 时，自相关函数为 0,也就是说MA (q )模型的自相关函数在 q 步以后是

偏自相关函数Lhj的计算方法。

设序列yt可由下回归方程估计：

y t =时y t i * 人2 y t 2 ■ IH ■ kk 丄y t 上 1. y t ±" a

根据回归方程的性质，式中估计系数翠慝即为偏自相关函数。为了估计回归系数，采用OLS方法，即

L = E( y t - k1y tl - ' k2 y t _2 「kk 丄y t _k _ ' kk y t Jk )达到取小。

对L关于各回归系数求偏导，可得到以下方程组：

：：1 ,2「HI 'kk 讥

：2 = k1 i 1 ：% k2 i 0 H I ' kk i k 2

―k

该方程组称为Yule- Wolker方程。根据自相关系数，求解Y-W方程即可得到偏自相关系数。

2、AR (p)的偏自相关函数

对于AR(p濮型，匕时，

P k

L ~ E((耳+-粥)斤可▼ g %%』尸

i=i j*

由于??与序列的滞后项无关，因此一：7 ,且当

:=打jP

"0,p 1：：jEk

由此，AR(p)模型的偏自相关函数如愛在k>p后等于0，即AR(p)模型的偏自相关函数具有截尾性。事实上，AR模型偏相关函数的截尾性也可直接从该模型的表达式看出。AR( p) 模型实质上假设序列至多只与滞后p阶的值相关，因此偏自相关函数至多在p阶处非0。

3、MA (q)和ARMA ( p, q)的偏自相关函数

由于MA (q )和ARMA ( p, q)相当于无穷阶的AR模型，因此这两个模型的偏自相

关函数均不具有截尾性，而是拖尾性。

221.3 ARMA 模型自相关系数与偏自相关系数的估计与检验

根据以上分析，不同 ARMA 模型自相关系数与偏自相关系数的表现存在明显的差异。

表2.1给出了三类模型 ACF 与PACF 的特征。

因此，我们可以通过观察偏自相关函数来识别并确定

AR 模型的滞后阶数，通过自相关

函数来识别并确定 MA 模型的滞后阶数q 。那么对于给定的样本数据，如何估计样本 ACF 与

PACF 并从统计角度检验两者是否为 0呢？下面分别介绍 ACF 与PACF 的估计与检验。

1、样本ACF 与PACF 的估计与实现

而言，样本自协方差表示为:

1 T_k

%二百》侨一 0如+厂刃

j=i

其中「表示样本均值。 那么 SACF=.。

对于PACF 主要是利用 Yule- Wolker 方程求解。 当滞后阶数较大时，Y-W 方程直接计算较

难，目前多采用递推算法来求解。 2、样本ACF 与PACF 的显著性检验

若序列满足独立性，则由统计渐进分布的有关定理可知，当样本个数充分大时, PACF 均满足均值为0，方差1/T 的正态分布，即

因此若區|：；［?匸，|轨「：1.筠\匸，则可认为样本数据是独立的，即自相关系数

和偏自相关系数均不显著异于

0。该检验法即为正态检验法。

Portmanteau 检验法是联合检验法，即检验直到 k 阶的自相关系数是否同时为

0。该检

验法使用Q 统计量进行检验。Q 统计量具体形式为：

Q=T(T + 2)

》去

j=i

对于给定样本，只需估计样本的自协方差与方差, 将两者相除即可得到样本

ACF 具体

常用的检验方法主要包括两类：正态检验法和

Portmanteau 检验法。

ACF 和

其中T为样本容量，k为设定的滞后阶数。Q统计量服从诜I；；分布。当Q统计量超过设定的临界值时，就拒绝原假设，即序列至少存在k阶以内的自相关性。

2.2.2 时间序列平稳性检验

建立ARMA的前提是序列是平稳的。检验平稳性常用的方法主要有三种：经验法、自/偏自相关系数法、单位根检验法。

1经验法

经验法是通过观察图形的方式来初步判断时间序列是否平稳的。首先画出时间序列的

图形，如果该图形围绕某一直线上下以较小的幅度波动，则该序列一般是平稳的，否则是不

平稳的。

2、自/偏自相关系数法

由于ARMA模型的自/偏自相关系数要么是截尾的，要么是拖尾的，因此可以观察时间

序列的自/偏自相关图，如果时间序列的自/偏自相关系数从某个滞后期开始均与0无差异，可以认

为该时间序列是平稳的；若自/偏自相关系数衰减很慢，且与0存在明显的差异，贝U 时间序列是非

平稳的。

3、单位检验法

常用的单位根检验法主要包括DF检验法和ADF检验法。

(1) DF检验法

DF检验包括三种形式：

% 二九? ；t

y t 八y tv c ；t

% 二L y t4 c t t

其中，c为常数项，t表示线性趋势，随机干扰项独立同分布，且服从N (0,二2)。根据平稳性的概念，若序列y t是不平稳的，则回归系数防辺。一般乍：〉1较易识别。因此判断序列y t是否平稳，主要是判断是否为1。如果］=1，则说明序列存在单位根，是不平稳的，否则是平稳的。

yi = yn a

L y t = y t c 亠冷

.：y t = y t二c 亠E t ；t

其中，，因此可以将DF检验的原假设和备择假设分别为：

H。：

已：：、：：:0

相应的统计量为：

P

DF=^—

std（p）

DF的形式与t统计量相似，但是该统计量并不服从t分布，Dickey和Fuller（1979）给出了利用蒙特卡罗模拟方法模拟的临界值，因此该检验称为DF检验。DF检验是左侧检验，且不同形式的方程临界值是不同的。注意DF检验只有当时间序列为AR（1）过程时才有效。如果存在高阶滞后相关，那么将违背随机干扰项独立同分布的假设。因此，Dickey-Fuller提出来ADF检验来弥补DF检验的不足。

（2）ADF检验

假设时间序列存在p阶自相关，那么用p阶自回归方程来判断单位根，形式为：

y^ i y t丄；yt2 申；t

上式两边同时减去y t ，通过整理可得：

p丄

闷=讪丄、:i "y t^ ■ *

p p丄

其中，i T， S「-' j

i=±j=t-H

上述检验形式是在DF检验方程中加入了力的高阶滞后项，因此可以看成是DF检验的增广形式，简称ADF检验。与DF检验类似，ADF检验也存在三种形式，

p_L

?y t二讪1亠二'：.汇?寸

i 士

p -L

"y t 二c f 、：i 3t 丄"

i 士

Pl

绍=c ?、：t ? Ty t 丄一二-i .-y t 丄? t

i .1

不难看出，当p=1时，ADF检验就是DF检验，因此DF检验是ADF检验的特例。

检验的原假设和备择假设为：

H o： 40

H i :匸：：：0

即原假设为序列至少存在一个单位根，备择假设为序列不存在单位根。使用ADF检验时，应该注意如下几个问题：

首先，要确定合理的滞后阶数。

其次，因为检验统计量的临界值依赖于方程的形式，因此选择检验的方程形式很重要。

再者，如果检验的结果是拒绝原假设，那么原序列就不存在单位根，即原序列是平稳的；如果接受原假设，则序列是不平稳的，需要进行若干次差分，直到拒绝原假设，从而确定序

列单整的阶数。

例2.1 以2003年2月到2010年4月上证国债（交易代码000012 ）月末收盘指数为原始数据进行分析。首先对指数序列取对数，然后对这个收益率序列进行单位根检验，结果如下：

表2.2上证综指收益率单位根检验表

Null Hypothesis: LSP has a unit root

Exogenous: Constant, Linear Trend

Lag Length: 0 (Automatic based on SIC, MAXLAG=11)

从检验结果来看，应该接受原假设，即认为该序列存在单位根。

第三节季节性ARMA模型

2.3.1 时间序列的季节性

GDP、消费支出等。它们在现实中，存在一些经济时间序列具有明显的季节性，比如

具有共同的特点：有规律，每年重复出现，表现为逐年同时期有相同的变化方向和大致相同

的变动幅度。一般我们认为季节性是一种对数据的干扰，在数据分析中通常将其过滤。那么

该如何识别时间序列的季节性。常用的分析方法是自相关图分析法。

观察时间序列的自相关图，如果是月度数据，可以观察滞后期为12、24、36……的自

相关系数；如果是季度数据，可以观察滞后期为4、8、12……的自相关系数。如果自相关

系数与0无差异，那么同月或者同季度之间不相关，即不存在季节性，否则该时间序列具有

季节性。值得注意的是，若时间序列的趋势性较强时，这种趋势性会掩盖季节性，因此分析季节性时，应该先消除时间序列的趋势性。

由于季节性是对数据的一种干扰，因此在数据分析中常常消除季节性。常用的处理季节性的方法主要包括三种：（1）移动平均法；（2）用X-11方法去掉时间序列的季节性；

（3）

采用逐期差分和季节差分，建立SARMA模型。下面先简单介绍一下移动平均法和X-11方法，再重点介绍SARMA模型。

1、移动平均法

移动平均法是将原时间序列的两个或多个时期的数据进行平均，用平均值代替原序列

值，以此来平滑周期性。假设周期长度为s,则在时点t上的移动平均值可以表示为：

s丄

' %丄

MA t=広

s

2、X-11方法

X-11方法最初是由美国人口普查局在1965年研究开发出来的方法。它的基本原理是滑动平均法法。这种方法能适应各经济指标的性质，并根据季节调整的目的进行相应的调整及进行算法的选择。

X-11方法包括两种模型：加法模型和乘法模型。

加法模型的基本形式为：

y t 二人S h U

乘法模型的基本形式为：

y t =T t S t I t D t

其中T t为趋势循环要素，S t为季节要素，I t为不规则变动要素，D t为周工作日变动要素。常用的是乘法模型，但注意乘法模型仅适用于时间序列数据为正的情况。

例2.2以江苏省2001年1月到2009年12月的每月的社会消费品零售总额为原始数

据，进行季节性分析。数据来源于中经网统计数据库。数据的时间序列图如下。

从图中可以看出，该序列存在一定的周期性季节性变动。采用X-11方法进行季节性调整，结果如下：SERIES0SA

从图中可以看出，经过调整，季节性因素与不规则因素得到了一定程度的消除。

2.3.2季节性ARMA模型

SARMA模型是在回归项中加入季节自回归项（SAR和季节移动平均项（SMA）来考虑季节性的。假设季节长度为S,则

SMA（q）（m）S的一般形式为：

y t -u =（1 nL ? PL2-Vq L q）（i ? Ps L S rs L2S - rs L mS）t

其中q为一般移动平均项的滞后阶数，m为季节性滞后阶数。

SAR（ p）（k）沖勺一般形式为：

（1 一丄一2L2-1|1一p L P）（1—!S L S—2S L2S -1|1一kS L kS）y t ? I

S

SARMA的一般形式只需将SAR和SMA合并即可，简记为：SARMA(p,q)(k,m)。

为了进一步理解季节模型的特点，以滞后两阶的AR模型为例，无季节性的AR (2)模

型如下：

y t =c 叽2%N；t

用滞后算子表示：

(1 - 丄 - 2L2)y t= c ；t

对于带有季节性的季度数据，若k=1，则模型就会变成下式：

(1 -丄- 2L2)(1 - - L4)y t二c ；t

它等价于：

y t =c「办」」2%丄一勺心」y t_6 ；t

因此该季节性AR模型实质上具有稀疏系数的更高阶AR模型。

SARMA模型的自相关系数与偏自相关系数与普通的ARMA模型类似，只是ACF和PACF 的计算需要计算滞后期分别为S 2S……mS处的大小，性质与ARMA模型相同，即季节自回

归模型的季节偏相关函数按季节周期的增加而截尾；季节移动平均模型的季节自相关函数按

季节周期增加而截尾；混合的季节自回归移动平均模型的季节自相关函数和偏相关函数均按指数衰减。

例2.3引用例22的数据进行季节性分析。

为了便于分析，首先将数据取自然对数，然后进行差分去除趋势性。下图是新序列的自相关图，结果如下:

Autocorrelation Partial Correialion AC PAC Q-Stat Prob

从图中可以看出，该序列在滞后6期、12期处的自相关函数和偏自相关函数显著异于0,因此可认为该序列存在一个周期为6的季节性变动。建立ARMA模型时，可在模型中加

入周期为6的sar或sma项。

第四节ARMA模型的构建与Eviews实现

2.4.1 ARMA模型的具体构建步骤

ARMA模型的具体运用主要包括以下几个步骤：

1、判断序列的平稳性

由于ARMA模型只适用于平稳时间序列，因此在建立ARMA模型前首先要检验序列的平稳性。若序列是平稳的，则可以直接建立ARMA模型，特别当存在季节性时，应建立SARMA 模型。若序列是非平稳的，则应做差分使得序列变成平稳序列后才能建立ARMA模型。

2、ARMA模型滞后阶数的选择

一般根据自相关图和偏自相关图初步确定ARMA模型的滞后阶数可能取值，再根据各

可能滞后阶数估计ARMA模型的AIC或SC信息准则做进一步判断。

3、ARMA模型的参数估计与检验

ARMA模型参数主要是利用极大似然法进行估计。该模型的检验主要是对残差是否是白

噪声进行检验。

4、模型的预测

模型的预测包括动态预测和静态预测，其中动态预测是在给定样本值时往前多步预测，

其中样本区间长度不变，静态预测是给定部分样本值时往前一步预测，其中样本区间会发生调整。

2.3.2 ARMA模型在Eviews中的窗口实现

1判断序列的平稳性

根据前文的介绍，平稳性的判断主要包括三种方法：经验法、自相关图法和单位根检

验。经验法主要是观察时间序列图形。以例 2.1的数据为例，计算上证综指收益率序列，单

击工具栏中的View/Graph，显示图形如下：

从图可以看出，该序列基本在

0上下发生小幅度波动，因此从图形上可以初步判断该

序列是平稳的。

自相关图检验法是在序列显示窗口中单击view/correlogram，出现如下窗口:

SERIES01

选择最大滞后阶数，即可该序列的 ACF 、PACF 置信区间以及 Q 统计量。具体显示结

构可见下节的例表 2.4。

单位根检验主要是利用 DF 检验法或者 ADF 检验法。在序列显示窗口中单击 view/unit

root test …，即出现单位根检验窗口:

左上方test type 是选择单位根检验的方法， 默认的选项是 ADF 方法，此外还包括phillips- perron 、Dickey-Fuller GLS 等，主要是用来处理序列的异方差和自相关的。

下方Test for unit root in 是选择检验的对象，其中 差分序列，2nd differenee 表示二阶差分序列。

In elude in test equation 是选择检验的形式，其中 项，trend and intereept 是指同时包括趋势项和截距项,

势项。检验形式的选择往往是根据该序列的图形特征来确定， 若该序列在图形上无明显的趋 势且明显异于 0时，应选择intereept ;若该序列既无趋势也与 0无差异时，应选择 none ;

若序列呈现出一定的趋势性，则选择

trend and intereept 。

右下方lag length 是选择滞后阶数。User speeified 是用户自己定义滞后阶数， automatie

seleetion 是在maximum lags （最大滞后阶数） 范围内根据信息准则选择最优的滞后阶数，

其

中信息准则选项包括 Akaike Info Criterion （AIC ,默认）、Sehwarz info eriterion （SQ 、Hannan- Quinn eriterion 等，前两种最为常见。

选择各选项后，点击确定即出现单位根检验结果，可见示例

2.1中表2.2。在表中给出

了检验统计量t 的值和相应的p 值，t 值下方是1%、5%和10%显著性水平下的临界值。 若 t 值小于临界值，则应拒绝原假设，即原序列是平稳的，否则存在单位根。

level 表示原序列， 1st differenee 表示

intereept 是指检验方程中只包括截距 none 是既不包括截距项也不包括趋

2、ARMA 模型滞后阶数的选择与模型的估计

由于AR(p 濮型的偏自相关函数在 p 阶处截尾，MA(q)模型的自相关函数在 q 阶处截尾, 因此可以通过观察样本自相关函数和偏自相关函数的截尾特征来确定 择不同p 和q 的各种可能，再根据信息准则来筛选最优的滞后阶数。

AIC 和SC 定义分别为:

从AIC 和SC 的定义可以看出，残差方差增加或者滞后阶数增加都会使得 AIC 和SC 增加。

而增加滞后阶数往往会使得残差方差减少，

因此最优滞后阶数的选择应使在减少方差和减少

自由度之间进行权衡， 使得AIC 和SC 达到最小。值得注意的是，AIC 和SC 的基本形式类似， 但是SC 对于自由度损失的惩罚要大于

AIC 。

选择一组p 和q 后，对ARMA 模型进行估计。 点击主菜单中的 Quick/Estimate Equation ， 弹出如下对话框：

确定 I [—取消

若估计AR(p)模型，在对话框中输入

y c ar(1)??….ar(p);若估计MA(q)模型，在对话框中

输入y c ma(1)…?ma(q);若估计 ARMA 模型，在 对话框 中输入y c ar(1) ??….ar(p) ma(1)??…ma(q)。再单击 OK 即可。也可以直接在命令行中输入

ls y c ar(1) ??….ar(p)

p 和q 的范围，从而选 常用的信息准则主要是

AK = log

2k

T

其中此表示残差， 表示样本容量， k 表示滞后阶数。

Specification tions

Eqpia.ti on speci ￡ication

Dependent variable followed by list of regressors and PEL terms, OR an e^plicit cn like