VAR模型与向量VECM模型(7)

《VAR模型与向量VECM模型(7)》由会员分享，可在线阅读，更多相关《VAR模型与向量VECM模型(7)（20页珍藏版）》请在装配图网上搜索。

1、 第7章 向量自回归模型VAR与向量误差修正模型VEC 7.1 向量自回归模型VAR(p) 传统的经济计量学联立方程模型建摸方法, 是以经济理论为根底来描述经济变量之间的结构关系，采用的是结构方法来建立模型，所建立的就是联立方程结构式模型。这种模型其优点是具有明显的经济理论含义。但是，从计量经济学建摸理论而言，也存在许多弊端而受到质疑。一是在模型建立之处，首先需要明确哪些是生变量，哪些是外生变量，尽管可以根据研究问题和目的来确定，但有时也并不容易；二是所设定的模型，每一结构方程都含有生多个生变量，当将某一生变量作为被解释变量出现在方程左边时，右边将会含有多个其余生变量，由于它们与扰动项相关，

2、从而使模型参数估计变得十分复杂，在未估计前，就需要讨论识别性；三是结构式模型不能很好地反映出变量间的动态联系。为了解决这一问题，经过一些现代计量经济学家门的研究，就给出了一种非结构性建立经济变量之间关系模型的方法，这就是所谓向量自回归模型VAR模型不是以经济理论为根底描述经济变量之间的结构关系来建立模型的，它是以数据统计性质为根底，把某一经济系统中的每一变量作为所有变量的滞后变量的函数来构造模型的。它是一种处理具有相关关系的多变量的分析和预测、随机扰动对系统的动态冲击的最方便的方法。而且在一定条件下，多元MA模型、ARMA模型，也可化为VAR模型来处理，这为研究具有相关关系的多变量的分析和预测

3、带来很大方便。 VAR模型的一般形式 1、非限制性VAR模型高斯VAR模型,或简化式非限制性VAR模型设为一维随机时间序列，为滞后阶数，为一维随机扰动的时间序列，且有结构关系 711 假设引入矩阵符号，记 可写成 ， 712 进一步，假设引入滞后算子，那么又可表示成 7. 1. 3其中: ,为滞后算子多项式. 如果模型满足的条件: 参数阵特征方程 的根全在单位园外；，即 相互独立，同服从以为期望向量、为方差协方差阵的维正态分布。这时，是维白噪声向量序列，由于没有结构性经济含义，也被称为冲击向量；，即与与各滞后期不相关。那么称上述模型为非限制性VAR模型高斯VAR模型,或简化式非限制性VAR模型

4、。 2、受限制性VAR模型，或简化式受限制性VAR模型 如果将做为一维生的随机时间序列，受维外生的时间序列影响限制，那么VAR模型为 ， 714或利用滞后算子表示成 7. 1. 5 其中： 此时称该模型为受限制性VAR模型，简化式受限制性VAR模型。对于受限制性VAR模型，可通过对作OLS回归，得到残差估计，从而将变换成 ， 716 7. 1. 7这说明受限制性VAR模型可化为非限制性VAR模型。 简化式非限制、受限制VAR模型，皆简记为。 3、结构式非限制性VAR模型 如果中的每一分量受其它分量当期影响, 无维外生的时间序列影响限制,那么模型化为 ， 718或利用滞后算子表示成 7. 1.

5、9其中： ,这时的此时称该模型为结构式非限制性VAR模型。 如果可逆，既逆阵存在，那么结构式非限制性VAR模型可化为简化式非限制性VAR模型， 7110或利用滞后算子表示成 7. 1. 11 这时，其中的 4、结构式受限制性VAR模型 如果将做为一维生的随机时间序列，其中每一分量受其它分量当期影响,且还受维外生的时间序列影响限制，那么VAR模型为 ， 7112或利用滞后算子表示成 7. 1. 13此时称该模型为结构式受限制性VAR模型。如果可逆，既逆阵存在，那么结构式受限制性VAR模型可化为简化式受限制性VAR模型， 7114或利用滞后算子表示成 7. 1. 15这时，其中的结构式非限制、受限

