计量经济学课件随机时间序列分析模型-

1、§9.2 随机时间序列分析模型一、时间序列模型的基本概念及其适用性二、随机时间序列模型的平稳性条件三、随机时间序列模型的识别四、随机时间序列模型的估计五、随机时间序列模型的检验经典计量经济学模型与时间序列模型确定性时间序列模型与随机性时间序列模型一、时间序列模型的基本概念及其适用性 1、时间序列模型的基本概念随机时间序列模型(time series modeling是指仅用它的过去值及随机扰动项所建立起来的模型,其一般形式为 Xt=F(Xt-1, Xt-2, , ?t 建立具体的时间序列模型,需解决如下三个问题: (1模型的具体形式 (2时序变量的滞后期 (3随机扰动项的结构例如,取

2、线性方程、一期滞后以及白噪声随机扰动项( ?t =?t,模型将是一个1阶自回归过程AR(1: Xt=?Xt-1+ ?t 这里, ?t特指一白噪声。一般的p阶自回归过程AR(p是Xt=?1Xt-1+ ?2Xt-2 + + ?pXt-p + ?t (* 将纯AR(p与纯MA(q结合,得到一个一般的自回归移动平均(autoregressive moving average过程ARMA(p,q:经典回归模型的问题:迄今为止,对一个时间序列Xt的变动进行解释或预测,是通过某个单方程回归模型或联立方程回归模型进行的,由于它们以因果关系为基础,且具有一定的模型结构,因此也常称为结构式模型(structura

3、l model。然而,如果Xt波动的主要原因可能是我们无法解释的因素,如气候、消费者偏好的变化等,则利用结构式模型来解释Xt的变动就比较困难或不可能,因为要取得相应的量化数据,并建立令人满意的回归模型是很困难的。有时,即使能估计出一个较为满意的因果关系回归方程,但由于对某些解释变量未来值的预测本身就非常困难,甚至比预测被解释变量的未来值更困难,这时因果关系的回归模型及其预测技术就不适用了。例如,时间序列过去是否有明显的增长趋势,如果增长趋势在过去的行为中占主导地位,能否认为它也会在未来的行为里占主导地位呢?或者时间序列显示出循环周期性行为,我们能否利用过去的这种行为来外推它的未来走向?随机时间

4、序列分析模型,就是要通过序列过去的变化特征来预测未来的变化趋势。使用时间序列分析模型的另一个原因在于:如果经济理论正确地阐释了现实经济结构,则这一结构可以写成类似于ARMA(p,q式的时间序列分析模型的形式。例如,对于如下最简单的宏观经济模型:这里,Ct、It、Yt分别表示消费、投资与国民收入。 Ct与Yt作为内生变量,它们的运动是由作为外生变量的投资It的运动及随机扰动项?t的变化决定的。上述模型可作变形如下:两个方程等式右边除去第一项外的剩余部分可看成一个综合性的随机扰动项,其特征依赖于投资项It的行为。如果It是一个白噪声,则消费序列Ct就成为一个1阶自回归过程AR(1,而收入序列Yt就

5、成为一个(1,1阶的自回归移动平均过程ARMA(1,1。二、随机时间序列模型的平稳性条件自回归移动平均模型(ARMA是随机时间序列分析模型的普遍形式,自回归模型(AR和移动平均模型(MA是它的特殊情况。关于这几类模型的研究,是时间序列分析的重点内容:主要包括模型的平稳性分析、模型的识别和模型的估计。考虑p阶自回归模型AR(pmoving average时间序列,记为ARIMA(p,d,q。例如,一个ARIMA(2,1,2时间序列在它成为平稳序列之前先得差分一次,然后用一个ARMA(2,2模型作为它的生成模型的。当然,一个ARIMA(p,0,0过程表示了一个纯AR(p平稳过程;一个ARIMA(0

6、,0,q表示一个纯MA(q平稳过程。三、随机时间序列模型的识别所谓随机时间序列模型的识别,就是对于一个平稳的随机时间序列,找出生成它的合适的随机过程或模型,即判断该时间序列是遵循一纯AR过程、还是遵循一纯MA过程或ARMA过程。所使用的工具主要是时间序列的自相关函数(autocorrelation function,ACF及偏自相关函数(partial autocorrelation function, PACF 。 1、AR(p过程 (1自相关函数ACF 1阶自回归模型AR(1 Xt=?Xt-1+ ?t 的k阶滞后自协方差为:Xt=?1Xt-1+ ?2Xt-2 + ?t 该模型的方差?0以及

7、滞后1期与2期的自协方差?1, ?2分别为一般地,p阶自回归模型AR(p Xt=?1Xt-1+ ?2Xt-2 + ?pXt-p + ?t 其中:1/zi是AR(p特征方程?(z=0的特征根,由AR(p平稳的条件知,|zi| 1; 因此,当1/zi均为实数根时,?k呈几何型衰减(单调或振荡;当存在虚数根时,则一对共扼复根构成通解中的一个阻尼正弦波项, ?k呈正弦波衰减。(2偏自相关函数从Xt中去掉Xt-1的影响,则只剩下随机扰动项?t,显然它与Xt-2无关,因此我们说Xt与Xt-2的偏自相关系数为零,记为在实际识别时,由于样本偏自相关函数rk*是总体偏自相关函数?k*的一个估计,由于样本的随机性

8、,当k p时,rk*不会全为0,而是在0的上下波动。但可以证明,当k p时,rk*服从如下渐近正态分布: rk*N(0,1/n 式中n表示样本容量。因此,如果计算的rk*满足对MA(1过程 MA(1过程可以等价地写成?t关于无穷序列Xt,Xt-1,的线性组合的形式:其自协方差系数为与MA(1相仿,可以验证MA(q过程的偏自相关函数是非截尾但趋于零的。ARMA(p,q的自相关函数,可以看作MA(q的自相关函数和AR(p的自相关函数的混合物。当p=0时,它具有截尾性质;当q=0时,它具有拖尾性质;当p、q都不为0时,它具有拖尾性质从识别上看,通常: ARMA(p,q过程的偏自相关函数(PACF可能

9、在p阶滞后前有几项明显的尖柱(spikes,但从p阶滞后项开始逐渐趋向于零;而它的自相关函数(ACF则是在q阶滞后前有几项明显的尖柱,从q阶滞后项开始逐渐趋向于零。四、随机时间序列模型的估计 AR(p、MA(q、ARMA(p,q模型的估计方法较多,大体上分为3类:(1最小二乘估计;(2矩估计;(3利用自相关函数的直接估计。下面有选择地加以介绍。 AR(p模型的Yule Walker方程估计利用实际时间序列提供的信息,首先求得自相关函数的估计值 MA(q模型的矩估计(1MA(1模型的直接算法对于MA(1模型,(*式相应地写成(2MA(q模型的迭代算法第二步,将第一次迭代值代入(*式,计算出第二次

10、迭代值 ARMA(p,q模型的矩估计第二步,改写模型,求?1,?2,?,?q以及?2的估计值 AR(p的最小二乘估计 为了与 AR(p模型的 Yule Walker 方程估计进行比较， 需要说明的是，在上述 将(*改写成： 解该方程组，得到： 模型的平稳性、识别与估计的讨论中，ARMA(p,q模型中均未包含常 数项。 五、模型的检验 由于 ARMA(p,q模型的识别与估计 是在假设随机扰动项是一白噪声的基础上进行的，因此，如果估计的 模型确认正确的话，残差应代表一白噪声序列。 如果通过所 估计的模型计算的样本残差不代表一白噪声， 则说明模型的识别与估 计有误，需重新识别与估计。 序列是否存在自

11、相关。 在实际检验时，主要检验残差 2 、 AIC 与 SBC 模 型 选 择 标 准 另外一个遇到的问题是，在实际识别 ARMA(p,q模型时，需多次反复 偿试，有可能存在不止一组（p,q）值都能通过识别检验。 显 然，增加 p 与 q 的阶数，可增加拟合优度，但却同时降低了自由度。 因此，对可能的适当的模型，存在着模型的“简洁性”与模型的拟合 优度的权衡选择问题。 其中，n 为待估参数个数（p+q+可能 存在的常数项） 为可使用的观测值，RSS 为残差平方和（Residual ，T sum of squares） 。 越好 在选择可能的模型时，AIC 与 SBC 越小 显然，如果添加的滞后项没有解释能力，则对 RSS 值的 减小没有多大帮助，却增加待估参数的个数，因此使得 AIC 或 SBC 的 值增加。 相同的时间段。 需注意的是：在不同模型间进行比较时，必须选取 由第一节知：中国支出法 GDP 是非平稳的， 可 记 但它的一阶差分是平稳的， 即支出法 GDP 是 I(1时间序列。 以对经过一阶差分后的 GDP 建立适当的 ARMA(p,q模型