Zaman Serileri İçin Durağanlık Analizi – 3

Z

 

4.  Sabit ve Trend İçeren Rassal Yürüyüş Modeli Çerçevesinde Durağanlık Analizi

Aşağıda gösterilen sabit terimi ve trendi bulunan model, \displaystyle {Y_t}’nin açıklanan, \displaystyle {Y_t}’nin bir dönem gecikmeli değerinin ve kronolojik dönemlerin açıklayıcı değişken olduğu rassal yürüyüş modelini temsil etmektedir.

\displaystyle {Y_t} = \beta  + \alpha T + \rho {Y_{t - 1}} + {\varepsilon _t}

T ifadesi gözlem dönemlerinin kronolojik sırasını ifade eder. Söz gelimi 2011’den 2016’ya kadar 6 gözlemden oluşan bir Y serisi olsun. Yukarıdaki modeli EKK ile tahmin etmek için kronolojik olmak koşuluyla T açıklayıcı değişkeni olarak 2011’den 2016’ya kadar olan yıllar kullanılabileceği gibi 1’den 6’ya kadar olan gözlem sırası da kullanılabilir. T açıklayıcı değişkeni için kronolojik dönemleri ya da gözlem sırasını kullanmak T’nin katsayısı olan \displaystyle \alpha  ’nın tahmin edilen değerini değiştirmeyecektir.  Dolayısıyla yukarıdaki modelde açıklayıcı değişken olarak yer alan T’nin aynı zamanda gözlem sırasını ifade ettiğini söylemek yanlış olmayacaktır.

n=3 basitleştirici varsayımı altında yukarıdaki modele göre \displaystyle {Y_t}’nin her alternatif t gözlemi için izlediği patika aşağıda gibi olacaktır.

t=1:

\displaystyle {Y_1} = \beta  + \alpha T + \rho {Y_0} + {\varepsilon _1}

t=2:

\displaystyle {Y_2} = \beta  + \alpha T + \rho {Y_1} + {\varepsilon _2}

\displaystyle {Y_2} = \beta  + \alpha T + \rho \underbrace {\left( {\beta  + \alpha T + \rho {Y_0} + {\varepsilon _1}} \right)}_{{Y_1}} + {\varepsilon _2}

\displaystyle {Y_2} = \beta  + \rho \beta  + \alpha T + \rho \alpha T + {\rho ^2}{Y_0} + \rho {\varepsilon _1} + {\varepsilon _2}

t=3:

\displaystyle {Y_3} = \beta  + \alpha T + \rho {Y_2} + {\varepsilon _3}

\displaystyle {Y_3} = \beta  + \alpha T + \rho \underbrace {\left( {\beta  + \rho \beta  + \alpha T + \rho \alpha T + {\rho ^2}{Y_0} + \rho {\varepsilon _1} + {\varepsilon _2}} \right)}_{{Y_2}} + {\varepsilon _3}

\displaystyle {Y_3} = \beta  + \rho \beta  + {\rho ^2}\beta  + \alpha T + \rho \alpha T + {\rho ^2}\alpha T + {\rho ^3}{Y_0} + {\rho ^2}{\varepsilon _1} + \rho {\varepsilon _2} + {\varepsilon _3}

Eğer \displaystyle \rho  = 1 ise serinin birim kökü olduğu anlamına gelir. Böyle bir durumda da \displaystyle {Y_t} aşağıdaki gibi tahmin edilemez yapıya sahip stokastik (olasılıksal) geçmiş dönem hata terimlerinin birikimli bir fonksiyonu olacaktır.

\displaystyle {Y_1} = \beta  + \alpha T + \underbrace \rho _1{Y_0} + {\varepsilon _1}

\displaystyle {Y_1} = \beta  + \alpha T + {Y_0} + {\varepsilon _1}

\displaystyle {Y_2} = \beta  + \underbrace \rho _1\beta  + \alpha T + \underbrace \rho _1\alpha T + \underbrace {{\rho ^2}}_1{Y_0} + \underbrace \rho _1{\varepsilon _1} + {\varepsilon _2}

\displaystyle {Y_2} = 2\beta  + 2\alpha T + {Y_0} + {\varepsilon _1} + {\varepsilon _2}

\displaystyle {Y_3} = \beta  + \underbrace \rho _1\beta  + {\underbrace \rho _1^2}\beta  + \alpha T + \underbrace \rho _1\alpha T + {\underbrace \rho _1^2}\alpha T + {\underbrace \rho _1^3}{Y_0} + {\underbrace \rho _1^2}{\varepsilon _1} + \underbrace \rho _1{\varepsilon _2} + {\varepsilon _3}

\displaystyle {Y_3} = 3\beta  + 3\alpha T + {Y_0} + {\varepsilon _1} + {\varepsilon _2} + {\varepsilon _3}

\displaystyle {Y_t}’nin \displaystyle \beta  gibi bir sabit, trend değişkeni T ve \displaystyle {Y_0} başlangıç değeri ile birlikte geçmiş dönem hata terimlerinin birikimli fonksiyonu olduğu yukarıdaki durum “sabit ve trend içeren rassal yürüyüş modeli” (random walk with drift around a trend model) olarak adlandırılır.

4.1. Ortalama için Durağanlık Analizi

Sabit ve trend içeren rassal yürüyüş modeli çerçevesinde \displaystyle {Y_t}’nin beklenen değeri (ortalaması) için durağanlık analizi aşağıdaki model yardımıyla yapılabilir.

