Ein reines saisonales Modell anpassen
Wie bei anderen Modellen kannst du saisonale Modelle in R mit dem Befehl sarima() aus dem Paket astsa schätzen.
Um ein Gefühl dafür zu bekommen, wie reine saisonale Modelle funktionieren, betrachtest du am besten simulierte Daten. Wir haben 250 Beobachtungen aus einem reinen saisonalen Modell erzeugt, gegeben durch $$X_t = .9 X_{t-12} + W_t + .5 W_{t-12}\,,$$ was wir als SARMA(P = 1, Q = 1)S = 12 bezeichnen würden. Drei Jahre an Daten sowie die Modell-ACF und -PACF sind für dich geplottet.
Du vergleichst die Stichprobenwerte von ACF und PACF der generierten Daten mit den angegebenen wahren Werten.
Das Paket astsa ist bereits geladen und die generierten Daten liegen in x.
Diese Übung ist Teil des Kurses
ARIMA-Modelle in R
Anleitung zur Übung
- Verwende
acf2(), um die Stichproben-ACF und -PACF der generierten Daten bis Lag 60 zu plotten, und vergleiche sie mit den tatsächlichen Werten. Um bis Lag 60 zu schätzen, setze das Argumentmax.lagauf60. - Passe das Modell an die generierten Daten mit
sarima()an. Gib in deinemsarima()-Befehl zusätzlich zu den Argumentenp,dundqauchP,D,QundSan (beachte die Groß- und Kleinschreibung in R).
Interaktive Übung
Vervollständige den Beispielcode, um diese Übung erfolgreich abzuschließen.
# Plot sample P/ACF to lag 60 and compare to the true values
acf2(___, max.lag = ___)
# Fit the seasonal model to x
sarima(x, p = 0, d = 0, q = 0, P = ___, D = 0, Q = ___, S = ___)