第9章平稳时间序列分析

Download Report

Transcript 第9章平稳时间序列分析

第9章
平稳时间序列分析
平稳时间序列分析
9.1 时间序列的概念
9.2 时间序列模型
型
9.2.1
9.2.2
9.2.3
9.2.4
白噪声序列
自回归模型
移动平均模型
自回归模型转化为移动平均模
9.3 自回归模型的平稳性和相关函数
9.3.1 自回归模型的平稳性
9.3.2 自回归模型的自相关函数
平稳时间序列分析
9.4 自回归模型定阶和估计
9.4.1 自回归模型定阶
9.4.2 自回归模型估计
9.4.3 自回归模型再定阶—信息准则
9.5 自回归分布滞后模型
9.5.1 自回归分布滞后模型
9.5.2 格兰杰因果关系检验
9.6 ARCH模型
9.6.1 ARCH模型的定义
9.6.2 ARCH模型估计
重要概念
9.1 时间序列的概念
• 设时点 t  1,2,,T 处的观测为随机变量 y1, y2 ,, yT ,
T
{
y
}
这些随机变量形成一个时间序列,记为 t t 1
或者 {yt , t  1,2,,T} ,y1, y2 ,, yT 的一组具体取
值称为时间序列的实现值(realization)。
• 自相关函数(ACF:AutoCorrelation function )
cov(ys , yt )
 (s, t ) 
  (s, s  k ),其中k  t  s
Var ( ys ) Var ( yt )
9.1 时间序列的概念
定义1(平稳性):如果时间序列 {yt , t  1,2,,T}
的数学期望、方差和协方差不随时间变化,
即
E( yt )   , Var ( yt )   2 , Cov( yt , yt  k )  C (k ),
t  1,2,,T ; k  1,2,T
{称
yt }Tt1
为宽平稳(wide-sense
stationary)时间序列。宽平稳也称为协方
差平稳或者二阶矩平稳。
9.1 时间序列的概念
• 严平稳:时间序列中任意一组随机变量的
联合分布不随时间发生变化,即对任意一
组时间点 t1  t2    tn 和时间间隔 s ,
{ yt , yt ,, yt } 的联合分布与 { yt , yt ,, yt }
T
{
y
}
的联合分布相同,称 t t 1 严平稳。
1
2
n
1 s
2 s
ns
• 二阶矩存在的严平稳时间序列一定宽平稳,
宽平稳的时间序列不一定严平稳,本书只
讨论宽平稳,将宽平稳时间序列简称为平
稳时间序列。
9.1 时间序列的概念
T
{
y
}
• 若 t t 1 为平稳时间序列,则:
cov( ys , ys  k )
C (k )

(1) (k )   ( s, s  k ) 
2
2


T
{
y
}
(2) t t 1 满足大数定律,因此
ˆ 
1
T
t 1 yt  y, ˆ 
T

T k
1
ˆ (k ) 
2
T
t 1
1
T 1
1
T
  (k )
2
(
y

y
)
t 1 t
T
( yt  y )( yt  k  y )
T
2
(
y

y
)
 t
t 1
2


分别是
、 和 (k ) 的一致估计。
9.1 时间序列的概念
• yt  k 表示 yt 的 k 阶滞后,用滞后算符 L 表示
为
Lk yt  yt k ,
k  0,1,2,
例如 2 yt 3  5 yt 2  3 yt 1  2 yt
用滞后算符多项式表示为:
(2L3  5L2  3L  2) yt
9.2 时间序列模型
9.2.1
9.2.2
9.2.3
9.2.4
白噪声序列
自回归模型
移动平均模型
自回归模型转化为移动平均模型
9.2 时间序列模型
9.2.1 白噪声序列
定义2(白噪声):如果时间序列 { t , t  1,,T}
满足:
2
E
(

)

0
Var
(

)


(1) t
t, 
 s 和 不相关,即