\displaystyle {Y_t} = \beta  + \alpha T + \rho {Y_{t - 1}} + {\varepsilon _t}

Yukarıdaki ifadeye beklenen değer operatörü ilave edilir.

\displaystyle E\left[ {{Y_t}} \right] = E\left[ \beta  \right] + E\left[ {\alpha T} \right] + E\left[ {\rho {Y_{t - 1}}} \right] + \underbrace {E\left[ {{\varepsilon _t}} \right]}_0

Bilindiği gibi hata terimlerinin beklenen değeri 0’dır.

\displaystyle E\left[ {{Y_t}} \right] = E\left[ \beta  \right] + E\left[ {\alpha t} \right] + E\left[ {\rho {Y_{t - 1}}} \right]

Bilindiği gibi \displaystyle \beta  ve \displaystyle \rho  birer sabittir. Ayrıca yukarıdaki modelde bulunan T açıklayıcı değişkenin değeri modelden bağımsız belirlendiği için T bağımsız değişkeni de \displaystyle {Y_t} için bir sabittir.

\displaystyle E\left[ {{Y_t}} \right] = \beta  + \alpha T + \rho E\left[ {{Y_{t - 1}}} \right]

Yukarıdaki sonuçtan yola çıkılarak gözlem sayısı n=3 basitleştirici varsayımı altında her alternatif t gözlemi için \displaystyle {Y_t}’nin beklenen değeri aşağıda analiz edilmektedir.

t=1:

\displaystyle E\left[ {{Y_1}} \right] = \beta  + \alpha T + \rho E\left[ {{Y_0}} \right]

\displaystyle {Y_0} bir sabit olduğundan beklenen değer operatöründen çıkartılabilir.

\displaystyle E\left[ {{Y_1}} \right] = \beta  + \alpha T + \rho {Y_0}

t=2:

\displaystyle E\left[ {{Y_2}} \right] = \beta  + \alpha T + \rho E\left[ {{Y_1}} \right]

\displaystyle E\left[ {{Y_1}} \right] = \beta  + \alpha T + \rho {Y_0} olduğundan denklemde yerine yazılır.

\displaystyle E\left[ {{Y_2}} \right] = \beta  + \alpha T + \rho \left( {\beta  + \alpha T + \rho {Y_0}} \right)

\displaystyle E\left[ {{Y_2}} \right] = \beta  + \rho \beta  + \alpha T + \rho \alpha T + {\rho ^2}{Y_0}

t=3:

\displaystyle E\left[ {{Y_3}} \right] = \beta  + \alpha T + \rho E\left[ {{Y_2}} \right]

\displaystyle E\left[ {{Y_2}} \right] = \beta  + \rho \beta  + \alpha T + \rho \alpha Tt + {\rho ^2}{Y_0} olduğundan denklemde yerine yazılır.

\displaystyle E\left[ {{Y_3}} \right] = \beta  + \alpha T + \rho \left( {\beta  + \rho \beta  + \alpha T + \rho \alpha T + {\rho ^2}{Y_0}} \right)

\displaystyle E\left[ {{Y_3}} \right] = \beta  + \rho \beta  + {\rho ^2}\beta  + \alpha T + \rho \alpha T + {\rho ^2}\alpha T + {\rho ^3}{Y_0}

Yukarıdaki sonuçlar her alternatif t gözlemi için aşağıda özetlenmiştir.

\displaystyle t = 1 \to E\left[ {{Y_1}} \right] = \beta  + \alpha T + \rho {Y_0}

\displaystyle t = 2 \to E\left[ {{Y_2}} \right] = \beta  + \rho \beta  + \alpha T + \rho \alpha T + {\rho ^2}{Y_0}

\displaystyle t = 3 \to E\left[ {{Y_3}} \right] = \beta  + \rho \beta  + {\rho ^2}\beta  + \alpha T + \rho \alpha T + {\rho ^2}\alpha T + {\rho ^3}{Y_0}

Yukarıda \displaystyle {Y_t}’nin beklenen değerleri için gösterilen süreç matematiksel olarak aşağıdaki gibi daha sade ifade edilebilir.

\displaystyle E\left[ {\sum\limits_{t = 1}^{n = 3} {{Y_t}} } \right] = {\rho ^t}{Y_0} + \sum\limits_{t = 1}^{n = 3} {{\rho ^{t - 1}}} \beta  + \sum\limits_{t = 1}^{n = 3} {{\rho ^{t - 1}}\alpha T}

Gözlem sayısı n iken yukarıdaki ifade aşağıdaki gibi genel bir forma dönüştürülebilir.

\displaystyle E\left[ {\sum\limits_{t = 1}^n {{Y_t}} } \right] = {\rho ^t}{Y_0} + \sum\limits_{t = 1}^n {{\rho ^{t - 1}}} \beta  + \sum\limits_{t = 1}^n {{\rho ^{t - 1}}\alpha T}

\displaystyle E\left[ {{Y_t}} \right] = {\rho ^t}{Y_0} + \sum\limits_{t = 1}^n {{\rho ^{t - 1}}} \beta  + \sum\limits_{t = 1}^n {{\rho ^{t - 1}}\alpha T}

