9干预分析模型预测法9.1干预分析模型概述9.2单变量干预分析模型的识别与估计9.3干预分析模型的应用实例9.1干预分析模型概述一、干预模型简介干预的含义:时间序列经常会受到特殊事件及态势的影响,称这类外部事件为干预。研究干预分析的目的:从定量分析的角度来评估政策干预或突发事件对经济环境和经济过程的具体影响。二、干预分析模型的基本形式干预变量的形式:干预分析模型的基本变量是干预变量,有两种常见的干预变量。一种是持续性的干预变量,表示T时刻发生以后,一直有影响,这时可以用阶跃函数表示,形式是:)干预事件发生之后()干预事件发生之前(TtTtSTt,1,0第二种是短暂性的干预变量,表示在某时刻发生,仅对该时刻有影响,用单位脉冲函数表示,形式是:1,0,TttTPtT干预事件发生时()其它时间()干预事件的形式:干预事件虽然多种多样,但按其影响的形式,归纳起来基本上有四种类型:1.干预事件的影响突然开始,长期持续下去设干预对因变量的影响是固定的,从某一时刻T开始,但影响的程度是未知的,即因变量的大小是未知的。这种影响的干预模型可写为:TttSYω表示干预影响强度的未知参数。Yt不平稳时可以通过差分化为平稳序列,则干预模型可调整为:其中B为后移算子。如果干预事件要滞后若干个时期才产生影响,如b个时期,那么干预模型可进一步调整为:TttSYB)1(TtbtSBY2.干预事件的影响逐渐开始,长期持续下去有时候干预事件突然发生,并不能立刻产生完全的影响,而是随着时间的推移,逐渐地感到这种影响的存在。这种形式的最简单情形的模型方程为:更一般的模型是:10,1TttSBBY10,11TtrrbtSBBBY3.干预事件突然开始,产生暂时的影响这类干预现象可以用数学模型描述如下:当时,干预的影响只存在一个时期,当时,干预的影响将长期存在。10,1TtbtPBBY014.干预事件逐渐开始,产生暂时的影响干预的影响逐渐增加,在某个时刻到达高峰,然后又逐渐减弱以至消失。这类干预现象可用以下模型描绘:TtrrtPBBY1019.2单变量干预分析模型的识别与估计一、干预模型的构造与干预效应的识别单变量时间序列的干预模型,就是在时间序列模型中加进各种干预变量的影响。设平稳化后的单变量序列满足下述模型(ARMA):()()ttByaB又设干预事件的影响为:其中为干预变量,它等于或,则单变量序列的干预模型为:TttIBBZ)()(TtITtSTtPtTttaBBIBBy)()()()(tTtIB)()()()(BBBttaBB)()(,这里:二、干预效应的识别在对实际数据进行干预分析的过程中,一个主要的困难是,观察到的序列现实值是受到了干预变量影响的数据,不能保证自相关函数与偏自相关函数所反映的ARMA模型是真实的。下面我们介绍两种应对方法。(1)根据序列的具体情况和干预变量的性质进行识别确定干预变量的影响是短暂的还是长期的,需要进行具体的识别工作。它是利用干预变量产生影响之前或干预影响过后,也就是消除了干预影响或没有干预影响的净化数据,计算出自相关函数与偏自相关函数。首先识别ARMA模型中的p和q,然后估计出,中的参数。)(B)(B假定假定干预模型的模式为:)(1)(1BB)(1)(1BB0()()1TTttBBISBB组合这两个模型,便得到单变量序列的干预分析模型:或:01111()11()TtttBBxSaBB01111(1())(1)(1())1TtttBBBxSBaB(2)已知干预影响的情形假定在模型识别之前,对干预的影响已很清楚,以至于通过数据分析,能够确定干预变量的影响部分并估计出这部分的参数,然后计算出残差序列:这个序列是一个消除了干预变量影响的序列,可计算出它的自相关与偏自相关函数,从而识别出ARMA模型的阶数。)()(BBTtttIBBx)(ˆ)(ˆt三、干预模型的建模步骤1.利用干预影响产生前的数据,建立单变量的时间序列模型。然后利用此模型进行外推预测,得到的预测值,作为不受干预影响的数值。2.将实际值减去预测值,得到受干预影响的具体结果,利用这些结果求估干预影响的参数。3.利用排除干预影响后的全部数据,识别与估计出一个单变量的时间序列模型。4.求出总的干预分析模型。9.3干预分析模型的应用实例•例1我国国民收入增长的政策干预分析:现在采用按可比价格计算的国民收入指数来反映国民收入,研究其在1952~1993年间的增长模型。由于国民收入的增长一方面源于政策干预调节的影响,另一方面又包含自然增长的趋势,因此,把干预分析模型和一般的时间序列增长模型结合起来进行研究。已知1978年是我国一系列改革开放政策措施出台的开始,之后中国经济出现了呈加快增长的新形势,可以确定1978年为干预事件发生的开始时间,在建模中纳入政策变化等干预变量的影响。试确定干预分析模型。t123456789101112xt100114.0120.6128.3146.4153.0186.7202.0199.1140.0130.9144.9t131415161718192021222324xt168.8197.4231.0214.3200.3239.0294.6315.3324.3351.2355.2384.7t2526272829303132333435xt374.5403.7453.4485.1516.3541.5585.8644.2731.9830.6894.5t36373839404142xt985.71097.21133.41191.71283.41480.91704.6解答:(1)根据1952~1993年的数据建立一个干预模型如下:其中,t为自变量,xt为因变量,Zt表示干预事件对因变量的影响,它的确定是整个模型的关键。由于改革的影响是逐渐加强的,其作用又是长期深远的,因而干预变量可选取如下的形式:其中:tttZtbtbbx3210TttSBz1年及其后年前1978,11978,0TtS先对1952~1977年的国民收入指数建立时间序列模型,结果如下:该模型拟合度较好,可以通过参数的显著性检验和整个回归方程的显著性检验。301788.08744.72702.94ttxt20.972,278.084RF(2)在此基础上分离出干预影响的具体数值,求估干预模型的参数。用刚才的模型进行1978~1993年的国民收入指数的预测,然后用实际值减去预测值得到的差值就是改革所产生的干预值,记为Zt。求得具体数值见下表:t19781979198019811982198319841985Zt3.805.153.73-6.040.8319.2364.25117.49t19861987198819891990199119921993Zt133.04172.89229.94212.28209.60237.50354.96404.24利用上表数据可以估计出干预模型:的参数与实际上是自回归方程:的参数:TttSBz11ttzzˆˆ13.4012,1.11341ˆ1.113413.4012ttzz(3)计算净化序列净化序列是指消除了干预影响的序列即:T=27t=1,2,3…,42(4)对净化序列建立拟合模型ˆˆ=ˆ1TtttttyxzxSB396.59567.59250.0182tytt20.99322274.88RF该模型拟合度较好,可以通过参数的显著性检验和整个回归方程的显著性检验,因此模型是合理的。(5)组建干预分析模型经过以上各步的参数估计,可以组建最终的干预分析模型如下:(6)预测313.401296.59567.59250.018211.1134TttxttSB年及其后年前1978,11978,0TtS