根據 Engle and Granger(1987),所謂共整合(cointegration)的定義,是指一 組非定態的時間序列變數,透過線性組合而得到定態的時間序列,我們稱這些變 數之間具有共整合關係;一序列要進行共整合檢定的前提是原始序列必須為非定 態,且需以同階做整合(即以同階差分使成定態)檢定;而本文中共整合檢定的目
的在於確保(3.12)式所估出來的消費偏向θˆ 具有經濟涵義。
至於如何檢定共整合的時間序列,茲分別介紹如下:
(1) (1) (1)
(1) Engle and Granger(1987)所提出的兩階段估計法(two-stage estimation),主 要是將兩個變數以普通最小平方法(OLS)做迴歸,並利用所獲得的殘差項來檢 定兩變數的共整合關係,若殘差項呈現 I
( )
0 序列,則表示兩變數有共整合關係 存在。兩階段共整合檢定法用來檢定非同條件程序間是否存在著長期的趨勢均衡 關係。若序列X 須經過 d 次差分過後成為同條件程序,則稱其整合級次為 d ,t 表示為Xt~I
( )
d 。若符合兩個條件,則X 與t Y 共整合關係存在:(i)序列t X 與t Yt 的整合級次相同,皆為 d 。(ii)存在一X 與t Y 的線性組合t Zt =Yt −AXt,其整合 級次為 b ,且b< 。則此時序列d X 與t Y 的共整合級次為t(
d,b)
,表示為CI(
d,b)
,其中(1,−A)為共整合向量(cointegration vector)。此方法是以最小平方法估計共 整合向量。當 n 個變數組成的向量都具有同階的整合級次時,最多可存在n−1條 共整合向量。但在兩階段共整合檢定中只能估計出至多一條共整合向量存在,且 為n−1條共整合向量中的線性組合,因此最小平方法得到的共整合向量只能代表 變數間確實存在共整合關係,並不代表只存在一條共整合關係,亦即二階段共整 合檢定法只適用於尋找兩變數之間的共整合向量,對於多個變數則不適用。
(2) (2) (2)
(2) Johansen(1988)則是將 Engle and Granger 觀念做延伸,提出軌跡檢定(trace test)以及最大特性根檢定(max-eigen test) 兩種概似比檢定統計量,來檢定共 整合向量的個數,以判斷變數間有無存在共整合關係。分述如下:
(i) (i) (i)
(i) 軌軌軌軌跡跡跡跡((((tracetracetracetrace))))統計量統計量統計量統計量::::
( ) ∑ ( )
+
=
−
−
=
n
r i
i
trace r T
1
1 ln λ
λ ) ; r= rank
( )
Π 其中,T 代表樣本總數,而λ)i是第 i 個特性根的估計值。用以檢定 Ho:rank
( )
Π ≤r(至多有 r 個共整合向量)H1:rank
( )
Π >r(最少有r+1個共整合向量)這個檢定的原理是,如果沒有任何共整合向量存在(即r=0),隠含
2 0
1 =λ =⋅⋅⋅=λn =
λ ,這將使得λtrace =0;而如果有 r 個共整合向量存 在,λ1 ≠0,λ2 ≠0,⋅⋅⋅,λr ≠0但是λr+1 =λr+2 =⋅⋅⋅=λn =0,則檢定統計量 之值會接近零。
(ii) (ii) (ii) (ii) 最大特性根最大特性根最大特性根最大特性根((((max eigenmax eigenmax eigenmax eigen))))統計量統計量統計量統計量::::λmax
(
r,r+1)
=−Tln(
1−λ)r+1)
用以檢定 Ho:rank
( )
Π =r則 H1:rank( )
Π =r+1同理,如果有 r 個共整合向量存在,則 λ1 ≠0,λ2 ≠0,⋅⋅⋅,λr ≠0 但是
2 0
1 = + =⋅⋅⋅= =
+ r n
r λ λ
λ ,因為第r+1個特性根λr+1 =0,故檢定統計量之 值也應該接近零。
上述兩種統計量的漸進分配(asymptotic distribution)皆為布朗移動(Brownian motion)的函數。其臨界值可參閱 Johansen–Juselius(1990)以及 Osterwald–Lenum
(1992)所建立之模型分配表。
這兩種檢定方法事先假設變數間不存在共整合關係。若拒絕該假設則依次增 加向量個數再行檢定,直到完全無法拒絕假設為止;檢定結果若存在有一個或多 個顯著的特性根,則表示相關變數之間具有長期均衡的關係。
(3) (3) (3)
(3) Philip and Loretan(1990)則是提出在估計長期均衡關係時,單一方程式的誤 差修正模型並不是一個漸進的有效估計法,在所有單一方程式模型 single-equation method)估計法中,唯有完全修正最小平方統計量(Fully Modified OLS estimator)
才能符合完全效率性(fully efficient)與漸進中位數不偏性(asymptotically median unbiased)。
FM-OLS 估計法在時間序列資料為非恆定時估計較為準確。因為完全修正 估計法修正了共整合關係中模型誤差項與自變數誤差項的的長期回饋效果
(feedback effect)。而且完全修正最小平方法可直接利用非恆定資料做迴歸分 析,並自動修正解釋變數與模型誤差項的序列相關偏誤及內生性偏誤問題。在本 文所用到的變數裡,變數之間的確存在相互影響的問題,例如:當暫時性所得增 加時∆Yt =1,消費也會增加0<∆Ct <1,故經常帳也會增加0<∆CAt <1,也因 此我們必須注意序列相關偏誤及內生性偏誤的問題。
本文採用了 Philips and Hansen(1990)提出的完全修正最小平方估計法(Fully Modified OLS,FM-OLS)來估計(3.12)式單一方程式的共整合關係,同時也採用 Hansen(1992)提出的參數穩定性檢定法來檢定期共整合關係是否存在。
第三節 向量自我迴歸模型
傳統上計量模型的建立,皆事先判定變數是外生或是內生,然後再建立模 型,但是實際上變數之間會相互影響,若不分辨因果關係,則模型參數估計之後 所做的推論,其正確性將大受影響。
Sims(1980)認為根據先驗理論建立的模型,估計得到的實證結果並無法證 明是來自所有經濟變數的聯合過程(joint process)之特性,亦即無法確定是否表 現真正資料的特性,因此便提出向量自我迴歸模型(Vector Autoregressive Model,
VAR),其主要特點在於 Sims 假設各個變數皆為內生的,並建立 VAR 型來觀察 變數之間互動關係,來反應模型參數內真實的互動效果,避免結構設定上的錯誤。
但是 Sims 的方法限定所有的變數皆為穩定的或是在小樣本之下才可以使 用,然而我們所使用的資料透過檢定不一定能符合穩定的性質,因此若貿然的將 變數直接以傳統的 VAR 來估計,可能會誤導所估計的結果,因此本文中我們改 用 Phillips(1995)所發展的完全修正向量自我迴歸模型(FM-VAR),我們以矩 陣表達 VAR(p)如下:
t t p t p t
t
t y y y x
y = Λ1 −1 +Λ2 −2 +L+Λ − +ψ +ε (4.12) 式(4.12)中y 為 k 個內生變數所組成的向量;t Λ ,1
21
Λ ,……..,Λ 為p
(
k×k)
的係數矩陣;x 為外生變數所組成的向量,也可包含時間趨勢項,季節或是結構虛擬t 變數;ψ為x 所對應的係數矩陣;t εt為誤差項,誤差項之間可同期相關,但是不 可跨期相關,並且與解釋變數獨立。至於落後期數 p 的選取標準,則以 SBC
(Schwarz’s Bayesian Criterion)作為判斷的準則。
我們將(4.12)式以誤差修正模型(Vector Error Correction Model ,VECM)來 表示如下:
t t t p t p t
t
t y y y y x
y =Φ1∆ −1+Φ2∆ −2 +⋅⋅⋅+Φ −1∆ − +1+ς −1+ψ +ε (4.13)
其中,
第四節 因果關係檢定(Granger Causality Test)
早期經濟學家對因果關係的定義爭論不已,但一般而言,至少同意因果需要 有三個條件:第一是變數間的相關;第二為現象之間(即時間先後順序)暫時不 對稱;第三則是消除另一個現象會產生觀察到的共變數之機率。
直至 Granger(1969)從「預測能力(Predictability)」的角度來定義變數的 因果關係,亦即欲預測Xt時,除考慮變數X 過去的所有資料以外,若將變數Y的
式(4.15)中每一列共有 k 個變數,具有 n 期;γi為常數項,βij為係數參數,其中 k
j
i, =1,2,⋅⋅⋅, , εi為誤差項其必須符合白噪音的假設。我們可由式中發現到yit不 但可以由自己的過去值來解釋,也可由其他變數的過去值來解釋,假如第j個變 數之落後期可以解釋現在第 i 個變數,那麼參數βij則會明顯的異於零。因此,倘 若y 可以解釋jt yit,那麼則稱yjt Granger-Causes yit。其假設檢定以及檢定統計
量如下:
虛無假設 H0:有p個βij為0
對立假設 H1:上面p個βij至少有一個不為零
當H0成立時,表示變數yjt不會對變數yit造成 Granger 因果關係;假如拒絕掉虛
無假設H0,則表示變數yjt Granger-Causes 變數yit。對此所使用的 F 檢定統計量 如下:
) 1 ( /
/ ) (
+
−
= −
k k n SSE
p SSE F SSE
U
U
R ~F
(
p,n−k(k+1))
(4.16)(4.16)式中SSE 表示受限制模型(restricted model)的誤差平方和,R SSE 表示未受u 限制模型(unrestricted model)之迴歸的殘差平方和。
若單一βij的檢定,H0:βij =0,因為 FM-VAR 之估計參數巳達效率性
(efficiency),因此, Phillips and Hansen(1990)以及 Phillips(1995)建議可 以一般的 t 檢定來做參數係數值的檢定;是以,本文即以簡單個別估計參數之 t 值
來做為Granger因果關係的檢定,其檢定結果與(4.16)式之 F 檢定結果是一致的。
第五節 Wald 檢定
Wald 檢定(Wald test)主要的目的在於檢視是否將模型適當的限制其參數,
由於本文主要是檢驗根據模型所估計出來的經常帳的預測值與實際的預測值是
⇒
[ ] [ ]
Φ第六節 預測模型的檢定
依循 Diebold and Mariano(1995)中所介紹的六種主要檢定方法,來比較不 同模型間的預測能力,以期希望能找出預測誤差相對比較小的一組模型,來做為
= 為預測誤差(forecast error),而損失函數(loss function)為預測誤 差的函數;
( ) ( ) ∑ ( )
k 為向前預測值之期數(k-step-ahead forecast)。
若檢定出的統計量S 太大而拒絕了虛無假設,就表示1
{ }
y)it Tt=1這一組所對應的預測 誤差e 相對比較大,因此i{ }
y)it Tt=1的預測能力比較差;反之,若檢定出的統計量太 小落入了左尾的拒絕域,就表示{ }
yjt Tt=1
) 這一組所對應的預測誤差e 相對比較大,j
因此
{ }
yjt Tt=1
) 的預測能力比較差;但若檢定統計量無法拒絕虛無假設,則表示
{ }
y)it Tt=1與
{ }
y)jt Tt=1之預測能力不分軒輊。故藉由符號檢定可以檢定出哪一個模型的 預測值比較好。
(2) (2) (2)
(2) 符號檢定符號檢定符號檢定符號檢定(((The Sign Test(The Sign TestThe Sign Test)The Sign Test)))
用以檢定兩組預測模型是否一樣好(即med
[
g( )
eit −g( )
ejt]
=0)的虛無假設 是否成立。若檢定出的統計量S2a其值太大而拒絕了虛無假設,就表示{ }
y)it Tt=1這 一組所對應的預測誤差e 相對比較大,因此i{ }
y)it Tt=1的預測能力比較差;反之若檢 定出的統計量其值太小落入了左尾的拒絕域,就表示
{ }
yjt Tt=1
) 這一組所對應的預測
誤差e 相對比較大,因此j
{ }
yjt Tt=1
) 的預測能力比較差;但若檢定統計量無法拒絕虛
無假設,則表示
{ }
y)it Tt=1與{ }
yjt Tt 1=
) 之預測能力不分軒輊。故藉由符號檢定可以檢定
出哪一個模型的預測值比較好。其檢定統計量介紹如下:
T
s
Ts
2a= 20−.250.5 ~N( )
0,1 (4.23)其中, T
( )
tt
d I
S
∑
=
= + 1
2 且 dt ≡
[
g( )
eit −g( )
ejt]
( )
>+ =
otherwise d d if
I t t
0
0 1
((((3333)))) 魏克森符號排序檢定魏克森符號排序檢定(魏克森符號排序檢定魏克森符號排序檢定(((Wilcoxon's SignedWilcoxon's SignedWilcoxon's SignedWilcoxon's Signed----Rank TestRank TestRank Test)Rank Test)))
預測能力比較好。
=0
若檢定出的統計量 MR 顯著大於 0,拒絕了虛無假設,就表示
{ }
y)it Tt=1這一組 所對應的預測誤差e 相對比較大,因此i{ }
y)it Tt=1的預測能力比較差,{ }
yjt Tt1
=
) 的預測
能力表現較佳;反之若檢定出的統計量太小落入左尾的拒絕域,就表示
{ }
yjt Tt=1
) 這
一組所對應的預測誤差e 相對比較大,因此j
{ }
yjt Tt=1
) 的預測能力比較差;但若檢定
統計量無法拒絕虛無假設,則表示
{ }
y)it Tt=1與
{ }
yjt Tt=1
) 之預測能力不分軒輊。
第伍章 實證分析與結果
第一節 實證資料的選取與處理
資料來源 資料來源 資料來源 資料來源::: :
本研究所採用的資料來源取自於國際貨幣基金(International Monetary Fund,
IMF)之國際金融統計資料庫(International Financial Statistics),IFS CD-ROM,
Version 1.1.53。所擷取的時間序列(time series data)資料詳列如下:
Version 1.1.53。所擷取的時間序列(time series data)資料詳列如下: