Eğilim ve değişen varyansla başa çıkma
Burada, durağan olmayan verileri, aşağıdaki gibi getiri ya da büyüme oranını hesaplayarak durağanlığa zorlayacağız.
Zaman serileri çoğunlukla $$X_t = (1 + p_t) X_{t-1}$$ şeklinde üretilir; bu da \(t\) anında gözlenen zaman serisi değerinin, \(t-1\) anındaki değere ve \(t\) anındaki küçük yüzde değişim \(p_t\)'ye eşit olduğu anlamına gelir.
Basit deterministik bir örnek, sabit faiz oranı \(p\) ile bankaya para yatırmaktır. Bu durumda \(X_t\), başlangıç yatırımı \(X_0\) olan bir hesabın \(t\) dönemindeki değeridir.
Genellikle \(p_t\), bir zaman serisinin getirisi veya büyüme oranı olarak adlandırılır ve bu süreç çoğu zaman kararlıdır.
Bu dersin kapsamı dışında kalan nedenlerle, büyüme oranı \(p_t\)'nin $$Y_t = \log X_t - \log X_{t-1} \approx p_t$$ ile yaklaşıklanabileceği gösterilebilir.
R'da \(p_t\) sıklıkla diff(log(x)) olarak hesaplanır ve bunu çizdirmek tek satırda plot(diff(log(x))) ile yapılabilir.
Bu egzersiz
R ile ARIMA Modelleri
kursunun bir parçasıdırEgzersiz talimatları
- Daha önce olduğu gibi, astsa ve xts paketleri önceden yüklüdür.
- Çoklu şekilli bir grafik üret: (1) üç aylık ABD GSYİH (
gnp) verilerini çiz ve durağan olmadığını fark et; (2)diff()velog()kullanarak ABD GSYİH'nin yaklaşık büyüme oranını çiz. - Çoklu şekilli bir grafik kullanarak (1) günlük DJIA kapanışlarını (
djia$Close) çiz ve bunun durağan olmadığını fark et. Veriler birxtsnesnesidir. Ardından (2)diff()velog()kullanarak yaklaşık DJIA getirilerini çiz. Bu, GSYİH'nin büyüme oranıyla nasıl karşılaştırılır?
Uygulamalı interaktif egzersiz
Bu örnek kodu tamamlayarak bu egzersizi bitirin.
# astsa and xts are preloaded
# Plot GNP series (gnp) and its growth rate
par(mfrow = c(2,1))
plot(gnp)
# Plot DJIA closings (djia$Close) and its returns
par(mfrow = c(2,1))
plot(djia$Close)