Daha önce belirtildiği gibi gözlem sayısı n sonsuza giderken (\displaystyle n \to \infty  ) yukarıdaki denklemin sağındaki ifadeler sonucu aşağıda verilen birer geometrik seri oluştururlar.

\displaystyle E\left[ {{Y_t}} \right] = {\rho ^t}{Y_0} + \underbrace {\sum\limits_{t = 1}^n {{\rho ^{t - 1}}} \beta }_{\frac{\beta }{{1 - \rho }}} + \underbrace {\sum\limits_{t = 1}^n {{\rho ^{t - 1}}\alpha T} }_{\frac{{\alpha T}}{{1 - \rho }}}

\displaystyle E\left[ {{Y_t}} \right] = {\rho ^t}{Y_0} + \frac{\beta }{{1 - \rho }} + \frac{{\alpha T}}{{1 - \rho }}

\displaystyle E\left[ {{Y_t}} \right] = {\rho ^t}{Y_0} + \frac{{\beta  + \alpha T}}{{1 - \rho }}

\displaystyle {Y_t}’nin birim kökü varken, diğer bir ifadeyle \displaystyle \rho  = 1 iken sabit ve trend içeren rassal yürüyüş modeli için \displaystyle {Y_t}’nin beklenen değeri asimptotik açıdan (yani \displaystyle n \to \infty  ) tanımsız olmaktadır.

\displaystyle E\left[ {{Y_t}} \right] = {\rho ^t}{Y_0} + \frac{{\beta  + \alpha T}}{{1 - \rho }}

\displaystyle E\left[ {{Y_t}} \right] = \frac{{{\rho ^t}{Y_0}\left( {1 - \rho } \right) + \beta  + \alpha T}}{{1 - \rho }}

\displaystyle E\left[ {{Y_t}} \right] = \frac{{{1^t}{Y_0}\left( {1 - 1} \right) + \beta  + \alpha T}}{{1 - 1}}

\displaystyle E\left[ {{Y_t}} \right] = \frac{{\beta  + \alpha T}}{0} =  Tanımsız

Diğer taraftan n sonlu ve \displaystyle {Y_t}’nin birim kökü var (\displaystyle \rho  = 1) iken  sabit ve trend içeren rassal yürüyüş modeli için \displaystyle {Y_t}’nin beklenen değeri, gözlem sırasını temsil eden t ile \displaystyle \beta  ve \displaystyle \alpha T’nin çarpımı ve başlangıç değerini temsil eden \displaystyle {Y_0}’ın toplamına eşit olacaktır. Diğer bir ifadeyle \displaystyle \rho  = 1 iken \displaystyle {Y_t}’nin ortalaması zamana bağımlı bir şekilde artış göstermektedir.

\displaystyle E\left[ {{Y_t}} \right] = {\rho ^t}{Y_0} + \sum\limits_{t = 1}^n {{\rho ^{t - 1}}} \beta  + \sum\limits_{t = 1}^n {{\rho ^{t - 1}}\alpha T}

\displaystyle E\left[ {{Y_t}} \right] = {1^t}{Y_0} + \sum\limits_{t = 1}^n {{1^{t - 1}}} \beta  + \sum\limits_{t = 1}^n {{1^{t - 1}}\alpha T}

\displaystyle E\left[ {{Y_t}} \right] = {Y_0} + \underbrace {\sum\limits_{t = 1}^n \beta  }_{t\beta } + \underbrace {\sum\limits_{t = 1}^n {\alpha T} }_{t\alpha T}

\displaystyle E\left[ {{Y_t}} \right] = {Y_0} + t\beta  + t\alpha T

\displaystyle E\left[ {{Y_t}} \right] = {Y_0} + t\left( {\beta  + \alpha T} \right)

\displaystyle {Y_t}’nin birim kökü yok (\displaystyle \rho  \ne 1) ve \displaystyle \rho  = 0 iken sabit ve trend içeren rassal yürüyüş modeli için \displaystyle {Y_t}’nin beklenen değeri zamandan bağımsız şekilde \displaystyle \beta  + \alpha T ifadesine eşit olur.

\displaystyle E\left[ {{Y_t}} \right] = {\rho ^t}{Y_0} + \frac{{\beta  + \alpha T}}{{1 - \rho }}

\displaystyle E\left[ {{Y_t}} \right] = {0^t}{Y_0} + \frac{{\beta  + \alpha T}}{{1 - 0}}

\displaystyle E\left[ {{Y_t}} \right] = \beta  + \alpha T

\displaystyle {Y_t}’nin birim kökü yok (\displaystyle \rho  \ne 1) ve \displaystyle \rho  \ne 0, yani \displaystyle \rho  0 ile 1 arasında bir değere sahip iken, sabit ve trend içeren rassal yürüyüş modeli için \displaystyle {Y_t}’nin beklenen değeri zamandan bağımsız şekilde \displaystyle \beta  + \alpha T ifadesinin \displaystyle 1 - \rho  ’ya bölümüne eşit olacaktır.

\displaystyle 0 < \rho  < 1

\displaystyle \mathop {\lim }\limits_{t \to \infty } {\text{ }}f\left( t \right) = {\rho ^t} = 0

\displaystyle E\left[ {{Y_t}} \right] = {\rho ^t}{Y_0} + \frac{{\beta  + \alpha T}}{{1 - \rho }}

\displaystyle E\left[ {{Y_t}} \right] = \underbrace {{\rho ^t}}_0{Y_0} + \frac{{\beta  + \alpha T}}{{1 - \rho }}

