Bir ARMA modelini uydurma
Şimdi AR modelini ve MA modelini ARMA modelinde birleştirmeye hazırsın. ARMA(2,1) modelinden veri ürettik, $$X_t = X_{t-1} - .9 X_{t-2} + W_t + .8 W_{t-1}, $$ x <- arima.sim(model = list(order = c(2, 0, 1), ar = c(1, -.9), ma = .8), n = 250). Olası bir modeli belirlemek için simüle edilen veriye ve örnek ACF ile PACF çiftine bak.
ARMA(\(p, q\)) modellerinde hem teorik ACF hem de PACF’in yavaşça sönümlendiğini (tail off) unutma. Bu durumda dereceleri veriden ayırt etmek zor olabilir ve örnek ACF ya da örnek PACF’in kesin olarak kesilip kesilmediği veya sönümlenip sönümlenmediği net olmayabilir. Bu örnekte gerçek model derecelerini bildiğin için üretilen veriye bir ARMA(2,1) uydur. Genel modelleme stratejileri kursun ilerleyen kısımlarında daha ayrıntılı tartışılacak.
Bu egzersiz
R ile ARIMA Modelleri
kursunun bir parçasıdırEgzersiz talimatları
- astsa paketi önceden yüklendi.
xiçinde 250 adet ARMA(2,1) gözlemi var. - Önceki egzersizlerde olduğu gibi,
xiçindeki üretilmiş veriyi çizmek içinplot()kullan ve örnek ACF ile PACF çiftlerini görmek içinacf2()kullan. - Üretilen veriye ARMA(2,1) uydurmak için
sarima()kullan. t-tablosunu incele ve tahminleri gerçek değerlerle karşılaştır.
Uygulamalı interaktif egzersiz
Bu örnek kodu tamamlayarak bu egzersizi bitirin.
# astsa is preloaded
# Plot x
# Plot the sample P/ACF of x
# Fit an ARMA(2,1) to the data and examine the t-table