• Sonuç bulunamadı

şeklinde verilmiş olsun. Buradan, ρ nun EKK tahmin edicisi,

N/A
N/A
Protected

Academic year: 2021

Share "şeklinde verilmiş olsun. Buradan, ρ nun EKK tahmin edicisi,"

Copied!
16
0
0

Yükleniyor.... (view fulltext now)

Tam metin

(1)

5.5. Asimptotik Özellikler

Bu kısımda, birim köklü zaman serilerinde parametrelerin EKK tahmin edicilerinin bazı asimptotik özellikleri incelenecektir. Bu özelliklerin çoğu Hamilton (1994) den özetlenerek aktarılmıştır.

Birinci dereceden otoregresif zaman serisi modeli,

X t = ρ X t−1 +e t , t=1,2,3, ...,n

şeklinde verilmiş olsun. Buradan, ρ nun EKK tahmin edicisi,

ρ ^ n = [ t=1 n X t−1 2 ] −1 t=1 n X t X t−1 = ρ + [ t=1 n X t−1 2 ] −1 [ t=1 n e t X t−1 ] şeklinde olup, ρ∈[ 0 , 1]

ve n→∞ iken ρ ^ nP ρ ve | ρ |<1 için,

n ( ρ ^ nρ ) D N (0 , 1−ρ 2 )

olduğunu biliyoruz. ρ=1 için asimptotik dağılım normal değildir. Bilindiği gibi, ρ=1 için n ( ρ ^ n −1 ) = [ n 1 2 t=1 n X t−1 2 ] −1 [ 1 n t=1 n X t−1 e t ] =O P ( 1)

dir. Bu terimin payındaki istatistiğin asimptotik dağılımını bulmak kolaydır. X 0 =0 olmak üzere,

X n =( e 1 +e 2 +.. .+e n ) olacağından,

t=1 n

X t−1 e t = ∑

t=1

n [ t−1 j=1 e j ] e t = 1 2 [ X n 2 t=1 n e t 2 ]

yazılır. Ayrıca, e t ~ WN (0, σ 2 ) olduğundan zayıf büyük sayılar yasasına göre n→∞ iken,

1

n

t =1 n

e t 2P σ 2

ve merkezi limit teoremine göre de n→∞ iken, (0, 2 )

D

n e n  N

olduğunu biliyoruz. Buradan da n→∞ iken

  2 2 2 2 2 1

1

1 n n D

n t

t

n e e X

n n  

 

    

  

asimptotik dağılımı elde edilir. Bu iki sonuç Slutsky teoremi ile beraber, ρ nun EKK tahmin edicisinin payının asimptotik dağılımı n→∞ iken,

1

n

t=1 n

X t−1 e t = 1

2 [ X n n 2 1 n t=1 n e t 2 ] D σ 2 2 ( χ 1 2 −1 )

olarak bulunur. n ( ρ ^ n −1 ) nin paydasındaki istatistiğin asimptotik dağılımı bu kadar kolay değildir.

1 ( 1 2 ... 1 )

t t

X e    e e olduğundan,

E( X t−1 )=0 , Var( X t−1 )= E ( X t−1 2 )=σ 2 ( t−1)

(2)

dir. O halde, paydadaki istatistiğin beklenen değeri

E [ t=1 n X t−1 2 ] 2 t=1 n ( t−1)= σ 2 n ( n−1) 2

olduğundan,

n −2

t=1 n

X t−1 2 =O P (1)

dir. Asimptotik dağılımı elde etmek için Brownian Hareketini (Brownian Motion) hatırlayalım.

Brownian Hareketi: e t ~ (0,1) N normal dağılımlı beyaz gürültü serisi olsun. Rasgele yürüyüş süreci (random walk), X 0 =0 olmak üzere X t = X t−1 +e t şeklindedir. Buna göre,

X t ~ N ( 0 , t ) olup s >t için

X sX t =e t+1 + e t+2 + .. .+e s ~ N ( 0 , s−t )

dir. Ayrıca, 0<t <s <r<q için X sX t ile X qX r bağımsızdır.

Tanım 5.5.1 Standart Brownian Hareketi W (t ) , aşağıdaki üç özelliği her bir t∈[ 0 , 1 ] için

sağlayan sürekli zamanlı stokastik süreçtir. Bunlar;

a) W (0) =0

b) 0     t 1 ... t k 1 için W t ( ) 2W t ( ) , 1 W t ( ) 3W t ( ) ,..., ( ) 2 W t kW t ( k 1 ) rasgele değişkenleri bağımsız normal dağılıma sahiptir.

c) W (t ) hemen hemen her yerde t nin sürekli bir fonksiyonudur

Örnek 5.5.1 Aşağıda bazı asimptotik özellikler verilmiştir (Levin ve Lin (1992)). Bu örnekte aksi söylenmedikçe, W (r ) standart Brownian hareketidir. Aşağıdaki özelliklerden bazılarının çözümüne girilmemiştir.

1) W (0)=0 dır.

2) E [ W (r ) ] =0 , ∀ r , E [ W (1) ] =0 dır.

3) Var [ W (r) ] =r , ∀ r

4) W (r ) ~ N (0,r ) dır.

5) Cov W s W t( ), ( )   E W s W t( ) ( )   min{ , } s t dir.

İspat: s≤t için iddia

   

 

( ) , ( ) ( ), ( ) ( ) ( ) ( ), ( ) ( ) ( ( )

( ( )) ( uncorrelated increments) Cov W s W t Cov W s W t W s W s

Cov W s W t W s Var W s Var W s s

  

  

 

eşitliğinden açıktır.

6) E [ W 2 (r ) ] =r

7) E W r

4

( ) 3 r

2

. Not: (6) ve (7) den Var W r ( 2 ( ) ) 3  r 2r 2  2 r 2 bulunur.

8) E W s W t ( 2 ( ) 2 ( ))   s t 2 min{ , } s t 2 2 dir.

İspat: Yine, s≤t için W t 2 ( ) nin

(3)

W 2 ( t ) =W (t ) W (t )= [ W (t )−W (s)+W ( s) ] [ W (t )−W ( s)+W (s) ]

şeklinde yazılması ile,

 

 

 

2 2

2 2 2

2 2 3 4

2 2 2 2

( ) ( ) [ ( ) ( ) ( )][ ( ) ( ) ( )] ( ) ( ) ( ) ( )] 2 ( )[ ( ) ( )] ( ) ( ) [ ( ) ( )] ( ) 2 ( )[ ( ) ( )] ( )

( ) 0 3 2 min{ , }

[

E W s W t E W t W s W s W t W s W s W s W s E W t W s W s W t W s W s W s E W t W s W s W s W t W s W s

t s s s t s s t s t s

    

    

    

       

 

 

bulunur.

9) Cov W s W t ( 2 ( ), 2 ( )) 2 min{ , }  t s 2 2 dir.

İspat: Yine, s≤t için kovaryansın tanımından

   

 

2 2 2 2 2 2

2 2 2 2

2 2 2 2

( ) , ( ) [ ( ) ( ( ))][ ( ) ( ( ))]

[ ( ) ( )] [ ( )] [ ( )]

2 min{ , } 2 min{ , }

Cov W s W t E W t E W t W s E W s

E W t W s E W s E W t

s t t s s t t s

  

 

   

bulunur

Brownian hareketi, merkezi limit teoreminin genel bir kullanılmına olanak sağlar. Merkezi limit teoremi basit hali ile söyledir: Z 1 , Z 2 , ..., Z n ,... beklenen değeri μ varyansı σ 2 olan

bağımsız aynı dağılımlı rasgele değişkenler ve Z ¯ n de örneklem ortalaması olmak üzere, n→∞

iken

n ( ¯ Z nμ ) D N ( 0 , σ 2 )

dir.

Beklenen değeri sıfır varyansı σ 2 olan aynı dağılımlı rasgele değişkenlerin dizisi , 1, 2,3,...

e t t  olsun. Örneklemin ilk yarısını (  n / 2  kadarını) kullanmak isteyelim. Burada, n çift ise  n / 2  *n / 2 , n tek olduğunda ise  n / 2  *  ( n  1) / 2 dir. Bu örnekleme merkezi limit teoremi uygulandığında,

/ 2 *

/ 2 *

1

1 / 2

n

n t

t

e e

n

  

olmak üzere n→∞ iken,

* 2

/ 2 n / 2 D ( 0, ) n e   N

 

yazılır. Aynı şekilde örneklemin ikinci yarısına merkezi limit teoremi uygulandığında aynı asimptotik özellik elde edilir. Bu iki dağılım birbirinden bağımsızdır. Bunun daha genel hali, r∈[ 0 , 1 ] olmak

üzere X n (r) ,

/ 2 *

1

( ) 1

n

n t

t

X r e

n

(4)

şeklinde verilebilir. Burada, X n ( r) , r ye göre bir basamak fonksiyonudur. Daha açık bir ifade ile X n ( r) ,

1

1 2

1 2

0 , 0 (1/ )

/ , (1/ ) (2 / )

( ) / , (2 / ) (3/ )

( ) . . .

. . .

( ... ) / , 1

n

n

r n

e n n r n

e e n n r n

X r

e e e n r

  

  

    

  

 

    



olarak yazılabilir. Buradan,

* * *

|

1 * 1

1 1

( )

n r n r

n t t

t t

n X r e n r e

n n n r

 

 

 

 

 

 

     

  dir. Diğer taraftan,

*

2

* 1

1 ( 0, )

n r

D t t

e N

n r

 

 

  ve

n r *

n r

 

  

 

 

 

olduğundan, n→∞ iken √ n X n ( r) ⃗ D N ( 0 , r σ 2 ) veya

n ( X n σ ( r ) ) D N ( 0 , r )

elde edilir. Benzer şekilde, r 2 > r 1 ve n r 1  * ve n r 2  * tam değer kısımlarını göstermek üzere n→∞ iken,

n ( X n ( r 2 )− σ X n ( r 1 ) ) D N ( 0 , r 2 −r 1 )

dir. O halde, W (.) standart normal rasgele değişkeni göstermek üzere,

n ( X n σ ( . ) ) D W (. )

dir. Bu gösterim, fonksiyonel merkezi limit teoremi (functional central limit theorem) olarak bilinir.

Burada, X n ( .) rasgele bir fonksiyonu göstermekte olup, X n (r) bu rasgele fonksiyonun r zamanındaki değeridir. Dolayısı ile, X n ( .) bir fonksiyon olup X n (r) de bir rasgele değişkendir. X n ( 1) de bilinen örneklem ortalamasıdır. Yani,

X n ( 1) =¯e n =n −1

t=1 n

e t

dir. Fonksiyonel merkezi limit teoreminden, bilinen merkezi limit teoremine geçiş r=1 ile yapılır.

Yani,

1

(1) / 1 n n D (1) ~ ( 0 , 1)

n t

t

n X e n e W N

n

   

(5)

dir. Bu sonuçlar ve sürekli dönüşüm teoremi (continuous mapping theorem) yardımı ile rasgele değişken dizilerinin bir çok özelliği elde edilir. Önce süreki dönüşüm teoremini hatırlayalım:

Sürekli Dönüşüm Theoremi (Continuous Mapping Theorem, Hamilton (1994) s.482) X n rasgele

değişkenlerin bir dizisi, g :    sürekli bir fonksiyon ve X de bir rasgele değişken olsun.

Bu durumda n→∞ iken,

X nD X ise g( X n ) ⃗ D g( X )

dir 

Buna göre n→∞ iken, n X n (.)  DW (.) , W r ( ) ~ ( 0 , ) N r ve

( ) ( )

n

n X r  DW r olup, ( ) ( ) 2

n n

S r    n X r   ise n X n ( .)⃗ D σ W (. )

dir. Buradan da

S n ( r ) ⃗ D σ 2 W 2 ( r )

elde edilmiş olur. Şimdi, sürekli dönüşüm teoremi kullanılarak durağan olmayan zaman serilerinde bu sonuçların uygulamalarını görelim.

Birim Köklü Zaman Serileri İçin Asimptotik Özellikler: Birim köklü zaman serilerinde parametrelerin EKK tahmin edicilerinin asimptotik dağılımlarının normal olmadığını biliyoruz. Önce, AR(1) modeli X 0 =0 olmak üzere, X t   X t 1e t olarak verilmiş olsun. ρ=1 için,

X t = e 1 + e 2 +.. .+e t olup, X n (r) rasgele değişkenini,

1

1 2

1 2

0 , 0 (1/ )

/ , (1/ ) (2 / )

( ) / , (2 / ) (3/ )

( ) . . .

. . .

( ... ) / , 1

n

n

r n

e n n r n

e e n n r n

X r

e e e n r

  

  

    

  

 

    



şeklinde yazalım. Burada X n ( r) , r nin bir fonksiyonudur. Buradaki n tane dikdörtgenin alanları toplamı tüm eğri altındaki alanı verir. Örneğin, t. inci dikdörtgenin alanı, 1/n uzunluğu

ile X t−1 / n yüksekliğinin çarpımıdır. O zaman, X n (r) fonksiyonunun integrali,

0 1

X n ( r ) d r =n −2

t=1 n

X t−1

şeklinde yazılabilir. Bu aşağıdaki grafikte gösterilmeye çalışılmıştır.

(6)