\displaystyle E\left[ {{Y_t}} \right] = \frac{{\beta  + \alpha T}}{{1 - \rho }}

4.2. Varyans için Durağanlık Analizi

Varyans için durağanlık analizine geçmeden önce varyansın durağanlık analizinde yararlanılacak olan \displaystyle {Y_t}’nin \displaystyle {Y_1}, \displaystyle {Y_2} ve \displaystyle {Y_3}’ün beklenen değerleri, beklenen değerlerinin kareleri ve karelerinin beklenen değerleri aşağıda sunulmuştur.

\displaystyle E\left[ {{Y_1}} \right] = {A_1} = \beta  + \alpha T + \rho {Y_0}

\displaystyle {\left( {E\left[ {{Y_1}} \right]} \right)^2} = A_1^2 = {\left( {\beta  + \alpha T + \rho {Y_0}} \right)^2}

\displaystyle E\left[ {Y_1^2} \right] = E\left[ {{{\left( {{A_1} + {\varepsilon _1}} \right)}^2}} \right]

\displaystyle E\left[ {Y_1^2} \right] = E\left[ {A_1^2 + 2{A_1}{\varepsilon _1} + \varepsilon _1^2} \right]

\displaystyle E\left[ {Y_1^2} \right] = A_1^2 + 2{A_1}\underbrace {E\left[ {{\varepsilon _1}} \right]}_0 + E\left[ {\varepsilon _1^2} \right]

\displaystyle E\left[ {Y_1^2} \right] = A_1^2 + E\left[ {\varepsilon _1^2} \right]

Yukarıda \displaystyle {Y_1}’in beklenen değerinin tüm elemanları sabit olup \displaystyle {A_1}’e eşit olduğu kabul edilmiştir.

Hata terimlerinin beklenen değeri 0 olduğu için diğer elemanlar ile çapraz çarpımlarının beklenen değerleri de yukarıda gösterildiği gibi 0’a eşit olacaktır.

Benzer yaklaşımlar diğer \displaystyle {Y_t}’ler için de aynı gerekçelerle kabul edilmiş olup aşağıda gösterilmiştir.

\displaystyle E\left[ {{Y_2}} \right] = {A_2} = \beta  + \rho \beta  + \alpha T + \rho \alpha T + {\rho ^2}{Y_0}

\displaystyle {\left( {E\left[ {{Y_2}} \right]} \right)^2} = A_2^2 = {\left( {\beta  + \rho \beta  + \alpha T + \rho \alpha T + {\rho ^2}{Y_0}} \right)^2}

\displaystyle E\left[ {Y_2^2} \right] = E\left[ {{{\left( {{A_2} + \rho {\varepsilon _1} + {\varepsilon _2}} \right)}^2}} \right]

\displaystyle E\left[ {Y_2^2} \right] = A_2^2 + {\rho ^2}E\left[ {\varepsilon _1^2} \right] + E\left[ {\varepsilon _2^2} \right]

 

\displaystyle E\left[ {{Y_3}} \right] = {A_3} = \beta  + \rho \beta  + {\rho ^2}\beta  + \alpha T + \rho \alpha T + {\rho ^2}\alpha T + {\rho ^3}{Y_0}

\displaystyle {\left( {E\left[ {{Y_3}} \right]} \right)^2} = A_3^2 = {\left( {\beta  + \rho \beta  + {\rho ^2}\beta  + \alpha T + \rho \alpha T + {\rho ^2}\alpha T + {\rho ^3}{Y_0}} \right)^2}

\displaystyle E\left[ {Y_3^2} \right] = E\left[ {{{\left( {{A_3} + {\rho ^2}{\varepsilon _1} + \rho {\varepsilon _2} + {\varepsilon _3}} \right)}^2}} \right]

\displaystyle E\left[ {Y_3^2} \right] = A_3^2 + {\rho ^4}E\left[ {\varepsilon _1^2} \right] + {\rho ^2}E\left[ {\varepsilon _2^2} \right] + E\left[ {\varepsilon _3^2} \right]

t dönemi temsil etmek üzere \displaystyle {Y_t}’nin varyansı aşağıdaki gibi hesaplanmaktaydı.

\displaystyle Var\left( {{Y_t}} \right) = E\left[ {{{\left( {{Y_t} - E\left( {{Y_t}} \right)} \right)}^2}} \right]

\displaystyle Var\left( {{Y_t}} \right) = E\left[ {Y_t^2} \right] - {\left( {E\left( {{Y_t}} \right)} \right)^2}

Yukarıdaki sonuçtan yola çıkılarak gözlem sayısı n=3 basitleştirici varsayımı altında her alternatif t gözlemi için \displaystyle {Y_t}’nin varyansı aşağıda analiz edilmektedir.

t=1:

\displaystyle Var\left( {{Y_1}} \right) = E\left[ {Y_1^2} \right] - {\left( {E\left( {{Y_1}} \right)} \right)^2}

\displaystyle E\left[ {Y_1^2} \right] = A_1^2 + E\left[ {\varepsilon _1^2} \right] ve \displaystyle {\left( {E\left[ {{Y_1}} \right]} \right)^2} = A_1^2 olduğundan denklemde yerine yazılır.

\displaystyle Var\left( {{Y_1}} \right) = \cancel{{A_1^2}} + E\left[ {\varepsilon _1^2} \right] - \cancel{{A_1^2}}