E( t s )  0
(2)对任意 s  t  t ,
{称
为白噪声序
t , t  1,, T }
列,简称白噪声 (white noise)。
是平稳时间序列的极端例子。
9.2 时间序列模型
9.2.2 自回归模型
一阶自回归模型AR(1)
{ t } 为白噪声
yt  c  1 yt 1   t ,  t ~ N (0, 2 ) | 1 | 1 ,
除了常数项以外,y 在 t 时刻的值由前定项
(predetermined term)1 yt 1和与前期值不
相关的新息(innovation) t 组成。
k 阶自回归模型AR(k)
yt  c  1 yt 1  2 yt 2  ... p yt k   t ,  t ~ N (0,2 )
| 1  2  ...  k | 1
9.2 时间序列模型
9.2.2 自回归模型
t ~ N (0,1)
模型yt  1.5  0.7 yt 1   t ,
150个样本的时序图:
一阶自回归时序图
9.00
8.00
7.00
6.00
5.00
4.00
3.00
2.00
1.00
0.00
1
10 19 28 37 46 55 64 73 82 91 100 109 118 127 136 145
9.2 时间序列模型
9.2.3 移动平均模型
对一阶自回归模型进行递推:
yt  c  1 yt 1   t
 c  1 (c  1 yt  2   t 1 )  c  1c  12 yt  2  1 t 1   t

 c(1  1  12    1k 1 )   t  1 t 1    1k 1 t  k  1k yt  k
当
k 
时,1k  0 ,
c
yt 
  t  1 t 1    1k 1 t  k  
1  1
9.2 时间序列模型
9.2.3 移动平均模型
取有限项,
yt  c1   t  1 t 1    k t  k
上式即为 k 阶移动平均模型。
9.2 时间序列模型
9.2.4 自回归模型转化为移动平均模型
( L) yt  (1  1L) yt  c   t
由于
(1  1L   L  )1  1L)  1
2 2
1
在上述一阶自回归模型两边同乘  1 ( L)
c
yt  (1  1L) yt 
 (1  1L  12 L2  ) t
1  1
1
c

  t  1 t 1  12 t  2  
1  1
9.2 时间序列模型
9.2.4 自回归模型转化为移动平均模型
• 可以转化为移动平均模型的自回归模型称
为可逆的(invertible)。
• 从上面推导可以看出,一阶自回归模型可
逆的条件是 | 1 | 1 。实际上,自回归模型的
可逆条件,是滞后多项式的根在单位圆外。
滞后多项式即:
(L)  (1  1L  2 L2    k Lk )
9.3 自回归模型的平稳性和相关函数
9.3.1 自回归模型的平稳性
9.3.2 自回归模型的自相关函数
9.3 自回归模型的平稳性和相关函数
9.3.1 自回归模型的平稳性
k 阶自回归模型:
yt  c  1 yt 1    k yt k   t
 t ~ N (0,  2 ) 且{ t }为白噪声序列
用滞后算子表达上式为:
(1  1L    k Lk ) yt  (L) yt  c   t
9.3 自回归模型的平稳性和相关函数
9.3.1 自回归模型的平稳性
结论1:自回归模型平稳的充分必要条件为:
滞后多项式的根都在单位圆之外,即方程
( L)  0
的根 L 满足 | L | 1 。其中 L 为实根时 | L |表
示绝对值,L 为虚根时 | L | 表示虚数的模。
9.3 自回归模型的平稳性和相关函数
9.3.1 自回归模型的平稳性
• 滞后多项式的单位根均在单位圆内,则时
间序列平稳;若有根为1,则不平稳,此时
称存在单位根。
• 把模型是否平稳的检验称为单位根检验。
9.3 自回归模型的平稳性和相关函数
9.3.1 自回归模型的平稳性
例子9.1
yt  1.5  0.7 yt 1   t
为平稳时间序列
yt  yt 1  yt  2   t
为
非平稳时间序列
9.3 自回归模型的平稳性和相关函数
9.3.2 自回归模型的自相关函数
对 k 阶自回归模型:
yt  c  1 yt 1    k yt k   t
 t ~ N (0,  2 )