Burada,

∫ 0 1

X n (r ) d r =n −2

t=1 n

X t−1

ifadesinin her iki tarafı √ n ile çarpılırsa,

∫ 0 1

n X n (r) d r =n −3 /2

t=1 n

X t−1

eşitliği elde edilir. Buradan, sürekli dönüşüm teoremine göre n→∞ iken

1 1

0 0

( ) D ( )

n X r d r n    W r dr

 

veya

1 3 / 2

0

1 1

( )

n D

t t

n X W r dr

 

 

asimptotik sonucu elde edilmiş olur. Daha açık olarak,

 

 

3/ 2 3/ 2

1 1 1 2 1 2 3 1 2 1

1

3/ 2 1 2 3 2 1

( ) ( ) ... ( ... )

( 1) ( 2) ( 3) ... 2

n

t n

t

n n

n X n e e e e e e e e e

n n e n e n e e e

 

 

  

          

        

3/ 2 1/ 2 3/ 2

1 1 1

( )

n n n

t t t

t t t

n n t e n e n t e

  

      

yazılabilir. O halde,

2 1 1/ 2

1/ 2 1/ 3

V   

  

 

olmak üzere n→∞ iken,

 

1/ 2 3/ 2

1 1

, ,

n n

D

t t

t t

n e n t e N V

 

  

  

    0

olup,

3 / 2 1 1 n

t t

n

X

 terimi beklenen değeri sıfır varyansı,

 

2 1 2 1/ 2 1/ 3 2 / 3

        

olan bir rasgele değişken olup asimptotik dağılımı normaldir. Buna göre, n→∞ iken

3/ 2 2

1 1

(0, / 3)

n D

t t

n X N

 

(7)

elde edilir. Burada N ( 0 , σ 2 / 3) dağılımı

σ

0 1

W (r ) dr

ile ifade edilmiştir. Drift (yığılma)

olmayan X t = X t−1 +e t şeklindeki bir model için n→∞ iken

n −1

t=1 n

X t

ıraksak,

n −3 /2

t=1 n

X t

yakınsaktır. Ayrıca, n→∞ iken

n −3 /2

t=1 n

t e t = n −1/2

t=1 n

e t −n −3 /2

t=1 n

X t−1D σ W (1) −σ

0 1

W (r ) dr

dir.

S n ( r)=n [ X n ( r) ] 2 denirse S , n sonlu sayıdaki dikdörtgenin alanları toplamı olup,

1 n 2

t=1 n

X t−1 2 = ∫

0 1

S n (r ) d r

dir. Sürekli dönüşüm teoremine göre n→∞ iken,

1 n 2

t=1 n

X t−1 2D σ 2

0 1

W 2 ( r ) d r

asimptotik yaklaşımı yazılır. Bu integral, S nin n

1 2 2 2

2

0 , 0 (1/ )

/ , (1/ ) (2 / ) / , (2 / ) (3/ )

( ) . . .

. . .

/ , 1

n

n

r n

X n n r n

X n n r n

S r

X n r

  

  

   

  

 

  

şeklindeki ifadesinden açıktır. Bilindiği gibi n→∞ iken

1 n

t=1 n

X t−1 e tD σ 2

2 [ W 2 ( 1 ) − 1 ]

dir. Bu iki sonuç birleştirilirse X t   X t 1e t modeline göre H 0 : ρ=1 hipotezini test etmek için kullanılan n ( ^ρ n −1 ) istatistiğinin asimptotik dağılımı (Dickey-Fuller dağılımı) n→∞

iken,

n ( ρ ^ n −1 ) =

1 n

t=1 n

X t−1 e t

1 n 2

t=1 n

X t−1 2

D

σ 2

2 [ W 2 ( 1)−1 ]

σ 2

0 1

W 2 ( r ) d r

= [ W 2 ( 1 )−1 ]

2 ∫

0 1

W 2 ( r ) d r

şeklinde ifade edilir. Diğer taraftan n→∞ iken,

n −5/2

t=1 n

t X t−1 =n −3/2

t=1

n ( n t ) X t−1 D σ 0 1 r W (r) d r

ve

(8)

n −3

t=1 n

t X t−1 2 =n −2

t=1

n ( n t ) X t−1 2 D σ 2 0 1 r W 2 ( r) d r

dir. Bununla birlikte,

1 n

t=1 n

X t−1 e t = 1

2 n [ X n 2 t=1 n e t 2 ] = S n 2 ( 1 ) 2 n 1 t=1 n e t 2

olup n→∞ iken,