\displaystyle Var\left( {{Y_1}} \right) = E\left[ {\varepsilon _1^2} \right]

\displaystyle Var\left( {{Y_1}} \right) = \delta _\varepsilon ^2

t=2:

\displaystyle Var\left( {{Y_2}} \right) = E\left[ {Y_2^2} \right] - {\left( {E\left( {{Y_2}} \right)} \right)^2}

\displaystyle E\left[ {Y_2^2} \right] = A_2^2 + {\rho ^2}E\left[ {\varepsilon _1^2} \right] + E\left[ {\varepsilon _2^2} \right] ve \displaystyle {\left( {E\left[ {{Y_2}} \right]} \right)^2} = A_2^2 olduğundan denklemde yerine yazılır.

\displaystyle Var\left( {{Y_2}} \right) = \cancel{{A_2^2}} + {\rho ^2}E\left[ {\varepsilon _1^2} \right] + E\left[ {\varepsilon _2^2} \right] - \cancel{{A_2^2}}

\displaystyle Var\left( {{Y_2}} \right) = {\rho ^2}E\left[ {\varepsilon _1^2} \right] + E\left[ {\varepsilon _2^2} \right]

\displaystyle Var\left( {{Y_2}} \right) = {\rho ^2}\delta _\varepsilon ^2 + \delta _\varepsilon ^2

t=3:

\displaystyle Var\left( {{Y_3}} \right) = E\left[ {Y_3^2} \right] - {\left( {E\left( {{Y_3}} \right)} \right)^2}

\displaystyle E\left[ {Y_3^2} \right] = A_3^2 + {\rho ^4}E\left[ {\varepsilon _1^2} \right] + {\rho ^2}E\left[ {\varepsilon _2^2} \right] + E\left[ {\varepsilon _3^2} \right] ve \displaystyle {\left( {E\left[ {{Y_3}} \right]} \right)^2} = A_3^2 olduğundan denklemde yerine yazılır.

\displaystyle Var\left( {{Y_3}} \right) = \cancel{{A_3^2}} + {\rho ^4}E\left[ {\varepsilon _1^2} \right] + {\rho ^2}E\left[ {\varepsilon _2^2} \right] + E\left[ {\varepsilon _3^2} \right] - \cancel{{A_3^2}}

\displaystyle Var\left( {{Y_3}} \right) = {\rho ^4}E\left[ {\varepsilon _1^2} \right] + {\rho ^2}E\left[ {\varepsilon _2^2} \right] + E\left[ {\varepsilon _3^2} \right]

\displaystyle Var\left( {{Y_3}} \right) = {\rho ^4}\delta _\varepsilon ^2 + {\rho ^2}\delta _\varepsilon ^2 + \delta _\varepsilon ^2

Yukarıdaki sonuçlar her alternatif t gözlemi için aşağıda özetlenmiştir.

\displaystyle t = 1 \to Var\left( {{Y_1}} \right) = \delta _\varepsilon ^2

\displaystyle t = 2 \to Var\left( {{Y_2}} \right) = {\rho ^2}\delta _\varepsilon ^2 + \delta _\varepsilon ^2

\displaystyle t = 3 \to Var\left( {{Y_3}} \right) = {\rho ^4}\delta _\varepsilon ^2 + {\rho ^2}\delta _\varepsilon ^2 + \delta _\varepsilon ^2

Yukarıda \displaystyle {Y_t}’nin varyansları için gösterilen süreç matematiksel olarak aşağıdaki gibi daha sade ifade edilebilir.

\displaystyle Var\left( {{Y_t}} \right) = \sum\limits_{t = 1}^{n = 3} {{\rho ^{2t - 2}}\delta _\varepsilon ^2}

Gözlem sayısı n iken yukarıdaki ifade aşağıdaki gibi genel bir forma dönüştürülebilir.

\displaystyle Var\left( {{Y_t}} \right) = \sum\limits_{t = 1}^n {{\rho ^{2t - 2}}\delta _\varepsilon ^2}

Yukarıda \displaystyle {Y_t}’nin varyansı için gösterilen süreç sabit içermeyen rassal yürüyüş modelinde elde edilen \displaystyle {Y_t}’nin varyansı ile aynıdır. Dolayısıyla sabit ve trend içeren rassal yürüyüş modeli için \displaystyle {Y_t}’nin varyansı da sabit içermeyen rassal yürüyüş modeli için \displaystyle {Y_t}’nin varyansı ile aynı değerlendirmelere tabi olacaktır.

\displaystyle {Y_t}’nin varyansı birim kök yokken aşağıdaki gibi olacaktır.

\displaystyle Var\left( {{Y_t}} \right) = \sum\limits_{t = 1}^{n = \infty } {{\rho ^{2t - 2}}\delta _\varepsilon ^2}

\displaystyle Var\left( {{Y_t}} \right) = \frac{{\delta _\varepsilon ^2}}{{1 - {\rho ^2}}}

\displaystyle {Y_t}’nin varyansı birim kök varken aşağıdaki gibi olacaktır.

\displaystyle Var\left( {{Y_t}} \right) = \sum\limits_{t = 1}^n {{\rho ^{2t - 2}}\delta _\varepsilon ^2}  = \sum\limits_{t = 1}^n {{1^{2t - 2}}\delta _\varepsilon ^2}

