第9章平稳时间序分析ppt课件

1、第第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

2、.5.2 格兰杰因果关系检验9.6 arch模型模型 9.6.1 arch模型的定义 9.6.2 arch模型估计重要概念重要概念9.1 时间序列的概念时间序列的概念 设时点 处的观测为随机变量 ，这些随机变量形成一个时间序列，记为 或者 ， 的一组具体取值称为时间序列的实现值（realization）。 自相关函数（acf:autocorrelation function ）tt, 2 , 1 tyyy,21ttty1, 2 , 1,ttyttyyy,21stkkssyyyytststs其中),()(var)(var),cov(),(9.1 时间序列的概念时间序列的概念定义1（平稳性）：如果

3、时间序列 的数学期望、方差和协方差不随时间变化，即 称 为宽平稳（wide-sense stationary）时间序列。宽平稳也称为协方差平稳或者二阶矩平稳。, 2 , 1,ttyttkttkcyyyyktttt, 2 , 1;, 2 , 1),(),(cov,)(var,)(e2ttty19.1 时间序列的概念时间序列的概念 严平稳：时间序列中任意一组随机变量的联合分布不随时间发生变化，即对任意一组时间点 和时间间隔 ， 的联合分布与 的联合分布相同，称 严平稳。 二阶矩存在的严平稳时间序列一定宽平稳，宽平稳的时间序列不一定严平稳，本书只讨论宽平稳，将宽平稳时间序列简称为平稳时间序列。ntt

4、t21s,21ntttyyy,21snsstttyyyttty19.1 时间序列的概念时间序列的概念 若 为平稳时间序列，则：（1）（2） 满足大数定律，因此 分别是 、 和 的一致估计。ttty1)()(),cov(),()(22kkcyyksskksstttttkttttttttyyyyyykyyyytkttt1211221)()()()(,11111ttty12)(k9.1 时间序列的概念时间序列的概念 表示 的 阶滞后，用滞后算符 表示为 例如 用滞后算符多项式表示为：ktytykl, 2 , 1 , 0,kyylkttkttttyyyy2352123tylll) 2352(239.2

5、 时间序列模型时间序列模型9.2.1 白噪声序列9.2.2 自回归模型9.2.3 移动平均模型9.2.4 自回归模型转化为移动平均模型9.2 时间序列模型时间序列模型9.2.1 白噪声序列定义2（白噪声）：如果时间序列 满足：（1） ， （2）对任意 ， 和 不相关，即 称 为白噪声序列，简称白噪声 (white noise)。 是平稳时间序列的极端例子。, 1,ttt0)(et2)(varttsts0)(est, 1,ttt9.2 时间序列模型时间序列模型9.2.2 自回归模型一阶自回归模型ar(1) ， 为白噪声除了常数项以外， 在 时刻的值由前定项 （predetermined term

6、） 和与前期值不相关的新息（innovation） 组成。 阶自回归模型ar(k), 0(,211nycyttttyt11tyt), 0 (,.22211nyyycyttktpttt1|11|.|21kkt9.2 时间序列模型时间序列模型9.2.2 自回归模型 模型 ， 150个样本的时序图：tttyy17 . 05 . 1) 1 , 0( nt9.2 时间序列模型时间序列模型9.2.3 移动平均模型对一阶自回归模型进行递推：当 时， ，ktkktkttkttttttttycyccyccycy1111111211112211121111)1 ()(k01kktktttcy1111119.2 时

7、间序列模型时间序列模型9.2.3 移动平均模型取有限项，上式即为 阶移动平均模型。ktktttcy111k9.2 时间序列模型时间序列模型9.2.4 自回归模型转化为移动平均模型 由于 在上述一阶自回归模型两边同乘tttcylyl)1 ()(11)1)1 (12211lll)(1l22111122111111)1 (1)1 (ttttttcllcyly9.2 时间序列模型时间序列模型9.2.4 自回归模型转化为移动平均模型 可以转化为移动平均模型的自回归模型称为可逆的（invertible）。 从上面推导可以看出，一阶自回归模型可逆的条件是 。实际上，自回归模型的可逆条件，是滞后多项式的根在单

8、位圆外。滞后多项式即：1|1)1 ()(221kkllll 9.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.1 自回归模型的平稳性9.3.2 自回归模型的自相关函数9.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.1 自回归模型的平稳性 阶自回归模型： 且 为白噪声序列 用滞后算子表达上式为：k), 0(211nyycyttktkttttttkkcylyll)()1 (19.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.1 自回归模型的平稳性u结论1：自回归模型平稳的充分必要条件为：滞后多项式的根都在单位圆之外，即方程 的根

9、满足 。其中 为实根时 表示绝对值， 为虚根时 表示虚数的模。0)( lll1|l|ll|l9.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.1 自回归模型的平稳性 滞后多项式的单位根均在单位圆内，则时间序列平稳；若有根为1，则不平稳，此时称存在单位根。 把模型是否平稳的检验称为单位根检验。9.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.1 自回归模型的平稳性例子9.1 为平稳时间序列 为非平稳时间序列 tttyy17 . 05 . 1ttttyyy219.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.2 自回归模型的自相关函

10、数 对 阶自回归模型： 且 为白噪声序列 若 为平稳时间序列，则 两边取期望，得 代入原模型，整理可得零均值化的 阶自回归模型：k), 0(211nyycyttktkttttykttyyee)1/(1kctktkttyyy)()(11k9.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.2 自回归模型的自相关函数 重新将记为 ，即用 表示原模型中的 。此时， 。后面的自回归模型都将采用这种零均值化后的模型 。 tktkttyyy)()(11tktkttyyy11tyty0ety9.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.2 自回归模型的自相关函数

11、 ar(1)模型的自相关函数 用 表示变量的方差 得出 表示 与 的协方差 020212212121120)(e)(e)(e)(etttttyyy21201)(varty212101111111)(e)(covtttttyyyy1ty1ty9.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.2 自回归模型的自相关函数 ar(1)模型的自相关函数 同理， 则ar(1)的自相关函数为 。, 2 , 1,1),(cov212101kyykkkttk0/)(kk 9.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.2 自回归模型的自相关函数u结论2：ar(1)模

12、型的自相关函数（acf）为 平稳性要求 ，当 时， ，即自相关系数随时间间隔增加指数递减到0，但不等于0。这种现象称为自回归模型自相关函数的拖尾性。, 2 , 1,)(1kkk1|1k0)(k9.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.2 自回归模型的自相关函数例子9.2tttyy17 . 05 . 19.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.2 自回归模型的自相关函数 ar(2)模型的自相关函数 当 时ttttyyy2211021122211221201122111122211221120)(e)(e)(e)(e)(e)(ettttt

13、ttttttttttttyyyyyyyyyyyyyy22112211)()(kkkttttkttkyyyeyye2k9.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.2 自回归模型的自相关函数 ar(2)模型的自相关函数 由上式可解得 )1)(1 ()1 (212222209.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.2 自回归模型的自相关函数 ar(2)模型的自相关函数u结论2：ar(2)模型的自相关函数（acf）为 下面一行等式称为尤勒-沃尔克方程（yule-walker equations）。2),2() 1()(1)2(,1) 1 (21

14、221221kkkk9.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.2 自回归模型的自相关函数 ar(2)模型的自相关函数例子 9.3ttttyyy217 . 09 . 09.3 自回归模型的平稳性和相关函数自回归模型的平稳性和相关函数9.3.2 自回归模型的自相关函数 ar(2)模型的偏自相关函数 称 为 阶自回归模型的偏自相关系数（pac：partial auto-correlation）tktkktktktttttttttttttyyyyyyyyyyyyy2211333232131222121111iiiki,2, 19.3 自回归模型的平稳性和相关函数自回归模型

15、的平稳性和相关函数9.3.2 自回归模型的自相关函数 ar(2)模型的偏自相关函数 ar(k)模型的偏相关函数为u 阶数大于 时，偏自相关系数为0，这种现象称为ar模型偏相关函数的截尾性。u偏自相关系数 是剔除 对 的影响后， 和 的相关系数。, 2 , 1, 0,)(ikikiiii，ii11,ittyytytyityk9.4自回归模型的定阶和估计自回归模型的定阶和估计9.4.1 自回归模型定阶9.4.2 自回归模型估计9.4.3 自回归模型再定阶信息准则9.4自回归模型的定阶和估计自回归模型的定阶和估计9.4.1 自回归模型定阶 自回归模型的确立：确定阶数 估计 再次确定阶数的循环 自相关

16、函数用来确定采用自回归模型是否合适。如果自相关函数具有拖尾性，则ar模型为合适模型。 偏自相关函数用来确定模型的阶数。如果从某个阶数之后，偏自相关函数的值都很接近0，则取相应的阶数作为模型阶数。9.4自回归模型的定阶和估计自回归模型的定阶和估计9.4.1 自回归模型定阶例子9.3 库存投资模型定阶 打开包含库存投资变量invent的工作文件，在主菜单中点击quick series statistics correlogram9.4自回归模型的定阶和估计自回归模型的定阶和估计9.4.1 自回归模型定阶例子9.3 库存投资模型定阶 在出现的对话框中输入序列（变量）名称，点击ok按钮，弹出的对话框（

17、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 库存投资模型定阶 看自相关（au

18、tocorrelation），发现有拖尾性，故可以选择自回归模型，看偏自相关（partial correlation），发现有截尾现象：高于四阶的偏自相关均为0，故可以建立4阶自回归模型ar(4)。9.4自回归模型的定阶和估计自回归模型的定阶和估计9.4.2 自回归模型估计自回归模型估计 最小二乘估计 ar(k)模型仍为线性模型，且误差项 满足基本第4章的假设1假设4，故得出的估计仍然就有一致性和马尔科夫性。 当样本量较大时，采用滞后变量导致的回归样本减少对估计精度的影响不大。t自回归模型估计 极大似然估计 假设误差项服从正态分布，可以用极大似然估计。 故211,(1ktktyytyycnyk

19、tt)2/()(5 . 02/12,|22111)2()(ktkttktttyycytyyyeyftktyycytttyyyyyyyyytktktttteyyyfyfyfyfyyl1)2/()( 5 . 02/211,|3,|2|11221111213121)2(),|()()()(),(自回归模型估计 极大似然估计 对数似然函数为 其最大化得出的估计与最小二乘估计一致。 tktktktttyycycyyl12111221)()2()ln(5 . 0),(自回归模型估计 eviews操作 第一种方法：点击主菜单的quick estimate equation，若按一般线性回归模型进行设定，则输

20、入 自回归模型估计 eviews操作 输出结果自回归模型估计 eviews操作 若按自回归模型进行设定，则输入注意：估计ar(4)时要将ar(1)到ar(4)全部写出。自回归模型估计 eviews操作 输出结果自回归模型估计 eviews操作 上述两种估计方法仅截距项的估计不一样，原因是前者采用普通最小二乘法，后者采用约束最小二乘法。在用eviews估计时间序列模型时，应采用第二种设定方法。自回归模型估计 eviews操作 第二种方法：采用向量自回归var估计方法。点击主菜单的quickestimate var，点选var type中的unrestricted var，在endogenous

21、variables栏中输入变量名invent，在lag intervals for endogenous栏中输入1 2 3 4，在exogenous variables栏输入c：自回归模型估计 eviews操作 自回归模型估计 eviews操作 输出结果与普通最小二乘法一致自回归模型估计 脉冲响应函数 脉冲响应是指 处新息 的变化对后续 的影响。 平稳ar模型可以转换为无穷阶的移动平均模型，设转换后的ma模型为 的脉冲响应函数（irf：impulse response function）定义为:ttyjtjttttcy2211ty, 2 , 1,)irf(jyjjtjt自回归模型估计 脉冲响应

22、函数 例子9.3（续） 库存投资模型脉冲响应 分析 首先用var估计法估计ar模型，在输出结果界面，点击viewimpulse response，弹出对话框自回归模型估计 脉冲响应函数 例子9.3（续） 库存投资模型脉冲响应 分析 点选display format下的combined graph，在impulse和response中均填入invent，在periods（最大滞后期）中输入数字，点击确定输出结果如右：9.4自回归模型的定阶和估计自回归模型的定阶和估计9.4.3 自回归模型的再定阶信息准则 前面自回归模型的定阶是通过计算样本的自相关函数和偏自相关函数得出的，存在一定的偏差。 信息准

23、则也可以被用来确定阶数，常用的有赤池信息准则（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.032800sc ：-5.864837、-5.893780、-5.872679、-5.820648、-5.77

24、5950故aic建议选择3阶滞后，sc建议选择2阶滞后。9.5 自回归分布滞后模型自回归分布滞后模型9.5.1 自回归分布滞后模型9.5.2 格兰杰因果关系检验9.5 自回归分布滞后模型自回归分布滞后模型9.5.1 自回归分布滞后模型 ardl(autoregression distribution lag) 适应性预期模型（消费与收入的关系） tptpttktkttxxxyycy11011tettttxxxy110)(11ettetetxxxxtettttxxxy1110)()(ttttttuyxxxcy1221109.5 自回归分布滞后模型自回归分布滞后模型9.5.1 自回归分布滞后模型

25、ardl(autoregression distribution lag) 部分调整模型（宏观经济变量）tptpttktkttxxxyycy11011)(1*1ttttyyyytttxy*ttttuxycy0119.5 自回归分布滞后模型自回归分布滞后模型9.5.1 自回归分布滞后模型例子9.5 货币需求 弗里德曼货币需求函数： 取对数， 货币需求部分调整 货币需求的时间序列模型tuttteraym*ttttryam*)(1*1ttttmmmmtttttvmrym132109.5 自回归分布滞后模型自回归分布滞后模型9.5.2. 格兰杰因果关系检验 检验变量之间领先关系的方法:如果序列 的信息

26、 对 有显著解释能力，表明变量 是 的原因，这种关系称为格兰杰因果关系。 检验 （没有格兰杰因果关系）xltttxxx,21tyxytltltktkttxxyycy11110:h210l9.5 自回归分布滞后模型自回归分布滞后模型9.5.2. 格兰杰因果关系检验例子9.6 石油与经济 表示价格变化， 为国民生产总值变化的百分比。tltltmtmttxbxbzazaaz11110ttpz)gdp/gdpln(1tttx9.5 自回归分布滞后模型自回归分布滞后模型9.5.2. 格兰杰因果关系检验例子9.6 石油与经济 在eviews中将两个变量打开（openas group），在数据表格界面点击v

27、iewgranger 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模型的定义 股票价格 随机游走，从而股票收益为白噪声， 之间没有关系，但 之间则可能存在关系，能为资产定价提供信息。tptrtr2tr9.6 arch模型模型9.6.1 arch模型的定义定义1(arch模型)：设

28、时间序列 满足 阶自回归模型，误差项序列 为白噪声。如果误差项平方形成的序列 服从 阶自回归模型，称时间序列 为带 误差项的自回归模型，表示为 其中 和 为相互独立的白噪声序列。上式中第一个方程为均值方程，第二个方程称为方差方程或者波动方程。, 2 , 1,tytk, 2 , 1,tt, 2 , 1,2ttq, 2 , 1,tyttmtmtttktkttyycy22110211, 2 , 1,tt, 2 , 1,tt)(qarch9.6 arch模型模型9.6.1 arch模型的定义 上述定义中中均值和方差模型也可写作： 若序列 的前后相关性持续时间太长（方差聚集效应），则需要选取较大的 ，而

29、garch 模型 也能描述持续的相关性，因此通常使用garch(1,1)来描述条件方差的行为。 22110211mtmtttktkttyycy, 2 , 1,2ttm21121102ttt9.6 arch模型模型9.6.1 arch模型的定义 arch-m模型（风险溢价） 方差还可以以 或 的形式引入 21121102211tttttktkttyycyt)log(2t9.6 arch模型模型9.6.2 arch模型估计 不可观测，故只能用极大似然估计2t)2/()(2/12|22111)2()(tttttycyttyyeyf)2/()(122/,|2|12122111112)2()()()()

30、,(tttttycytttttyyyyyteyfyfyfyyyl2111211211022111121221)()2()ln(),(tttttttttttttttycyycycyyyl9.6 arch模型模型9.6.2 arch模型估计 给出初值 和 ，即可以估计参数。初值 和 可以采用如下两种方法得到：（1）用ols方法估计均值模型，回归残差作为对初值的估计 （2）回代法 ， 为平滑参数，通常取 。 误差项的分布还可以选择t-分布或者ged分布。20202020tttt1212020ttttjttttttt12111212020)1 (7 . 09.6 arch模型模型9.6.2 arch模型估计例子9.6 沪深300指数（日收益率建模） 打开包含日收益率的eviews文件，quick estimat