S n (1 )⃗ D σ 2 W 2 ( 1 ) ve 1

n

t =1 n

e t 2P σ 2

olduğundan, Slutsky teoremine göre n→∞ iken,

1 n

t=1 n

X t−1 e tD σ 2

2 [ W 2 ( 1)−1 ]

asimptotik dağılım elde edilir. Burada, W 2 (1) ~ 1 2 dir. Özet olarak, e t bağımsız aynı dağılımlı rasgele değişkenler olmak üzere, AR(1) modeli, X

t

X

t1

e

t

olarak verildiğinde, aşağıdaki sonuçlar yazılabilir (Hamilton, 1994). Bunlar n→∞ iken,

a)

1

n

t=1 n

e tD σ W (1)

b)

1 n

t=1 n

X t−1 e tD σ 2

2 [ W 2 ( 1)−1 ]

c)

3/ 2 1

1 0

(1) ( )

n D

t t

n t eWW r dr

 

 

d)

1 n 3/2

t=1 n

X t−1D σ

0 1

W (r) d r

e)

1 n 2

t=1 n

X t−1 2D σ 2

0 1

W 2 ( r ) d r

f)

n −5/2

t=1 n

t X t−1D σ

0 1

r W ( r) d r

g)

1 n 3

t=1 n

t X t−1 2D σ 2

0 1

r W 2 ( r ) d r

h)

1 n υ+1

t=1 n

t υ → 1

υ+1 , υ=0 , 1 , 2 , 3, ...

dır.

Şimdi bu sonuçları, birim kök testleri için neredeyse standart hale gelen Dickey-Fuller test

istatistiklerine uygulayalım. Üç farklı durum ayrı ayrı incelenecektir. Aşağıdaki her durumda beyaz

gürültü serisi, birbirinden bağımsız aynı dağılımlı rasgele değişkenler olarak alınacaktır.

(9)

Durum I: AR(1) modeli X t = ρ X t−1 +e t olarak verildiğinde, H 0 : ρ=1 hipotezi altında ρ nun EKK tahmin edicisi

1 2

1 1

1 1

ˆ n n t n t t

t t

X X X

 

 

  

 

   

ve   2 2 1 1 1

1 1

1 1

ˆ n 1 n t n t t

t t

n X X e

n n

 

 

   

 

   

şeklinde olup n→∞ iken,

2 2

1 1

1 (1) 1

2

n D

t t

t

X e W

n

        

 ve

2 2 1 2

2 0

1 1

1 n t D ( )

t X W r d r

n

  

 

olduğundan,

 

1 1

2 2

0

ˆ n 1 D 2 ( ) (1) 1

nW r d r W

  

 

           

olduğunu biliyoruz. AR(1) modeli regresyon modeline benzediğinden, H 0 : ρ=1 hipotezini test etmek için t− türü istatistik daha yaygın şekilde kullanılmaktadır. Bu istatistiğin asimptotik dağılımı da bilinen t− dağılımı değildir. Dağılım farklı olduğundan, istatistik t yerine τ ^ n

kullanılmıştır. τ ^ n istatistiği de basit aritmetik işlemlerden sonra

2 1

1 1

1

1 2

ˆ 2 2 1

1 1

1

ˆ 1 ˆ 1

ˆ

n n t t t

n n t t

n n

n n

n t t t

t

X e X

S S X X S

 

  

 

  

 

 

 

 

 

1 1

1 1

2 2

1 2 1

1 1

1

1

n n

t t t t

t t

n n

n t n t

t t

X e X e

n

S X S X

n

 

 

 

 

   

 

olarak yazılabilir. Burada,

S n 2 = 1 n−1

t=1 n

( X t − ^ ρ n X t−1 ) 2

ve

1

2 2 2

ˆ 1

1 n

t n

t

S X S

 

 

      dir. Ayrıca n→∞ iken,

S n 2 = 1

n−1

t =1 n

( X t − ^ ρ n X t −1 ) 2P σ 2

olup τ ^ n nin asimptotik dağılımı n→∞ iken,

τ ^ nD 1

2 [ 0 1 W 2 ( r ) d r ] −1/2 [ W 2 (1) −1 ]

dir.

Durum II: AR(1) modeli ( X tμ )=ρ ( X t−1μ )+e t olarak verilmiş olsun. Model

α=μ(1− ρ) için X t = α+ ρ X t−1 + e t şeklinde de yazılabilir. H 0 : ρ=1 hipotezi altında

(10)

