金融計量VARVEC模型講義

1、1 傳統(tǒng)的經(jīng)濟計量方法是以經(jīng)濟理論為基礎來描述變傳統(tǒng)的經(jīng)濟計量方法是以經(jīng)濟理論為基礎來描述變量關系的模型。但是，經(jīng)濟理論通常并不足以對變量之間量關系的模型。但是，經(jīng)濟理論通常并不足以對變量之間的動態(tài)聯(lián)系提供一個嚴密的說明，而且內(nèi)生變量既可以出的動態(tài)聯(lián)系提供一個嚴密的說明，而且內(nèi)生變量既可以出現(xiàn)在方程的左端又可以出現(xiàn)在方程的右端使得估計和推斷現(xiàn)在方程的左端又可以出現(xiàn)在方程的右端使得估計和推斷變得更加復雜。變得更加復雜。為了解決這些問題而出現(xiàn)了一種用非結(jié)構性方法來為了解決這些問題而出現(xiàn)了一種用非結(jié)構性方法來建立各個變量之間關系的模型。本章所要介紹的向量自回建立各個變量之間關系的模型。本章所要介紹的

2、向量自回歸模型歸模型(vector autoregression，VAR)和向量誤差修正模型和向量誤差修正模型(vector error correction model，VEC)就是非結(jié)構化的多就是非結(jié)構化的多方程模型。方程模型。 2 向量自回歸向量自回歸(VAR)是基于數(shù)據(jù)的統(tǒng)計性質(zhì)建立模型，是基于數(shù)據(jù)的統(tǒng)計性質(zhì)建立模型，VAR模型是處理多個相關經(jīng)濟指標的分析與預測模型是處理多個相關經(jīng)濟指標的分析與預測最容易操作的模型之一，并且在一定的條件下，多元最容易操作的模型之一，并且在一定的條件下，多元MA和和ARMA模型也可轉(zhuǎn)化成模型也可轉(zhuǎn)化成VAR模型，因此近年來模型，因此近年來VAR模型受到越

3、來越多的經(jīng)濟工作者的重視。模型受到越來越多的經(jīng)濟工作者的重視。 3 VAR(p) 模型的數(shù)學表達式是模型的數(shù)學表達式是 (9.1.1)其中：其中：yt是是 k 維內(nèi)生變量列向量，維內(nèi)生變量列向量，xt 是是d 維外生變量列向量，維外生變量列向量，p是滯后階數(shù)，是滯后階數(shù)，T是樣本個數(shù)。是樣本個數(shù)。k k 維矩陣維矩陣 1， p 和和 k d 維維矩陣矩陣 H 是待估計的系數(shù)矩陣。是待估計的系數(shù)矩陣。 t 是是 k 維擾動列向量，它們相互維擾動列向量，它們相互之間可以同期相關，但不與自己的滯后值相關且不與等式右邊之間可以同期相關，但不與自己的滯后值相關且不與等式右邊的變量相關，假設的變量相關，假

4、設 是是 t 的協(xié)方差矩陣，是一個的協(xié)方差矩陣，是一個(k k)的正定的正定矩陣。式矩陣。式(9.1.1)可以展開表示為可以展開表示為 Tt,2, 1ttptpttHxyyy 114(9.1.2) 即含有即含有 k 個時間序列變量的個時間序列變量的VAR(p)模型由模型由 k 個方程個方程組成。組成。Tt,2, 1ktttdtttptkptptptkttktttxxxyyyyyyyyy21212111211121H5ttttttMbIPbMaIPacM, 222222112211 , 22111其中其中， ci , aij , bij 是要被估計的參數(shù)。也可表示成：是要被估計的參數(shù)。也可表示成

5、：ttttttttMIPbbbbMIPaaaaccMIP, 2, 12222211211112221121121111作為作為VAR的一個例子，假設工業(yè)產(chǎn)量（的一個例子，假設工業(yè)產(chǎn)量（IP）和貨和貨幣供應量（幣供應量（M1）聯(lián)合地由一個雙變量的聯(lián)合地由一個雙變量的VAR模型決定。內(nèi)生模型決定。內(nèi)生變量滯后二階的變量滯后二階的VAR(2)模型是：模型是： ttttttMbIPbMaIPacIP, 12122111121111116 一般稱式一般稱式(9.1.1)為為(unrestricted VAR)。沖擊向量沖擊向量 t 是白噪聲向量，因為是白噪聲向量，因為 t 沒有結(jié)構性的含沒有結(jié)構性的含義