且 { t } 为白噪声序列
若 { yt }为平稳时间序列,则 Eyt    Eyt k  
两边取期望,得   c /(1  1    k )
代入原模型,整理可得零均值化的 k 阶自回归模型:
yt    1 ( yt 1   )    k ( yt k   )   t
9.3 自回归模型的平稳性和相关函数
9.3.2 自回归模型的自相关函数
重新将 yt    1( yt 1  )   k ( yt k  )  t
记为 yt  1 yt 1    k yt k   t ,即用 yt 表示原
模型中的 yt   。此时,Eyt  0。后面的自回
归模型都将采用这种零均值化后的模型 。
9.3 自回归模型的平稳性和相关函数
9.3.2 自回归模型的自相关函数
• AR(1)模型的自相关函数
用  0 表示变量的方差
 0  E( yt2 )  E(1 yt 1   t )2  12E( yt21)  E( t2 )  12 0   2
得出
 2
 0  Var ( yt ) 
1  12
 1 表示 yt 与 yt 1 的协方差
1 2
 1  Cov( yt yt 1 )  E[(1 yt 1   t ) yt 1 ]  1 0 
1  12
9.3 自回归模型的平稳性和相关函数
9.3.2 自回归模型的自相关函数
• AR(1)模型的自相关函数
同理,
k 2


 k  Cov( yt , yt  k )  1k  0  1 2 , k  1,2,
1  1
则AR(1)的自相关函数为  (k )   k /  0 。
9.3 自回归模型的平稳性和相关函数
9.3.2 自回归模型的自相关函数
结论2:AR(1)模型的自相关函数(ACF)为
 (k )   , k  1,2,
平稳性要求 | 1 | 1,当k   时, (k )  0 ,
k
1
即自相关系数随时间间隔增加指数递减到0,
但不等于0。这种现象称为自回归模型自相
关函数的拖尾性。
9.3 自回归模型的平稳性和相关函数
9.3.2 自回归模型的自相关函数
例子9.2
yt  1.5  0.7 yt 1   t
9.3 自回归模型的平稳性和相关函数
9.3.2 自回归模型的自相关函数
• AR(2)模型的自相关函数
yt  1 yt 1  2 yt  2   t
 0  E( yt2 )  E( yt (1 yt 1  2 yt  2   t ))  1 1  2 2   2
 1  E( yt yt 1 )  E((1 yt 1  2 yt  2   t ) yt 1 )  1 0  2 1
 2  E( yt yt  2 )  E((1 yt 1  2 yt  2   t ) yt  2 )  1 1  2 0
当k  2
时
 k  E( yt yt k )  E((1 yt 1  2 yt 2   t ) yt k )  1 k 1  2 k 2
9.3 自回归模型的平稳性和相关函数
9.3.2 自回归模型的自相关函数
• AR(2)模型的自相关函数
由上式可解得
(1  2 ) 2
0 
(1  2 )[(1  2 )2  12 ]
9.3 自回归模型的平稳性和相关函数
9.3.2 自回归模型的自相关函数
• AR(2)模型的自相关函数
结论2:AR(2)模型的自相关函数(ACF)为
1
12
 (1) 
,  (2)  2 
1  2
1  2
 (k )  1 (k  1)  2  (k  2), k  2
下面一行等式称为尤勒-沃尔克方程(YuleWalker equations)。
9.3 自回归模型的平稳性和相关函数
9.3.2 自回归模型的自相关函数
• AR(2)模型的自相关函数
例子 9.3
yt  0.9 yt 1  0.7 yt 2   t
9.3 自回归模型的平稳性和相关函数
9.3.2 自回归模型的自相关函数
• AR(2)模型的偏自相关函数
yt  11 yt 1   t
yt  21 yt 1  22 yt  2   t
yt  31 yt 1  32 yt  2  33 yt  3   t