X t = X t−1 + e t şekline dönüşür. X t = α+ ρ X t−1 +e t modeline göre ρ nun EKK tahmin edicisi Durum I den farklıdır. Parametrelerin EKK tahmin edicileri, 1

1 1 1

1 2 1 1

1 1 1

ˆ ˆ

n n

t t

n t t

n n n

n t t t t

t t t

n X X

X X X X

  

  

  

   

   

      

     

     

   

   

 

  

denklem sisteminin çözümüleridir. H 0 : ρ=1 hipotezi altında, X t yerine X t 1e t yazıldığında denklem sistemi, 1

1 1 1

1 2 1 1

1 1 1

ˆ ˆ 1

n n

t t

n t t

n n n

n t t t t

t t t

n X e

X X X e

  

  

  

   

   

      

      

     

   

   

 

  

şekline dönüşür. Buradaki ifadelerin yakınsama hızları,

t=1 n

X t−1 =O P ( n 3/2 )

, ∑

t=1 n

X t−1 e t =O P ( n)

,

t=1 n

X t−1 2 =O P ( n 2 )

, ∑

t=1 n

e t =O P ( n 1/2 )

şeklinde olup n yeterince büyük ise bu ifadeleri,

3/2 1 1/2

3/2 2

ˆ ( ) ( ) ( )

ˆ 1 ( ) ( ) ( )

n P P P

n P P P

O n O n O n

O n O n O n

    

 

 

  

      

     

şeklinde yazalım. Son eşitliğin her iki tarafı A n matrisi ile çarpıldığında,

1

1/ 2 1 1/ 2

1

1 1

2

1 1

1 1

1/ 2

1 1

1 1

ˆ 0 0

0

ˆ 1

0 0 0

* 0 0

n

n t t

n n

n

t t

t t

n t t

n

t t

t

n X

n n

n

n n X X n

n e

n X e

  

 

 

 

 

 

 

   

       

             

           

       

 

   

 

 

 

 

 

 

 

 

 

 

 

 

veya

 

1

3/ 2 1

1 1

3/ 2 2 2

1 1 1

1 1 1

1 1 ˆ

ˆ 1 1

n n

t t

t t

n

n n n

n t t t t

t t t

n X e

n n

n n X n X X e

n

 

 

 

  

  

   

   

 

   

  

   

  

         

 

  

eşitliği elde edilir. Burada

(11)

0

n 0 A n

n

 

  

  ve

1 0 A 0

 

  

 

dir. Buradan n→∞ iken,

3/ 2 1

1 1 0

1 1

3/ 2 2 2 2 2

1 1

1 1 0 0

1 ( )

1

( ) ( )

D n

t t

n n

t t

t t

W r d r

n X

n X n X W r d r W r d r

 

 

 

 

 

 

 

 

 

 

    

 

 

 

     

 

   

veya,

1 1

0 0

1 1 1 1

2 2 2

0 0 0 0

1 ( ) 1 ( )

( ) ( ) ( ) ( )

A A

W r d r W r d r

W r d r W r d r W r d r W r d r

 

   

   

   

   

   

   

   

 

   

elde edilir. Buna göre n→∞ iken,

2 2 2

1

1 1

(1) (1)

1 (1) 1

(1) 1 2 2

1

1 0 1 0

D n

t t n

t t

t

W W

W W n e

n X e

  

 

 

 

   

       

          

               

 

 

asimptotik dağılımı elde edilir. Bu sonuçlardan parametrelerin EKK tahmin edicilerinin asimptotik dağılımı n→∞ iken,

 

1 1

2 1 1

0

1 1

2

0 0

ˆ (1)

1 (1) 1

ˆ 1

2

1 ( )

( ) ( )

n D

n

n W

A A A

n W

W r dr

W r dr W r dr

 

 

 

 

   

       

       

 

       

 

 

 

2 1 1

0

1 1

2

0 0

0 (1)

0 1 1

(1) 1 2

1 ( )

( ) ( )

W

W W r dr

W r dr W r dr

 

   

   

 

   

 

   

 

   

 

             

 

olarak bulunur. Diğer taraftan,

Δ=

0 1

W 2 ( r) dr− [ 0 1 W (r ) dr ] 2

denirse yukarıdaki matrisin tersi,

1

1 1 1 1

2

0 0 0

1 1 1

2

0 0 0

1 ( ) ( ) ( )

( ) ( ) ( ) 1

W r d r W r d r W r d r

W r d r W r d r W r d r

 

   

    

   

   

   

   

   

  

  

(12)