6、，被稱為簡化形式的沖擊向量。義，被稱為簡化形式的沖擊向量。 為了敘述方便，下面考慮的為了敘述方便，下面考慮的VAR模型都是不含外生變量模型都是不含外生變量的非限制向量自回歸模型，用下式表示的非限制向量自回歸模型，用下式表示 或或 其中其中: :tptpttyyy 11ttLy)(ppkLLLLI221)(9.1.5)7 如果行列式如果行列式det (L)的根都在單位圓外，則式的根都在單位圓外，則式(9.1.5)滿足穩(wěn)定性條件，可以將其表示為無窮階的向量動平均滿足穩(wěn)定性條件，可以將其表示為無窮階的向量動平均(VMA()形式形式 (9.1.6)其中其中 ttL Ay)(1)()(LLA2210)(

7、LLLAAAAkIA 08 對對VAR模型的估計可以通過最小二乘法來進行和極大模型的估計可以通過最小二乘法來進行和極大似然估計法，在特定條件下，兩者的估計系數(shù)是完全相同似然估計法，在特定條件下，兩者的估計系數(shù)是完全相同的。當隨機擾動項滿足獨立同分布（的。當隨機擾動項滿足獨立同分布（iid.N(0,），可以證），可以證明最小二乘法和極大似然法估計的系數(shù)結(jié)果一致。明最小二乘法和極大似然法估計的系數(shù)結(jié)果一致。所以如果所以如果VAR模型的隨機擾動項滿足上述條件，則可以對模型的隨機擾動項滿足上述條件，則可以對VAR模型中的每個方程分別進行模型中的每個方程分別進行OLS估計，獲得的系數(shù)是估計，獲得的系數(shù)是

8、有效的一致估計值。另外即使跨等式的隨機擾動項存在相有效的一致估計值。另外即使跨等式的隨機擾動項存在相關性，但只要各個等式擾動項自身不存在序列相關，那么關性，但只要各個等式擾動項自身不存在序列相關，那么OLS也可以給出有效的估計值。也可以給出有效的估計值。 9注意，由于任何序列相關都可以通過增加更多的注意，由于任何序列相關都可以通過增加更多的yt 的滯的滯后而被消除后而被消除，所以擾動項序列不相關的假設并不要求非所以擾動項序列不相關的假設并不要求非常嚴格。常嚴格。 10 為了研究貨幣供應量和利率的變動對經(jīng)濟波動的長為了研究貨幣供應量和利率的變動對經(jīng)濟波動的長期影響和短期影響及其貢獻度，采用我國期

9、影響和短期影響及其貢獻度，采用我國1995年年1季度季度2007年年4季度的季度數(shù)據(jù)，并對變量進行了季節(jié)調(diào)整。設季度的季度數(shù)據(jù)，并對變量進行了季節(jié)調(diào)整。設居民消費價格指數(shù)為居民消費價格指數(shù)為CPI_90 (1990年年1季度季度=1)、居民消費、居民消費價格指數(shù)增長率為價格指數(shù)增長率為CPI 、實際、實際GDP的對數(shù)的對數(shù)ln(GDP/CPI_90) 為為ln(gdp) 、實際實際M1的對數(shù)的對數(shù)ln(M1/CPI_90) 為為ln(m1) 和實和實際利率際利率rr (一年期存款利率一年期存款利率R-CPI ）。）。 11tttptptptpttttttgdpmrrgdpmrrcccgdpmr

10、r3211111321)ln() 1ln()ln() 1ln()ln() 1ln( 利用利用VAR(p)模型對模型對 ln(gdp) ， ln(m1) 和和 rr，3個變量之個變量之間的關系進行實證研究，其中實際間的關系進行實證研究，其中實際GDP和實際和實際M1以對數(shù)差分以對數(shù)差分的形式出現(xiàn)在模型中，而實際利率沒有取對數(shù)。的形式出現(xiàn)在模型中，而實際利率沒有取對數(shù)。 12 為了創(chuàng)建一個為了創(chuàng)建一個VAR對象，應選擇對象，應選擇Quick/Estimate VAR或者選擇或者選擇Objects/New object/VAR或者在命令窗口中鍵入或者在命令窗口中鍵入var。便會出現(xiàn)下圖的對話框便會出

11、現(xiàn)下圖的對話框(以例以例9.1為例為例)： 13 無約束向量自回歸（無約束向量自回歸（Unrestricted VAR）或者向量或者向量誤差修正（誤差修正（Vector Error Correction）。）。無約束無約束VAR模模型是指型是指VAR模型的簡化式。模型的簡化式。 14 在在Lag Intervals for Endogenous編輯框中輸入滯后信編輯框中輸入滯后信息，表明哪些滯后變量應該被包括在每個等式的右端。息，表明哪些滯后變量應該被包括在每個等式的右端。例例如，滯后對如，滯后對 1 4表示用系統(tǒng)中所有內(nèi)生變量的表示用系統(tǒng)中所有內(nèi)生變量的1階到階到4階滯后變量作為等式階滯后變

12、量作為等式右端的變量。右端的變量。 也可以添加代表滯后區(qū)間的任意數(shù)字，但都要成對輸也可以添加代表滯后區(qū)間的任意數(shù)字，但都要成對輸入。例如：入。例如： 2 4 6 9 12 12即為用即為用24階，階，69階及第階及第12階滯后變量。階滯后變量。 15 EViews允許允許VAR模型中包含外生變量，模型中包含外生變量，其中其中 xt 是是 d 維外生變量向量維外生變量向量 , k d 維矩陣維矩陣 H 是要被估計的系數(shù)是要被估計的系數(shù)矩陣。可以在矩陣?？梢栽贓xogenous Variables編輯欄中輸入相應的外生變編輯欄中輸入相應的外生變量。系統(tǒng)通常會自動給出常數(shù)量。系統(tǒng)通常會自動給出常數(shù)

13、c 作為外生變量。作為外生變量。 其余兩個菜單（其余兩個菜單（Cointegration 和和 Restrictions）僅與僅與VEC模型有關，將在下面介紹。模型有關，將在下面介紹。 ttptpttHxyyy 1116 VAR對象的設定框填寫完畢，單擊對象的設定框填寫完畢，單擊OK按紐，按紐，EViews將會在將會在VAR對象窗口顯示如下估計結(jié)果：對象窗口顯示如下估計結(jié)果： 17 表中的每一列對應表中的每一列對應VAR模型中一個內(nèi)生變量的方模型中一個內(nèi)生變量的方程。對方程右端每一個變量，程。對方程右端每一個變量，EViews會給出會給出、估計、估計及及。例如，在例如，在D(log(M1_SA

14、_P)的方程中的方程中RR_SA(-1)的系數(shù)是的系數(shù)是-0.002187。 同時，有兩類回歸統(tǒng)計量出現(xiàn)在同時，有兩類回歸統(tǒng)計量出現(xiàn)在VAR對象估計輸對象估計輸出的底部：出的底部： 18 輸出的第一部分顯示的是每個方程的標準輸出的第一部分顯示的是每個方程的標準OLS回歸統(tǒng)回歸統(tǒng)計量。根據(jù)各自的殘差分別計算每個方程的結(jié)果，并顯示計量。根據(jù)各自的殘差分別計算每個方程的結(jié)果，并顯示在對應的列中。在對應的列中。 輸出的第二部分顯示的是輸出的第二部分顯示的是VAR模型的回歸統(tǒng)計量。模型的回歸統(tǒng)計量。19 殘差的協(xié)方差的行列式值殘差的協(xié)方差的行列式值(自由度調(diào)整自由度調(diào)整)由下式得出：由下式得出： 其中其

15、中 m 是是VAR模型每一方程中待估參數(shù)的個數(shù)，不做自由模型每一方程中待估參數(shù)的個數(shù)，不做自由度調(diào)整的殘差協(xié)方差行列式計算中不減度調(diào)整的殘差協(xié)方差行列式計算中不減 m。 是是 k 維殘差維殘差列向量。通過假定服從多元正態(tài)（高斯）分布計算對數(shù)似列向量。通過假定服從多元正態(tài)（高斯）分布計算對數(shù)似然值：然值： AIC和和SC兩個信息準則的計算將在后文詳細說明。兩個信息準則的計算將在后文詳細說明。 tttmT 1dett ln22ln12TTnl20 例例9.1結(jié)果如下：結(jié)果如下： 盡管有一些系數(shù)不是很顯著，我們?nèi)匀贿x擇滯后階數(shù)為盡管有一些系數(shù)不是很顯著，我們?nèi)匀贿x擇滯后階數(shù)為2。3個方程擬合優(yōu)度分別

16、為：個方程擬合優(yōu)度分別為： 可以利用這個模型進行預測及下一步的分析。可以利用這個模型進行預測及下一步的分析。 37. 0,17. 0,85. 02212GDPMRRRRtttttttttttteeegdpmrrgdpmrrgdpmrr321222111)ln() 1ln(0.034-0.0150.0040.002-0.124-0.00317.5511.2-0.387-)ln() 1ln(0.495-0.004-0.005-0.404-0.1780.002-4.0-1.51-1.320.0390.040.17)ln() 1ln(21 同時，為了檢驗擾動項之間是否存在同期相關關系，同時，為了檢驗擾

17、動項之間是否存在同期相關關系，可用殘差的同期相關矩陣來描述。用可用殘差的同期相關矩陣來描述。用ei 表示第表示第 i 個方程的殘個方程的殘差，差，i =1，2，3。其結(jié)果如表其結(jié)果如表9.1所示。所示。 e1e 2e 3e 110.36-0.4e 20.3610.15 e 3-0.40.15 122 從表中可以看到實際利率從表中可以看到實際利率rr、實際實際M1的的 ln(m1) 方程和實際方程和實際GDP的的 ln(gdp)方程的殘差項之間存在的方程的殘差項之間存在的同期相關系數(shù)比較高，進一步表明實際利率同期相關系數(shù)比較高，進一步表明實際利率、實際貨實際貨幣供給量幣供給量(M1)和實際和實際

18、GDP之間存在著同期的影響關系，之間存在著同期的影響關系，盡管得到的估計量是一致估計量，但是在本例中卻無盡管得到的估計量是一致估計量，但是在本例中卻無法刻畫它們之間的這種同期影響關系。法刻畫它們之間的這種同期影響關系。 23 在式在式(9.1.1)或式或式(9.1.3)中，可以看出，中，可以看出，VAR模型并模型并沒有給出變量之間當期相關關系的確切形式，即在模沒有給出變量之間當期相關關系的確切形式，即在模型的右端不含有當期的內(nèi)生變量，而這些當期相關關型的右端不含有當期的內(nèi)生變量，而這些當期相關關系隱藏在誤差項的相關結(jié)構之中，是無法解釋的，所系隱藏在誤差項的相關結(jié)構之中，是無法解釋的，所以將式以

19、將式(9.1.1)和式和式(9.1.3)稱為稱為VAR模型的簡化形式模型的簡化形式。本。本節(jié)要介紹的結(jié)構節(jié)要介紹的結(jié)構VAR模型模型(Structural VAR，SVAR)，實際是指實際是指VAR模型的結(jié)構式，即模型的結(jié)構式，即。 24 為了明確變量間的當期關系，首先來研究兩變量的為了明確變量間的當期關系，首先來研究兩變量的VAR模型結(jié)構式和簡化式之間的轉(zhuǎn)化關系。如含有兩個模型結(jié)構式和簡化式之間的轉(zhuǎn)化關系。如含有兩個變量變量(k=2)、滯后一階滯后一階(p=1)的的VAR模型結(jié)構式可以表示模型結(jié)構式可以表示為下式為下式 (9.1.8)Tt,2, 1ztttttxtttttuzxxczuzxz

20、cx1221212120112111121025 在模型在模型(9.1.8)中假設：中假設： （1）隨機誤差）隨機誤差 uxt 和和 uzt 是白噪聲序列，不失一般性，是白噪聲序列，不失一般性，假設方差假設方差 x2 = z2 =1 ； （2）隨機誤差）隨機誤差 uxt 和和 uzt 之間不相關，之間不相關，cov(uxt , uzt )=0 。 式式(9.1.8)一般稱為一般稱為。 26 它是一種結(jié)構式經(jīng)濟模型，引入了變量之間的作它是一種結(jié)構式經(jīng)濟模型，引入了變量之間的作用與反饋作用，其中系數(shù)用與反饋作用，其中系數(shù) c12 表示變量表示變量 zt 的單位變化對的單位變化對變量變量 xt 的的

21、， 21表示表示 xt-1的單位變化對的單位變化對 zt 的的。雖然。雖然 uxt 和和 uzt 是單純出現(xiàn)在是單純出現(xiàn)在 xt 和和 zt 中的隨機沖擊，中的隨機沖擊，但如果但如果 c21 0，則作用在則作用在 xt 上的隨機沖擊上的隨機沖擊 uxt 通過對通過對 xt 的影響，能夠即時傳到變量的影響，能夠即時傳到變量 zt 上，這是一種上，這是一種；同樣，如果；同樣，如果 c12 0，則作用在則作用在 zt 上的隨機沖擊上的隨機沖擊 uzt 也可以對也可以對 xt 產(chǎn)生間接的即時影響。沖擊的交互影響體現(xiàn)產(chǎn)生間接的即時影響。沖擊的交互影響體現(xiàn)了變量作用的雙向和反饋關系。了變量作用的雙向和反饋

22、關系。 ztttttxtttttuzxxczuzxzcx1221212120112111121027 為了導出為了導出VAR模型的簡化式方程，將上述模型表示為模型的簡化式方程，將上述模型表示為矩陣形式矩陣形式 該模型可以簡單地表示為該模型可以簡單地表示為 (9.1.9)Tt,2, 1ztxtttttuuzxzxcc11222112112010211211tttuyyC110028 假設假設 C0可逆，可導出簡化式方程為可逆，可導出簡化式方程為 其中其中 (9.1.10)tttuCyCCy101110010tty11020100100C222112111101Ctttt2110uC29 從而可以

23、看到，簡化式擾動項從而可以看到，簡化式擾動項 t 是結(jié)構式擾動項是結(jié)構式擾動項 ut 的線性組合，因此代表一種復合沖擊。因為的線性組合，因此代表一種復合沖擊。因為 uxt 和和 uzt 是不是不相關的白噪聲序列，則可以斷定上述相關的白噪聲序列，則可以斷定上述 1t 和和 2t 也是白噪聲也是白噪聲序列，并且均值和方差為序列，并且均值和方差為 2211221222112221221111111)var(, 0)(, 0)(cccccctsEEzxttst2211222122112222122222111)var(, 0)(, 0)(cccccctsEExzttst30 同期的同期的 1t 和和

24、2t 之間的協(xié)方差為之間的協(xié)方差為 從式從式(9.1.11)可以看出當可以看出當 c12 0 或或 c21 0 時，時，VAR模模型簡化式中的擾動項不再像結(jié)構式中那樣不相關，正如型簡化式中的擾動項不再像結(jié)構式中那樣不相關，正如例例9.1中的表中的表9.1所顯示的情況。所顯示的情況。(9.1.11)22112122122112212221212111)(),cov(ccccccccEzxtttt31 下面考慮下面考慮k個變量的情形，個變量的情形，p階結(jié)構向量自回歸模型階結(jié)構向量自回歸模型SVAR(p)為為 (9.1.13)其中其中： , , , , piikkikikikiiikiii,2, 1

25、,)()(2)(1)(2)(22)(21)(1)(12)(11kttttuuu21utptptttuyyyyC22110111212211120kkkkccccccC32 可以將式可以將式(9.1.13)寫成滯后算子形式寫成滯后算子形式 (9.1.14)其中：其中：C(L) = C0 1L 2L2 pLp ，C(L)是滯后算是滯后算子子L的的 k k 的參數(shù)矩陣，的參數(shù)矩陣，C0 Ik。需要注意的是，需要注意的是，如果如果 C0 是一個下三角矩陣，則是一個下三角矩陣，則SVAR模型稱為遞歸模型稱為遞歸的的SVAR模型。模型。 kttttELIuuuyC)(,)(33 用時間序列模型來分析影響關

26、系的一種思路，是考慮擾用時間序列模型來分析影響關系的一種思路，是考慮擾動項的影響是如何傳播到各變量的。下面先根據(jù)兩變量的動項的影響是如何傳播到各變量的。下面先根據(jù)兩變量的VAR(2)模型來說明脈沖響應函數(shù)的基本思想。模型來說明脈沖響應函數(shù)的基本思想。 ttttttttttttzdzdxcxczzbzbxaxax222112211122112211(9.4.1)其中，其中，ai，bi，ci，di 是參數(shù)，是參數(shù)， t = ( 1t , 2t ) 是擾動項，假是擾動項，假定是具有下面這樣性質(zhì)的白噪聲向量：定是具有下面這樣性質(zhì)的白噪聲向量： Tt,2, 134(9.4.2) 假定上述系統(tǒng)從期開始活動

27、，且設假定上述系統(tǒng)從期開始活動，且設 x-1=x-2= z-1=z-2= 0，又設于第期給定了擾動項又設于第期給定了擾動項 10 =1， 20 =0，并且其后均為，并且其后均為，即即 1t = 2t =0 (t 1，2，)，稱此為第期給稱此為第期給 x 以脈沖。以脈沖。stttsttttt,0)(,)()var(,0)( 35 下面討論下面討論 xt 與與 zt 的響應，的響應，t = 0 時：時： 將其結(jié)果代入式將其結(jié)果代入式(9.4.1) ，當，當t = 1時時1111,czax 再把此結(jié)果代入式再把此結(jié)果代入式(9.4.1) ，當，當t =2時時 ,112212cbaax112112cd

28、cacz繼續(xù)這樣計算下去，設求得結(jié)果為繼續(xù)這樣計算下去，設求得結(jié)果為稱為稱為。同時所求得。同時所求得 ,43210 xxxxx0, 100zxttttttttttttzdzdxcxczzbzbxaxax22211221112211221136,43210zzzzz。 當然，第期的脈沖反過來，從當然，第期的脈沖反過來，從 10 =0， 20 =1 出發(fā)，出發(fā)，可以求出由可以求出由 z 的脈沖引起的的脈沖引起的 x 的響應函數(shù)和的響應函數(shù)和 z 的響應函數(shù)。的響應函數(shù)。因為以上這樣的脈沖響應函數(shù)明顯地捕捉對沖擊的效果，因為以上這樣的脈沖響應函數(shù)明顯地捕捉對沖擊的效果，所以同用于計量經(jīng)濟模型的沖擊乘

29、數(shù)分析是類似的。所以同用于計量經(jīng)濟模型的沖擊乘數(shù)分析是類似的。 37 本例選擇鋼鐵行業(yè)及其主要的下游行業(yè)的銷售收入本例選擇鋼鐵行業(yè)及其主要的下游行業(yè)的銷售收入數(shù)據(jù)做為各行業(yè)的需求變量，利用脈沖響應函數(shù)分析各下數(shù)據(jù)做為各行業(yè)的需求變量，利用脈沖響應函數(shù)分析各下游行業(yè)自身需求的變動對鋼鐵行業(yè)需求的影響。游行業(yè)自身需求的變動對鋼鐵行業(yè)需求的影響。 分別用分別用 表示表示； 表示表示 表示表示； 表示表示；表示表示。樣本區(qū)間為。樣本區(qū)間為1999年年1月月2002年年12月，所采用月，所采用數(shù)據(jù)均作了季節(jié)調(diào)整，指標名后加上后綴數(shù)據(jù)均作了季節(jié)調(diào)整，指標名后加上后綴sa，并進行了協(xié)并進行了協(xié)整檢驗，存在協(xié)

30、整關系，這表明，所選的各下游行業(yè)的銷整檢驗，存在協(xié)整關系，這表明，所選的各下游行業(yè)的銷售收入與鋼鐵工業(yè)的銷售收入之間具有長期的均衡關系。售收入與鋼鐵工業(yè)的銷售收入之間具有長期的均衡關系。38 為了得到脈沖響應函數(shù)，先建立一個為了得到脈沖響應函數(shù)，先建立一個VAR模型，然后在模型，然后在VAR工具欄中選擇工具欄中選擇View/Impulse Response或者在工具欄或者在工具欄選擇選擇Impulse，并得到下面的對話框，有兩個菜單：并得到下面的對話框，有兩個菜單：Display 和和 Impulse Definition。39 選擇以圖或表來顯示結(jié)果。如果選擇選擇以圖或表來顯示結(jié)果。如果選擇

31、Combined Graphs 則則Response Standard Error選項是灰色，不顯示標準誤差。選項是灰色，不顯示標準誤差。而且應注意：而且應注意： 輸入產(chǎn)生沖擊的變量（輸入產(chǎn)生沖擊的變量（Impulses）和希望觀察其脈沖響和希望觀察其脈沖響應的變量（應的變量（Responses）。）?？梢暂斎雰?nèi)生變量的名稱，也可可以輸入內(nèi)生變量的名稱，也可以輸入變量的對應的序數(shù)。以輸入變量的對應的序數(shù)。40 例如，如果例如，如果VAR模型以模型以GDP、M1、CPI的形式定義，的形式定義，則既可以以：則既可以以： GDP CPI M1 的形式輸入，也可以以的形式輸入，也可以以 1 3 2 的

32、形式輸入。輸入變量的順序僅僅影響結(jié)果的顯示。的形式輸入。輸入變量的順序僅僅影響結(jié)果的顯示。 還應定義一個確定響應函數(shù)軌跡的期間的正整數(shù)。如還應定義一個確定響應函數(shù)軌跡的期間的正整數(shù)。如果想顯示累計的響應，則需要單擊果想顯示累計的響應，則需要單擊Accumulate Response選選項。對于穩(wěn)定的項。對于穩(wěn)定的VAR模型，脈沖響應函數(shù)應趨向于模型，脈沖響應函數(shù)應趨向于0，且累，且累計響應應趨向于某些非計響應應趨向于某些非0常數(shù)。常數(shù)。 41 提供計算脈沖響應標準誤差的選項。解析的或提供計算脈沖響應標準誤差的選項。解析的或Monte Carlo標準誤差對一些標準誤差對一些Impulse選項和誤

33、差修正模型（選項和誤差修正模型（VEC）一般不一定有效。若選擇了一般不一定有效。若選擇了Monte Carlo，還需在下面的編還需在下面的編輯框確定合適的迭代次數(shù)。輯框確定合適的迭代次數(shù)。 如果選擇表的格式，被估計的標準誤差將在響應函數(shù)值如果選擇表的格式，被估計的標準誤差將在響應函數(shù)值下面的括號內(nèi)顯示。如果選擇以多圖來顯示結(jié)果，曲線圖將下面的括號內(nèi)顯示。如果選擇以多圖來顯示結(jié)果，曲線圖將包括關于脈沖相應的正負（包括關于脈沖相應的正負（+/-）兩個標準偏離帶。在）兩個標準偏離帶。在Combined Graphs中將不顯示標準誤差偏離帶。中將不顯示標準誤差偏離帶。 42 設置脈沖為殘差的一個單位的

34、沖擊。這個選項忽略了設置脈沖為殘差的一個單位的沖擊。這個選項忽略了VAR模型殘差的單位度量和相關性，所以不需要轉(zhuǎn)換矩陣的模型殘差的單位度量和相關性，所以不需要轉(zhuǎn)換矩陣的選擇。這個選項所產(chǎn)生的響應函數(shù)是選擇。這個選項所產(chǎn)生的響應函數(shù)是VAR模型相對應模型相對應VMA()模型的系數(shù)。模型的系數(shù)。 設置脈沖為殘差的一個標準偏差的沖擊。這個選項忽略設置脈沖為殘差的一個標準偏差的沖擊。這個選項忽略了了VAR模型殘差的相關性。模型殘差的相關性。 43 用殘差協(xié)方差矩陣的用殘差協(xié)方差矩陣的Cholesky因子的逆來正交化脈沖。因子的逆來正交化脈沖。這個選項為這個選項為VAR模型的變量強加一個次序，并將所有影

35、響變模型的變量強加一個次序，并將所有影響變量的公共因素歸結(jié)到在量的公共因素歸結(jié)到在VAR模型中第一次出現(xiàn)的變量上。注模型中第一次出現(xiàn)的變量上。注意：如果改變變量的次序，將會明顯地改變響應結(jié)果?？梢砸猓喝绻淖冏兞康拇涡?，將會明顯地改變響應結(jié)果?？梢栽谠贑holesky Ordering 的編輯框中重新定義的編輯框中重新定義VAR模型中變量模型中變量的次序。的次序。 44 Cholesky分解有分解有2種選擇：種選擇： ：在估計的殘差協(xié)方在估計的殘差協(xié)方差矩陣利用差矩陣利用Cholesky 因子時進行小樣本的自由度修正。具有因子時進行小樣本的自由度修正。具有自由度修正的殘差協(xié)方差矩陣的第自由度修

36、正的殘差協(xié)方差矩陣的第(i, j)元素的計算是按下列公元素的計算是按下列公式計算的：式計算的： 其中其中m是是VAR模型中每一個方程中待估計參數(shù)的個數(shù)。模型中每一個方程中待估計參數(shù)的個數(shù)。 ：估計殘差協(xié)估計殘差協(xié)方差矩陣的第方差矩陣的第(i, j) 元素的計算是按下列公式計算的：元素的計算是按下列公式計算的：tt jt imT1tt jt iT145 用結(jié)構因子分解矩陣估計的正交轉(zhuǎn)換矩陣。如果沒有先用結(jié)構因子分解矩陣估計的正交轉(zhuǎn)換矩陣。如果沒有先估計一個結(jié)構因子分解矩陣，或者沒有對模型施加約束，這估計一個結(jié)構因子分解矩陣，或者沒有對模型施加約束，這個選項不能用。個選項不能用。 描述描述Pesa

37、ran和和Shin(1998)構建的不依賴于構建的不依賴于VAR模型中模型中變量次序的正交的殘差矩陣。應用按上面的變量次序的正交的殘差矩陣。應用按上面的Cholesky順序計順序計算的第算的第j個變量的個變量的Cholesky因子得到第因子得到第j個變量的擾動項的廣個變量的擾動項的廣義脈沖響應。義脈沖響應。 46 這個選項允許用戶定義脈沖。建立一個包含脈沖的矩這個選項允許用戶定義脈沖。建立一個包含脈沖的矩陣（或向量），并在編輯框中輸入矩陣的名字。如果陣（或向量），并在編輯框中輸入矩陣的名字。如果VAR模型中有模型中有k個內(nèi)生變量，則脈沖矩陣必須是個內(nèi)生變量，則脈沖矩陣必須是k行和行和1列或列或

38、k列列的矩陣，每一列代表一個脈沖向量。的矩陣，每一列代表一個脈沖向量。 例如：一個有例如：一個有k（= 3）個變量的個變量的VAR模型，希望同模型，希望同步對第一個變量有一個正的一個單位的沖擊，給第二個變步對第一個變量有一個正的一個單位的沖擊，給第二個變量一個負的一個單位的沖擊，可以建立一個量一個負的一個單位的沖擊，可以建立一個3 1的脈沖矩的脈沖矩陣陣 S ，其值分別為：其值分別為：1， 1，0。在編輯框中鍵入矩陣的。在編輯框中鍵入矩陣的名字名字: S。 47 例例9.4建立建立5變量的變量的VAR(3)模型，下面分別給各下游模型，下面分別給各下游行業(yè)銷售收入一個沖擊（選擇廣義脈沖）行業(yè)銷售

39、收入一個沖擊（選擇廣義脈沖） ，得到關于鋼，得到關于鋼材銷售收入的脈沖響應函數(shù)圖。在下列各圖中，橫軸表材銷售收入的脈沖響應函數(shù)圖。在下列各圖中，橫軸表示沖擊作用的響應期間數(shù)示沖擊作用的響應期間數(shù)(單位：月度單位：月度)，縱軸表示鋼材，縱軸表示鋼材銷售收入銷售收入(億元億元)，實線表示脈沖響應函數(shù)，代表了鋼材，實線表示脈沖響應函數(shù)，代表了鋼材銷售收入對相應的行業(yè)銷售收入的沖擊的反應，虛線表銷售收入對相應的行業(yè)銷售收入的沖擊的反應，虛線表示正負兩倍標準差偏離帶示正負兩倍標準差偏離帶 。 48y1：鋼材鋼材； y2：建材建材； y3：汽車汽車； y4：機械機械； y5：家電家電49 從第一個圖中可以

40、看出，當在本期給建材行業(yè)銷售收入從第一個圖中可以看出，當在本期給建材行業(yè)銷售收入一個正沖擊后，鋼材銷售收入在前一個正沖擊后，鋼材銷售收入在前4期內(nèi)小幅上下波動之后期內(nèi)小幅上下波動之后在第在第6期達到最高點期達到最高點( =12.03，即在第，即在第6期期 y1 對對 y2 的響應是的響應是12.03) ；從第；從第9期以后開始穩(wěn)定增長。這表明建材行業(yè)受外期以后開始穩(wěn)定增長。這表明建材行業(yè)受外部條件的某一沖擊后，經(jīng)市場傳遞給鋼鐵行業(yè)，給鋼鐵行業(yè)部條件的某一沖擊后，經(jīng)市場傳遞給鋼鐵行業(yè)，給鋼鐵行業(yè)帶來同向的沖擊，而且這一沖擊具有顯著的促進作用和較長帶來同向的沖擊，而且這一沖擊具有顯著的促進作用和較

41、長的持續(xù)效應。的持續(xù)效應。 從第二幅圖中可以看出，當在本期給汽車行業(yè)銷售收入從第二幅圖中可以看出，當在本期給汽車行業(yè)銷售收入一個正沖擊后，鋼材銷售收入在前一個正沖擊后，鋼材銷售收入在前5期內(nèi)會上下波動；從第期內(nèi)會上下波動；從第5期以后開始穩(wěn)定增長期以后開始穩(wěn)定增長( =1.76)。這表明汽車行業(yè)的某一沖。這表明汽車行業(yè)的某一沖擊也會給鋼鐵行業(yè)帶來同向的沖擊，即汽車行業(yè)銷售收入增擊也會給鋼鐵行業(yè)帶來同向的沖擊，即汽車行業(yè)銷售收入增加會在加會在5個月后對鋼材的銷售收入產(chǎn)生穩(wěn)定的拉動作用。個月后對鋼材的銷售收入產(chǎn)生穩(wěn)定的拉動作用。 )6(12a)5(13a50 從第三幅圖中可以看出，機械行業(yè)銷售收入

42、的正沖擊從第三幅圖中可以看出，機械行業(yè)銷售收入的正沖擊經(jīng)市場傳遞也會給鋼材銷售收入帶來正面的影響，并且此經(jīng)市場傳遞也會給鋼材銷售收入帶來正面的影響，并且此影響具有較長的持續(xù)效應。從第四幅圖中可以看出當在本影響具有較長的持續(xù)效應。從第四幅圖中可以看出當在本期給家電行業(yè)銷售收入一個正沖擊后，也會給鋼材銷售收期給家電行業(yè)銷售收入一個正沖擊后，也會給鋼材銷售收入帶來正面的沖擊，但是沖擊幅度不是很大。入帶來正面的沖擊，但是沖擊幅度不是很大。 綜上所述，由于市場化程度、政府保護政策等各方面綜上所述，由于市場化程度、政府保護政策等各方面的原因，使得各下游相關行業(yè)的外部沖擊會通過市場給鋼的原因，使得各下游相關行業(yè)的外部沖擊會通過市場給鋼鐵行業(yè)帶來不同程度的影響，但是都是同向的影響。政府鐵行業(yè)帶來不同程度的影響，但是都是同向的影響。政府可以利用這種現(xiàn)象，對市場進行有區(qū)別、有重點的調(diào)整，可以利用這種現(xiàn)象，對市場進行有區(qū)別、有重點的調(diào)整，減少盲目的重復建設項目。減少盲目的重復建設項目。 回歸分析中考慮一個變量依賴于另一個變量，但這不一定意味著因果關系。 格蘭杰因果檢驗的思想：時間不會倒退，即，如果事件A在事件B之前發(fā)生，那么可能是A導致了B，但不可能是B導致了A，換言之，過去的事件