LoslegenKostenlos starten

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

<Kurs>ARIMA-Modelle in R</Kurs>
Kurs ansehen

Übungsanweisungen

  • 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 Argument max.lag auf 60.
  • Passe das Modell an die generierten Daten mit sarima() an. Gib in deinem sarima()-Befehl zusätzlich zu den Argumenten p, d und q auch P, D, Q und S an (beachte die Groß- und Kleinschreibung in R).

Interaktive praktische Übung

Versuche dich an dieser Übung, indem du diesen Beispielcode vervollständigst.

# 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 = ___)
Code bearbeiten und ausführen