olur. H 0 : ρ=1 hipotezini test etmek için n ( ^ρ n −1 ) test istatistiğinin kullanılabileceğini biliyoruz. Bunun asimptotik dağılımı ise n→∞ iken,

2 1

0 2

1 1

2

0 0

1 (1) 1 (1) ( )

ˆ 2 ( 1)

( ) ( )

n D

W W W r dr

n

W r dr W r dr

   

 

 

 

 

    

 

2 1

0

1 (1) 1 (1) ( )

2   W     W W r dr

 

şeklindedir. Denklem sisteminin çözümünden hem α ^ n hem de ρ ^ n EKK tahmin edicilerinin asimptotik dağılımları normal değildir. Ayrıca, ortak dağılım da asimptotik normal değildir.

Model ( X tμ )=ρ ( X t−1μ )+e t şeklinde ise ρ nun EKK tahmin edicisi ρ ^ μ ile

gösterilir (örneklem hacmine bağlı olarak bazen ρ ^ n , μ kullanılır).

Benzer şekilde, H 0 : ρ=1 hipotezini test etmek için, t− türü istatistiğin asimptotik dağılımını elde etmeye çalışalım. Önce,

S n 2 = 1 n−2

t=1 n

( X t − ^ α n −^ ρ μ X t−1 ) 2

olmak üzere,

 

1 1 1

2 2

1 2 1

1 1

1 0

( ˆ ) 0 1

1

n t t

n n n n

t t

t t

X

S S

X X

 

 

 

 

 

   

 

      

 

 

 

den t− türü istatistik de

τ ^ n , μ = [ S ^ ρ n

2 ] −1 ( ^ ρ μ −1 )

şeklinde yazılabilir. Buradan,

 

1 1 1

2 2 2

1 2 1

1 1

1 0

( ˆ ) 0

n t t

n n n n

t t

t t

X

n S S n

X X n

 

 

 

 

 

 

 

 

        

 

 

 

1 1 1

1 2 1

1 1

1 0

1

n t t

n n n

t t

t t

X

A

X X

 

 

 

 

   

 

        

 

 

 

olup,

(13)

1 1

1 1 1

2

1 1

1 1

1 1 1

1 2 1

1 1

1 1 n t

t

n n n n n n

t t

t t

n t t

n n

t t

t t

X

A A A A

X X

X

X X

 

 

 

 

 

 

 

     

       

     

     

         

 

 

 

3/ 2 1

1

3/ 2 2 2

1 1

1 1

1 n t

t

n n

t t

t t

n X

n X n X

 

 

 

 

 

 

 

 

 

 

 

eşitliği yukarıdaki sonuçlar ile kullanıldığında n→∞ iken S n 2P σ 2 olduğu görülür.

Buradan da,

 

1 1

2 2 0

1 1

2

0 0

1 ( )

ˆ 0

( ) 0 1

( ) ( ) 1

n P

W r d r n S

W r d r W r d r

 

  

 

   

        

 

 

 

elde edilir. Ayrıca,

  2

1 1

2

0 0

1 1

0

1 1

2

0 0

0 1

0 1 1

( ) ( )

1 ( )

( ) ( ) W r d r W r d r

W r d r

W r d r W r d r

 

 

   

       

   

 

     

 

olduğundan τ ^ n , μ istatistiğinin asimptotik dağılımı da

 

1 2

, , 0

, , 2 2 , 1 2 1 2

0 0

1 (1) 1 (1) ( )

ˆ

ˆ 1 1 2

ˆ ( ˆ ) ( ˆ )

( ) ( )

n n D

n

n n n n

W W W r d r

n

S n S

W r d r W r d r

 

 

 

  

  

  

 

 

 

 

 

 

 

şeklinde elde edilmiş olur.

Durum III (Random Walk with Drift) Burada H 0 : ρ=1 hipotezi altında model,

X t = α + X t−1 +e t şeklinde yazılır. Daha açık olarak X t için zaman serisi modeli,

X t = X 0 + α t +( e 1 + e 2 + .. .+e t )= X 0 + α t +ξ t

olarak yazılabilir. Buna göre,

t=1 n

X t−1 = ∑

t=1 n

( X 0 + α (t−1) +ξ t−1 )

olup n→∞ iken,

1 n 2

t=1 n

α (t−1) = α n 2

n(n−1) 2 = α

2 n−1

nα 2

ve

(14)

n −3 /2

t=1 n

ξ t−1D σ

0 1

W (r ) dr

yakınsamaları elde edilir. Buradan n örneklem hacmi yeterince büyük ise

