AR(1) modelini uydurma
ACF ve PACF çiftini, bir ARMA modelinin \(p\) ve \(q\) mertebelerini belirlemeye yardımcı olmak için kullandığını unutma. Aşağıdaki tablo sonuçların özetidir:
| AR(\(p\)) | MA(\(q\)) | ARMA(\(p,q\)) | |
|---|---|---|---|
| ACF | Yavaşça sönümlenir | Gecikme \(q\)'den sonra kesilir |
Yavaşça sönümlenir |
| PACF | Gecikme \(p\)'den sonra kesilir |
Yavaşça sönümlenir | Yavaşça sönümlenir |
Bu egzersizde, AR(1) modelinden veri üreteceksin: $$X_t = .9 X_{t-1} + W_t,$$ simüle edilen veriye ve örnek ACF-PACF çiftine bakarak mertebeyi belirleyeceksin. Sonrasında modeli uydurup, tahmin edilen parametreleri gerçek parametrelerle karşılaştıracaksın.
Kurs boyunca, modellere kolayca uyum sağlamak için astsa paketindeki sarima() fonksiyonunu kullanacaksın. Bu komut, artık daha sonra bölümde ele alacağımız tanı grafiğini üretir; şimdilik bunu görmezden gelebilirsin.
Bu egzersiz
R ile ARIMA Modelleri
kursunun bir parçasıdırEgzersiz talimatları
- astsa paketi önceden yüklüdür.
- AR parametresi .9 olan bir AR(1) modelinden 100 gözlem üretmek için önceden yazılmış
arima.sim()komutunu kullan. Bunuxolarak kaydet. - Üretilen veriyi
plot()ile görselleştir. - Örnek ACF ve PACF çiftlerini
astsapaketindekiacf2()komutuyla çiz. - Daha önce üretilen veriye AR(1) uydurmak için
astsaiçindensarima()yı kullan. t-tablosunu incele ve tahminleri gerçek değerlerle karşılaştır. Örneğin, zaman serisixiçindeyse, verilere AR(1) uydurmak içinsarima(x, p = 1, d = 0, q = 0)ya da kısacasarima(x, 1, 0, 0)kullan.
Uygulamalı interaktif egzersiz
Bu örnek kodu tamamlayarak bu egzersizi bitirin.
# Generate 100 observations from the AR(1) model
x <- arima.sim(model = list(order = c(1, 0, 0), ar = .9), n = 100)
# Plot the generated data
# Plot the sample P/ACF pair
# Fit an AR(1) to the data and examine the t-table