\displaystyle Var\left( {{Y_t}} \right) = \sum\limits_{t = 1}^n {\delta _\varepsilon ^2}

\displaystyle Var\left( {{Y_t}} \right) = t \cdot \delta _\varepsilon ^2

4.3. Kovaryans için Durağanlık Analizi

Kovaryans için durağanlık analizine geçmeden önce kovaryansın durağanlık analizinde yararlanılacak olan \displaystyle {Y_t}’nin \displaystyle {Y_1}, \displaystyle {Y_2}, \displaystyle {Y_3}, ve \displaystyle {Y_4} için beklenen değerleri, beklenen değerlerinin çapraz çarpımları ve çapraz çarpımlarının beklenen değerleri aşağıda sunulmuştur.

\displaystyle E\left[ {{Y_1}} \right] = {A_1} = \rho {Y_0} + \beta  + \alpha T

\displaystyle E\left[ {{Y_2}} \right] = {A_2} = {\rho ^2}{Y_0} + \rho \beta  + \beta  + \rho \alpha T + \alpha T

\displaystyle E\left[ {{Y_1}} \right]E\left[ {{Y_2}} \right] = {A_1}{A_2}

\displaystyle E\left[ {{Y_1}{Y_2}} \right] = E\left[ {\left( {{A_1} + {\varepsilon _1}} \right)\left( {{A_2} + \rho {\varepsilon _1} + {\varepsilon _2}} \right)} \right]

\displaystyle E\left[ {{Y_1}{Y_2}} \right] = E\left[ {{A_1}{A_2} + {A_1}\rho {\varepsilon _1} + {A_1}{\varepsilon _2} + {A_2}{\varepsilon _1} + \rho \varepsilon _1^2} \right]

\displaystyle E\left[ {{Y_1}{Y_2}} \right] = {A_1}{A_2} + {A_1}\rho \overbrace {E\left[ {{\varepsilon _1}} \right]}^0 + {A_1}\overbrace {E\left[ {{\varepsilon _2}} \right]}^0 + {A_2}\underbrace {E\left[ {{\varepsilon _1}} \right]}_0 + \rho E\left[ {\varepsilon _1^2} \right]

\displaystyle E\left[ {{Y_1}{Y_2}} \right] = {A_1}{A_2} + \rho E\left[ {\varepsilon _1^2} \right]

Yukarıda \displaystyle {Y_1} ve \displaystyle {Y_2}’nin beklenen değerlerinin tüm elemanları sabit olup \displaystyle {A_1} ve \displaystyle {A_2}’ye eşit olduğu kabul edilmiştir.

Hata terimlerinin beklenen değeri 0 olduğu için diğer elemanlar ile çapraz çarpımlarının beklenen değerleri de yukarıda gösterildiği gibi 0’a eşit olacaktır.

Benzer yaklaşımlar diğer \displaystyle {Y_t}’ler için de aynı gerekçelerle kabul edilmiş olup aşağıda gösterilmiştir.

\displaystyle E\left[ {{Y_2}} \right] = {A_2} = {\rho ^2}{Y_0} + \rho \beta  + \beta  + \rho \alpha T + \alpha T

\displaystyle E\left[ {{Y_3}} \right] = {A_3} = {\rho ^3}{Y_0} + \beta  + \rho \beta  + {\rho ^2}\beta  + \alpha T + \rho \alpha T + {\rho ^2}\alpha T

\displaystyle E\left[ {{Y_2}} \right]E\left[ {{Y_3}} \right] = {A_2}{A_3}

\displaystyle E\left[ {{Y_2}{Y_3}} \right] = E\left[ {\left( {{A_2} + \rho {\varepsilon _1} + {\varepsilon _2}} \right)\left( {{A_3} + {\rho ^2}{\varepsilon _1} + \rho {\varepsilon _2} + {\varepsilon _3}} \right)} \right]

\displaystyle E\left[ {{Y_2}{Y_3}} \right] = {A_2}{A_3} + {\rho ^3}E\left[ {\varepsilon _1^2} \right] + \rho E\left[ {\varepsilon _2^2} \right]

 

\displaystyle E\left[ {{Y_3}} \right] = {A_3} = {\rho ^3}{Y_0} + \beta  + \rho \beta  + {\rho ^2}\beta  + \alpha T + \rho \alpha T + {\rho ^2}\alpha T

\displaystyle E\left[ {{Y_4}} \right] = {A_4} = {\rho ^4}{Y_0} + \beta  + \rho \beta  + {\rho ^2}\beta  + {\rho ^3}\beta  + \alpha T + \rho \alpha T + {\rho ^2}\alpha T + {\rho ^3}\alpha T

\displaystyle E\left[ {{Y_3}} \right]E\left[ {{Y_4}} \right] = {A_3}{A_4}

\displaystyle E\left[ {{Y_3}{Y_4}} \right] = E\left[ {\left( {{A_3} + {\rho ^2}{\varepsilon _1} + \rho {\varepsilon _2} + {\varepsilon _3}} \right)\left( {{A_4} + {\rho ^3}{\varepsilon _1} + {\rho ^2}{\varepsilon _2} + \rho {\varepsilon _3} + {\varepsilon _4}} \right)} \right]

\displaystyle E\left[ {{Y_3}{Y_4}} \right] = {A_3}{A_4} + {\rho ^5}E\left[ {\varepsilon _1^2} \right] + {\rho ^3}E\left[ {\varepsilon _2^2} \right] + \rho E\left[ {\varepsilon _3^2} \right]

