拒绝访问 |
| 百度云加速
请打开cookies.
此网站 () 的管理员禁止了您的访问。原因是您的访问包含了非浏览器特征(3388-ua98).
重新安装浏览器,或使用别的浏览器王燕-应用时间序列分析
中国人民大学出版社 中国人民大学音像出版社《应用时间序列分析》目? ? ? ? ? ?录第一章 第二章 第三章 第四章 第五章 第六章时间序列分析简介 时间序列的预处理 平稳时间序列分析 非平稳序列的确定性分析 非平稳序列的随机分析 多元时间序列分析第一章时间序列分析简介本章结构? ? ? ?引言 时间序列的定义 时间序列分析方法简介 时间序列分析软件1.1 引言?最早的时间序列分析可以追溯到7000年前的古 埃及。?古埃及人把尼罗河涨落的情况逐天记录下来,就构 成所谓的时间序列。对这个时间序列长期的观察使 他们发现尼罗河的涨落非常有规律。由于掌握了尼 罗河泛滥的规律,使得古埃及的农业迅速发展,从 而创建了埃及灿烂的史前文明。?按照时间的顺序把随机事件变化发展的过程记 录下来就构成了一个时间序列。对时间序列进 行观察、研究,找寻它变化发展的规律,预测 它将来的走势就是时间序列分析。1.2 时间序列的定义?随机序列:按时间顺序排列的一组随机变量 观察值序列:随机序列的 n 个有序观察值,称之 为序列长度为 n 的观察值序列 x1 , x 2 , ? , xt 随机序列和观察值序列的关系? ? ?? , X 1 , X 2 ,? , X t ,???观察值序列是随机序列的一个实现 我们研究的目的是想揭示随机时序的性质 实现的手段都是通过观察值序列的性质进行推断1.3 时间序列分析方法?描述性时序分析 统计时序分析?描述性时序分析??通过直观的数据比较或绘图观测,寻找 序列中蕴含的发展规律,这种分析方法 就称为描述性时序分析 描述性时序分析方法具有操作简单、直 观有效的特点,它通常是人们进行统计 时序分析的第一步。描述性时序分析案例?德国业余天文学家施瓦尔发现太阳黑子的活动具有11年左右的周期统计时序分析?频域分析方法 时域分析方法?频域分析方法?原理?假设任何一种无趋势的时间序列都可以分解成若干不同频率 的周期波动 早期的频域分析方法借助富里埃分析从频率的角度揭示时间 序列的规律 后来借助了傅里叶变换,用正弦、余弦项之和来逼近某个函 数 20世纪60年代,引入最大熵谱估计理论,进入现代谱分析阶 段 非常有用的动态数据分析方法,但是由于分析方法复杂,结 果抽象,有一定的使用局限性?发展过程? ???特点?时域分析方法?原理?事件的发展通常都具有一定的惯性,这种惯性用统 计的语言来描述就是序列值之间存在着一定的相关 关系,这种相关关系通常具有某种统计规律。 寻找出序列值之间相关关系的统计规律,并拟合出 适当的数学模型来描述这种规律,进而利用这个拟 合模型预测序列未来的走势 理论基础扎实,操作步骤规范,分析结果易于解 释,是时间序列分析的主流方法?目的??特点?时域分析方法的分析步骤? ? ? ? ?考察观察值序列的特征 根据序列的特征选择适当的拟合模型 根据序列的观察数据确定模型的口径 检验模型,优化模型 利用拟合好的模型来推断序列其它的统 计性质或预测序列将来的发展时域分析方法的发展过程? ? ?基础阶段 核心阶段 完善阶段基础阶段?G.U.Yule?1927年,AR模型 1931年,MA模型,ARMA模型?G.T.Walker?核心阶段?G.E.P.Box和 G.M.Jenkins?? ?1970年,出版《Time Series Analysis Forecasting and Control》 提出ARIMA模型(Box―Jenkins 模型) Box―Jenkins模型实际上是主要运用于单变 量、同方差场合的线性模型完善阶段???异方差场合 ? Robert F.Engle,1982年,ARCH模型 ? Bollerslov,1985年GARCH模型 多变量场合 ? C.Granger ,1987年,提出了协整(cointegration)理论 非线性场合 ? 汤家豪等,1980年,门限自回归模型1.4 时间序列分析软件?常用软件?S-plus,Matlab,Gauss,TSP,Eviews 和SAS 在SAS系统中有一个专门进行计量经济与时间序列 分析的模块:SAS/ETS。SAS/ETS编程语言简洁, 输出功能强大,分析结果精确,是进行时间序列分 析与预测的理想的软件 由于SAS系统具有全球一流的数据仓库功能,因此 在进行海量数据的时间序列分析时它具有其它统计 软件无可比拟的优势?推荐软件――SAS??第二章时间序列的预处理本章结构? ?平稳性检验 纯随机性检验2.1平稳性检验? ? ? ? ?特征统计量 平稳时间序列的定义 平稳时间序列的统计性质 平稳时间序列的意义 平稳性的检验概率分布?概率分布的意义?随机变量族的统计特性完全由它们的联合分布函数 或联合密度函数决定?时间序列概率分布族的定义{Ft1,t2 ,? ,tm (x1, x2,? , xm)} ?m?(1,2,? , m),?t1, t2,? ,tm ?T?实际应用的局限性特征统计量?均值 方差 自协方差 自相关系数? t ? EX t ? ? xdFt ( x)????DX t ? E ( X t ? ? t ) ? ? ( x ? ? t ) dFt ( x )2 ???2?? (t , s ) ? E ( X t ? ? t )( X s ? ? s )? (t , s) ? ? (t , s)DX t ? DX s?平稳时间序列的定义?严平稳?严平稳是一种条件比较苛刻的平稳性定义,它认为 只有当序列所有的统计性质都不会随着时间的推移 而发生变化时,该序列才能被认为平稳。 宽平稳是使用序列的特征统计量来定义的一种平稳 性。它认为序列的统计性质主要由它的低阶矩决 定,所以只要保证序列低阶矩平稳(二阶),就能 保证序列的主要性质近似稳定。?宽平稳?平稳时间序列的统计定义?满足如下条件的序列称为严平稳序列?正整数m, ?t1 , t 2 , ? , t m ? T, ?正整数?, 有Ft1 ,t 2 ? t m ( x1 , x 2 , ? , x m ) ? Ft1?? ,t 2 ?? ? t m ?? ( x1 , x 2 , ? , x m )?满足如下条件的序列称为宽平稳序列1) EX t2 ? ?, ?t ? T 2) EX t ? ? , ?为常数, ?t ? T 3) ? (t , s ) ? ? ( k , k ? s ? t ),?t , s, k且k ? s ? t ? T严平稳与宽平稳的关系?一般关系?严平稳条件比宽平稳条件苛刻,通常情况下,严平 稳(低阶矩存在)能推出宽平稳成立,而宽平稳序 列不能反推严平稳成立 不存在低阶矩的严平稳序列不满足宽平稳条件,例 如服从柯西分布的严平稳序列就不是宽平稳序列 当序列服从多元正态分布时,宽平稳可以推出严平 稳?特例??平稳时间序列的统计性质? ?常数均值 自协方差函数和自相关函数只依赖于时 间的平移长度而与时间的起止点无关?延迟k自协方差函数? (k ) ? ? (t , t ? k ), ?k为整数?延迟k自相关系数? k ?? ? (k ) ? (0)自相关系数的性质? ? ? ?规范性 对称性 非负定性 非唯一性平稳时间序列的意义?时间序列数据结构的特殊性?可列多个随机变量,而每个变量只有一个样 本观察值 极大地减少了随机变量的个数,并增加了待 估变量的样本容量 极大地简化了时序分析的难度,同时也提高 了对特征统计量的估计精度?平稳性的重大意义??平稳性的检验(图检验方法)?时序图检验?根据平稳时间序列均值、方差为常数的性 质,平稳序列的时序图应该显示出该序列始 终在一个常数值附近随机波动,而且波动的 范围有界、无明显趋势及周期特征 平稳序列通常具有短期相关性。该性质用自 相关系数来描述就是随着延迟期数的增加, 平稳序列的自相关系数会很快地衰减向零?自相关图检验?例题?例2.1?检验1964年――1999年中国纱年产量序列的平稳性 检验1962年1月――1975年12月平均每头奶牛月产 奶量序列的平稳性 检验1949年――1998年北京市每年最高气温序列的 平稳性?例2.2??例2.3?例2.1时序图例2.1自相关图例2.2时序图例2.2 自相关图例2.3时序图例2.3自相关图2.2 纯随机性检验? ? ?纯随机序列的定义 纯随机性的性质 纯随机性检验纯随机序列的定义?纯随机序列也称为白噪声序列,它满足 如下两条性质(1) EX t ? ? , ?t ? T ?? 2 , t ? s (2)? (t , s ) ? ? , ?t , s ? T ? 0, t ? s标准正态白噪声序列时序图白噪声序列的性质?纯随机性??(k) ? 0,?k ? 0各序列值之间没有任何相关关系,即为 “没有记忆” 的序列?方差齐性DX t ? ? (0) ? ? 2?根据马尔可夫定理,只有方差齐性假定成立时,用 最小二乘法得到的未知参数估计值才是准确的、有 效的纯随机性检验? ? ? ?检验原理 假设条件 检验统计量 判别原则Barlett定理?如果一个时间序列是纯随机的,得到一 个观察期数为n 的观察序列,那么该序列 的延迟非零期的样本自相关系数将近似 服从均值为零,方差为序列观察期数倒 数的正态分布1 ?k ~ ? ?N (0, ) n , ?k ? 0假设条件?原假设:延迟期数小于或等于 m 期的序列 值之间相互独立H 0:?1 ? ? 2 ? ? ? ? m ? 0, ?m ? 1?备择假设:延迟期数小于或等于 m 期的 序列值之间有相关性H 1:至少存在某个? k ? 0, ?m ? 1,k ? m检验统计量?Q统计量? k2 ~ ? 2 (m) Q ? n? ?k ?1 m?LB统计量LB ? n(n ? 2)? (k ?1 m? k2 ? n?k) ~ ? 2 ( m)判别原则?拒绝原假设?当检验统计量大于 ?12?? (m)分位点,或该统计 量的P值小于 ? 时,则可以以 1 ? ? 的置信水 平拒绝原假设,认为该序列为非白噪声序列 量的P值大于 ? 时,则认为在 1 ? ? 的置信水 平下无法拒绝原假设,即不能显著拒绝序列 为纯随机序列的假定2 ? 1 当检验统计量小于 ?? (m)分位点,或该统计?接受原假设?例2.4:标准正态白噪声序列纯随机性检验样本自相关图检验结果延迟 延迟6期 延迟12期QLB统计量检验QLB 统计量值P值 0.42.36 5.35由于P值显著大于显著性水平? ,所以该序列不能拒 绝纯随机的原假设。例2.5?对1950年――1998年北京市城乡居民定 期储蓄所占比例序列的平稳性与纯随机 性进行检验例2.5时序图例2.5自相关图例2.5白噪声检验结果LB统计量检验 延迟阶数 6 12 LB检验统计 量的值 75.46 82.57 P值 &0.1第三章平稳时间序列分析本章结构? ? ? ?方法性工具 ARMA模型 平稳序列建模 序列预测3.1 方法性工具? ? ?差分运算 延迟算子 线性差分方程差分运算?一阶差分?xt ? xt ? xt ?1?p 阶差分? p xt ? ? p ?1 xt ? ? p ?1 xt ?1?k 步差分? k ? xt ? xt ? k延迟算子??延迟算子类似于一个时间指针,当前序 列值乘以一个延迟算子,就相当于把当 前序列值的时间向过去拨了一个时刻 记B为延迟算子,有xt ? p ? B p xt , ?p ? 1延迟算子的性质? ? ? ?B0 ? 1B (c ? xt ) ? c ? B( xt ) ? c ? xt ?1 , c为任意常数B ( xt ? y t ) ? xt ?1 ? y t ?1B n xt ? xt ? nn n i ?0i n i i ? (1 ? B ) ? ? ( ?1) C n B ,其中 C n ?n! i!(n ? i )!用延迟算子表示差分运算?p 阶差分? p xt ? (1 ? B) p xt ? ? (?1) p C ip xt ?ii ?0 p?k 步差分? k ? xt ? xt ? k ? (1 ? B ) xtk线性差分方程?线性差分方程z t ? a1 z t ?1 ? a 2 z t ? 2 ? ? ? a p z t ? p ? h(t )?齐次线性差分方程z t ? a1 z t ?1 ? a 2 z t ? 2 ? ? ? a p z t ? p ? 0齐次线性差分方程的解?特征方程 p ? ? a1? p ?1 ? a2 ? p ?2 ? ? ? a p ? 0 特征方程的根称为特征根,记作?1 , ? 2 ,? , ? p 齐次线性差分方程的通解?? ??不相等实数根场合 t z t ? c1?1 ? c 2 ?t2 ? ? ? c p ?tp 有相等实根场合 t z t ? (c1 ? c2 t ? ? ? c d t d ?1 )?1 ? cd ?1?td ?1 ? ? ? c p ?tp 复根场合?zt ? r t (c1e it? ? c2 e ? it? ) ? c3?t3 ? ? ? c p ?tp非齐次线性差分方程的解?非齐次线性差分方程的特解?使得非齐次线性差分方程成立的任意一个解z t??z t?? ? a1 z t???1 ? a 2 z t??? 2 ? ? ? a p z t??? p ? h(t )?非齐次线性差分方程的通解?齐次线性差分方程的通解和非齐次线性差分方程的 特解之和 ztz t ? z t? ? z t??3.2 ARMA模型的性质? ? ?AR模型(Auto Regression Model) MA模型(Moving Average Model) ARMA模型(Auto Regression Moving Average model)AR模型的定义?具有如下结构的模型称为 p 阶自回归模 型,简记为 AR ( p ) ? xt ? ? 0 ? ?1 xt ?1 ? ? 2 xt ? 2 ? ? ? ? p xt ? p ? ? t?? ?? p ? 0 ? 2 ? E (? t ) ? 0,Var (? t ) ? ? ? , E (? t ? s ) ? 0, s ? t ? Ex ? ? 0, ?s ? t ? s t 特别当? 0 ? 0 时,称为中心化 AR ( p ) 模型AR(P)序列中心化变换?称 { y t }为 {xt } 的中心化序列 ,令???01 ? ?1 ? ? ? ? py t ? xt ? ?自回归系数多项式?引进延迟算子,中心化 AR ( p ) 模型又可 以简记为? ( B ) xt ? ? t?自回归系数多项式?( B) ? 1 ? ?1 B ? ? 2 B ? ? ? ? p B2pAR模型平稳性判别?判别原因?AR模型是常用的平稳序列的拟合模型之一, 但并非所有的AR模型都是平稳的 单位根判别法 平稳域判别法?判别方法? ?例3.1:考察如下四个模型的平稳性(1) xt ? 0.8 xt ?1 ? ? t(2) xt ? ?1.1xt ?1 ? ? t(3) xt ? xt ?1 ? 0.5 xt ? 2 ? ? t(4) xt ? xt ?1 ? 0.5 xt ?1 ? ? t例3.1平稳序列时序图(1) xt ? 0.8 xt ?1 ? ? t(3) xt ? xt ?1 ? 0.5 xt ? 2 ? ? t例3.1非平稳序列时序图(2) xt ? ?1.1xt ?1 ? ? t(4) xt ? xt ?1 ? 0.5 xt ?1 ? ? tAR模型平稳性判别方法?特征根判别??AR(p)模型平稳的充要条件是它的p个特征根都在单 位圆内 根据特征根和自回归系数多项式的根成倒数的性 质,等价判别条件是该模型的自回归系数多项式的 根都在单位圆外?平稳域判别?平稳域 {?1 ,?2 ,? ,? p 单位根都在单位圆内 }AR(1)模型平稳条件?特征根? ???平稳域? 〈1AR(2)模型平稳条件?特征根?1 ? ?12 ? 4?22?平稳域?1 ??2 ??1 ? ?12 ? 4?22{?1 , ? 2 ? 2 ? 1,且? 2 ? ?1 ? 1}例3.1平稳性判别模 型(1) (2) (3) (4)?1 ?特征根判别?1 ? 0.8?1 ? ?1.11? i 2平稳域判别结 论平稳 非 平稳 平稳 非 平稳? ? 0.8? ? ?1.1?2 ? 0.5, ?2 ? ?1 ? 0.5, ?2 ? ?1 ? ?1.51? i ?2 ? 21? 3 ? ? 1? 3 ?1 ? 2 2 2?2 ? 0.5, ?2 ? ?1 ? 1.5, ?2 ? ?1 ? ?0.5平稳AR模型的统计性质? ? ? ? ?均值 方差 协方差 自相关系数 偏自相关系数均值?如果AR(p)模型满足平稳性条件,则有Ext ? E (?0 ? ?1 xt ?1 ? ? ? ? p xt ? p ? ? t )?根据平稳序列均值为常数,且 {? t }为白噪声序 列,有 Ext ? ? , E (? t ) ? 0 , ?t ? T 推导出?? ?0?1 ? ?1 ? ? ? ? pGreen函数定义?AR模型的传递形式p ? ki xt ? ?? ? t ? ?? ki (?i B) j ? t ? ( B ) i ?1 1 ? ?i B i ?1 j ? 0 p?t? ?? ki ?ij ? t ? j ? ? G j ? t ? jj ? 0 i ?1 j ?0?p??其中系数 {G j , j ? 1,2,? } 称为Green函数Green函数递推公式?原理?? ( B) xt ? ? t ? ? ( B)G ( B)? t ? ? t ? ? xt ? G ( B)? t?方法?待定系数法?递推公式?G0 ? 1 ??k , k ? p ? j , 其中?k? ? ? ? ?0, k ? p ?G j ? ? ?k?G j ? k,j ? 1,2,? k ?1 ?方差?平稳AR模型的传递形式xt ? ? G j ? t ? jj ?0 ??两边求方差得2 Var ( xt ) ? ? G 2 ? j ? , G j为Green函数 j ?0 ?例3.2:求平稳AR(1)模型的方差?平稳AR(1)模型的传递形式为? ? ?t i i xt ? ? ? (?1 B) ? t ? ? ?1 ? t ?i 1 ? ?1 B i ?0 i ?0?Green函数为G j ? ?1 , j ? 0,1, ?j?平稳AR(1)模型的方差2 ? 2j 2 2 Var ( xt ) ? ? G j Var (? t ) ? ? ?1 ? ? ? ? 2 1 ? ?1 j ?0 j ?0 ? ?协方差函数?在平稳AR(p)模型两边同乘 xt ? k , ?k ? 1,再求期望E(xt xt?k ) ? ?1E(xt?1xt?k ) ?? ??p E(xt? p xt?k ) ? E(?t xt?k )?根据E (? t xt ? k ) ? 0?, ?k ? 1得协方差函数的递推公式? k ? ?1? k ?1 ? ? 2? k ? 2 ? ? ? ? p ? k ? p例3.3:求平稳AR(1)模型的协方差?递推公式? k ? ?1? k ?1 ? ? ? 0k 1?平稳AR(1)模型的方差为? ?2 ?0 ? 1 ? ?12?协方差函数的递推公式为2 ? ? k ? ?1k ? 2 , ?k ? 1 1 ? ?1例3.4:求平稳AR(2)模型的协方差?平稳AR(2)模型的协方差函数递推公式为1 ? ?2 ? 2 ? ? ? ? 0 (1 ? ? )(1 ? ? ? ? )(1 ? ? ? ? ) ? 2 1 2 1 2 ? ? ?1? 0 ?? 1 ? 1 ? ?2 ? ?? k ? ?1? k ?1 ? ? 2 ? k ? 2,k ? 2 ? ?自相关系数?自相关系数的定义?k ?k ? ?0?平稳AR(P)模型的自相关系数递推公式? k ? ?1 ? k ?1 ? ?2 ? k ? 2 ? L ? ? p ? k ? p常用AR模型自相关系数递推公式?AR(1)模型? k ? ?1k , k ? 0?AR(2)模型?1, ? ? ?1 ?k ? ? ?1 ? ? 2 ? ??1 ? k ?1 ? ? 2 ? k ? 2k ?0 k ?1 k?2AR模型自相关系数的性质?拖尾性? (k ) ? ? ci ?iki ?1 pc1 , c2 , ? , c p不能恒等于零?呈复指数衰减? (k ) ? ? ci ?iki ?1 p?0例3.5:考察如下AR模型的自相关图(1) xt ? 0.8 xt ?1 ? ? t (2) xt ? ?0.8 xt ?1 ? ? t (3) xt ? xt ?1 ? 0.5 xt ? 2 ? ? t (4) xt ? ? xt ?1 ? 0.5 xt ? 2 ? ? t例3.5― (1) xt ? 0.8 xt ?1 ? ? t?自相关系数按复指数单调收敛到零例3.5:― (2) xt ? ?0.8 xt ?1 ? ? t例3.5:― (3) xt ? xt ?1 ? 0.5 xt ?2 ? ? t?自相关系数呈现出“伪周期”性例3.5:―?(4) xt ? ? xt ?1 ? 0.5 xt ? 2 ? ? t自相关系数不规则衰减偏自相关系数?定义对于平稳AR(p)序列,所谓滞后k偏自相关系数就 xt ? k ?1 是指在给定中间k-1个随机变量 xt ?1 , xt ?2 ,? ,的 条件下,或者说,在剔除了中间k-1个随机变 x 量的干扰之后, xt对 影响的相关度量。用数 t ?k 学语言描述就是? x ,xt t ?kxt ?1 ,? , xt ? k ?1? x )( x ? E ? x )] E[( xt ? E t t ?k t ?k ? ?x ) 2 E[( x ? Et ?k t ?k偏自相关系数的计算?滞后k偏自相关系数实际上就等于k阶自 回归模型第个k回归系数的值。? ?1 ? ? k1 ? 0 ? ? k 2 ?1 ? ? ? ? kk ? k ?1 ? ? ? ? ? ?? ? ?? ?? ? ? 2 k1 1 k2 0 kk k ? 2 ? ???????? ? ? ? ? k ? ? k1 ? k ?1 ? ? k 2 ? k ? 2 ? ? ? ? kk ? 0? x )( x ? E ? x )] E[( xt ? E t t ?k kt ? kk ? ?x ) 2 ] E[( xt ? k ? E t ?k偏自相关系数的截尾性?AR(p)模型偏自相关系数P阶截尾?kk ? 0 , k ? p例3.5续:考察如下AR模型的偏自相关图(1) xt ? 0.8 xt ?1 ? ? t (2) xt ? ?0.8 xt ?1 ? ? t (3) xt ? xt ?1 ? 0.5 xt ? 2 ? ? t (4) xt ? ? xt ?1 ? 0.5 xt ? 2 ? ? t例3.5― (1) xt ? 0.8 xt ?1 ? ? t?理论偏自相关系数?样本偏自相关图?0.8 , k ? 1 ?kk ? ? ,k ? 2 ?0例3.5:― (2) xt ? ?0.8 xt ?1 ? ? t?理论偏自相关系数?样本偏自相关图??0.8 , k ? 1 ?kk ? ? ,k ? 2 ?0例3.5:― (3) xt ? xt ?1 ? 0.5 xt ?2 ? ? t?理论偏自相关系数?2 ,k ?1 ?3 ? ?kk ? ??0.5 , k ? 2 ?0 ,k ? 3 ? ??样本偏自相关图例3.5:―?(4) xt ? ? xt ?1 ? 0.5 xt ? 2 ? ? t?理论偏自相关系数,k ?1 ,k ? 2 ,k ? 3样本偏自相关系数图? 2 ?? 3 ? ?kk ? ??0.5 ?0 ? ?MA模型的定义?具有如下结构的模型称为 q 阶自回归模 型,简记为 MA(q )? xt ? ? ? ? t ? ?1? t ?1 ? ? 2? t ? 2 ? L ? ? q? t ? q ? ?? q ? 0 ? 2 ? ? ( ) 0 ( ) E ? Var ? ? , ? , E (? t ? s ) ? 0, s ? t t t ??特别当 ? ? 0 时,称为中心化 MA(q ) 模型移动平均系数多项式?引进延迟算子,中心化 MA( q ) 模型又可 以简记为xt ? ?( B)? t?q 阶移动平均系数多项式?( B ) ? 1 ? ? 1 B ? ? 2 B 2 ? ? ? ? q B qMA模型的统计性质?常数均值Ext ? E(? ? ? t ? ?1? t ?1 ? ? 2? t ? 2 ? ? ? ? q? t ? q) ???常数方差Var ( xt ) ? Var ( ? ? ? t ? ? 1? t ?1 ? ? 2 ? t ? 2 ? ? ? ? q ? t ? q ) ? (1 ? ? 12 ? ? ? ? q2 )? ?2MA模型的统计性质?自协方差函数P阶截 尾?自相关系数P阶截尾?(1 ? ?12 ? ? ? ?q2 )? ?2 , k ? 0 ? q ?k ? ? k ? ?(??k ? ??i?k ?i )? ?2 , 1 ? k ? q i ?1 ? ?0, k ?q ?k ?0 ?1, ? q?k ???k ? ??i?k ?i ? i ?1 ?k ? ? , 1? k ? q 2 2 ?1??1 ?? ??q ?0, k ?q ? ?常用MA模型的自相关系数?MA(1)模型,k ? 0 ,k ?1 ,k ? 2?MA(2)模型,k ? 0 ,k ?1 ,k ? 2 ,k ? 3?1 ? ? ? ?1 ?k ? ? 2 ? ? 1 1 ? ? ?0?1 ??? ?? ? 1 2 ? 1 ?1 ? ? 12 ? ? 22 ?k ? ? ? ?? 2 ?1 ? ? 12 ? ? 22 ? ?0MA模型的统计性质?偏自相关系数拖尾?kk ? (??1? t ?1 ? ? ? ? q? t ? q )(??1? t ? k ?1 ? ? ? ? q? t ? k ? q ?1 )?1 ,? ? q不恒为零 ? ?kk 不会在有限阶之后恒为零例3.6:考察如下MA模型的相关性质 (1) xt ? ? t ? 2? t ?1 (2) xt ? ? t ? 0.5? t ?1 4 16 (3) xt ? ? t ? ? t ?1 ? ? t ? 2 5 25 5 25 (4) xt ? ? t ? ? t ?1 ? ? t ? 2 4 16MA模型的自相关系数截尾?() 1 xt ? ? t ? 2? t ?1(2)xt ?? ? t ? 0.5? t ?1MA模型的自相关系数截尾4 16 ? 3)xt ? ? t ? ? t ?1 ? ( ? t ?2 5 255 25 ( ? t ?2 ? 4)xt ? ? t ? ? t ?1 ? 4 16MA模型的偏自相关系数拖尾?() 1 xt ? ? t ? 2? t ?1(2)xt ?? ? t ? 0.5? t ?1MA模型的偏自相关系数拖尾4 16 ? 3)xt ? ? t ? ? t ?1 ? ( ? t ?2 5 255 25 ( ? t ?2 ? 4)xt ? ? t ? ? t ?1 ? 4 16MA模型的可逆性?MA模型自相关系数的不唯一性?例3.6中不同的MA模型具有完全相同的自相 关系数和偏自相关系数(1) xt ? ? t ? 2? t ?1 ? (2) xt ? ? t ? 0.5? t ?1 4 16 5 25 (3) xt ? ? t ? ? t ?1 ? ? t ? 2 ? (4) xt ? ? t ? ? t ?1 ? ? t ? 2 5 25 4 16可逆的定义?可逆MA模型定义?若一个MA模型能够表示称为收敛的AR模 型形式,那么该MA模型称为可逆MA模型 一个自相关系数列唯一对应一个可逆MA 模型。?可逆概念的重要性?可逆MA(1)模型?xt ? ? t ? ?? t ?1?? ?? 1?? 2?xt ? ? t ? ? t ?11?xt ? ?t 1 ? ?Bxt ? ?t 1 1? B?? ? 1, 可逆? ? 1, 可逆MA模型的可逆条件?MA(q)模型的可逆条件是:?MA(q)模型的特征根都在单位圆内?i ? 1?等价条件是移动平滑系数多项式的根都在单 位圆外1?i?1逆函数的递推公式?原理? xt ? ?( B )? t ? ?( B ) I ( B ) xt ? xt ? ?? t ? I ( B ) xt?方法?待定系数法?递推公式?I 0 ? 1 ?? k , k ? q ? j , 其中? k? ? ? ? ?0, k ? q ? I j ? ? ? k? I j ? k,j ? 1,2, ? k ?1 ?例3.6续:考察如下MA模型的可逆性 (1) xt ? ? t ? 2? t ?1 (2) xt ? ? t ? 0.5? t ?1 4 16 (3) xt ? ? t ? ? t ?1 ? ? t ? 2 5 25 5 25 (4) xt ? ? t ? ? t ?1 ? ? t ? 2 4 16(1)―(2)? ? ?xt ? ? t ? 2? t ?1 ? ? ? 2 ? 1 ? 不可逆 xt ? ? t ? 0.5? t ?1 ? ? ? 0.5 ? 1 ? 可逆逆函数?1 Ik ? ? k ?0.5 , k ? 1?逆转形式? t ? ? 0.5k xt ? kk ?0?(3)―(4)?4 16 xt ? ? t ? ? t ?1 ? ? t ? 2 ? ? 2 ? 1, ? 2 ? ?1 ? 1 ? 可逆 5 25 5 25 25 xt ? ? t ? ? t ?1 ? ? t ? 2 ? ? 2 ? ? 1 ? 不可逆 4 16 16??逆函数?(?1) n ?1k , k ? 3n 或 3n ? 1 , n ? 0,1, ? Ik ? ? ?0, k ? 3n ? 2? n 3n ? n 3 n ?1?逆转形式n ?0? t ? ? (?1) 0.8 xt ?3n ? ? (?1) 0.8n ?0xt ?3n ?1ARMA模型的定义?具有如下结构的模型称为自回归移动平 均模型,简记为 ARMA( p, q )? xt ? ?0 ? ?1 xt ?1 ? ? ? ? p xt ? p ? ? t ? ?1? t ?1 ? ? ? ? q ? t ?q ? ?? p ? 0,? q ? 0 ? 2 , E ( ? ) 0 Var ( ? ) ? ? ? ? t t ? , E (? t ? s ) ? 0, s ? t ? Ex ? ? 0, ?s ? t ? s t?特别当 ?0 ? 0 时,称为中心化 ARMA( p, q ) 模型系数多项式?引进延迟算子,中心化 ARMA( p, q )模型 又可以简记为?( B) xt ? ?( B)? t?p 阶自回归系数多项式?( B) ? 1 ? ?1 B ? ?2 B 2 ? ? ? ? p B p?q 阶移动平均系数多项式?( B ) ? 1 ? ? 1 B ? ? 2 B 2 ? ? ? ? q B q平稳条件与可逆条件?ARMA(p,q)模型的平稳条件? ?P阶自回归系数多项式 ? ( B) ? 0 的根都在单位圆外 即ARMA(p,q)模型的平稳性完全由其自回归部分的 平稳性决定 q阶移动平均系数多项式?( B) ? 0的根都在单位圆外 即ARMA(p,q)模型的可逆性完全由其移动平滑部分 的可逆性决定?ARMA(p,q)模型的可逆条件? ?传递形式与逆转形式?传递形式? ? t ? ? G j? t ? jj ?1 ??逆转形式? t ? ? ?1 ( B ) ? ( B ) xt? xt ? ? I j xt ? jj ?1 ?xt ? ? ?1 ( B)?( B)? t?G0 ? 1 ? k ? ?Gk ? ? ? ?j Gk ? j ? ? ?j , j ?1 ?k ?1?I 0 ? 1 ? k ? , j ? I k ? ?? ?j I k ? j ? ? ? j ?1 ?k ?1ARMA(p,q)模型的统计性质?均值Ext ?1 ? ?1 ? ? ? ? p??0?协方差 自相关系数? (k ) ? ? ?2 ? Gi Gi ? ki ?0?? (k ) ? ? (k ) ? ? ( 0)?G Gj ?0 j ? j ?0?j?k2 G ? jARMA模型的相关性? ?自相关系数拖尾 偏自相关系数拖尾例3.7:考察ARMA模型的相关性?拟合模型ARMA(1,1): xt ? 0.5 xt ?1 ? ? t ? 0.8? t 并直观地考察该模型自相关系数和偏自 相关系数的性质。自相关系数和偏自相关系数拖尾性?样本自相关图?样本偏自相关图ARMA模型相关性特征模型AR(P) MA(q) ARMA(p,q)自相关系数 拖尾q阶截尾偏自相关系数P阶截尾拖尾 拖尾拖尾3.3平稳序列建模? ? ? ? ? ?建模步骤 模型识别 参数估计 模型检验 模型优化 序列预测建模步骤平 稳 非 白 噪 声 序 列 计 算 样 本 相 关 系 数 模型 识别 参数 估计N模型 检验Y模 型 优 化序 列 预 测计算样本相关系数?样本自相关系数?样本偏自相关系数?k ? ?? (xt ?1n?kt? x )( xt ? k ? x )2 ( x ? x ) ? t t ?1n? D k ? ?kk ? ? D模型识别?基本原则k? ?? ?选择模型AR(P) MA(q) ARMA(p,q)kk拖尾q阶截尾P阶截尾拖尾 拖尾拖尾模型定阶的困难因为由于样本的随机性,样本的相关系数不会 ?k 呈现出理论截尾的完美情况,本应截尾的 ? ? 仍会呈现出小值振荡的情况 或? kk ? 由于平稳时间序列通常都具有短期相关性,随 ? 都会衰减至零值附 ?k 与 ? ? 着延迟阶数 k ? ? , kk 近作小值波动 ? 在延迟若干阶之后衰减为小值波动 ? k 或? ?当 ? kk 时,什么情况下该看作为相关系数截尾,什么 情况下该看作为相关系数在延迟若干阶之后正 常衰减到零值附近作拖尾波动呢??样本相关系数的近似分布?Barlett1 ? k ~ N (0, ) , n ? ? ? n?Quenouille1 ? ?kk ~ N (0, ) , n ? ? n模型定阶经验方法?95%的置信区间? Pr ? ? ? ? Pr ? ? ? 2 2 ? ?k ? ?? ? ? 0.95 n n? 2 2 ? ? ? ?kk ? ? ? 0.95 n n??模型定阶的经验方法?如果样本(偏)自相关系数在最初的d阶明显大于两倍 标准差范围,而后几乎 95 %的自相关系数都落在 2 倍标准差的范围以内,而且通常由非零自相关系数 衰减为小值波动的过程非常突然。这时,通常视为 (偏)自相关系数截尾。截尾阶数为d。例2.5续?选择合适的模型 ARMA 拟合 1950 年 ―― 1998 年北京市城乡居民定期储蓄比例序 列。序列自相关图序列偏自相关图拟合模型识别???自相关图显示延迟3阶之后,自相关系数全部衰减 到2倍标准差范围内波动,这表明序列明显地短期 相关。但序列由显著非零的相关系数衰减为小值 波动的过程相当连续,相当缓慢,该自相关系数 可视为不截尾 偏自相关图显示除了延迟1阶的偏自相关系数显著 大于2倍标准差之外,其它的偏自相关系数都在2 倍标准差范围内作小值随机波动,而且由非零相 关系数衰减为小值波动的过程非常突然,所以该 偏自相关系数可视为一阶截尾 所以可以考虑拟合模型为AR(1)例3.8美国科罗拉多州某一加油站连续57天的OVERSHORT序列序列自相关图序列偏自相关图拟合模型识别?? ?自相关图显示除了延迟1阶的自相关系数在2倍 标准差范围之外,其它阶数的自相关系数都在2 倍标准差范围内波动。根据这个特点可以判断 该序列具有短期相关性,进一步确定序列平稳。 同时,可以认为该序列自相关系数1阶截尾 偏自相关系数显示出典型非截尾的性质。 综合该序列自相关系数和偏自相关系数的性 质,为拟合模型定阶为MA(1)例3.9?全球气表平均温度改变值差分序列序列自相关图序列偏自相关图拟合模型识别? ? ?自相关系数显示出不截尾的性质 偏自相关系数也显示出不截尾的性质 综合该序列自相关系数和偏自相关系数 的性质,可以尝试使用ARMA(1,1)模型拟 合该序列参数估计?待估参数?p ? q ? 2个未知参数?1 ,L , ? p ,?1 ,L ,? q , ? , ? ?2?常用估计方法? ? ?矩估计 极大似然估计 最小二乘估计矩估计?原理?样本自相关系数估计总体自相关系数?1 ? ?1 (?1 ,L , ? p ,?1 ,L ,? q ) ? ? ? M ? ? ? (? ,L , ? ,? ,L ,? ) ? ? ? p?q 1 p q ? p?q 1?样本一阶均值估计总体均值,样本方差估计总 体方差 n??x? ??xi ?1in?2 ? ? ? ? ?2 ? ? 1 1 p 2 ? ? ?2 ? ? ? ?2 ? ? ? ? ?2 x 1?? 1 q例3.10:求AR(2)模型系数的矩估计? ?AR(2)模型 xt ? ?1 xt ?1 ? ? 2 xt ? 2 ? ? t Yule-Walker方程? ?1 ? ?1 ? ? 2 ?1 ? ? ? 2 ? ?1 ?1 ? ? 2?矩估计(Yule-Walker方程的解) 2 ? ? ? 1 ? ? ? 2 ? ?1 2 ? ? ?1 ?1 ? ? ?2 ? 2 ?12 ?1 1? ? 1? ?例3.11:求MA(1)模型系数的矩估计? ?MA(1)模型 xt ? ? t ? ?1? t ?1 方程2 2 ? ? (1 ? ) ? ? ? ?1 ??1 ? 0 1 ? ? ?1 ? ? ? 2 2 ? ? ? 1 ? 0 1 ?? 1 ? ??1? ??矩估计?12 ? ? ? ? 1 1 4 ??1 ? ?1 2?例3.12:求ARMA(1,1)模型系数的矩估计? ?ARMA(1,1)模型 xt ? ?1 xt ?1 ? ? t ? ?1? t ?1 方程? 1 (?1 ? ?1 )(1 ? ?1?1 ) ? ? ? ?1 ? ?0 1 ? ?12 ? 2?1?1 ? ?? ? ? ? 1 1 ? 2?矩估计?c ? ? ? ? ? 2 ? ? ?1 ? , ?1 ? ? ?1 ? ?c ? ? ?c2 ? 4 , c ? ?2 2 c2 ? 4 ,c ? 2 2?2 1 ? ?12 ? 2 ? , c? ?1 ?1 ? ?对矩估计的评价?优点? ? ?估计思想简单直观 不需要假设总体分布 计算量小(低阶模型场合) 信息浪费严重??缺点?只用到了p+q个样本自相关系数信息,其他信息都被忽略?估计精度差?通常矩估计方法被用作极大似然估计和最小二 乘估计迭代计算的初始值极大似然估计?原理?在极大似然准则下,认为样本来自使该样本 出现概率最大的总体。因此未知参数的极大 似然估计就是使得似然函数(即联合密度函 数)达到最大的参数值~ ? ,? ? ,? , ? ? ;x ,~ L( ? x ) ? max{ p ( x ); ?1 , ? 2 ,? , ? k } 1 2 k 1似然方程~ ? ? n S (? ) ~ ~ ? ?0 ? 2 l (? ; x ) ? 2 4 2? ? 2? ? ? ?? ? ? ~ 1 ?S ( ? ) ~ ~ 1 ? ln ? ? ? l (? ;x) ? ~ ? ~ ?0 ~ 2 ? ?? 2 ?? 2? ? ?? ??% %的显式表达式。 ) 和 ln ? 都不是 ? 由于S ( ? 因而似然方程组实际上是由p+q+1个超越 方程构成,通常需要经过复杂的迭代算 法才能求出未知参数的极大似然估计值对极大似然估计的评价?优点??极大似然估计充分应用了每一个观察值所提 供的信息,因而它的估计精度高 同时还具有估计的一致性、渐近正态性和渐 近有效性等许多优良的统计性质 需要假定总体分布?缺点?最小二乘估计?原理?使残差平方和达到最小的那组参数值即为最 小二乘估计值~ ? Q( ? ) ? min Q( ? ) ? min ? ( xt ? ?1 xt ?1 ? ? ? ? p xt ? p ? ?1? t ?1 ? ? ? ? q? t ? q )t ?1 n 2条件最小二乘估计? ?实际中最常用的参数估计方法 假设条件xt ? 0 , t ? 0n 2 t n?残差平方和方程Q ( ? ) ? ? ? ? ? [ xt ? ? ? i xt ?1 ]i ?1 i ?1 i ?1?~t2解法?迭代法对最小二乘估计的评价?优点??最小二乘估计充分应用了每一个观察值所提 供的信息,因而它的估计精度高 条件最小二乘估计方法使用率最高 需要假定总体分布?缺点?例2.5续?确定1950年――1998年北京市城乡居民 定期储蓄比例序列拟合模型的口径? ? ?拟合模型:AR(1) 估计方法:极大似然估计 模型口径xt ? 25.17 ? 0.69 xt ?1 ? ? t? ?2 ) ? 16.17 Var (?例3.8续?确定美国科罗拉多州某一加油站连续57 天的OVERSHORTS序列拟合模型的口径? ? ?拟合模型:MA(1) 估计方法:条件最小二乘估计 模型口径xt ? ?4.40351 ? (1 ? 0.82303B)? t2 ? Var (? ? ) ?
例3.9续?确定全球气表平均温度改变值 差分序列拟合模型的口径? ? ?拟合模型:ARMA(1,1) 估计方法:条件最小二乘估计 模型口径xt ? 0.003 ? 0.407 xt ?1 ? ? t ? 0.9? t ?1? ?2 ) ? 0.016 Var (?模型检验?模型的显著性检验?整个模型对信息的提取是否充分 模型结构是否最简?参数的显著性检验?模型的显著性检验?目的?检验模型的有效性(对信息的提取是否充分) 残差序列 一个好的拟合模型应该能够提取观察值序列中几乎所 有的样本相关信息,即残差序列应该为白噪声序列 反之,如果残差序列为非白噪声序列,那就意味着残 差序列中还残留着相关信息未被提取,这就说明拟合 模型不够有效?检验对象??判定原则??假设条件?原假设:残差序列为白噪声序列H 0:?1 ? ? 2 ? L ? ? m ? 0, ?m ? 1?备择假设:残差序列为非白噪声序列H 1:至少存在某个? k ? 0, ?m ? 1,k ? m检验统计量?LB统计量2 ? ?k 2 LB ? n(n ? 2)? ( ) ~ ? ( m) k ?1 n ? k m例2.5续??检验1950年――1998年北京市城乡居民 定期储蓄比例序列拟合模型的显著性 残差白噪声序列检验结果LB统计量 5.83 10.28 11.38 P值 0.0 0.8361延迟阶数 612 18检验结论 拟合模型 显著有效参数显著性检验?目的?检验每一个未知参数是否显著非零。删除不 显著参数使模型结构最精简H0 : ? j ? 0 ? H1 : ? j ? 0 ?1 ? j ? m?假设条件 检验统计量T ? n?m?? ?? ? j j~ ~ t ( n ? m) a jj Q( ? )例2.5续?检验1950年――1998年北京市城乡居民定期储 蓄比例序列极大似然估计模型的参数是否显著 参数检验结果 检验参数 t统计量 均值46.12 6.72 P值 &0.1?结论 显著 显著?1例3.8续:对OVERSHORTS序列的拟合模型进行检验?残差白噪声检验 延迟阶数 6 12 LB统计量 P值3.15 9.05t统计量 -3.75 10.600.1P值 &0.1结论 模型显著 有效?参数显著性检验 检验参数 均值 结论 显著 显著?1例3.9续:对全球气表平均温度改变值差分 序列拟合模型进行检验?残差白噪声检验 延迟阶数 6 12 LB统计量 P值5.28 10.30t统计量 16.34 3.50.7P值 &0.7结论 模型显著 有效?参数显著性检验 检验参数?1 ?1结论 显著 显著模型优化?问题提出?当一个拟合模型通过了检验,说明在一定的 置信水平下,该模型能有效地拟合观察值序 列的波动,但这种有效模型并不是唯一的。 选择相对最优模型?优化的目的?例3.13:拟合某一化学序列序列自相关图序列偏自相关图拟合模型一? ?根据自相关系数2阶截尾,拟合MA(2)模型 参数估计yield t ? 51.17301 ? (1 ? 0.32286 B ? 0.31009 B 2 )? t?模型检验? ?模型显著有效 三参数均显著拟合模型二? ?根据偏自相关系数1阶截尾,拟合MA(1)模型 参数估计yield t ? 51.26169 ???t1 ? 0.42481B模型检验 ? 模型显著有效 ? 两参数均显著问题??同一个序列可以构造两个拟合模型,两 个模型都显著有效,那么到底该选择哪 个模型用于统计推断呢? 解决办法?确定适当的比较准则,构造适当的统计量, 确定相对最优AIC准则??最小信息量准则(An Information Criterion) 指导思想? ?似然函数值越大越好 未知参数的个数越少越好2 ? AIC ? n ln(? ? ) ? 2(未知参数个数 )?AIC统计量SBC准则?AIC准则的缺陷?在样本容量趋于无穷大时,由AIC准则选择 的模型不收敛于真实模型,它通常比真实模 型所含的未知参数个数要多?SBC统计量? ?2 ) ? ln(n)(未知参数 ) SBC ? n ln(?例3.13续?用AIC准则和SBC准则评判例3.13中两个 拟合模型的相对优劣模型 MA(2)AR(1) AIC 536.6 SBC 543.6?结果?AR(1)优于MA(2)序列预测?线性预测函数xt ? ? Ci xt ?1?ii ?0??预测方差最小原则Varx ?t ( l ) ? et (l ) ? ? min ?Var ? et (l ) ??序列分解xt ?l ? ? ? t ?l ? G1? t ?l ?1 ? L ? Gl ?1? t ?1 ? ? ? Gl ? t ? Gl ?1? t ?1 ? L ? et (l )预测误差???t (l ) x预测值? (l ) E ( xt ?l xt , xt ?1 , ? ) ? x Var ( xt ? l xt , xt ?1 , ? ) ? Var[et (l )]误差分析?估计误差et (l ) ? ? t ? l ? G1? t ? l ?1 ? ? Gl ?1? t ?1?期望 方差E[et (l )] ? 0Var[et (l )] ? ? Gi2 ? ? ?2i ?0 l ?1?AR(p)序列的预测?预测值 ? (l ) ? ?1 x ?t (l ? 1) ? ? ? ? p x ?t (l ? p ) x 预测方差?Var[et (l )] ? (1 ? G ? ? ? G )? ?2 1 2 l ?1?295%置信区间1 ? ? 2 2 ?t (l ) mz ? ? ?1 ? G1 ? L ? Gl ?1 ? 2 ? ? ? ? ?x 1? ? ? 2例3.14?已知某超市月销售额近似服从AR(2)模型 (单位:万元/每月)xt ? 10 ? 0.6 xt ?1 ? 0.3 xt ? 2 ? ? t , ? t ~ N (0,36)??今年第一季度该超市月销售额分别为: 101,96,97.2万元 请确定该超市第二季度每月销售额的95 %的置信区间例3.14解:预测值计算?四月份? 3 (1) ? 10 ? 0.6 x3 ? 0.3x 2 ? 97.12 x?五月份? 3 (2) ? 10 ? 0.6 x ? 3 (1) ? 0.3 x3 ? 97.432 x?六月份? 3 (3) ? 10 ? 0.6 x ? 3 (2) ? 0.3x ? 3 (1) ? 97.5952 x例3.14解:预测方差的计算?GREEN函数G0 ? 1 G1 ? ?1G0 ? 0.6 G2 ? ?1G1 ? ?2G0 ? 0.36 ? 0.3 ? 0.66?方差Var[e3 (1)] ? G02? ?2 ? 36 Var[e3 (2)] ? (G02 ? G12 )? ?2 ? 48.962 Var[e3 (3)] ? (G02 ? G12 ? G2 )? ?2 ? 64.6416例3.14解:置信区间?公式? 3 (l ) ? 1.96 Var[e3 (l )] , x ? 3 (l ) ? 1.96 Var[e3 (l )] ) (x?估计结果预测时期 四月份 五月份 六月份95%置信区间 (85.36,108.88) (83.72,111.15) (81.84,113.35)例2.5:北京市城乡居民定期储蓄比例序列 拟合与预测图MA(q)序列的预测?预测值q ? ? ? ? ? ? i ? t ? l ?i ?t (l ) ? ? x i ?l ?? ?,l ? q ,l ? q?预测方差2 ? (1 ? ?12 ? ? ? ? l2 ) ? ?1 ? ? Var[et (l )] ? ? (1 ? ? 2 ? ? ? ? 2 )? 2 q ? 1 ?,l ? q ,l ? q例3.15?已知某地区每年常驻人口数量近似服从MA(3)模型(单 位:万):xt ? 100 ? ? t ? 0.8? t ?1 ? 0.6? t ?2 ? 0.2? t ?1?最近3年的常驻人口数量及一步预测数量如下:年份 04?统计人数 104 108 105预测人数 110 100 109预测未来5年该地区常住人口的95%置信区间例3.15解:随机扰动项的计算? 2001 (1) ? 104 ? 110 ? ?6 ? t ?2 ? x2002 ? x ? 2002 (1) ? 108 ? 100 ? 8 ? t ?1 ? x2003 ? x ? 2003 (1) ? 105 ? 109 ? ?4 ? t ? x2004 ? x例3.15解:估计值的计算? t (1) ? 100 ? 0.8? t ? 0.6? t ?1 ? 0.2? t ?2 ? 109.2 x ? t (2) ? 100 ? 0.6? t ? 0.2? t ?1 ? 96 x ? t (3) ? 100 ? 0.2? t ? 100.8 x ? t (4) ? 100 x ? t (5) ? 100 x例3.15解:预测方差的计算Var[et (1)] ? ? ? ? 252Var[et (2)] ? (1 ? ? )? ? ? 412 1 2Var[et (3)] ? (1 ? ? ? ? )? ? ? 502 1 2 2 2Var[et (4)] ? (1 ? ? ? ? ? ? )? ? ? 512 1 2 2 2 3 2Var[et (5)] ? (1 ? ? ? ? ? ? )? ? ? 512 1 2 2 2 3 2例3.15解:置信区间的计算预测年份 07 95%置信区间 (99,119) (83,109) (87,115) (86,114) (86,114)ARMA(p,q)序列预测?预测值? t (k ) , k ? 1 ?x ? t (k ) ? ? x ? xt ? k , k ? 0?预测方差Var[et (l )] ? (G ? G ? ? ? G )? ?1 0 2 1 2 l ?1 2例3.16?已知模型为:xt ? 0.8 xt ?1 ? ? t ? 0.6? t ?1 ? ?2 ? 0.0025?且x100 ? 0.3?100 ? 0.01?预测未来3期序列值的95%的置信区间。例3.16解:估计值的计算?100 (1) ? 0.8 x100 ? 0.6? 100 ? 0.234 x ?100 (2) ? 0.8 x ?100 (1) ? 0.1872 x ?100 (3) ? 0.8 x ?100 (2) ? 0.14976 x例3.16解:预测方差的计算?Green函数G0 ? 1 G1 ? ?1G0 ? ?1 ? 0.2 G2 ? ?1G1 ? 0.16?方差Var[e100 (1)] ? G02? ?2 ? 0.0025 Var[e100 (2)] ? (G02 ? G12 )? ?2 ? 0.00262 Var[e100 (3)] ? (G02 ? G12 ? G2 )? ?2 ? 0.002664例3.16解:置信区间的计算时期101 102 103 95%置信区间(0.136,0.332) (0.087,0.287) (-0.049,0.251)修正预测?定义?所谓的修正预测就是研究如何利用新的信息 去获得精度更高的预测值 在新的信息量比较大时――把新信息加入到 旧的信息中,重新拟合模型 在新的信息量很小时――不重新拟合模型, 只是将新的信息加入以修正预测值,提高预 测精度?方法??修正预测原理?xt ?l 的预测值为 在旧信息的基础上, ? t (l ) ? Gl ? t ? G1? t ?1 ? ? x假设新获得一个观察值 xt ?1 ,则 ? xt ?l 的修正预测值为?t ?1 (l ? 1) ? Gl ?1? t ?1 ? Gl ? t ? Gl ?1? t ?1 ? ? ? Gl ?1? t ?1 ? x ?t (l ) x?修正预测误差为 et ?1 (l ? 1) ? G0? t ?l ? ? ? Gl ?2? t ? 2 ? 预测方差为?Var[et ?1 (l ? 1)] ? (G02 ? ? ? Gl2? 2 )? ?2一般情况?假设新获得p个观察值 xt ?1 ,? , xt ? p ,则 ? xt ? l 的修正预测值为?t ? p (l ? p) ? Gl ? p? t ? p ? ? ? Gl ?1? t ?1 ? x ?t (l ) x?修正预测误差为 et ? p (l ? p ) ? G0? t ?l ? ? ? Gl ? p ?1? t ? p ?1 预测方差为Var[et ? p (l ? p)] ? (G02 ? ? ? Gl2? p ?1 )? ?2?例3.14续:假如四月份的真实销售额为100万元, 求二季度后两个月销售额的修正预测值?计算四月份的预测误差?3 (1) ? 100 ? 97.12 ? 2.88 ? 4 ? x4 ? x?计算修正预测值?4 (1) ? G1? 4 ? x ?3 (2) ? 99.16 x ?4 (2) ? G2? 4 ? x ?3 (3) ? 99.50 x?计算修正方差Var[e4 (1)] ? Var[e3 (1)] ? G02? ?2 ? 36 Var[e4 (2)] ? Var[e3 (2)] ? (G02 ? G12 )? ?2 ? 48.96修正置信区间预测时期 修正前置信区间 修正后置信区间 四月份 五月份 六月份(85.36,108.88) (83.72,111.15)(87.40,110.92) (81.84,113.35)(85.79,113.21)第四章非平稳序列的确定性分析本章结构? ? ? ? ? ?时间序列的分解 确定性因素分解 趋势分析 季节效应分析 综合分析 X-11过程4.1 时间序列的分解? ?Wold分解定理 Cramer分解定理Wold分解定理(1938)?{xt }它都可以分解为 对于任何一个离散平稳过程 两个不相关的平稳序列之和,其中一个为确定 性的,另一个为随机性的,不妨记作 xt ? Vt ? ? t ? ? t ? ?? j? t? j ?? t ? 为随机序列, {Vt } 其中: 为确定性序列, j ?0 它们需要满足如下条件 ? 2 ? ? ? ~ WN ( 0 , ? (1)? 0 ? 1, ? ? 2 ( 2 ) ? ? t ? ) jj ?0E (Vt , ? s ) ? 0, ?t ? s (3)确定性序列与随机序列的定义?对任意序列 ?yt ? 而言,令 yt 关于q期之前 的序列值作线性回归 y t ? ? 0 ? ? 1 y t ? q ? ? 2 y t ? q ?1 ? ? ? ? t2 其中 {? t }为回归残差序列, 。 Var (? t ) ? ? q?2 lim ? 确定性序列,若 q ?? q ? 0?2 lim ? 随机序列,若 q ?? q ? Var ( yt )ARMA模型分解?( B ) xt ? ? ? ?t ?( B)确定性序列随机序列Cramer分解定理(1961)?任何一个时间序列 {xt } 都可以分解为两部分的叠 加:其中一部分是由多项式决定的确定性趋势成 分,另一部分是平稳的零均值误差成分,即xt ? ?t ? ? tj ? t ? j j ?0 d? ( B )at随机性影响确定性影响对两个分解定理的理解??Wold分解定理说明任何平稳序列都可以分解为 确定性序列和随机序列之和。它是现代时间序 列分析理论的灵魂,是构造ARMA模型拟合平 稳序列的理论基础。 Cramer 分解定理是Wold分解定理的理论推 广,它说明任何一个序列的波动都可以视为同 时受到了确定性影响和随机性影响的综合作用。 平稳序列要求这两方面的影响都是稳定的,而 非平稳序列产生的机理就在于它所受到的这两 方面的影响至少有一方面是不稳定的。4.2确定性因素分解?传统的因素分解? ? ? ??现在的因素分解? ? ?长期趋势 循环波动 季节性变化 随机波动长期趋势波动 季节性变化 随机波动确定性时序分析的目的??克服其它因素的影响,单纯测度出某一 个确定性因素对序列的影响 推断出各种确定性因素彼此之间的相互 作用关系及它们对序列的综合影响4.3趋势分析?目的?有些时间序列具有非常显著的趋势,我们分 析的目的就是要找到序列中的这种趋势,并 利用这种趋势对序列的发展作出合理的预测 趋势拟合法 平滑法?常用方法? ?趋势拟合法??趋势拟合法就是把时间作为自变量,相 应的序列观察值作为因变量,建立序列 值随时间变化的回归模型的方法 分类? ?线性拟合 非线性拟合线性拟合?使用场合?长期趋势呈现出线形特征?模型结构? xt ? a ? bt ? I t ? ? E ( I t ) ? 0,Var ( I t )例4.1:拟合澳大利亚政府1981―― 1990年每季度的消费支出序列线性拟合?模型, t ? 1,2 ? ,40 ? xt ? a ? bt ? I t ? 2 E ( I ) 0 , Var ( I ) ? ? ? t t ??参数估计方法?最小二乘估计?参数估计值 ? ? 8498.69 a? ? 89.12 , b拟合效果图非线性拟合?使用场合?长期趋势呈现出非线形特征 能转换成线性模型的都转换成线性模型,用 线性最小二乘法进行参数估计 实在不能转换成线性的,就用迭代法进行参 数估计?参数估计指导思想??常用非线性模型模型Tt ? a ? bt ? ct 2变换t2 ? t 2? Tt ? ln Tt变换后模型Tt ? a ? bt ? ct 2参数估计方法线性最小二乘估计Tt ? ab tTt ? a ? bc t1 Tt ? a ? bc ta ? ? ln a b? ? ln b? Tt ? a ? ? b?t线性最小二乘估计- - -- - -迭代法 迭代法 迭代法Tt ? ea ? bc t例4.2: 对上海证券交易所每月末上 证指数序列进行模型拟合非线性拟合?模型 变换Tt ? a ? bt ? ct 2t2 ? t 2???参数估计方法 ? 线性最小二乘估计 拟合模型口径Tt ? 502.2517 ? 0.0952t 2拟合效果图平滑法??平滑法是进行趋势分析和预测时常用的 一种方法。它是利用修匀技术,削弱短 期随机波动对序列的影响,使序列平滑 化,从而显示出长期趋势变化的规律 常用平滑方法? ?移动平均法 指数平滑法移动平均法?基本思想?假定在一个比较短的时间间隔里,序列值之 间的差异主要是由随机波动造成的。根据这 种假定,我们可以用一定时间间隔内的平均 值作为某一期的估计值 n期中心移动平均 n期移动平均?分类? ?n期中心移动平均?1 ? n ( xt ? n ?1 ? xt ? n ?1 ?1 ? ? ? xt ? ? ? xt ? n ?1 ?1 ? xt ? n ?1 ),n为奇数 ? 2 2 2 2 ~ xt ? ? 1 ?1 (1 x ? x ? ? ? xt ? ? ? x n ? x n ),n为偶数 n n 1 ? ? ? t t t ? ?1 ? 2 t? 2 2 2 ?n 2 25xt ? 2xt ?1xtxt ?1xt ? 2期中心移动平均xt ? 2 ? xt ?1 ? xt ? xt ?1 ? xt ? 2 ~ xt ? 5n期移动平均1 ~ xt ? ( xt ? xt ?1 ? ? ? xt ? n ?1 ) n5 xt ? 4 xt ? 3 xt ? 2 xt ?1 xt 期移动平均x ?x ?x ?x ?x ~ xt ? t ? 4 t ?3 t ? 2 t ?1 t 5移动平均期数确定的原则?事件的发展有无周期性?以周期长度作为移动平均的间隔长度 ,以 消除周期效应的影响 移动平均的期数越多,拟合趋势越平滑 移动平均的期数越少,拟合趋势越敏感?对趋势平滑的要求??对趋势反映近期变化敏感程度的要求?移动平均预测?T ?l x 1 ? ?l ?1 ? xT ? ?l ? 2 ? ? ? xT ? ?l ?n ) ? ( xT n? ? l ?i xT?T ?l ?i , l ? i ?x ?? ? xT ?l ?i , l ? i例4.3?某一观察值序列最后4期的观察值为: 5,5.5,5.8,6.2 ?T ? 2。 (1)使用4期移动平均法预测 x ?T ? 2 中 xT 前面的系数 (2)求在二期预测值 x 等于多少?例4.3解1 5 ? 5 .4 ? 5 . 8 ? 6 . 2 ?T ?1 ? ? xT ? xT ?1 ? xT ? 2 ? xT ?3 ? ? ? 5.6 (1)x 4 4 1 5.6 ? 5 ? 5.4 ? 5.8 ?T ? 2 ? ?x ?T ?1 ? xT ? xT ?1 ? xT ? 2 ? ? x ? 5.45 4 4?T ? 2 (2)x 1 ?T ?1 ? xT ? xT ?1 ? xT ? 2 ? ?x 4 1 ?1 ? ? ? ? xT ? xT ?1 ? xT ? 2 ? xT ? 2 ? ? xT ? xT ?1 ? xT ? 2 ? 4 ?4 ? 5 1 ? ? xT ? xT ?1 ? xT ? 2 ? ? xT ?3 16 16 5 在二期预测值中 xT 前面的系数等于 16 ?指数平滑法?指数平滑方法的基本思想?在实际生活中,我们会发现对大多数随机事件而 言,一般都是近期的结果对现在的影响会大些,远 期的结果对现在的影响会小些。为了更好地反映这 种影响作用,我们将考虑到时间间隔对事件发展的 影响,各期权重随时间间隔的增大而呈指数衰减。 这就是指数平滑法的基本思想 简单指数平滑 Holt两参数指数平滑?分类? ?简单指数平滑?基本公式~ xt ? ?xt ? ? (1 ? ? ) xt ?1 ? ? (1 ? ? ) 2 xt ? 2 ? ??等价公式 ~ ~ xt ? ?xt ? (1 ? ? ) xt ?1经验确定?初始值的确定~ x0 ? x1?平滑系数的确定? ? ??常取较小的值 一般对于变化缓慢的序列, ? 常取较大的值 对于变化迅速的序列, 经验表明 ? 的值介于0.05至0.3之间,修匀 效果比较好。简单指数平滑预测?一期预测值?T ?1 ? ~ x xT ? ?xT ? ? (1 ? ? ) xT ?1 ? ? (1 ? ? ) 2 xT ? 2 ? ??二期预测值?T ? 2 ? ?x ?T ?1 ? ? (1 ? ? ) xT ? ? (1 ? ? ) 2 xT ?1 ? ? x ?T ?1 ? (1 ? ? ) x ?T ?1 ? x ?T ?1 ? ?x?l 期预测值 ?T ?l ? x ?T ?1 , l ? 2 x例4.4?对某一观察值序列?xt ?使用指数平滑法。 ~ xT ?1 ? 10.5 ,平滑系数 ? ? 0.25 已知 xT ? 10 , ?T ? 2 。 (1) 求二期预测值 x ?T ? 2中 xT 前面的系数等于 (2)求在二期预测值 x 多少?例4.4解~ ~ ? ? ? 0 . 25 ? 0 . 75 x x x xT ?1 ? 10.3 (1) T ?1 T T ?T ? 2 ? x ?T ?1 ? 10.3 x?T ? 2 ? x ?T ?1 ? ?xT ? ? (1 ? ? ) xT ?1 ? ? (2)x所以使用简单指数平滑法二期预测值中 xT 前面的 系数就等于平滑系数 ? ? 0.25Holt两参数指数平滑?使用场合?适用于对含有线性趋势的序列进行修匀 假定序列有一个比较固定的线性趋势? t ? xt ?1 ? r x?构造思想??两参数修匀xt ? ?xt ? (1 ? ? )( ~ xt ?1 ? rt ?1 ) ?~ ? ~ ~ r ? ( x xt ?1 ) ? (1 ? ? )rt ?1 ? ? t ?t初始值的确定?平滑序列的初始值~ x0 ? x1?趋势序列的初始值x n ?1 ? x1 r0 ? nHolt两参数指数平滑预测?l 期预测值?T ?l x~ ? xT ? l ? rT例4.5?对北京市1978 ―― 2000年报纸发行量序 列进行Holt两参数指数平滑。指定~ x0 ? x1 ? 51259x 23 ? x1 r0 ? ? 4325 23? ? 0.15? ? 0 .1例4.5平滑效果图4.3 季节效应分析【例4.6】以北京市1995年――2000年月平 均气温序列为例,介绍季节效应分析的 基本思想和具体操作步骤。时序图季节指数?季节指数的概念?所谓季节指数就是用简单平均法计算的周期 内各时期季节性影响的相对数?季节模型xij ? x ? S j ? I ij季节指数的计算?计算周期内各期平均数xk ??xi ?1nik?计算总平均数x?n, k ? 1,2, ? , m?? xi ?1 k ?1nmik?计算季节指数nmxk Sk ? x, k ? 1,2,? , m季节指数的理解????季节指数反映了该季度与总平均值之间 的一种比较稳定的关系 如果这个比值大于1,就说明该季度的值 常常会高于总平均值 如果这个比值小于1,就说明该季度的值 常常低于总平均值 如果序列的季节指数都近似等于1,那就 说明该序列没有明显的季节效应例4.6季节指数的计算例4.6季节指数图综合分析?常用综合分析模型??加法模型 xt ? Tt ? S t ? I t 乘法模型xt ? Tt ? S t ? I t?混合模型a ) xt ? S t ? Tt ? I t b) xt ? S t ? (Tt ? I t )例4.7?对1993年 ―― 2000年中国社会消费品零 售总额序列(数据见附录1.11)进行确 定性时序分析。(1)绘制时序图(2)选择拟合模型?长期递增趋势和以年为固定周期的季节 波动同时作用于该序列,因而尝试使用 混合模型(b)拟合该序列的发展xt ? S t ? (Tt ? I t )(3)计算季节指数月份 1 2 3 4 5 6 季节指数 0.982 0.943 0.920 0.911 0.925 0.951 月份 7 8 9 10 11 12 季节指数 0.929 0.940 1.001 1.054 1.100 1.335季节指数图季节调整后的序列图xt ? Tt ? I t ? St(4)拟合长期趋势? ?
? 20.93178t T t(5)残差检验xt ? ? Tt ? I t ? St(6)短期预测? ?T ? ?t (l ) ? S x t ?l t ?lX-11过程?简介?X-11过程是美国国情调查局编制的时间序列季节调整过 程。它的基本原理就是时间序列的确定性因素分解方法 长期趋势起伏 季节波动 不规则波动 交易日影响 加法模型 乘法模型?因素分解? ? ? ??模型? ?方法特色?普遍采用移动平均的方法? ? ?用多次短期中心移动平均消除随机波动 用周期移动平均消除趋势 用交易周期移动平均消除交易日影响例4.7续??对1993年――2000年中国社会消费品零 售总额序列使用X-11过程进行季节调整 选择模型(无交易日影响)xt ? Tt St I tX11过程获得的季节指数图季节调整后的序列图趋势拟合图随机波动序列图第五章非平稳序列的随机分析本章结构? ? ? ? ? ?差分运算 ARIMA模型 Auto-Regressive模型 异方差的性质 方差齐性变化 条件异方差模型5.1 差分运算? ? ?差分运算的实质 差分方式的选择 过差分差分运算的实质???差分方法是一种非常简便、有效的确定 性信息提取方法 Cramer分解定理在理论上保证了适当阶 数的差分一定可以充分提取确定性信息 差分运算的实质是使用自回归的方式提 取确定性信息i ? d xt ? (1 ? B ) d xt ? ? (?1) i C d xt ?i i ?0 d差分方式的选择???序列蕴含着显著的线性趋势,一阶差分 就可以实现趋势平稳 序列蕴含着曲线趋势,通常低阶(二阶 或三阶)差分就可以提取出曲线趋势的 影响 对于蕴含着固定周期的序列进行步长为 周期长度的差分运算,通常可以较好地 提取周期信息例5.1【例1.1】1964年――1999年中国纱年产 量序列蕴含着一个近似线性的递增趋势。 对该序列进行一阶差分运算?xt ? xt ? xt ?1考察差分运算对该序列线性趋势信息的提 取作用差分前后时序图?原序列时序图?差分后序列时序图例5.2?尝试提取1950年――1999年北京市民用 车辆拥有量序列的确定性信息差分后序列时序图?一阶差分?二阶差分例5.3?差分运算提取1962年1月――1975年12月平均 每头奶牛的月产奶量序列中的确定性信息差分后序列时序图?一阶差分?1阶-12步差分过差分??足够多次的差分运算可以充分地提取原 序列中的非平稳确定性信息 但过度的差分会造成有用信息的浪费例5.4?假设序列如下 x t ? ? 0 ? ? 1t ? a t 考察一阶差分后序列和二阶差分序列 的平稳性与方差?比较?一阶差分??二阶差分(过差分)?平稳?xt ? xt ? xt ?1 ? ?1 ? at ? at ?1平稳? 2 xt ? ?xt ? ?xt ?1 ? at ? 2at ?1 ? at ? 2?方差小?方差大? 6? 2Var (?xt ) ? Var (at ? at ?1 ) ? 2?2Var(?2 xt ) ? Var(at ? 2at ?1 ? at ?2 )5.2 ARIMA模型? ? ? ? ? ?ARIMA模型结构 ARIMA模型性质 ARIMA模型建模 ARIMA模型预测 疏系数模型 季节模型ARIMA模型结构?使用场合?差分平稳序列拟合?模型结构??( B)? d xt ? ?( B)? t ? 2 , E ? Var ? ( ? ) 0 ( ? ) ? ? t t ? , E (? t ? s ) ? 0, s ? t ? Ex ? ? 0, ?s ? t ? s tARIMA 模型族?d=0 ARIMA(p,d,q)=ARMA(p,q) P=0 ARIMA(P,d,q)=IMA(d,q) q=0 ARIMA(P,d,q)=ARI(p,d) d=1,P=q=0 ARIMA(P,d,q)=random walk model???随机游走模型( random walk)?模型结构? xt ? xt ?1 ? ? t ? 2 ? ? E ( ? ) 0 Var ( ? ) ? , ? ? , E (? t ? s ) ? 0, s ? t t t ? Ex ? ? 0, ?s ? t ? s t?模型产生典故?Karl Pearson(1905)在《自然》杂志上提问:假如有个 醉汉醉得非常严重,完全丧失方向感,把他放在荒郊 野外,一段时间之后再去找他,在什么地方找到他的 概率最大呢?ARIMA模型的平稳性?ARIMA(p,d,q) 模 型 共 有 p+d 个 特 征 根,其中p个在单 位圆内,d个在单 位圆上。所 以 d? 0当 时 ARIMA(p,d,q) 模 型非平稳。?例5.5 ARIMA(0,1,0)时序图ARIMA模型的方差齐性?d ? 0时,原序列方差非齐性ARIMA(0,1,0)模型 Var ( xt ) ? Var ( x0 ? ? t ? ? t ?1 ? ? ? 1 ) ? t? ?2?d阶差分后,差分后序列方差齐性ARIMA(0,1,0)模型 Var (?xt ) ? Var (? t ) ? ? ?2ARIMA模型建模步骤获 得 观 察 值 序 列 Y Y平稳性 检验 N 差分 运算白噪声 检验 N 拟合 ARMA 模型分 析 结 束例5.6?对1952年――1988年中国农业实际国民 收入指数序列建模一阶差分序列时序图一阶差分序列自相关图一阶差分后序列白噪声检验延迟阶数 6 12 18? 2 统计量P值 0.0 0.134415.33 18.33 24.66拟合ARMA模型?偏自相关图建模?定阶?ARIMA(0,1,1)?参数估计(1 ? B) xt ? 4.99661 ? (1 ? 0.70766 B)? t Var (? t ) ? 56.48763?模型检验? ?模型显著 参数显著ARIMA模型预测?原则?最小均方误差预测原理?Green函数递推公式?? 1 ? ?1 ? ?1 ?? ? ? ? ? ? ? ? ? 2 1 1 2 2 ? ?? ?? j ? ?1? j ?1 ? ? ? ? p ? d? j ? p ? d ? ? j ?预测值xt?l ? (?t?l ??1?t?l?1 ?? ??l?1?t?1) ?(?l?t ??l?1?t?1 ?? )et (l )? t (l ) xE[et (l )] ? 0 Var[et (l )] ? (1 ? ? ? ? ? ?2 1 2 l ?1)? ?2例5.7?已知ARIMA(1,1,1)模型为(1 ? 0.8 B)(1 ? B) xt ? (1 ? 0.6 B)? t且 xt ?1 ? 4.5?xt ? 5.3? t ? 0.8 ? ?2 ? 1求 xt ? 3 的95%的置信区间预测值?等价形式(1?1.8B ? 0.8B2 )xt ? (1? 0.6B)?t xt ?1.8xt?1 ? 0.8xt?2 ? ?t ? 0.6?t?1?计算预测值? t (1) ? 1.8 xt ? 0.8 xt ?1 ? 0.6? t ? 5.46 x ? t (2) ? 1.8 x ? t (1) ? 0.8 xt ? 5.59 x ? t (3) ? 1.8 x ? t (2) ? 0.8 x ? t (1) ? 5.69 x计算置信区间?Green函数值 ?? 1 ? 1.8 ? 0.6 ? 1.2 ? ?? 2 ? 1.8? 1 ? 0.8 ? 1.36 方差2 Var[e(3)] ? (1 ? ? 12 ? ? 2 )? ?2 ? 4.2896??95%置信区间?t (3) ? 1.96 Var (e(3)) , x ?t (3) ? 1.96 Var (e(3)) ) (x ? (1.63,9.75)例5.6续:对中国农业实际国民收入指 数序列做为期10年的预测疏系数模型??ARIMA(p,d,q)模型是指d阶差分后自相关 最高阶数为p,移动平均最高阶数为q的 模型,通常它包含p+q个独立的未知系 ?1 , ? , ? p , ? 1 , ? , ? q 数: 如果该模型中有部分自相关系数 ? ,1 ? j ? p 或部分移动平滑系数 ? ,1 ? k ? q 为零,即原 模型中有部分系数省缺了,那么该模型 称为疏系数模型。jk疏系数模型类型?如果只是自相关部分有省缺系数,那么该疏系 数模型可以简记为 ARIMA(( p1 ,? , p m ), d , q)?p1 ,? , p m 为非零自相关系数的阶数??如果只是移动平滑部分有省缺系数,那么该疏 系数模型可以简记为 ARIMA( p, d , (q1 ,? , q n )) ? q1 ,? , q n 为非零移动平均系数的阶数 如果自相关和移动平滑部分都有省缺,可以简 记为 ARIMA(( p ,? , p ), d , (q ,? , q ))1 m 1 n例5.8?对1917年-1975年美国23岁妇女每万人 生育率序列建模一阶差分自相关图偏自相关图建模?定阶?ARIMA((1,4),1,0)1 ? (1 ? B ) xt ? 4 t 1 ? 0.26633 B ? 0.33597 B?参数估计?模型检验? ?模型显著 参数显著季节模型? ?简单季节模型 乘积季节模型简单季节模型?简单季节模型是指序列中的季节效应和 其它效应之间是加法关系xt ? S t ? Tt ? I t?简单季节模型通过简单的趋势差分、季 节差分之后序列即可转化为平稳,它的 模型结构通常如下?( B ) ? D ? xt ? ?t ?( B)d例5.9?拟合1962――1991年德国工人季度失业率序列差分平稳?对原序列作一阶差分消除趋势,再作4步差分消除季节 效应的影响,差分后序列的时序图如下白噪声检验延迟阶数 6 12 18? 2 统计量P值 &0.1 &0.000143.84 51.71 54.48差分后序列自相关图差分后序列偏自相关图模型拟合?定阶?ARIMA((1,4),(1,4),0)?参数估计1 (1 ? B )(1 ? B ) xt ? ? 4 t 1 ? 0.44746 B ? 0.28132 B4模型检验残差白噪声检验 延迟 阶数 6 12? 2统 计量参数显著性检验 待估 参数P值 0.7191t统计量 5.48P值 &0.00012.09?1?410.99 0.3584-3.41 &0.0001拟合效果图乘积季节模型?使用场合?序列的季节效应、长期趋势效应和随机波动之间有着复 杂地相互关联性,简单的季节模型不能充分地提取其中 的相关关系 短期相关性用低阶ARMA(p,q)模型提取 季节相关性用以周期步长S为单位的ARMA(P,Q)模型提取 假设短期相关和季节效应之间具有乘积关系,模型结构 如下?构造原理? ? ??( B ) ? S ( B ) ? ? xt ? ?t ?( B) ? S ( B)d D S例5.10 :拟合1948――1981年美国女性月度失业率序列差分平稳?一阶、12步差分差分后序列自相关图差分后序列偏自相关图简单季节模型拟合结果拟合模型残差白噪声检验 延迟 阶数6 12 结果 AR(1,12) MA(1,2,12)ARMA((1,12),(1,12)? 2值P值? 2值9.5P值? 2值P值14.58 0.00570. 0.000416.42 0. 0. 0.0213 拟合模型均不显著乘积季节模型拟合?模型定阶?ARIMA(1,1,1)×(0,1,1)12?参数估计1 ? 0.66137 B ??12 xt ? (1 ? 0.77394 B12 )? t 1 ? 0.78978 B模型检验残差白噪声检验 延迟 阶数 6 12 18 结果?2统 计量参数显著性检验 待估 参数? 2统 计量P值 0.2P值4.50 9.42?1?12 ?1-4.66 &0. &0. &0.0001 参数均显著20.58 0.1507 模型显著乘积季节模型拟合效果图5.3 Auto-Regressive模型?构造思想?首先通过确定性因素分解方法提取序列中主 要的确定性信息xt ? Tt ? S t ? ? t?然后对残差序列拟合自回归模型,以便充分 提取相关信息? t ? ?1? t ?1 ? ? ? ? p ? t ? p ? a tAuto-Regressive模型结构? xt ? Tt ? S t ? ? t ? ?? t ? ?1? t ?1 ? ? ? ? p ? t ? p ? at ? 2 ? E (at ) ? 0,Var (at ) ? ? , Cov(at , at ?i ) ? 0, ?i ? 1对趋势效应的常用拟合方法?自变量为时间t的幂函数Tt ? ? 0 ? ? 1 ? t ? ? ? ? k ? t k ? ? t?自变量为历史观察值Tt ? ? 0 ? ?1 ? xt ?1 ? ? ? ? k ? xt ? k ? ? t对季节效应的常用拟合方法?给定季节指数S t ? S t??建立季节自回归模型Tt ? ? 0 ? ? 1 ? xt ? m ? ? ? ? l ? xt ?lm例5.6续? ?使用Auto-Regressive模型分析1952年-1988年中国农业 实际国民收入指数序列。 时序图显示该序列有显著的线性递增趋势,但没有季节 效应,所以考虑建立如下结构的Auto-Regressive模型? xt ? Tt ? ? t , t ? 1,2,3,? ? ?? t ? ?1? t ?1 ? ? ? ? p ? t ? p ? at ? 2 E ( a ) 0 , Var ( a ) ? , Cov(at , at ?i ) ? 0, ?i ? 1 ? ? t t ?趋势拟合?方法一:变量为时间t的幂函数Tt ? 66.1491 ? 4.5158t , t ? 1,2,3, ??方法二:变量为一阶延迟序列值 xt ?1? t ? 1.0365 xt ?1 x, t ? 1,2,3, ?趋势拟合效果图残差自相关检验?检验原理?回归模型拟合充分,残差的性质E (? t , ? t ? j ) ? 0 , ?j ? 1?回归模型拟合得不充分,残差的性质E (? t , ? t ? j ) ? 0 , ? j ? 1Durbin-Waston检验(DW检验)?假设条件?原假设:残差序列不存在一阶自相关性H 0 : E (? t , ? t ?1 ) ? 0 ? H 0 : ? ? 0?备择假设:残差序列存在一阶自相关性H 0 : E (? t , ? t ?1 ) ? 0 ? H 0 : ? ? 0DW统计量?构造统计量DW ?2 ( ? ) ? ? ? t t ?1 t ?2 2 ? ? t t ?1 n n?DW统计量和自相关系数的关系DW ? 2?1 ? ? ?DW统计量的判定结果正 相 关相 关 性 待 定不相关相 关 性 待 定负 相 关0dLdU2 4 ? dU 4 ? d L4例5.6续?检验第一个确定性趋势模型xt ? 66.1491 ? 4.5158t ? ? t , t ? 1,2,3, ?残差序列的自相关性。DW检验结果?检验结果DW统计 量的值 0.1378dL1.42dU1.53P值 0.0001?检验结论?检验结果显示残差序列高度正自相关。Durbin h检验?DW统计量的缺陷?当回归因子包含延迟因变量时,残差序列的 DW 统计量是一个有偏统计量。在这种场合 下使用 DW 统计量容易产生残差序列正自相 关性不显著的误判?Durbin h检验n Dh ? DW 2 1 ? n? ?例5.6续?检验第二个确定性趋势模型xt ? 1.0365 xt ?1 ? ? t , t ? 1,2,3, ?残差序列的自相关性。Dh检验结果?检验结果Dh统计量的值 2.8038 P值 0.0025?检验结论?检验结果显示残差序列高度正自相关。残差序列拟合? ? ?确定自回归模型的阶数 参数估计 模型检验例5.6续?对第一个确定性趋势模型的残差序列? t ? xt ? Tt ? xt ? 66.1491 ? 4.5158t , t ? 1,2,?进行拟合残差序列自相关图残差序列偏自相关图模型拟合?定阶?AR(2) 极大似然估计?参数估计方法??最终拟合模型口径? xt ? 69.1491 ? 4.5158t ? ? t ? ?? t ? 1.4859? t ?1 ? 0.5848? t ?2 ? at例5.6?第二个Auto-Regressive模型的拟 合结果? xt ? 1.033xt ?1 ? ? t ? ?? t ? 0.4615? t ?1 ? at三个拟合模型的比较模型ARIMA(0,1,1)模型:(1 ? B) xt ? 4.99661 ? (1 ? 0.70766 B)? tAIC249.4 250.6317SBC252.1 253.7987Auto-Regressive模型一: ? xt ? 69.1491 ? 4.5158t ? ? t ? ?? t ? 1.4859? t ?1 ? 0.5848? t ? 2 ? at? xt ? 1.033xt ?1 ? ? t ? ?? t ? 0.4615? t ?1 ? atAuto-Regressive模型二:5.4 异方差的性质?异方差的定义?如果随机误差序列的方差会随着时间的变化 而变化,这种情况被称作为异方差 Var (? t ) ? h(t ) 忽视异方差的存在会导致残差的方差会被严 重低估,继而参数显著性检验容易犯纳伪错 误,这使得参数的显著性检验失去意义,最 终导致模型的拟合精度受影响。?异方差的影响?异方差直观诊断? ?残差图 残差平方图残差图?方差齐性残差图?递增型异方差残差图残差平方图?原理?残差序列的方差实际上就是它平方的期望。Var (? t ) ? E (? t2 )?所以考察残差序列是否方差齐性,主要是考 察残差平方序列是否平稳例5.11?直观考察美国1963年4月――1971年7月短 期国库券的月度收益率序列的方差齐性。一阶差分后残差图一阶差分后残差平方图异方差处理方法??假如已知异方差函数具体形式,进行方 差齐性变化 假如不知异方差函数的具体形式,拟合 条件异方差模型5.5 方差齐性变换?使用场合?序列显示出显著的异方差性,且方差与均值之间具 有某种函数关系h(?) 是某个已知函数 其中:?? t2 ? h( ? t )处理思路?尝试寻找一个转换函数 g (?),使得经转换后的变量 满足方差齐性 Var[ g ( xt )] ? ? 2转换函数的确定原理?转换函数g ( xt )在 ? t 附近作一阶泰勒展开 g ( xt ) ? g ( ? t ) ? ( xt ? ? t ) g ?( ? t ) 求转换函数的方差Var[ g ( xt )] ? Var[ g ( ?t ) ? ( xt ? ?t ) g ?( ?t )] ? [ g ?( ?t )]2 h( ?t )??转换函数的确定g ?( ? t ) ?1 h( ? t )常用转换函数的确定?假定? t ? ? t ? h( ? t ) ? ??2 t转换函数的确定? g ( ?t ) ? log(?t )1 1 g ?( ?t ) ? ? h( ? t ) ? t例5.11续??对美国1963年4月――1971年7月短期国 库券的月度收益率序列使用方差齐性变 换方法进行分析 假定? t ? xt?函数变换y t ? log( xt )对数序列时序图一阶差分后序列图白噪声检验延迟阶数 6 12 18 LB统计量 3.58 10.82 21.71 P值 0.1 0.2452拟合模型口径及拟合效果图? log( xt ) ? ? t5.6 条件异方差模型? ? ?ARCH模型 GARCH模型 GARCH模型的变体? ? ? ?EGARCH模型 IGARCH模型 GARCH-M模型 AR-GARCH模型ARCH模型?假定?ARCH(q)模型结构? ? x ? f (t , x , x , ? ) ? ? t ?1 t ?2 t ? t ? ?? t ? ht et ? q ?h ? ? ? ? ? ? 2 t j t? j ? j ? 1 ??t?ht ~ N (0,1)原理?通过构造残差平方序 列的自回归模型来拟 合异方差函数GARCH 模型结构?使用场合??模型结构?ARCH模型实际上适 用于异方差函数短期 自相关过程 GARCH模型实际上适 用于异方差函数长期 自相关过程? ? x ? f (t , x , x ,? ) ? ? t ?1 t ?2 t ? t ? ?? t ? ht et ? p q ?h ? ? ? ? ? h ? ? ? ? 2 t i t ?i j t? j ? i ?1 j ?1 ?GARCH模型的约束条件?参数非负? ? 0,? i ? 0, ? j ? 0?参数有界?? ? ? ?i ?1 i j ?1pqj?1EGARCH模型? xt ? f (t , xt ?1 , xt ? 2 ,? ) ? ? t ? ?? t ? ht et ? p q ? ?ln(ht ) ? ? ? ?? i ln(ht ?i ) ? ? ? j g (et ) i ?1 j ?1 ? ? g (et ) ? ?et ? ? [ et ? E et ] ?IGARCH模型? xt ? f (t , xt ?1 , xt ? 2 , ? ) ? ? t ? ?? t ? ht et ? q p ? 2 h h ? ? ? ? ? ? ? ? t ? ? j t? j i t ?i j ?1 i ?1 ? ? p q ???i ? ? ? j ? 1 ? j ?1 ? i ?1GARCH-M模型? ? x ? f (t , x , x , ? ) ? ? h ? ? t ?1 t ?2 t t ? t ? ?? t ? ht et ? p q ?h ? ? ? ?? h ? ? ? ? 2 t i t ?i j t? j ? i ?1 j ?1 ?AR-GARCH模型? xt ? ?? t ? ? ? ?? t ? ?ht ? ? ? f (t , xt ?1 , xt ? 2 ,? ) ? ? t ? ? ? k ? t ?k ? ? tk ?1 m? ht et ? ? ? ?? i ht ?i ? ? ? j?i ?1 j ?1 p q 2 t? jGARCH模型拟合步骤? ? ? ? ? ?回归拟合 残差自相关性检验 异方差自相关性检验 ARCH模型定阶 参数估计 正态性检验例5.12?使用条件异方差模型拟合某金融时间序列。回归拟合?拟合模型xt ? ?1 xt ?1 ? ? t?参数估计?1 ? 1.0053 ??参数显著性检验?P值&0.0001,参数高度显著残差自相关性检验?残差序列DW检验结果? ?Durbin h=-2.6011Pr( Dh ? ?2.6011) ? 0.0046?拟合残差自回归模型? ?方法:逐步回归 模型口径? t ? ?0.1559? t ?1 ? 0.407? t ? 2 ? ? t异方差自相关检验? ?Portmantea Q检验 拉格朗日乘子(LM)检验Portmantea Q检验?假设条件 检验统计量Q(q) ? n(n ? 2)?qH 0 : ?1 ? ? 2 ? ? ? ? q ? 0 ? H 1 : ?1 , ? 2 , ? , ? q不全为零??i2n?i~ ? 2 (q ? 1)?检验结果? ?i ?1拒绝原假设 接受原假设Q (q ) ? ? 12?? ( q ? 1) Q (q ) ? ?12?? (q ? 1)??LM检验?假设条件 检验统计量 检验结果? ?H 0 : ?1 ? ? 2 ? ? ? ? q ? 0 ? H 1 : ?1 , ? 2 , ? , ? q不全为零?2 2 ? ?12 ? 2 ? ?q LM ( q ) ? W ?W , W ? ? 2 , 2 , ? , 2 ? ?? ? ? ? ? ? ? ? ??拒绝原假设 接受原假设Q (q ) ? ? 12?? ( q ? 1) Q (q ) ? ?12?? (q ? 1)??例5.12残差序列异方差检验ARCH模型拟合? ? ?定阶:GARCH(1,1) 参数估计:极大似然估计 拟合模型口径:AR(2)-GARCH(1,1)? xt ?? ? t ? ?? t ?h ? t ? 1.0046 xt ?1 ? ? t ? ?0.1559? t ?1 ? 0.407? t ? 2 ? ? t ? ht et ? 0.0951 ? 0.8999ht ?1 ? 0.1053? t2?1模型检验? ? ?检验方法:正态性检验 假设条件: H : u ~ N (0,1) ? H : u not 0 0 t t ~ N (0,1) 检验统计量 n 2 n Tn ? b1 ? (b2 ? 3) 2 ~ ? 2 (2) 6 24 检验结果? ??2 拒绝原假设 Tn ? ? 1?? ( 2) 接受原假设 Tn ? ?12?? ( 2)例5.13正态性检验结果Tn ? 1.1585P值=0.5603?AR(2)-GARCH(1,1)模型显著成立拟合效果图第六章多元时间序列分析本章结构? ? ? ? ?平稳时间序列建模 虚假回归 单位根检验 协整 误差修正模型6.1 平稳时间序列建模?ARIMAX模型结构? i ( B ) li ? ? y t ? ? ? ? ? ( B) B xit ? ? t ? k ?1 i ? ?? ? ?( B) a t t ? ?( B) ?k例6.1?在天然气炉中,输入的是天然气,输出 CO2 的输出浓度与天然气的输入 的是 CO2 , 速率有关。现在以中心化后的天然气输 入速率为输入序列,建立 CO2 的输出百分 浓度模型。输入/输出序列时序图?输入序列?输出序列一元分析?拟合输入序列at xt ? ?0.1228 ? 1 ? 1.97607 B ? 1.37499B 2 ? 0.34336B 3?拟合输出序列y t ? 53.90176 ? at 1 ? 3.10703B ? 1.34005B 2 ? 0.21274B 4多元分析?协相关图拟合回归模型?模型结构? 0-?1 B-? 2 B 2 3 yt ? ? ? B xt ? ? t 1 ? ?1 B?模型口径? 0.5648 ? 0.42573B ? 0.29964 B 2 3 y t ? 53.32256 ? B xt ? ? t 1 ? 0.60057 B拟合残差序列?偏自相关图?残差拟合模型1 ?t ? at 2 1 ? 1.53B ? 0.64 B拟合模型模型结构at 一元 yt ? 53.9 ? 2 4 1 ? 3 . 1 B ? 1 . 3 B ? 0 . 2 B 模型比较 AIC=196.3 SBC=211.1 AIC=8.3 SBC=34.0多元 模型? ? 0.54 ? 0.38B ? 0.52 B 2 3 y ? 53.26 ? B xt ? ? t ? ? t 1 ? 0.55B ? 1 ?? ? a t 2 t ? 1 ? 1.53B ? 0.64 B ?ARIMAX模型拟合效果图6.2 虚假回归?假设条件H 0 : ?1 ? 0 ? H 1 : ?1 ? 0?检验统计量 虚假回归?1 t? ???Pr{ t ? t? 2 (n) 非平稳序列} ? ?6.3 单位根检验?定义?通过检验特征根是在单位圆内还是单位圆上 (外),来检验序列的平稳性 DF检验 ADF检验 PP检验?方法? ? ?DF检验?假设条件? ?原假设:序列非平稳 H 0:?1 备择假设:序列平稳 H 0:?1?1 ?1?检验统计量? ?? ?? ? ? ?? N (0,1) ?1 ? 1 时 t (?1 ) ? 1 ? 1 ?渐近 S (?1 ) ?1 ? 1 时??? ?1 ? 1 ?) S (? 1DF统计量??1 ? 1 时? ?? ? t (?1 ) ? 1 1 ?渐近 ? ?? N (0,1) ?) S (? 1??1 ? 1 时W (r )dW (r ) ? ?? ?? ?? ?) S (? ? ?W (r )? dr极限 0 1 1 2 0? ?1 ? 11DF检验的等价表达?等价假设H 0:? ? 0 ? H1:? ? 0 其中:? ? ?1 ? 1?检验统计量? ? ?? ?) S (?DF检验的三种类型?第一种类型xt ?? 1xt ?1 ? ? t?第二种类型xt ? ? ?? 1xt ?1 ? ? t?第三种类型xt ? ? ? ? t ?? 1xt ?1 ? ? t例6.2?对1978年-2002年中国农村居民家庭人 均纯收入对数序列{ln xt }和生活消费支出 对数序列{ln yt } 进行检验例6.2 时序图例6.2 输入序列的DF检验例6.2 输出序列的DF检验ADF检验?DF检验只适用于AR(1)过程的平稳性检 验 。为了使检验能适用于AR(p)过程的 平稳性检验,人们对检验进行了一定的 修正,得到增广检验(Augmented Dickey-Fuller),简记为ADF检验ADF检验的原理?若AR(p)序列有单位根存在,则自回归系 数之和恰好等于1? p ? ?1? p ?1 ? ? ? ? p ? 0? ?1?1 ? ? ??11?? ??p ? 0? ?2 ? ? ? ? p ? 1ADF检验?等价假设H 0:? ? 0 ? H1:? ? 0 其中:? ? ?1 ? ?2 ? ? ? ? p ? 1?检验统计量? ? ?? ?) S (?ADF检验的三种类型?第一种类型xt ?? 1xt ?1 ? ? ?? pxt ? p ? ? t?第二种类型xt ? ? ?? 1xt ?1 ? ? ?? pxt ? p ? ? t?第三种类型xt ? ? ? ?t ?? 1xt ?1 ? ? ?? pxt ? p ? ? t例6.2续?对1978年-2002年中国农村居民家庭人 均纯收入对数差分后序列 ?? ln xt ?和生活消 费支出对数差分后序列?? ln yt ?进行检验例6.2 ?? ln x ?序列的ADF检验t例6.2 ?? ln y ?序列的ADF检验tPP检验???ADF检验主要适用于方差齐性场合,它对 于异方差序列的平稳性检验效果不佳 Phillips和 Perron于1988年对ADF检验进 行了非参数修正,提出了PP检验统计量。 PP检验统计量适用于异方差场合的平稳 性检验,且服从相应的ADF检验统计量的 极限分布PP检验统计量2 2 2 2 ?2 ? ? Sl ? Sl ? 2 )T ? ? Sl Z (? ) ? ? (? ) ? (1 2)(? ?? ( x ? x ) ? t ?1 T ?1 t ?2 T?其中:? ?T (1) ?2 ?1 2 ? ? ? t t ?1TT? ?T (2) ?2 Sl?1? ??t ?12 t? 2T?1? w (l ) ? ?? ??j ?1 j t ? j ?1lTt t? j1 T ?1 xt (3) xT ?1 ? ? T ? 1 t ?1例6.2续?对1978年-2002年中国农村居民家庭人 均纯收入对数差分后序列 ?? ln xt ?和生活 消费支出对数差分后序列 ?? ln yt ?进行PP 检验例6.2 ?? ln x ?序列的pp检验t例6.2 ?? ln y ?序列的PP检验t例6.2 二阶差分后序列的PP检验6.4 协整?单整的概念???如果序列平稳,说明序列不存在单位根,这 时称序列为零阶单整序列,简记为 xt ~ I (0) 假如原序列一阶差分后平稳,说明序列存在 一个单位根,这时称序列为一阶单整序列, 简记为 xt ~ I (1) 假如原序列至少需要进行d阶差分才能实现 平稳,说明原序列存在d个单位根,这时称 原序列为阶单整序列,简记为xt ~ I (d )单整的性质?若xt~ I (0) ,对任意非零实数a,b,有a ? bxt ~ I (0)?若xt ~ I (d ) ,对任意非零实数a,b,有 a ? bxt ~ I (d ) 若xt~ I ( 0) , y t ~ I ( 0)?对任意非零实数a,b,有z t ? axt ? by t ~ I (0)?若xt ~ I (d ) ,yt ~ I (c) 对任意非零实数a,b,有z t ? axt ? by t ~ I (k )k ? max[ d , c ]协整的概念?假定自变量序列为 {x1},? ,{x k } ,响应变量 序列为 { y t },构造回归模型y t ? ? 0 ? ? ? i xit ? ? ti ?1 k假定回归残差序列 ?? t ? 平稳,我们称响应 序列 { y t }与自变量序列 {x1},? ,{xk } 之间具有 协整关系。协整检验?假设条件??原假设:多元非平稳序列之间不存在协整关 系 H 0 : ? t ~ I (k ), k ? 1 备择假设:多元非平稳序列之间存在协整关 H 1 : ? t ~ I (0) 系 建立响应序列与输入序列之间的回归模型 对回归残差序列进行平稳性检验?检验步骤? ?例6.2续?对1978年-2002年中国农村居民家庭人 均纯收入对数序列 ln{xt }和生活消费支出 对数序列 ln{ yt }进行EG检验。构造回归模型?拟合模型?一元线性模型 最小二乘估计?估计方法??拟合模型口径ln y t ? 0.96832 ln xt ? ? t残差序列单位根检验我们可以以91.55%(1-0.0845)的把握断定残差 序列平稳且具有一阶自相关性 ? t ? ?1? t ?1 ? ? t最终拟合模型ln y t ? 0.96821ln xt ?i .i . d?t1 ? 0.83713Bvt ~ N (0,0.000893)误差修正模型??误差修正模型(Error Correction Model) 简称为ECM,最初由Hendry和Anderson 于1977年提出,它常常作为协整回归模 型的补充模型出现 协整模型度量序列之间的长期均衡关 系,而ECM模型则解释序列的短期波动 关系短期影响因素分析?响应序列的当期波动 ?yt 主要会受到三方 面短期波动的影响? ? ?输入序列的当期波动 ?xt 上一期的误差 ECM t ?1 纯随机波动 ? tyt ? yt ?1 ? ?xt ? ?xt ?1 ? ? t ?1 ? ? t ?yt ? ??xt ? ECM t ?1 ? ? t误差修正模型?y t ? ? 0 ?xt ? ? 1 ECM t ?1 ? ? t例6.2续?对1978年-2002年中国农村居民家庭人 均纯收入对数序列 ln{x } 和生活消费支出 对数序列 ln{ y } 构造ECM模型t t例6.2 构造ECM模型?拟合长期协整关系ln y t ? 0.96832 ln xt ? ? t?拟合短期波动(ECM模型)? ln y t ? 0.9579? ln xt ? 0.1537 ECM t ?1 ? ? t中 国 人 民 大 学 出 版 社 中国人民大学音像出版社地址:中国北京市中关村大街31号(100080) 电话:(010) 网址:
更多相关文档