2 3 / 2 2

1 0 1

1 1 1 1

( 1) ( ) ( ) ( ) ( )

n n n n

t t P P P P

t t t t

X Xt O n O n O n O n

   

       

   

dir. Modeldeki α (t−1) terimi bütün asimptotik sonuçları etkilemektedir. Önce, n→∞ iken

1 1/ 2 3/ 2

2 1 0 2 1

1 1 1

1 1

( 1)

2

n n n P

t t

t t t

X n X t n n

n n

  

  

 

             

  

olduğu açıktır. Örneklem hacmi n yeterince büyük ise,

  2

2 2 2 2 2

0 0

0 0 1

3 2 2 3 / 2 5 / 2 3

1 1 1

1

1 1 1 1 1

1 1 1

( 1) ( 1)

2 ( 1) 2 2 ( 1)

( ) ( ) ( ) ( ) ( ) ( ) ( )

n t

P P P P P P P

n n n n n

t t t

n n

t

t t t t t

t t t

X X t X t

X t X t

O n O n O n O n O n O n O n

   

   

  

    

  

       

    

      

    

  

olup n→∞ iken

1 n 3

t =1 n

X t −1 2P α 2

3 dür. Diğer taraftan,

 

1 0 1 1 0 1

1 1 1 1 1

1/ 2 3/ 2 3 / 2

( 1) ( 1)

( ) ( ) ( ) ( )

n n n n n

t t t t t t t t

t t t t t

P P P P

X e X t e X e t e e

O n O n O n O n

   

   

    

       

   

    

olduğundan n yeterince büyük olduğunda,

3 / 2 3/ 2

1 1 1

( 1)

n n

t t t

t t

n X e n t e

 

 

 

yazılabilir. Başka bir ifade ile,

1 1/ 2

2

2 3 3/ 2

ˆ ( ) ( ) ( )

ˆ 1 ( ) ( ) ( )

n P P P

n P P P

O n O n O n

O n O n O n

 

  

 

  

    

      

      ve

1/ 2 3/ 2

0

n 0 A n

n

 

  

 

 

özellikleri dikkate alınarak elde edilen sonuçlar toparlandığında,

 

 

1

2 1/ 2

1 1 1

2 3 2 3/ 2

1 1 1

1 1 1

ˆ 1 ˆ 1

n n

t t

n t t

n n n

n t t t t

t t t

n X n e

n

n n X n X n X e

 

 

  

  

  

  

   

   

       

 

   

  

     

   

 

  

eşitliği yazılır. Ayrıca, n→∞ iken

2 1

1

2 3 2 2

1 1

1 1

1 1 / 2

/ 2 / 3 :

n

t t P

n n

t t

t t

n X

Q

n X n X

 

 

 

 

 

 

   

     

   

 

 

 

olduğu da dikkate alınırsa,

Referanslar

Benzer Belgeler

2011, çalışmalarında, Çin’in Beijing kentinde hava kirleticiler ile Alerjik Rinit şikâyetiyle hastanelere başvuranlar arasındaki yüksek derecede doğrusal

[65–67] zaman b¨ olgesini diskrize edebilmek i¸cin tamamen a¸cık ¸sema kullanarak farklı viskozite de˘ gerleri i¸cin kuadratik, k¨ ubik B-spline ve en k¨ u¸c¨ uk kareler

Beklenen değer ve otokovaryans fonksiyonu zamana bağlı olmadığından bu model de durağandır.. Otokorelasyonların grafiklerine bakıldığında, fonksiyon değerleri

Bu hata terimlerinin bazıları (modele bağlı olarak) kestirildikten sonra öngörüler AR modellerinde olduğu gibi ardışık olarak hesaplanabilir.. için X n+s

Parametrelerin tahmin edicilerinin özelliklerine göre, değişik birim kök testleri olmasına rağmen, bunlar arasında EKK tahmin edicisinin dağılımına bağlı

Durağan zaman serilerinde örneklem otokorelasyon fonksiyonunun serinin otokorelasyon fonksiyonu için tutarlılığı aşağıda incelenmeye çalışılmıştır... Böylece,

Bunun için eldeki veriler kullanılarak, tahmin edicinin asimptotik dağılımı bilindiğinde  1 için güven aralığı yazılabilir..  1 in EKK tahmin edicisinin

Sonlu Dağıtılmış Gecikme Modelleri (Finite Distributed Lag Models) - FDL Modelleri Yukarıda bahsedilen modeller, daha sonra göreceğimiz zaman serileri verisiyşe regresyon