k iki dönem arasındaki gecikme uzunluğunu, t dönemi temsil etmek üzere \displaystyle {Y_t}’nin kovaryansı aşağıda gösterildiği gibi hesaplanmaktaydı.

\displaystyle Cov\left( {{Y_t},{Y_{t + k}}} \right) = E\left[ {\left( {{Y_t} - E\left( {{Y_t}} \right)} \right)\left( {{Y_{t + k}} - E\left( {{Y_{t + k}}} \right)} \right)} \right]

\displaystyle  = E\left[ {{Y_t}{Y_{t + k}} - {Y_t}E\left( {{Y_{t + k}}} \right) - E\left( {{Y_t}} \right){Y_{t + k}} + E\left( {{Y_t}} \right)E\left( {{Y_{t + k}}} \right)} \right]

\displaystyle  = E\left[ {{Y_t}{Y_{t + k}}} \right] - E\left[ {{Y_t}} \right]E\left( {{Y_{t + k}}} \right) - \cancel{{E\left( {{Y_t}} \right)E\left[ {{Y_{t + k}}} \right]}} + \cancel{{E\left( {{Y_t}} \right)E\left( {{Y_{t + k}}} \right)}}

\displaystyle  = E\left[ {{Y_t}{Y_{t + k}}} \right] - E\left[ {{Y_t}} \right]E\left( {{Y_{t + k}}} \right)

\displaystyle Cov\left( {{Y_t},{Y_{t + k}}} \right) = E\left[ {{Y_t}{Y_{t + k}}} \right] - E\left[ {{Y_t}} \right]E\left( {{Y_{t + k}}} \right)

Yukarıdaki sonuçtan yola çıkılarak gecikme uzunluğu k=1 ve gözlem sayısı n=4 basitleştirici varsayımı altında her alternatif t gözlemi için \displaystyle {Y_t}’nin kovaryansı aşağıda analiz edilmektedir.

\displaystyle Cov\left( {{Y_t},{Y_{t + 1}}} \right) = E\left[ {{Y_t}{Y_{t + 1}}} \right] - E\left[ {{Y_t}} \right]E\left( {{Y_{t + 1}}} \right)

k=1 ve t=1:

\displaystyle Cov\left( {{Y_1},{Y_2}} \right) = E\left[ {{Y_1}{Y_2}} \right] - E\left[ {{Y_1}} \right]E\left( {{Y_2}} \right)

\displaystyle E\left[ {{Y_1}{Y_2}} \right] = {A_1}{A_2} + \rho E\left[ {\varepsilon _1^2} \right] ve \displaystyle E\left[ {{Y_1}} \right]E\left[ {{Y_2}} \right] = {A_1}{A_2} olduğundan denklemde yerine yazılır.

\displaystyle Cov\left( {{Y_1},{Y_2}} \right) = \cancel{{{A_1}{A_2}}} + \rho E\left[ {\varepsilon _1^2} \right] - \cancel{{{A_1}{A_2}}}

\displaystyle Cov\left( {{Y_1},{Y_2}} \right) = \rho E\left[ {\varepsilon _1^2} \right]

\displaystyle Cov\left( {{Y_1},{Y_2}} \right) = \rho \delta _\varepsilon ^2

k=1 ve t=2:

\displaystyle Cov\left( {{Y_2},{Y_3}} \right) = E\left[ {{Y_2}{Y_3}} \right] - E\left[ {{Y_2}} \right]E\left( {{Y_3}} \right)

\displaystyle E\left[ {{Y_2}{Y_3}} \right] = {A_2}{A_3} + {\rho ^3}E\left[ {\varepsilon _1^2} \right] + \rho E\left[ {\varepsilon _2^2} \right] ve \displaystyle E\left[ {{Y_2}} \right]E\left[ {{Y_3}} \right] = {A_2}{A_3} olduğundan denklemde yerine yazılır.

\displaystyle Cov\left( {{Y_2},{Y_3}} \right) = \cancel{{{A_2}{A_3}}} + {\rho ^3}E\left[ {\varepsilon _1^2} \right] + \rho E\left[ {\varepsilon _2^2} \right] - \cancel{{{A_2}{A_3}}}

\displaystyle Cov\left( {{Y_2},{Y_3}} \right) = {\rho ^3}E\left[ {\varepsilon _1^2} \right] + \rho E\left[ {\varepsilon _2^2} \right]

\displaystyle Cov\left( {{Y_2},{Y_3}} \right) = {\rho ^3}\delta _\varepsilon ^2 + \rho \delta _\varepsilon ^2

k=1 ve t=3:

\displaystyle Cov\left( {{Y_3},{Y_4}} \right) = E\left[ {{Y_3}{Y_4}} \right] - E\left[ {{Y_3}} \right]E\left( {{Y_4}} \right)

\displaystyle E\left[ {{Y_3}{Y_4}} \right] = {A_3}{A_4} + {\rho ^5}E\left[ {\varepsilon _1^2} \right] + {\rho ^3}E\left[ {\varepsilon _2^2} \right] + \rho E\left[ {\varepsilon _3^2} \right] ve \displaystyle E\left[ {{Y_3}} \right]E\left[ {{Y_4}} \right] = {A_3}{A_4} olduğundan denklemde yerine yazılır.