yt  k1 yt 1  k 2 yt  2    kk yt  k   t
称 ii 为 i 阶自回归模型的偏自相关系数
(PAC:Partial Auto-Correlation) i  1,2,, k
9.3 自回归模型的平稳性和相关函数
9.3.2 自回归模型的自相关函数
• AR(2)模型的偏自相关函数
AR(k)模型的偏相关函数为
ii , i  k
 (i)  
,i  1,2,
0, i  k
 阶数大于k 时,偏自相关系数为0,这种现象称
为AR模型偏相关函数的截尾性。
yt
偏自相关系数ii 是剔除yt 1,, yt i 1 对 yt 的影响后,
和 yt  i 的相关系数。
9.4自回归模型的定阶和估计
9.4.1 自回归模型定阶
9.4.2 自回归模型估计
9.4.3 自回归模型再定阶—信息准则
9.4自回归模型的定阶和估计
9.4.1 自回归模型定阶
• 自回归模型的确立:确定阶数 估计 再
次确定阶数的循环
• 自相关函数用来确定采用自回归模型是否
合适。如果自相关函数具有拖尾性,则AR
模型为合适模型。
• 偏自相关函数用来确定模型的阶数。如果
从某个阶数之后,偏自相关函数的值都很
接近0,则取相应的阶数作为模型阶数。
9.4自回归模型的定阶和估计
9.4.1 自回归模型定阶
例子9.3 库存投资模型—定阶
打开包含库存投资变量Invent的工作文件,
在主菜单中点击Quick →Series Statistics →
Correlogram
9.4自回归模型的定阶和估计
9.4.1 自回归模型定阶
例子9.3 库存投资模型—定阶
在出现的对话框中输入序列(变量)名称,
点击OK按钮,弹出的对话框(Correlogram
Specification)中有对原数据(level),一
阶差分后的数据(1st difference),二阶差
分后的数据(2nd difference)的选择,以
及自回归包含多少滞后项(Lags to
include)。
9.4自回归模型的定阶和估计
9.4.1 自回归模型定阶
例子9.3 库存投资模型—定阶
9.4自回归模型的定阶和估计
9.4.1 自回归模型定阶
例子9.3 库存投资模型—定阶
点击OK后,将显示结果
9.4自回归模型的定阶和估计
9.4.1 自回归模型定阶
例子9.3 库存投资模型—定阶
看自相关(Autocorrelation),发现有拖尾
性,故可以选择自回归模型,看偏自相关
(Partial Correlation),发现有截尾现象:
高于四阶的偏自相关均为0,故可以建立4
阶自回归模型AR(4)。
9.4自回归模型的定阶和估计
9.4.2 自回归模型估计
自回归模型估计
• 最小二乘估计
t
AR(k)模型仍为线性模型,且误差项
满足基本第4章的假设1~假设4,故得出的
估计仍然就有一致性和马尔科夫性。
当样本量较大时,采用滞后变量导致
的回归样本减少对估计精度的影响不大。
自回归模型估计
• 极大似然估计
假设误差项服从正态分布,可以用极
大似然估计。
yt
故
~ N (c  1 yt 1    k yt  k , 2
y t 1 ,, y t  k
2 1 / 2 0.5 ( y t  c 1 y t 1  k y t  k ) 2 /( 2 2 )

f yt | yt 1 ,, yt k ( yt )  (2 )
e
L( y1,, yT )
 f y1 ( y1 )  f y2 | y1 ( y2 )  f y3 | y1 , y 2 ( y3 )   f yT | y1 ,, yT 1 ( yT | y1,, yT 1 )
 (2 )
2 T / 2


T
t  k 1
e
 0.5( y t  c 1 yt 1  k y t k ) 2 /(2 2 )
自回归模型估计
• 极大似然估计
对数似然函数为
l ( y1,, yT )
 c  0.5 ln( )  (2 )
2

2 1

2
(
y

c


y




y
)
t k 1 t
1 t 1
k t k
T
其最大化得出的估计与最小二乘估计一致。
自回归模型估计
• EViews操作
第一种方法:点击主菜单的Quick→
Estimate Equation,若按一般线性回归模型
进行设定,则输入
自回归模型估计
• EViews操作
输出结果
自回归模型估计
• EViews操作
若按自回归模型进行设定,则输入
注意:估计AR(4)时要将AR(1)到AR(4)全部写出。
自回归模型估计
• EViews操作
输出结果
自回归模型估计
• EViews操作
上述两种估计方法仅截距项的估计不一样,
原因是前者采用普通最小二乘法,后者采
用约束最小二乘法。在用EViews估计时间序
列模型时,应采用第二种设定方法。
自回归模型估计
• EViews操作
第二种方法:采用向量自回归VAR估计方
法。点击主菜单的Quick→Estimate VAR,点
选VAR Type中的Unrestricted VAR,在
Endogenous Variables栏中输入变量名invent,
在Lag Intervals for Endogenous栏中输入1 2 3
4,在Exogenous Variables栏输入c:
自回归模型估计
• EViews操作
自回归模型估计
• EViews操作
输出结果与普通最小二乘法一致
自回归模型估计
• 脉冲响应函数
脉冲响应是指 t 处新息  t 的变化对后续 y
的影响。
• 平稳AR模型可以转换为无穷阶的移动平均模型,
设转换后的MA模型为
yt  c   t  1 t 1  2 t 2     j t  j  
yt 的脉冲响应函数(IRF:impulse response
function)定义为:
IRF( j ) 
yt  j
 t
  j , j  1,2,
自回归模型估计
• 脉冲响应函数
例子9.3(续) 库存投资模型—脉冲
响应
分析
首先用VAR估计法估计AR模型,在输出
结果界面,点击View→Impulse Response,
弹出对话框
自回归模型估计
• 脉冲响应函数
例子9.3(续) 库存投资模型—脉冲
响应
分析
点选Display Format
下的Combined Graph,在
Impulse和Response中均填
入invent,在Periods(最
大滞后期)中输入数字,
点击确定输出结果如右:
9.4自回归模型的定阶和估计
9.4.3 自回归模型的再定阶—信息准则
前面自回归模型的定阶是通过计
算样本的自相关函数和偏自相关函数得出
的,存在一定的偏差。
信息准则也可以被用来确定阶数,常
用的有赤池信息准则(Aikaike info Criterion)
AIC和施瓦兹信息准则(Schwarz info
Criterion)SC,最优阶数使得信息准则值最
小。
9.4自回归模型的定阶和估计
9.4.3 自回归模型的再定阶—信息准则
例子9.4 信息准则定阶
分别计算AR(1)到AR(5)的信息准则值,然
后
再比较,选出使得信息准则值最小的阶数
AIC:-6.003716、-6.061554、-6.069740、-6.047398、-6.032800
SC :-5.864837、-5.893780、-5.872679、-5.820648、-5.775950
故AIC建议选择3阶滞后,SC建议选择2阶滞后。
9.5 自回归分布滞后模型
9.5.1 自回归分布滞后模型
9.5.2 格兰杰因果关系检验
9.5 自回归分布滞后模型
9.5.1 自回归分布滞后模型
ARDL(AutoRegression Distribution Lag)
yt  c  1 yt 1    k yt k   0 xt  1xt 1     p xt  p  t
适应性预期模型(消费与收入的关系)
yt    0 xt  1xt 1  xte  t
xte  xte1  ( xt  xte1)
yt    (0  ) xt  1xt 1   (  ) xte1   t
yt  c  0 xt  1xt 1  2 xt  2  yt 1  ut
9.5 自回归分布滞后模型
9.5.1 自回归分布滞后模型
ARDL(AutoRegression Distribution Lag)
yt  c  1 yt 1    k yt k   0 xt  1xt 1     p xt  p  t
部分调整模型(宏观经济变量)
yt  yt 1   ( y  yt 1)
*
t
y    xt   t
*
t
yt  c  1 yt 1  0 xt  ut
9.5 自回归分布滞后模型
9.5.1 自回归分布滞后模型
例子9.5 货币需求
弗里德曼货币需求函数:
Mt*  AYt Rt eu
*
m
取对数,t  a  yt  rt   t
*
m

m


(
m
货币需求部分调整 t
t 1
t  mt 1 )
货币需求的时间序列模型
mt   0   1 yt   2rt   3mt 1  vt
t
9.5 自回归分布滞后模型
9.5.2. 格兰杰因果关系检验
检验变量之间领先关系的方法:如果序列 x
的信息 xt 1, xt  2 ,, xt l 对 yt 有显著解释能力,
表明变量 x 是 y 的原因,这种关系称为格
兰杰因果关系。
yt  c  1 yt 1    k yt k   1xt 1     l xt l   t
检验 H0 :  1   2     l  0 (没有格兰杰因果
关系)
9.5 自回归分布滞后模型
9.5.2. 格兰杰因果关系检验
例子9.6
石油与经济
zt  a0  a1zt 1    am zt m  b1xt 1    bl xt l   t
zt  Pt 表示价格变化,
xt  ln(GDPt / GDPt 1 )
为国民生产总值变化的百分比。
9.5 自回归分布滞后模型
9.5.2. 格兰杰因果关系检验
例子9.6 石油与经济
在Eviews中将两个变量打开(Open→as
Group),在数据表格界面点击
View→Granger Causality,在弹出对话框Lag
Specification中输入需要加入的滞后阶数。
9.5 自回归分布滞后模型
9.5.2. 格兰杰因果关系检验
例子9.6 石油与经济
输出结果
 格兰杰因果只是统计意义上的因果关系
9.6 ARCH模型
9.6.1 ARCH模型的定义
9.6.2 ARCH模型估计
9.6 ARCH模型
9.6.1 ARCH模型的定义
pt
股票价格
随机游走,从而股 rt
2
r
r
票收益为白噪声,
之间没有关系,但
t
t
之间则可能存在关系,能为资产定价提供
信息。
9.6 ARCH模型
9.6.1 ARCH模型的定义
k
 定义1(ARCH模型):设时间序列 { yt , t  1,2,}
{ t , t  1,2,}
满足
阶自回归模型,误差项序列
2
q
{

为白噪声。如果误差项平方形成的序列
t , t  1,2,}
{ yt , t  1,2,}
ARCH(q)
服从
阶自回归模型,称时间序列
为带
yt  c  1 yt 1    k yt  k   t
误差项的自回归模型,表示为
 t2   0  1 t21     m t2 m  t
{ t , t  1,2,}
{t , t  1,2,}
其中
和
为相互独立的白噪声序列。上式中第一个方程为
均值方程,第二个方程称为方差方程或者波动方
9.6 ARCH模型
9.6.1 ARCH模型的定义
上述定义中中均值和方差模型也可写
yt  c  1 yt 1    k yt  k   t
作:
 t2   0  1 t21     m t2 m
{ t2 , t  1,2,}
若序列
的前后相关性持续时间
2
2
2









m
太长(方差聚集效应),则需要选取较大
t
0
1 t 1
1 t 1
的 ,而GARCH 模型
也
能描述持续的相关性,因此通常使用
9.6 ARCH模型
9.6.1 ARCH模型的定义
ARCH-m模型(风险溢价)
yt  c  1 yt 1    k yt  k   t2   t
 t2   0  1 t21  1 t21
方差还可以以  t 或 log( t2 ) 的形式引入
9.6 ARCH模型
9.6.2 ARCH模型估计
 t2 不可观测,故只能用极大似然估计
f y t | y t 1 ( yt )  ( 2 )
2 1 / 2
t
e
 ( y t  c 1 y t 1 ) 2 /( 2 t2 )
L( y1, y2 ,, yT )  f ( y1 )  f y 2 | y1 ( y2 )   f yT | yt 1 ,, y1 ( yT )
 (2 )
T / 2

T
t 1
2  ( y t  c 1 y t 1 ) 2 /(2 t2 )
t
 e
l ( y1 , y2 ,, yT )  c  t 1 ln( )  t 1 (2 t2 )1 ( yt  c  1 yt 1 )2
T
 t2   0  1 t21  1 t21
 t 1  yt 1  c  1 yt  2
2
t
T
9.6 ARCH模型
9.6.2 ARCH模型估计
2
2

给出初值 0 和  0 ,即可以估计参数。初值  02
2

和 0 可以采用如下两种方法得到:(1)用
OLS方法估计均值模型,回归残差作为对初
T
2
2
1
ˆ
ˆ
值的估计  0   0  T t 1ˆt2 (2)回代法
T
T
2
2
T 1
2
1
ˆ 0  ˆ0   T t 1ˆt  (1   )T t 1 T  j 1ˆT2 t , 为平
滑参数,通常取   0.7 。
• 误差项的分布还可以选择t-分布或者GED分
布。
9.6 ARCH模型
9.6.2 ARCH模型估计
例子9.6 沪深300指数(日收益率建模)
打开包含日收益率的EViews文件,
Quick Estimation equation,在Estimation
settings的
Methods中选择ARCH-Autoregressive
Conditional Heteroskedasticity,弹出对话框
9.6 ARCH模型
9.6.2 ARCH模型估计
例子9.6
均值模型设
定,此处为
白噪声
对估计方法
的细节设定
是否在均值
沪深300指数(日收益率建模)
模型中加入
波动率为解
释变量
选择
GARCH的
类型
GARCH的
具体回归
阶数
误差项的分
布假设
9.6 ARCH模型
9.6.2 ARCH模型估计
对含移动平均
项的均值模型
进行设定
设定回代法
 的值
9.6 ARCH模型
9.6.2 ARCH模型估计
例子9.7 β系数
重要概念
1. 时间序列分析的重点,是建立合适的模型刻画不同时间
点上随机变量的相关性。常用的时间序列模型有自回归模
型和移动平均模型,平稳的自回归模型可以转化为无穷阶
的移动平均模型。时间序列的相关性,用自相关系数和偏
自相关系数表示。
2. 平稳性是时间序列的数学期望、方差和协方差不随时间
变化。自回归模型描述的时间序列的平稳性,可以用模型
的滞后多项式根来判断。如果滞后多项式的根都落在单位
圆外,则时间序列平稳。
3. 阶数的自回归模型的自相关函数和偏自相关函数具有不
同的性质,据此可以判断对给定的样本数据,多少阶的自
回归模型是合适的。
4. 自回归模型的估计可以采用OLS估计和极大似然估计,两
种方法得出的自回归系数估计相同。自回归模型的估计和
定阶是交替进行的。除了采用自相关函数和偏自相关函数
初步确定阶数之外,还可以采用信息准则确定模型阶数。
重要概念
5. 为了研究两个时间序列之间的关系,在自回归模型自
变量中引入另一个时间序列变量及其滞后项,形成自回
归分布滞后模型。适应性预期模型和部分调整模型,是
自回归分布滞后模型在经济理论中的应用。对自回归分
布滞后模型进行参数约束检验,可以对经济时间序列之
间的领先关系进行格兰杰因果关系检验。
6. 用自回归模型对时间序列的条件方差进行建模,与均
值模型一起形成带ARCH误差项的模型。GARCH模型不仅
能充分刻画条件方差相关性,还具有更为简洁的形式。
GARCH(1,1)模型应用最为广泛,其一般形式为
yt  c  1 yt 1    k yt  k  
 t2   0  1 t21  1 t21
GARCH模型采用极大似然估计方法进行估计,为保证
方差的非负性和方差模型的平稳性,需要对模型系数施
加复杂的约束。