6、制VAR模型，皆简记为。 简化式VAR模型的参数估计 VAR模型参数估计, 简化式VAR模型比拟简单可采用Yule-Walker估计、OLS估计、极大似然估计法等进展估计,且可获得具有良好统计性质的估计量。结构式VAR模型参数估计比拟复杂，可有两种途径：一种是化成简化式，直接估计简化式模型参数，然后再通过简化式模型参数与结构式模型参数的关系，求得结构式模型参数估计，但这存在一个问题是否可行，什么情况下可行，这与结构式模型的识别性有关。另一种途径是直接对结构式模型参数进展估计，但这也存在一个问题，上述方法不可应用，原因是每一方程含有众多生的与扰动项相关变量，那么，如何估计？这也与结构式模型的识别

7、性有关。对于简化式VAR模型，即 相互独立，同服从以为期望向量、为方差协方差阵的维正态分布。这时，是维白噪声向量序列的条件下，模型参数阵与也可采用Yule-Walker估计、OLS估计、极大似然估计。设,为长度为的样本向量 1、Yule-Walker估计 在充分大时, 首先估计自协方差阵 ()令 ,那么可得模型参数阵的Yule-Walker估计(矩估计)为 2、估计模型参数阵的OLS估计，即求使下的作为估计。 记 ()由此可推得 由此可见, 模型参数阵的OLS估计(但式中的的计算不同. 但是, 当充分大时,(在充分大时, 可直接采用Yule-Walker估计比拟简单方便。 而的估计为 其中：

8、3、极大似然估计可证明, 模型参数阵的极大似然估计与OLS估计完全等价。这里不在赘述。 简化式VAR模型的预测 在时，对的一步线性预测 其一步预测误差为 一步预测误差的方差阵为的估计为 在时，如果利用模型参数的估计量，对进展一步线性预测，那么的实际一步线性预测为 其一步预测误差为 一步预测误差的方差阵为的估计为 VAR模型阶数p确实定VAR模型的定阶是一个矛盾过程，阶数p确实定，既不能太大，又不能太小，必须兼顾。因为，一方面，希望滞后阶数p要大一些，以便使模型能更好地反映出动态特征，但另一方面，又不希望太大，否那么，阶数p太大，会造成需要估计的模型参数过多，而使模型自由度减少。因此，在定阶时需

9、要综合考虑，以既要有足够大的滞后项，又能有足够大的自由度为原那么确定阶数。VAR模型的定阶方法有多种：1、FPE准那么(最小最终预测误差准那么)FPE准那么(最小最终预测误差准那么)，即利用一步预测误差方差进展定阶。因为，如果模型阶数适宜，那么模型对实际数据拟合优度必然会高，其一步预测误差方差也必然会小；反之，那么相反。设给定时间序列向量长度为的样本向量为,，那么其一步预测误差方差阵的估计量为式,它是一个阶阵,因此可定义其最终预测误差为 显然, 是的函数。所谓最小最终预测误差准那么，就是分别取=1，2,，M, 来计算, 使值所对应的, 为模型适宜阶数。相应的模型参数估计为最正确模型参数估计。其

10、中，M为预先选定的阶数上界,一般取之间。 在实际计算过程中，可如下判断： 如果的值，随着从1开场逐渐增大就一直上升，那么可判定=1； 如果的值，随着从1开场逐渐增大就一直下降，那么可判定该随机时间序列不能用ARp模型来描述； 如果的值，在某一值下降很快，而后又缓慢下降，那么可判定该值为所确定的阶数； 如果的值，随着从1开场逐渐增大而上下剧烈跳动，难以找到最小值，这可能由于样本数据长度T太小造成的，应增大样本长度，重新进展定阶、估计模型参数，建立模型。利用FPE信息准那么还可以用来检验模型的建立是否可由局部分量，比方前个分量，来进展，方法如下：记式中的阶矩阵的左上角阶子方阵为, 那么前个分量，的

11、最终预测误差为 当时， 如果，那么可认为仅用前个分量，建立模型即可，没有必要采用维随机时间序列建立模型，因为从最小最终预测误差准那么角度，用维随机时间序列建立模型比仅采前个分量，建立模型，带来拟合优度的显著改善；反之，那么相反。2、AIC(Akaike Information Criterion)与SCBayes Information Criterion信息准那么AIC、SC信息准那么，也称最小信息准那么，定义 ， 其中：为模型需要估计参数个数，对,。所谓最小信息准那么，就是分别取=1，2, 来计算AIC或者SC, 使AIC或SC值所对应的, 为模型适宜阶数。相应的模型参数估计为最正确模型参

12、数估计。3、似然比检验法Likelihood Ratio,LR检验： 由于，即 相互独立，同服从以为期望向量、为方差协方差阵的维正态分布。因此，记，那么在给的条件下，的条件分布为 于是，在给的条件下，的联合分布密度，即似然函数为对数似然函数为 将参数估计代入，那么有， 又 因此，有 现在，欲检验假设 样本数据是由滞后阶数为的VAR模型生成；样本数据是由滞后阶数为的VAR模型生成 取似然比统计量为 分布 在给定的显著性水平下，当，那么拒绝，说明增加滞后阶数，可显著增大似然函数值；否那么，那么相反。检验在小样本下，可取似然比统计量为分布 其中,. VAR模型的Granger因果关系检验 VAR模型

13、的另一重要应用是可用来检验一个变量与另一变量间是否存在Granger因果关系，这也是建立VAR模型所需要的。1、 Granger因果关系的涵义 设为一维随机时间序列，如果在给定的滞后值下的条件分布与仅在给定的的滞后值下的条件分布一样，即 那么称对存在Granger非因果性关系，否那么，对存在Granger因果性关系。 Granger因果性关系涵义的另一表述：在其条件不变下，如果加上的滞后值，并不对只由的滞后值下对进展预测有显著改善，那么称对存在Granger非因果性关系，否那么，对存在Granger因果性关系。2、 Granger因果关系检验设为一维随机时间序列，为滞后阶数，为一维随机扰动的时

14、间序列，那么有2元VAR模型为 7131显然，欲检验对是否存在Granger非因果性关系，等价地，检验假设；中至少有一个不为0。其用于检验的统计量为 7132其中，为模型中第1方程残差平方和, 各期滞后项后拟合残差平方和。 在给定的显著性水平下，当时，拒绝。如果模型7131满足，即 相互独立，同服从以为期望向量、为方差协方差阵的维正态分布条件，那么 也可采用如下统计量进展检验 7133在给定的显著性水平下，当时，拒绝，上述Granger因果性关系检验，可推广到对任意维VAR模型以与SVAR模型中的某一或某几个随机时间序列包括生、外生变量是否对另一时间序列具有Granger因果性的检验上去。 7

15、.2 VARp模型的脉冲响应函数与方差分解在实际应用中，由于通常所设定的VAR模型都是非经济理论性的简化式模型，出它无需对变量作任何先验性约束，因此，在分析应用中，往往并不利用VAR模型去分析某一变量的变化对另一变量的影响如何，而是分析当某一扰动项发生变化，或者说模型受到某种冲击时，对系统的动态影响，这钟分析方法称为脉冲响应函数方法Impulse Response Function,IRF。 脉冲响应函数根本思想对VAR模型采用脉冲响应函数分析扰动项发生变化，或者说模型受到某种冲击时，对系统的动态影响，就是分析扰动项发生变化是如何传播到各变量的。设为一维随机时间序列，滞后阶数=2，为一维随机扰

16、动的时间序列，那么有2元VAR模型为 721扰动项满足白噪声假设条件，即； 现在假设上述VAR模型系统从时期开场运行，并设，在时给定扰动项 并且其后，即在时给定一脉冲，我们来讨论的响应。 由于 由721，在时，于是有，； 将上述结果再代入721，在时，于是有，； 再将上述结果代入1521，在2时，于是有，如此下去，可求得结果，称此结果为由的冲脉冲引起的的响应函数；所求得的，称为由的冲脉冲引起的的响应函数。 反过来，也可求得在时，给定扰动项并且其后，即在给定一脉冲时，由的冲脉冲引起的、的响应函数。 VAR模型的脉冲响应函数假设有VAR(p)模型 ， 722引入滞后算子，表示成 7.2. 3其中:

17、 ,为滞后算子多项式. 在满足特征方程 的根全在单位园外条件下，那么VAR(p)是可逆的，即可将表示成白噪声滑动和形式 7.2. 4其中：阶单位阵 7.2. 4中第方程为 7. 2. 5 当时, ()为 726现在假定在基期给一个单位脉冲, 即 而 那么可求得由的脉冲引起的响应函数为： 由此可看出，对于7.2. 4式的一般情形，由的脉冲引起的响应函数为：由的脉冲引起的累积响应函数为： 由7. 2. 4式, 其中的 中的第行、第列元素可表示为 7. 2. 7作为的函数，它描述了在时期，其他变量和早期变量不变的情况下，对的一个冲击的反响，称为脉冲响应函数。 用矩阵可表示为 = 7. 2. 8即 中

18、的第行、第列元素等于时期的第变量扰动项增加一个单位，其它时期扰动项为常数时，对时期的第个变量值的影响。723 方差分解VAR模型的脉冲响应函数是用来描述VAR模型中一个生变量的冲击给其它生变量所带来的影响的，它是随时间的推移，观察模型中各变量对于冲击是如何反响的。而方差分解是要通过分析每一结构冲击对生变量变化通常用方差来度量的奉献度，进一步评价不同结构冲击的重要性的，与脉冲响应函数相比，方差分解是一种比拟粗糙的把握变量间关系的方法，它给出的是对VAR模型中的变量产生影响的每个扰动项的相对重要信息。方差分解的根本思想是：由7. 2. 5式 7. 2. 9可知，左边括号为是第扰动项从过去无限远至现

19、在时点对第生变量影响的总和。在，无序列相关的假设下，对其求方差，可得 7. 2. 10它是把第扰动项从过去无限远至现在时点对第生变量影响总和，用方差加以评价的结果。如果为对角阵，那么的方差为 7. 2. 11由此可知，的方差可分解成个不相关的的影响。 由此，可测定出各个扰动项对方差的相对方差奉献率为 7. 2. 12 在实际应用计算中，不可能从过去无限远的来评价。在模型满足平稳性条件下，由于随着的增大是按几何级数衰减的，故只要取前有限项计算即可。其近似相对方差奉献率为 ， 7. 2. 13 有如下性质： 7. 2. 14 7. 2. 15如果大，那么意味着第变量第扰动项对第变量影响大，反之，那

20、么相反。7.3 Johansen协整检验与向量误差修正模型(VEC) 前面我们已经介绍了单方程的协整检验与误差修正模型。且其协整检验方法是以回归模型为根底的基于回归残差序列的ADF检验法进展检验的。现在我们把它推广到VAR模型上去，并给出以VAR模型为根底基于回归系数的协整检验方法。 在单方程协整检验中，由于是基于回归残差序列进展，故在第一阶段需要采用OLS进展回归分析，应用很不方便。为此，Johansen1988与Juselius(1990)提出了一个以VAR模型为根底的基于回归系数的特别适合于多变量的协整检验法。731 Johansen协整检验 1、协整定义：设为一维随机时间序列，如果 且

21、每一， 存在非零向量=，使那么称为协整，记为，为协整向量。假设为协整，那么最多存在个线性无关的协整向量。即假设记由的所有协整向量组成的矩阵为，那么秩，。例如，=2，假设有使，按照上述，最多存在个线性无关的协整向量，那么协整向量唯一。因为假设有 这与矛盾，故，即唯一。2、Johansen协整检验根本思想设为一维随机时间序列，且 即每一，受维外生的时间序列影响限制，那么首先可建立VAR模型 ， 731将上式进展差分变换，也称为协整变换，可写成 732其中， 733在732中，由于 所以、，因此，只要 那么，亦即之间具有协整关系，而之间是否具有协整关系取决于阶矩阵的秩。因为，与模型全部参数阵有关，故

22、称为压缩矩阵影响矩阵。设，那么有3种情况：如果，这意味着是一列满秩阵，那么只有当时，才能保证 但这与相矛盾，故如果，那么 由732，这时用不着讨论之间是否具有有协整关系。除上述两种极端情形外，一般情况是：如果，这意味着中一定存在个协整关系协整组合，其余个关系仍然为关系。在这种情况下，可将分解成两个阶阵的乘积 且、。 将其代入到742式中，有 734上式要求，向量，其每一行都是变量，即的每一列都是一协整向量， 所以决定了之间协整向量的个数和形式，故称称为协整向量阵，为协整向量个数。的每一行是出现在上述每一方程中的个协整组合的一组权数，故称为调整参数阵，或修正参数阵。显然，在假定条件下，最大可能，

23、这就是对于维向量最大可能存在个线性无关的协整向量的道理。根据上述分析，可知欲检验是否具有协整关系，就转化为对矩阵的秩数的检验，由于=的非零特征根的个数，因此，就可以通过检验的非零特征根的个数，来检验，从而来判定是否具有协整关系。这就是Johansen协整检验的根本思想。3、 Johansen协整检验现在假设的个特征根为。Johansen协整检验有两种方法： 1、特征根迹检验trace检验由于个最大特征根可得到个协整向量，而对于其余个非协整组合而言，应该有，因此，检验是否等于，等价地检验假设可用于检验的特征根迹统计量为 735具体显著性检验程序如下：当某一显著性水平下的Johansen分布临界值

24、，即不显著时，承受，说明有个特征根，0个协整向量，即不存在协整关系。当某一显著性水平下的Johansen分布临界值，即显著时，拒绝，说明至少有1协整向量。这时必须接着检验。当某一显著性水平下的Johansen分布临界值，即不显著时，承受，说明只有1个协整向量。依次进展下去，直到承受，说明存在个协整向量时为止。这时，这个协整向量就是最大的个特征根所对应的经过正规化的特征向量。显然整个检验过程应该是序贯进展的，整个序贯检验过程如下：当某一显著性水平下的Johansen分布临界值，即不显著时，承受，说明只有0个协整向量即不存在协整关系。当某一显著性水平下的Johansen分布临界值，即显著时，拒绝，

25、说明至少有1协整向量。这时必须接着检验。当某一显著性水平下的Johansen分布临界值，即不显著时，承受，说明只有1个协整向量。当某一显著性水平下的Johansen分布临界值，即显著时，拒绝，说明只少2个协整向量。当某一显著性水平下的Johansen分布临界值，即不显著时，承受，说明只有个协整向量。2、最大特征根检验由于个最大特征根可得到个协整向量，而对于其余个非协整组合而言，应该有，因此，最大特征根检验用于检验假设 用于检验的最大特征根检验的统计量为 736具体显著性检验程序如下：当临界值，不显著时，承受，说明最大特征根为0，无协整向量；当临界值，显著时，拒绝，承受，说明至少有1个最大特征根

26、不为0，至少有1个协整向量。须接着检验。当临界值，不显著时，承受，说明最大特征根不为0，其余特征根皆为0，只有1个协整向量；检验截止。当临界值，显著时，拒绝，承受，说明至少有两个最大特征根不为0，至少有2个协整向量。须接着检验。依次进展下去，直到承受，共有个协整向量时为止。4、协整方程形式732 向量误差修正模型(VEC)由式可知，设为一维随机时间序列，且 即每一，如果不受维外生的时间序列影响限制，VAR模型变为 ， 737将上式进展协整变换，可写成 738其中， 739 如果存在协整关系，那么的这时可写成 (7310其中, 即为误差修正项, 反映的是变量之间的长期均衡关系。即，上式可写成 (

27、7311()即为向量误差修正模型VEC,其中每一方程都是一个误差修正模型ECM。VEC模型中的参数向量，反映的是变量之间的均衡关系偏离长期均衡状态时，将其调整到均衡状态的调整速度，故称其为调整参数阵，或修正参数阵。所有作为解释变量的差分项的系数向量，反映的是各变量的短期波动对作为被解释变量的短期变化的影响。在实际应用中，对于影响不显著的那些短期波动的项可以从模型中剔除。上述只是讨论了简单的VEC模型，我们也可以象VAR模型那样构造结构式VEC模型，也可以对VEC模型讨论Granger因果关系检验、脉冲响应函数和方差分解等等。关于这些更详细的容，可参见Davidson和Mackinnon(199

28、3)以与汉蜜而顿1999的著作。詹姆。汉密尔顿：时间序列分析明志译，中国社会科学，1999，第19章。7.4 SVAR(p)模型 741 SVAR模型的识别与约束条件如果中的每一分量受其它分量当期影响, 无维外生的时间序列影响限制,那么由718式，结构式非限制性SVARp模型为 ， 741或利用滞后算子表示成 7. 4. 2其中： , 这时的此时称该模型为结构式非限制性SVAR模型。 结构式非限制性SVAR模型，即使在扰动项满足白噪声条件下也不能采用普通最小二乘法估计模型参数来建立模型，因为每一方程含有同期相关的变量。 如果可逆，既逆阵存在，那么结构式非限制性SVAR模型可化为简化式非限制性V

29、AR模型， 743或利用滞后算子表示成 7. 4. 4这时，其中的假设记 745那么7. 3. 4可写成， 746简化式非限制性模型VAR所含需要估计参数个数为 747其中，为扰动项的方差协方差阵所含未知待估计参数个数。在扰动项满足白噪声条件下，746式可采用普通最小二乘法估计上述模型参数，来建立其简化式非限制性VAR模型。我们知道，结构式非限制性SVAR模型741，即使在扰动项满足白噪声条件下也不能采用普通最小二乘法估计模型参数来建立模型，因为每一方程含有同期相关的变量。既然其简化式非限制性VAR模型746模型参数可以通过普通最小二乘法估计，那么，可否根据上述简化式非限制性VAR模型的模型参

30、数与结构式非限制性SVAR模型的模型参数之间的关系式745，通过已估计的简化式非限制性VAR模型参数，得到相应的结构式非限制性SVAR模型参数建立模型？这就涉与到结构式非限制性SVAR模型741的识别性关于识别性与其方法，可见14章联立方程容，或者说取决于对结构式非限制性SVAR模型所施加的约束条件。因为，由结构式非限制性SVAR模型741可知，其需要估计的模型参数个数共 748，所以，如果不对结构式非限制性SVAR模型741施加限制条件，其模型参数不可估计。那么，对结构式非限制性SVAR模型741需要施加多少限制或约束条件？需要施加的约束条件数恰好为 749即只要施加个约束条件，那么结构式非

31、限制性SVAR模型731的模型参数就可估计。所施加的约束条件既可以是短期同期的，也可以是长期的。1、 短期约束结构式非限制性SVAR模型741式 ， 其中： 在可逆，既逆阵存在时，可化成简化式非限制性VAR模型746， 进一步，在满足特征方程 的根全在单位园外条件下，那么VAR(p)可逆，从而又可将表示成白噪声滑动和形式 ， 其中， 根据Cholesky分解根本思想，短期约束可直接施加在矩阵上，只要使成为主对角线上元素为1的下三角形矩阵，即那么结构式非限制性SVAR模型741式就可变成一递归形式的结构式非限制性SVAR模型，从而为恰好识别，可直接采用OLS从第1方程开场估计该结构式模型的模型参

32、数，建立模型。在实际中，对结构式非限制性SVAR模型741式施加短期约束，也可以不呈下三角形，只要施加约束条件数，即可。例如，如果我们要建立一个以GDP、税收、政府支出为变量的的结构式非限制性SVAR模型，那么只需施加个约束条件：，当期GDP影响当期税收，不影响当期G政府支出；，当期税收影响当期G政府支出；，根据以往研究已得知，税收关于产出弹性为1.71, 那么所建结构式非限制性SVAR模型即可识别,从而可估计。2、 长期约束所谓长期约束, 通常是指施加在7410式上的约束, 也可是单独施加在某一上的约束而言。比拟简单的是一般都施加在上，与短期约束类似，也可将长期约束直接施加在上来进展。 74

33、2 SVAR模型的三种类型SVAR模型根据模型特点主要有三种类型：K型、C型和AB型。其中最常用的是AB型，K型和C型可视为是AB型的特殊形式。1、：K型SVAR模型设为一维随机时间序列，为滞后阶数，为一维随机扰动的时间序列，且其VAR模型结构关系为即 ， 7411或写成滞后算子形式 7. 4. 12其中: ,为滞后算子多项式.设为一个阶可逆阵，左乘7. 4. 12，那么 7. 4. 13如果，那么称满足上述条件的为K型SVAR模型。由于 在下，这意味着对已施加了个非线性约束条件，中还余下个自由参数，因此，只需给出个短期约束条件即可。3、 C型SVAR模型对于VAR模型7411 ， 或写成滞后

34、算子形式 7. 4. 14设为一个阶可逆阵，如果，且，那么称满足上述条件的7. 4. 14模型为C型SVAR模型。由于，在下，这意味着对已施加了个非线性约束条件，中还余下个自由参数。3、AB型SVAR模型设为阶可逆阵，左乘7. 4. 12，那么 7. 4. 15且满足条件：，那么称7. 4. 15为AB型SVAR模型。显然，当为单位阵时，AB型SVAR模型就化为C型SVAR模型；当为单位阵时，AB型SVAR模型就化为K型SVAR模型。由 可知 在下，它是对施加了个非线性约束条件，余下了个自由参数。 SVAR模型的脉冲响应函数设有VAR(p)模型参见722 ， 写成滞后算子形式为 其中: ,为滞

35、后算子多项式. 在满足特征方程 的根全在单位园外条件下，那么VAR(p)是可逆的，即可将表示成白噪声滑动和形式 其中：阶单位阵由于、称为非正交化，故需要进展正交化变换。因是正定对称阵，故可以将其分解成其中，是一下三角阵，是唯一的一个主对角线上元素大于零的对角阵。利用对作变换， ，即 那么 为正交这时有， 于是与7. 2. 7、7. 2. 8同理，可导出一个正交的脉冲响应函数为 7. 4. 16由的脉冲引起的累积响应函数为： 7. 4. 17其中，是的第行、第列元素。用矩阵可表示为 = 7. 4. 18即 中的第行、第列元素等于时期的第变量扰动项增加一个单位，其它时期扰动项为常数时，对时期的第个变量值的影响。对于SVAR模型741式 ， 由于在可逆，既逆阵存在时，可化成简化式非限制性VAR模型， 进一步，在满足特征方程 的根全在单位园外条件下，那么VAR(p)可逆，从而又可将表示成白噪声滑动和形式 因此，其脉冲响应函数用矩阵也可表示为 = 7. 4. 19 如果SCAR模型为AB型，那么其脉冲响应函数为 7. 4. 20累积脉冲响应函数为 7. 4. 2120 / 20