\displaystyle Cov\left( {{Y_3},{Y_4}} \right) = \cancel{{{A_3}{A_4}}} + {\rho ^5}E\left[ {\varepsilon _1^2} \right] + {\rho ^3}E\left[ {\varepsilon _2^2} \right] + \rho E\left[ {\varepsilon _3^2} \right] - \cancel{{{A_3}{A_4}}}

\displaystyle Cov\left( {{Y_3},{Y_4}} \right) = {\rho ^5}E\left[ {\varepsilon _1^2} \right] + {\rho ^3}E\left[ {\varepsilon _2^2} \right] + \rho E\left[ {\varepsilon _3^2} \right]

\displaystyle Cov\left( {{Y_3},{Y_4}} \right) = {\rho ^5}\delta _\varepsilon ^2 + {\rho ^3}\delta _\varepsilon ^2 + \rho \delta _\varepsilon ^2

Yukarıdaki sonuçlar k=1 iken her alternatif t gözlemi için aşağıda özetlenmiştir.

\displaystyle k = 1{\text{ ve }}t = 1 \to Cov\left( {{Y_1},{Y_2}} \right) = \rho \delta _\varepsilon ^2

\displaystyle k = 1{\text{ ve }}t = 2 \to Cov\left( {{Y_2},{Y_3}} \right) = {\rho ^3}\delta _\varepsilon ^2 + \rho \delta _\varepsilon ^2

\displaystyle k = 1{\text{ ve }}t = 3 \to Cov\left( {{Y_3},{Y_4}} \right) = {\rho ^5}\delta _\varepsilon ^2 + {\rho ^3}\delta _\varepsilon ^2 + \rho \delta _\varepsilon ^2

Yukarıda gecikme uzunluğu k=1 n=4 iken \displaystyle {Y_t}’nin kovaryansları için gösterilen süreç matematiksel olarak aşağıdaki gibi daha sade ifade edilebilir.

\displaystyle Cov\left( {{Y_t},{Y_{t + k}}} \right) = \sum\limits_{t = 1}^{n - k} {{\rho ^{2t + k - 2}}\delta _\varepsilon ^2}

\displaystyle Cov\left( {{Y_t},{Y_{t + 1}}} \right) = \sum\limits_{t = 1}^3 {{\rho ^{2t - 1}}\delta _\varepsilon ^2}

Gecikme uzunluğu k=1 ve gözlem sayısı n iken yukarıdaki ifade aşağıdaki gibi genel bir forma dönüştürülebilir.

\displaystyle Cov\left( {{Y_t},{Y_{t + k}}} \right) = \sum\limits_{t = 1}^{n - k} {{\rho ^{2t + k - 2}}\delta _\varepsilon ^2}

\displaystyle Cov\left( {{Y_t},{Y_{t + 1}}} \right) = \sum\limits_{t = 1}^{n - 1} {{\rho ^{2t - 1}}\delta _\varepsilon ^2}

Yukarıda \displaystyle {Y_t}’nin kovaryansı için gösterilen süreç sabit içermeyen rassal yürüyüş modelinde elde edilen \displaystyle {Y_t}’nin kovaryansı ile aynıdır. Dolayısıyla sabit ve trend içeren rassal yürüyüş modeli için \displaystyle {Y_t}’nin kovaryansı da sabit içermeyen rassal yürüyüş modeli için \displaystyle {Y_t}’nin kovaryansı ile aynı değerlendirmelere tabi olacaktır.

\displaystyle {Y_t}’nin kovaryansı birim kök yokken aşağıdaki gibi olacaktır.

\displaystyle Cov\left( {{Y_t},{Y_{t + k}}} \right) = \sum\limits_{t = 1}^{n - k} {{\rho ^{2t + k - 2}}\delta _\varepsilon ^2}

\displaystyle Cov\left( {{Y_t},{Y_{t + k}}} \right) = {\rho ^k}\frac{{\delta _\varepsilon ^2}}{{1 - {\rho ^2}}}

\displaystyle Cov\left( {{Y_t},{Y_{t + k}}} \right) = {\rho ^k} \cdot Var\left( {{Y_t}} \right)

\displaystyle {Y_t}’nin kovaryansı birim kök varken aşağıdaki gibi olacaktır.

\displaystyle Cov\left( {{Y_t},{Y_{t + k}}} \right) = \sum\limits_{t = 1}^{n - k} {{\rho ^{2t + k - 2}}\delta _\varepsilon ^2}  = \sum\limits_{t = 1}^{n - k} {{1^{2t + k - 2}}\delta _\varepsilon ^2}

\displaystyle Cov\left( {{Y_t},{Y_{t + k}}} \right) = \sum\limits_{t = 1}^{n - k} {\delta _\varepsilon ^2}

\displaystyle Cov\left( {{Y_t},{Y_{t + k}}} \right) = \left( {t - k} \right) \cdot \delta _\varepsilon ^2

 

 

İshak Kutlu, Ekonometri Notları

Yorum Yap

İshak Kutlu

Veri & Otomasyon Uzmanı

Kurumsal iş akışlarını düzenleyen veri odaklı otomasyon sistemleri ve makine öğrenmesi tabanlı çözümler geliştirir. Gerçek operasyonel süreçler için ölçeklenebilir ve izlenebilir yapılar tasarlar.

YouTube   YouTube   YouTube   YouTube