Kapitel 23Zeitreihenökonometrie

Einleitung

Die Ökonometrie aus Kapitel 10 ruhte auf einer einzigen stillen Annahme: Die Beobachtungen sind unabhängige Ziehungen aus einer festen Grundgesamtheit. Randomisieren Sie eine Behandlung, finden Sie ein Instrument, nutzen Sie eine Unstetigkeit, und die i.i.d.-Stichprobenziehung, die die Asymptotik rechtfertigt, erledigt den Rest. Die meisten Daten, mit denen ein Makroökonom oder ein Finanzmarktforscher tatsächlich arbeitet, kommen nicht so daher. Sie kommen in einer Reihenfolge. Vierteljährliches BIP, monatliche Inflation, der Tagesgeldsatz, tägliche Aktienrenditen: Jede Beobachtung steht zeitlich neben ihren Nachbarn, und was im letzten Quartal geschah, bestimmt dieses hier mit.

Diese Reihenfolge verändert alles Weitere. Serielle Abhängigkeit bedeutet, dass aufeinanderfolgende Beobachtungen einander überlappende Information tragen. Der effektive Stichprobenumfang ist deshalb kleiner, als die Zahl der Datenpunkte nahelegt, und die Standardfehlerformeln der Querschnitts-OLS unterschätzen die Unsicherheit. Schlimmer noch: Viele ökonomische Reihen weisen einen Trend auf und wandern über Jahrzehnte nach oben, ohne einen festen Mittelwert, zu dem sie zurückkehren könnten. Führen Sie eine gewöhnliche Regression über zwei solche Reihen aus, und Sie erhalten womöglich ein hohes $R^2$ und eine klare $t$-Statistik für Größen, die nichts miteinander zu tun haben. Auf geordneten Daten verliert das Querschnittsinstrumentarium an Effizienz, auf trendbehafteten Daten erzeugt es Beziehungen, die es gar nicht gibt.

Nötig ist ein anderes Instrumentarium, geordnet um eine Frage (welche Art von Prozess hat diese Reihe erzeugt?) und um eine zweite, multivariate Frage, sobald mehrere Reihen im Spiel sind: Welche Schocks haben das System bewegt, und lässt sich das überhaupt sagen? Die Antwort auf die erste Frage führt von der Stationarität über ARMA-Modelle zu Einheitswurzeltests. Die Antwort auf die zweite führt von der Vektorautoregression zu ihrer strukturellen Deutung, wo die Daten beschreiben, wie sich Variablen gemeinsam bewegen, aber nie sagen, welche sich zuerst bewegt hat.

Intuition

Warum das wichtig ist: Zeitlich geordnete Daten erinnern sich an ihre Vergangenheit, und die Verfahren, die diese Ordnung ignoriert haben, führen Sie in die Irre. Eine Regression, die das diesjährige BIP wie eine unabhängige Ziehung behandelt, liest immer wieder dieselbe langsame Drift und hält sie für einen Befund. Die Ordnung ernst zu nehmen heißt herauszufinden, ob sich eine Reihe wieder auf ein Zentrum einpendelt oder für immer davonwandert, und herauszufinden, welche von mehreren gemeinsam laufenden Reihen die anderen schiebt.

Am Ende dieses Kapitels werden Sie in der Lage sein:
  1. Kovarianzstationarität definieren und erklären, warum sie die Voraussetzung für die übliche Zeitreiheninferenz ist
  2. AR- und MA-Prozesse an ihren Autokorrelations- und partiellen Autokorrelationsmustern erkennen und aus einem ARMA-Modell prognostizieren
  3. Differenzstationäre von trendstationären Reihen unterscheiden, eine Einheitswurzel erkennen und die Falle der Scheinregression erklären
  4. Ein VAR in reduzierter Form spezifizieren, ein Ergebnis zur Granger-Kausalität lesen und eine Impulsantwort sowie eine Varianzzerlegung deuten
  5. Erklären, warum ein VAR in reduzierter Form nicht strukturell ist und wie eine Identifikationsrestriktion (Cholesky-Anordnung, Vorzeichenrestriktionen) die Antwort verändert
  6. Angeben, was Kointegration ist, das Fehlerkorrekturmodell zu einer Kointegrationsbeziehung aufschreiben und die Anpassungsgeschwindigkeit deuten
  7. Volatilitätscluster erkennen und erklären, wie ARCH- und GARCH-Modelle die bedingte Varianz prognostizierbar machen

Voraussetzungen: der Identifikationsrahmen, OLS und die Anmerkung zur seriellen Korrelation aus Kapitel 10 (Grundlagen der Ökonometrie); lineare Algebra (das VAR ist matrixwertig); Grundzüge der Wahrscheinlichkeitsrechnung. Die Intuition für stochastische Prozesse wird hier aufgebaut.

Die Verfahren unten gehörten nicht immer zur Volkswirtschaftslehre. Im Lauf der 1970er Jahre lernte das Fach, die Zeitdimension ernst zu nehmen, vor allem durch die Arbeiten von Clive Granger und Robert Engle zu langfristigen Beziehungen und zur Volatilität und durch Christopher Sims zur multivariaten Dynamik. Ihr Programm formte um, wie Makroökonomen mit Daten umgehen. Diese Denklinie, wer auf diese Ideen kam und gegen welchen Widerstand, gehört in den Band zur Ideengeschichte und wird dort im Kapitel zur Informationsökonomik und zur spieltheoretischen Revolution nachgezeichnet.

23.1 Stationarität und die Bausteine

Bevor sich irgendein Modell schätzen lässt, muss die Reihe die Art von Gegenstand sein, die ein Modell beschreiben kann. Der ordnende Gedanke ist die Stationarität. Ein stationärer Prozess ist einer, dessen wahrscheinlichkeitstheoretischer Charakter über die Zeit stabil bleibt: Am Anfang der Stichprobe erzeugt derselbe Mechanismus die Daten wie an ihrem Ende. Trifft das zu, dann sagt die Vergangenheit etwas über die Zukunft, und zwar auf eine Weise, die sich nicht selbst verschiebt. Trifft es nicht zu, dann zielen die geschätzten Parameter auf ein bewegliches Ziel.

Kovarianzstationarität (schwache Stationarität). Ein Prozess $\{y_t\}$ ist kovarianzstationär, wenn sein Mittelwert konstant ist, seine Varianz konstant und endlich ist und die Kovarianz zwischen $y_t$ und $y_{t-k}$ nur vom Lag $k$ abhängt und nicht vom Zeitpunkt $t$. Die ersten beiden Momente und die Autokovarianzstruktur verschieben sich im Zeitverlauf nicht.
Strikte Stationarität. Die stärkere Forderung, dass die gesamte gemeinsame Verteilung eines beliebigen Blocks von Beobachtungen gegenüber Zeitverschiebungen invariant ist: $(y_t, \ldots, y_{t+m})$ hat dieselbe Verteilung wie $(y_{t+h}, \ldots, y_{t+m+h})$ für alle $h$. Strikte Stationarität impliziert Kovarianzstationarität, sofern die zweiten Momente existieren. Die Umkehrung gilt nicht. Die Zeitreihenökonometrie arbeitet fast durchgehend mit der schwachen Variante.
$$E[y_t]=\mu, \qquad \mathrm{Var}(y_t)=\sigma^2, \qquad \mathrm{Cov}(y_t, y_{t-k})=\gamma_k$$ (Gl. 23.1)
Weißes Rauschen. Eine Folge $\{\varepsilon_t\}$ mit Erwartungswert null, konstanter Varianz $\sigma^2$ und ohne Autokorrelation bei jedem Lag ungleich null ($\mathrm{Cov}(\varepsilon_t, \varepsilon_{t-k})=0$ für $k\neq 0$). Weißes Rauschen ist die rohe, unvorhersagbare Innovation, der Überraschungsanteil, aus dem jeder Prozess dieses Kapitels aufgebaut ist.
Lag-Operator. Der Operator $L$, der eine Reihe um eine Periode zurückschiebt: $L y_t = y_{t-1}$, und $L^k y_t = y_{t-k}$. Ein Polynom in $L$ steht für ein ganzes Muster von Lags, und die Wurzeln dieses Polynoms entscheiden, ob ein Prozess stationär ist.

Der Lag-Operator macht aus Dynamik Algebra. Ein Polynom $\phi(L)=1-\phi_1 L-\cdots-\phi_p L^p$, angewandt auf $y_t$, verschlüsselt eine ganze Autoregression in einem einzigen Symbol, und der Prozess ist genau dann stationär, wenn die Wurzeln von $\phi(z)=0$ außerhalb des Einheitskreises liegen. Liegt eine Wurzel auf dem Einheitskreis, scheitert die Stationarität. Das ist die Einheitswurzel aus §23.3.

Wold-Zerlegung. Jeder kovarianzstationäre Prozess lässt sich als Summe einer deterministischen Komponente und eines unendlichen gleitenden Durchschnitts von Innovationen weißen Rauschens schreiben, also als gewichtete Summe aktueller und vergangener Schocks. Sie rechtfertigt es, stationäre Reihen als ARMA-Prozesse zu modellieren, denn jede solche Reihe ist ein gefilterter Strom von Überraschungen.
$$y_t = \sum_{j=0}^{\infty} \psi_j \, \varepsilon_{t-j} + \kappa_t, \qquad \sum_{j=0}^{\infty} \psi_j^2 < \infty$$ (Gl. 23.2)

Dabei ist $\{\varepsilon_t\}$ weißes Rauschen, die $\psi_j$ sind quadratisch summierbare Gewichte mit $\psi_0=1$, und $\kappa_t$ ist derjenige Teil von $y_t$, der sich aus seiner eigenen Vergangenheit vollkommen vorhersagen lässt (ein deterministischer Trend oder eine Saisonfigur). Die quadratische Summierbarkeit der Gewichte hält die Varianz endlich, und sie ist der formale Gehalt des Satzes „der Einfluss eines alten Schocks verblasst“.

Intuition

Warum das wichtig ist: Ein stationärer Prozess ist einer, dessen statistischer Charakter nicht wandert. Schieben Sie ein Fenster über die Reihe, und das Bild darin sieht immer gleich aus. Der Mittelwert, den Sie aus dem ersten Jahrzehnt schätzen würden, stimmt mit dem aus dem letzten überein, das übliche Auf und Ab hat durchgehend dieselbe Größe, und wie stark der heutige Wert mit dem des Vormonats zusammenhängt, hängt nicht davon ab, um welchen Monat es geht. Weißes Rauschen ist der reinste Fall: lauter Überraschung ohne Gedächtnis. Und Wolds Ergebnis besagt, dass jede andere gutartige Reihe nichts weiter ist als ein Strom solcher Überraschungen, in dem die alten verblassen, während neue eintreffen. In der Abbildung unten können Sie zusehen, wie eine Reihe an ihrem Mittelwert verankert bleibt, und Sie können sie an die Kante schieben, an der die Verankerung reißt.

Abbildung 23.1. Simulator für stationäre und nichtstationäre Prozesse. Ein simulierter AR(1)-Pfad $y_t = \delta + \phi y_{t-1} + \varepsilon_t$ mit seinem Stichprobenmittel. Bei kleinem $\phi$ bleibt die Reihe dicht an ihrem Mittelwert (stationär). Mit $\phi \to 1$ wandert sie weiter fort und kehrt langsamer zurück. Bei $\phi = 1$ wird sie zum Random Walk, ohne Rückkehr zum Mittelwert und mit dauerhaften Schocks (der Grenzfall aus §23.3). Ziehen Sie den Schieberegler über $\phi = 1$ hinweg, schalten Sie die Drift um und ziehen Sie die Schocks neu, um zu bestätigen, dass das Verhalten generisch ist.

Ist eine Reihe erst einmal als stationär erkannt, lautet die nächste Frage, wie man sie modelliert. Zwei elementare Arten, auf die sich eine Reihe an ihre Vergangenheit erinnern kann, verbinden sich zu einer einzigen Familie.

23.2 ARMA-Prozesse

Eine stationäre Reihe erinnert sich, und es gibt zwei saubere Arten aufzuschreiben, woran. Sie kann ihre eigenen vergangenen Werte fortschreiben, sodass eine hohe Produktion im letzten Quartal die erwartete Produktion in diesem Quartal anhebt. Oder sie kann das Echo vergangener Überraschungen fortschreiben, also einen Schock, der ein paar Perioden braucht, bis er durchgelaufen ist. Das Erste ist die Autoregression, das Zweite der gleitende Durchschnitt. Kombiniert man beide, hat man das Standardmodell der univariaten Zeitreihenanalyse.

Autoregressiver Prozess (AR). Ein Prozess, in dem der aktuelle Wert eine lineare Funktion seiner eigenen vergangenen Werte plus einer Innovation weißen Rauschens ist. Ein AR(p) verwendet $p$ Lags. Die Reihe erinnert sich, indem sie ihr eigenes Niveau fortschreibt.
$$y_t = c + \phi_1 y_{t-1} + \phi_2 y_{t-2} + \cdots + \phi_p y_{t-p} + \varepsilon_t$$ (Gl. 23.3)
Moving-Average-Prozess (MA). Ein Prozess, in dem der aktuelle Wert eine lineare Funktion der aktuellen und der vergangenen Schocks weißen Rauschens ist. Ein MA(q) verwendet $q$ verzögerte Schocks. Die Reihe erinnert sich, indem sie das Echo vergangener Überraschungen trägt, das nach jeweils $q$ Perioden verklungen ist.
$$y_t = \mu + \varepsilon_t + \theta_1 \varepsilon_{t-1} + \cdots + \theta_q \varepsilon_{t-q}$$ (Gl. 23.4)
ARMA-Prozess. Ein Prozess, der beide Darstellungen verbindet, $p$ autoregressive Lags und $q$ Moving-Average-Terme, im Lag-Operator kompakt geschrieben als $\phi(L)y_t = \theta(L)\varepsilon_t$. ARMA-Modelle sind sparsam: Ein ARMA niedriger Ordnung kann eine Autokorrelation abbilden, für die viele reine AR- oder reine MA-Terme nötig wären.
$$\phi(L)\, y_t = \theta(L)\, \varepsilon_t, \qquad \phi(L)=1-\phi_1 L-\cdots-\phi_p L^p, \quad \theta(L)=1+\theta_1 L+\cdots+\theta_q L^q$$ (Gl. 23.5)
Autokorrelationsfunktion (ACF). Die Korrelation der Reihe mit ihren eigenen Lags, $\rho_k = \gamma_k/\gamma_0$, aufgetragen gegen $k$. Ihre Gestalt diagnostiziert die Moving-Average-Ordnung.
Partielle Autokorrelationsfunktion (PACF). Die Korrelation der Reihe mit ihrem Lag $k$, nachdem der lineare Einfluss aller dazwischenliegenden Lags herausgerechnet wurde. Ihre Gestalt diagnostiziert die autoregressive Ordnung.

ACF und PACF sind die diagnostischen Instrumente des Box-Jenkins-Ansatzes zum Modellaufbau, und die Regel, die sie nützlich macht, liegt in einem Gegensatz: darin, wie sie abfallen. Ein reiner AR-Prozess hat eine PACF, die nach Lag $p$ scharf abbricht (kontrolliert man die ersten $p$ Lags, trägt nichts weiter zur Vorhersage bei), während seine ACF allmählich abklingt. Ein reiner MA-Prozess ist das Spiegelbild: Seine ACF bricht nach Lag $q$ ab, während seine PACF abklingt. An diesem Gegensatz von Abbruch und Abklingen liest ein Anwender die Ordnung eines Modells aus den Daten ab.

Invertierbarkeit. Die Bedingung, unter der sich ein MA-Prozess als konvergente unendliche Autoregression umschreiben lässt. Sie gilt, wenn die Wurzeln von $\theta(z)=0$ außerhalb des Einheitskreises liegen. Die Invertierbarkeit macht die MA-Darstellung eindeutig und die Parameter aus den Daten identifizierbar.

Die Prognose aus einem ARMA-Modell folgt aus der Algebra des Lag-Operators. Die Einschritt-Prognose setzt künftige Innovationen auf ihren Erwartungswert null und schreibt die geschätzten Koeffizienten fort. Die Mehrschritt-Prognose iteriert diese Rekursion, und weil der Prozess stationär ist, konvergiert die Prognose mit wachsendem Horizont gegen den unbedingten Mittelwert $\mu$, während die Varianz des Prognosefehlers auf die unbedingte Varianz steigt.

Intuition

Warum das wichtig ist: Eine AR-Reihe schreibt ihre eigene Vergangenheit fort, sodass das heutige Niveau eine verblasste Kopie des gestrigen ist. Eine MA-Reihe trägt das Echo vergangener Überraschungen, einen Schock, der ein paar Perioden nachklingt und dann verschwunden ist. Beide hinterlassen unterschiedliche Fingerabdrücke, und die beiden Diagnoseplots lesen genau diese Abdrücke: Bei einem AR-Prozess fällt der Plot der partiellen Autokorrelation abrupt auf null, während die einfache Autokorrelation ausläuft, bei einem MA-Prozess ist es umgekehrt. Für die Prognose braucht es danach keinen neuen Gedanken. Wenn keine frische Überraschung mehr zu erwarten ist, ist die beste Schätzung für die ferne Zukunft schlicht der langfristige Durchschnitt, und das Modell sagt Ihnen, wie schnell Sie dort ankommen. Stellen Sie die Schieberegler unten auf einen reinen AR-Prozess und sehen Sie zu, wie der eine Plot auf null springt. Wechseln Sie zu einem reinen MA-Prozess, und der andere tut es.

Abbildung 23.2. ARMA(1,1)-Simulator mit Stichproben-ACF und -PACF. Stellen Sie einen reinen AR(1) ein ($\theta_1=0$): Die PACF hat eine Spitze bei Lag 1 und bricht dann ab, während die ACF geometrisch abklingt. Stellen Sie einen reinen MA(1) ein ($\phi_1=0$): Die ACF hat eine Spitze bei Lag 1 und bricht dann ab, während die PACF abklingt. Der Gegensatz von Abbruch und Abklingen ist die Identifikationsregel von Box und Jenkins. Ziehen Sie die Schieberegler.

All das setzte voraus, dass die Reihe stationär ist. Die Diagnoseplots, die Prognosen, die Bedeutung der Koeffizienten: Jedes davon ruht auf einem Mittelwert, zu dem zurückgekehrt wird, und auf einer endlichen Varianz. Was geschieht, wenn die Reihe weder das eine noch das andere hat?

23.3 Einheitswurzeln und Scheinregression

Schieben Sie den AR(1)-Koeffizienten ganz bis auf eins, und die Maschinerie aus §23.2 zerbricht. Die Reihe hat dann keinen Mittelwert mehr, zu dem sie zurückkehren könnte. Abbildung 23.1 hat das direkt gezeigt, als sie über $\phi = 1$ hinweg gezogen wurde. Unterhalb von eins ist die Reihe verankert und ein Schock klingt ab. Genau bei eins reißt die Verankerung und ein Schock wird dauerhaft. Dieser Grenzfall ist der Random Walk, der die Makroökonomie beherrscht, weil sich so viele Reihen so verhalten.

Random Walk. Ein AR(1) mit Koeffizient eins, $y_t = y_{t-1} + \varepsilon_t$, in dem der Wert jeder Periode der Wert der Vorperiode plus einem frischen Schock ist. Mit einer hinzugefügten Konstanten, $y_t = \delta + y_{t-1} + \varepsilon_t$, ist es ein Random Walk mit Drift, der nach oben (oder nach unten) trendet und dabei weiter wandert.
$$y_t = y_{t-1} + \varepsilon_t \qquad (\text{mit Drift:} \;\; y_t = \delta + y_{t-1} + \varepsilon_t)$$ (Gl. 23.6)
Einheitswurzel. Eine Wurzel des autoregressiven charakteristischen Polynoms $\phi(z)=0$, die auf dem Einheitskreis liegt ($z=1$). Die Einheitswurzel ist die Quelle dieser Art von Nichtstationarität: Sie lässt Schocks sich anhäufen statt abklingen, sodass die Varianz mit der Zeit unbeschränkt wächst.
Integrationsordnung, I(d). Die Anzahl der Differenzenbildungen, die eine Reihe braucht, um stationär zu werden. Eine stationäre Reihe ist I(0); ein Random Walk ist I(1), weil seine erste Differenz $\Delta y_t = y_t - y_{t-1} = \varepsilon_t$ weißes Rauschen ist; eine Reihe, die zwei Differenzen braucht, ist I(2).
Differenzstationär vs. trendstationär. Zwei Reihen können beide nach oben trenden und doch entgegengesetzte Behandlungen verlangen. Eine differenzstationäre Reihe (mit stochastischem Trend), etwa ein Random Walk mit Drift, wird durch Differenzenbildung stationär. Eine trendstationäre Reihe (mit deterministischem Trend), also stationäre Schwankungen um eine feste Gerade, wird stationär, indem man den geschätzten Trend abzieht. Verwechselt man die eine mit der anderen, folgt daraus entweder eine überflüssige Differenzenbildung oder eine zurückgelassene Einheitswurzel.

Eine Einheitswurzel lässt sich deshalb nicht ignorieren, weil sie der gewöhnlichen Regression eine Falle stellt. Nehmen Sie zwei völlig unabhängig voneinander erzeugte Random Walks. Regressieren Sie den einen auf den anderen, und weit häufiger, als der Zufall es erlauben sollte, finden Sie ein großes $R^2$ und eine $t$-Statistik, die jede übliche Schwelle mühelos überspringt. Die Regression verkündet eine starke Beziehung zwischen Reihen, die nichts gemeinsam haben. Das ist das Ergebnis zur Scheinregression, vor dem Granger und Newbold 1974 gewarnt haben, und es kippt den Instinkt, ein hohes $R^2$ zusammen mit einem signifikanten $t$ bedeute etwas Wirkliches.

Scheinregression. Das Phänomen (Granger und Newbold, 1974), dass die Regression zweier zusammenhangloser I(1)-Reihen aufeinander ein aufgeblähtes $R^2$ und eine scheinbar signifikante $t$-Statistik hervorbringt. Die übliche Inferenz ist ungültig, weil die Residuen der Regression selbst nichtstationär sind und die gewohnte Verteilungstheorie deshalb nicht greift.
Dickey-Fuller-Test / erweiterter Dickey-Fuller-Test (ADF). Ein Test der Nullhypothese, dass eine Reihe eine Einheitswurzel enthält. Die erweiterte Fassung regressiert die erste Differenz auf das verzögerte Niveau plus so viele verzögerte Differenzen, dass die Residuen weiß werden. Dass der Koeffizient des verzögerten Niveaus null ist, bildet die Nullhypothese der Einheitswurzel. Die Teststatistik folgt nicht der gewohnten $t$-Verteilung. Unter einer Einheitswurzel ist die maßgebliche Verteilung nicht standardisiert, sodass die besonderen kritischen Werte von Dickey und Fuller zu verwenden sind.
$$\Delta y_t = \alpha + \gamma\, y_{t-1} + \sum_{i=1}^{k}\delta_i\, \Delta y_{t-i} + \varepsilon_t, \qquad H_0:\gamma = 0 \;\; (\text{Einheitswurzel})$$ (Gl. 23.7)

Unter $H_0:\gamma=0$ fällt das Niveau $y_{t-1}$ heraus, und $\Delta y_t$ wird allein von den eigenen verzögerten Differenzen und vom Rauschen getrieben: eine Einheitswurzel. $H_0$ zugunsten von $\gamma<0$ zu verwerfen heißt, dass eine Abweichung vom Niveau teilweise zurückgeholt wird, die Reihe also stationär ist. Die verzögerten Differenzen sind die „Erweiterung“, die serielle Korrelation in $\varepsilon_t$ aufnimmt und den Test dadurch gültig macht. Die Anzahl der Lags $k$ wird über ein Informationskriterium gewählt.

Intuition

Warum das wichtig ist: Ein Random Walk kennt keine Schwerkraft. Eine stationäre Reihe ist verankert: Ziehen Sie sie von ihrem Zentrum weg, und sie schnellt zurück. Einen Random Walk zieht nichts nach Hause, ein Schock von heute wäscht sich also nie aus, und die Reihe wandert einfach dorthin, wohin die aufgehäuften Schocks sie tragen. Eine Einheitswurzel bedeutet Dauer statt Abklingen. Zwei unabhängig voneinander losgelassene Reihen driften über eine lange Stichprobe beide irgendwohin, und zwei driftende Dinge sehen immer korreliert aus, weil eine Gerade durch die Punktwolke stets in die eine oder die andere Richtung steigt. Ihre Regression meldet dann eine starke Beziehung zwischen zwei Reihen, die einander nie begegnet sind. Der Ausweg besteht darin, zuerst die richtige Frage zu stellen, nämlich ob diese Reihe verankert ist oder wandert, und zwar mit einem formalen Einheitswurzeltest, und wandernde Reihen in ihren Veränderungen (Differenzen) statt in ihren Niveaus zu untersuchen. Erzeugen Sie in der Abbildung unten die beiden unabhängigen Random Walks neu und sehen Sie zu, wie eine „signifikante“ Beziehung immer wieder aus dem Nichts auftaucht.

Abbildung 23.3. Demonstration der Scheinregression. Zwei mit unabhängigen Schocks erzeugte Random Walks, gemeinsam abgetragen. Die Anzeige meldet die OLS-Regression des einen auf den anderen mit ihrem $R^2$ und ihrer $t$-Statistik. Ziehen Sie die Schocks neu und sehen Sie zu, wie hohes $R^2$ und „signifikantes“ $t$ über unabhängige Ziehungen hinweg immer wiederkehren. Die Falle ist systematisch. Stellen Sie die Regression auf erste Differenzen um, und die scheinbare Beziehung bricht zusammen. Ziehen Sie die Schocks neu und schalten Sie zwischen Niveaus und Differenzen um.

Ein durchgehendes Beispiel. Ist das reale BIP der USA ein Random Walk? Schätzt man einen AR(1) für den Logarithmus, kommt der geschätzte Koeffizient sehr nahe an eins heraus. Führt man einen ADF-Test durch, verwirft er die Nullhypothese der Einheitswurzel in aller Regel nicht. Praktisch folgt daraus, dass sich die Produktion besser in Wachstumsraten (ersten Differenzen) als in Niveaus modellieren lässt.

Bisher ging es um jeweils eine Reihe. Makroökonomische Fragen betreffen meist mehrere Reihen zugleich, etwa Produktion und Inflation oder Zins und Wechselkurs, die sich gemeinsam bewegen. Die Werkzeuge für eine einzelne Gleichung verallgemeinern sich zu einem System.

23.4 Vektorautoregressionen (VAR)

Die meisten interessanten makroökonomischen Fragen betreffen mehr als eine Variable. Inflation und Leitzins reagieren aufeinander. Produktion, Preise und Geldmenge bewegen sich als System. Die Vektorautoregression ist die natürliche Verallgemeinerung des AR-Modells auf diesen Fall. Man stapelt die Variablen zu einem Vektor und lässt jede Variable von der jüngsten Vergangenheit aller Variablen abhängen.

Vektorautoregression (VAR). Ein System, in dem ein Variablenvektor $\mathbf{y}_t$ auf $p$ eigene Lags regressiert wird, mit einem Vektor von Innovationen $\mathbf{u}_t$. Jede Gleichung sieht aus wie eine gewöhnliche Regression einer Variablen auf die Lags aller Variablen, sodass sich das ganze System Gleichung für Gleichung mit OLS schätzen lässt.
$$\mathbf{y}_t = \mathbf{c} + A_1\mathbf{y}_{t-1} + A_2\mathbf{y}_{t-2} + \cdots + A_p\mathbf{y}_{t-p} + \mathbf{u}_t, \qquad E[\mathbf{u}_t\mathbf{u}_t'] = \Sigma$$ (Gl. 23.8)

Weil jede Gleichung dieselben Regressoren hat, nämlich die Lags aller Variablen, ist Gleichung für Gleichung angewandtes OLS effizient, und für die reduzierte Form braucht es keinen Systemschätzer. Die geschätzten Matrizen $A_i$ und die Residuenkovarianz $\Sigma$ fassen die Dynamik zusammen.

Granger-Kausalität. $X$ ist Granger-kausal für $Y$, wenn die Vergangenheit von $X$ die Prognose von $Y$ über das hinaus verbessert, was die Vergangenheit von $Y$ allein leistet. Die Bezeichnung ist ebenso sehr Warnung wie Definition: Granger-Kausalität ist prognostisch, nicht strukturell.
Impulsantwortfunktion (VAR). Der dynamische Pfad, den jede Variable des Systems nach einem einmaligen Schock beschreibt, aufgetragen gegen den Horizont. Im VAR in reduzierter Form beschreibt die Impulsantwort, wie sich das System nach einem Schock bewegt. Welche Variable aber den Schock erhält und wie er isoliert wird, ist genau die Frage, der sich §23.5 stellen muss.
Zerlegung der Prognosefehlervarianz. Der Anteil der Prognosefehlervarianz jeder Variablen in einem gegebenen Horizont, der jedem Schock im System zuzurechnen ist. Sie beantwortet die Frage „Wie viel der unvorhersagbaren Bewegung der Produktion stammt aus wirtschaftspolitischen Schocks und wie viel aus ihren eigenen?“ und hängt, wie die Impulsantwort, davon ab, wie die Schocks definiert sind.

Ein VAR in reduzierter Form liefert zwei Dinge fast umsonst. Es prognostiziert das ganze System gemeinsam, oft besser als ein strukturelles Modell, weil es wenige Restriktionen auferlegt, und es prüft Granger-Kausalität, sagt Ihnen also, welche Reihen für welche anderen prognostischen Gehalt tragen. Impulsantwort und Varianzzerlegung scheinen noch ein Drittes zu versprechen, nämlich eine Erklärung dessen, was nach einem Schock geschieht. Doch die Innovationen $\mathbf{u}_t$ sind über die Gleichungen hinweg korreliert. Eine Bewegung im Residuum der Inflation tritt meist zusammen mit einer Bewegung im Residuum des Zinses auf. Wenn das System also „auf einen Schock reagiert“, wessen Schock war es dann? Die reduzierte Form kann es nicht sagen.

Intuition

Warum das wichtig ist: Ein VAR lässt jede Variable von der jüngsten Vergangenheit aller anderen abhängen. Das bringt sofort zweierlei: eine gemeinsame Prognose des ganzen Systems und einen Test darüber, welche Reihe hilft, welche andere vorherzusagen, die sogenannte „Granger-Kausalität“. Vorhersagen helfen ist nicht verursachen. Dass Eisverkäufe helfen, Ertrinkungsfälle vorherzusagen, heißt nicht, dass Eis das Ertrinken verursacht. Der Sommer treibt beides. Ein VAR kann Ihnen sagen, dass der Leitzins hilft, die Inflation vorherzusagen, und das ist nützlich, aber es geht nicht so weit zu sagen, der Zins habe die Inflation bewegt. Es geht deshalb nicht so weit, weil die Überraschungen in den verschiedenen Gleichungen gemeinsam eintreffen, ineinander verheddert, sodass „die Reaktion des Systems auf einen Schock“ mehrdeutig bleibt, solange Sie nicht entwirrt haben, welchen Schock Sie meinen. Dieses Entwirren ist das Identifikationsproblem.

Abbildung 23.4. VAR-/SVAR-Impulsantwort-Simulator. Ein System aus zwei Variablen: Inflation und Leitzins, das geldpolitische VAR des durchgehenden Beispiels. Wählen Sie, welche Variable den Schock erhält, und das Panel zeigt die dynamischen Reaktionen beider Variablen. Die Anzeigen zu Granger-Kausalität und Varianzzerlegung aktualisieren sich mit. In §23.5 kehren Sie die Cholesky-Anordnung um und sehen zu, wie sich die Impulsantworten ändern: dieselben Daten, eine andere Identifikationsannahme, eine andere ökonomische Erklärung. Wählen Sie den Schock, setzen Sie den Horizont und schalten Sie die Anordnung um.

Die Impulsantworten, die Sie eben erzeugt haben, setzten eine Anordnung voraus, also eine Entscheidung darüber, welche Variable die andere innerhalb der Periode bewegen kann. Diese Entscheidung fiel stillschweigend. Sie offenzulegen ist der strukturelle Schritt, und an dieser Stelle können die Daten den Streit nicht mehr entscheiden.

23.5 Strukturelles VAR und Identifikation

Die Residuen der reduzierten Form sind korreliert, und diese Korrelation ist das Problem. Ein ökonomischer Schock, etwa eine überraschende geldpolitische Straffung, sollte eine einzelne Störung mit klarer Deutung sein. Die Innovation der Zinsgleichung in reduzierter Form ist aber davon verunreinigt, was im selben Augenblick die Inflation bewegt hat, und umgekehrt. Die Residuen sind Gemische. Um die zugrunde liegenden strukturellen Schocks zurückzugewinnen, müssen Sie angeben, wie die Gemische entstanden sind, und diese Information steckt nicht in den Daten.

Strukturelles VAR (SVAR). Ein VAR, dessen Innovationen als orthogonale (unkorrelierte, auf Varianz eins normierte) strukturelle Schocks ökonomisch gedeutet werden. Die Residuen der reduzierten Form werden als Linearkombination dieser strukturellen Schocks geschrieben, $\mathbf{u}_t = B\boldsymbol{\varepsilon}_t$, und die Aufgabe der Identifikation besteht darin, die Matrix $B$ festzulegen.
$$\mathbf{u}_t = B\,\boldsymbol{\varepsilon}_t, \qquad E[\boldsymbol{\varepsilon}_t\boldsymbol{\varepsilon}_t'] = I, \qquad \text{also} \;\; \Sigma = BB'$$ (Gl. 23.9)
Identifikation (im Sinn der Zeitreihenanalyse). Das Problem, die strukturellen Schocks $\boldsymbol{\varepsilon}_t$ aus den Residuen der reduzierten Form $\mathbf{u}_t$ zurückzugewinnen, gleichbedeutend damit, $B$ festzulegen. Die geschätzte Kovarianz $\Sigma$ liefert nur $n(n+1)/2$ Gleichungen, $B$ hat aber $n^2$ Unbekannte. Die Lücke von $n(n-1)/2$ Restriktionen muss von außerhalb der Daten kommen.

Für ein System aus zwei Variablen liefert $\Sigma=BB'$ drei Gleichungen (zwei Varianzen und eine Kovarianz) für die vier Elemente von $B$. Eine Restriktion fehlt. Sie aufzuerlegen ist die Identifikation, und die Impulsantworten im Horizont $h$, gegeben durch $\Theta_h = A^h B$, erben, was auferlegt wurde.

$$\Theta_h = A^h B \qquad (\text{strukturelle Impulsantwort im Horizont } h)$$ (Gl. 23.10)
Cholesky-Identifikation / rekursive Identifikation. Ein Schema, das $B$ untere Dreiecksgestalt gibt, gleichbedeutend damit, die Variablen anzuordnen und anzunehmen, dass die erste nicht zeitgleich auf Schocks der nach ihr angeordneten Variablen reagieren kann. Setzt man im geldpolitischen VAR den Zins an die erste Stelle, so heißt das: Die Inflation reagiert innerhalb der Periode auf einen Zinsschock, der Zins aber nicht innerhalb der Periode auf einen Inflationsschock.
Vorzeichenrestriktionen. Ein Identifikationsansatz, der statt Nullrestriktionen verlangt, dass die strukturellen Impulsantworten über einen gewissen Horizont bestimmte Vorzeichen haben (ein geldpolitischer Straffungsschock hebt etwa den Zins und senkt die Inflation). Er liefert eine Menge von Impulsantworten, die mit den auferlegten Vorzeichen vereinbar sind, und damit eine ehrliche Darstellung dessen, was die Annahmen offenlassen.

Drei Familien von Restriktionen sind gebräuchlich, und Abbildung 23.4 lässt Sie die Folgen der häufigsten davon spüren. Schalten Sie die Cholesky-Anordnung zwischen „Zins zuerst“ und „Inflation zuerst“ um, und die Impulsantworten ändern sichtbar ihre Gestalt. An den Daten hat sich nichts geändert: dieselben Residuen, dieselbe geschätzte Dynamik. Geändert hat sich die Annahme darüber, welche Variable die andere innerhalb der Periode bewegen kann, und diese Annahme hat die ökonomische Erklärung neu geschrieben. Kurzfristige Nullrestriktionen wie die von Cholesky sind eine Möglichkeit; Langfristrestriktionen, die annehmen, dass ein Schock auf eine bestimmte Variable keine dauerhafte Wirkung hat (die Zerlegung von Nachfrage- und Angebotsschocks bei Blanchard und Quah ist das Standardbeispiel), sind eine zweite; Vorzeichenrestriktionen sind die moderne Alternative, die statt einer einzelnen Linie die identifizierte Menge ausweist.

Hier ist sich das Fach uneins. Christopher Sims, der das VAR in die Makroökonomie eingeführt hat, hielt die aufwendigen Identifikationsrestriktionen der älteren Strukturmodelle für „unglaubwürdig“, für Annahmen, die der Handhabbarkeit zuliebe auferlegt wurden, und meinte, eine ehrliche empirische Makroökonomie solle die Daten mit so wenigen Restriktionen wie möglich sprechen lassen. Die Gegenposition hält die Impulsantworten ohne ökonomische Struktur für nicht deutbar. Richtig sei deshalb, Restriktionen aufzuerlegen, die man verteidigen kann, und sie offen zu benennen. Beide Positionen sind ernst zu nehmen, und der Schalter, den Sie eben benutzt haben, ist genau der Streitpunkt: wie viel Identifikation die Daten ehrlich tragen können und was mit der Lücke zu tun ist.

Das Urteil darüber, ob die atheoretischen VARs oder die strukturelle Identifikation das Feld gewonnen haben, wird in der Großen Frage zur Glaubwürdigkeit der ökonometrischen Methodologie ausführlich verhandelt.

Intuition

Warum das wichtig ist: Die Daten reichen Ihnen korrelierte Schocks. Die Überraschungen in der Inflation und im Zins treffen ineinander verheddert ein. Um eine ökonomische Erklärung zu geben, müssen Sie sagen, wohin der Pfeil innerhalb derselben Periode zeigt: Hat sich der Zins zuerst bewegt und die Inflation folgte, oder umgekehrt? Diesen Pfeil muss die Ökonomie liefern, denn die Daten können es nicht. Sie zeigen nur, dass sich beide gemeinsam bewegen. Diese eine hinzugefügte Annahme, die die Daten nicht prüfen können, entscheidet die ganze Antwort. In der Abbildung blieb beim Umkehren der Anordnung jede Zahl im Datensatz gleich, und die Erklärung kippte trotzdem. Sims sorgte sich, dass Ökonomen Pfeile einschmuggelten, die sie nicht rechtfertigen konnten, und das Ergebnis dann einen Befund nannten. Das andere Lager sagt, um die Wahl eines Pfeils komme man nicht herum, also wähle man einen, den man verteidigen kann, und sage es laut.

Die bisherigen Systeme waren aus stationären oder differenzierten Reihen gebaut. Doch die Differenzenbildung kann etwas Wirkliches wegwerfen. Wenn zwei Reihen gemeinsam wandern, löscht das Bilden ihrer Veränderungen genau die Beziehung, die sie verbindet. Die Kointegration holt sie zurück.

23.6 Kointegration und Fehlerkorrektur

Erinnern Sie sich an die Warnung vor der Scheinregression. Zwei unabhängig wandernde Reihen sehen verbunden aus, und der Ausweg besteht darin, sie in Differenzen zu untersuchen. Manchmal aber sind zwei wandernde Reihen tatsächlich aneinander gebunden, nicht jede an einen festen Mittelwert, sondern aneinander. Kurz- und langfristige Zinssätze driften über die Jahrzehnte, und doch bleibt der Abstand zwischen ihnen in einem Band. Konsum und Einkommen trenden beide nach oben, und doch ist ihr Verhältnis stabil. Ein solches Paar zu differenzieren hieße, genau die langfristige Beziehung wegzudifferenzieren, auf die es ankommt. Die Kointegration ist das Instrumentarium für den Fall, dass die Beziehung wirklich besteht.

Kointegration. Zwei (oder mehr) I(1)-Reihen sind kointegriert, wenn eine Linearkombination von ihnen I(0), also stationär, ist. Die Reihen wandern einzeln, die Kombination aber nicht: Sie wird zurückgeholt, sobald sie ausschert. Diese stationäre Kombination ist eine langfristige Gleichgewichtsbeziehung.
$$z_t = y_t - \beta x_t \sim I(0) \quad \text{wenn} \quad y_t, x_t \sim I(1)$$ (Gl. 23.11)
Zweistufiges Verfahren von Engle und Granger. Ein Verfahren in zwei Schritten. (1) Es schätzt die Kointegrationsbeziehung mit OLS auf den Niveaus, was gerade deshalb zulässig ist, weil die Variablen kointegriert sind und die Regression somit keine Scheinregression ist, und prüft, ob das Residuum $z_t$ stationär ist. (2) Es verwendet das geschätzte Residuum als Fehlerkorrekturterm in einem dynamischen Modell.
Fehlerkorrekturmodell (ECM). Ein Modell der kurzfristigen Dynamik, das einen Term enthält, der proportional zur Abweichung der Vorperiode von der langfristigen Beziehung ist. Der Koeffizient dieses Terms, $\lambda$, ist negativ und misst, wie schnell das System korrigiert: Pro Periode wird ein Anteil $|\lambda|$ der Lücke geschlossen.
$$\Delta y_t = \lambda\,(y_{t-1}-\beta x_{t-1}) + \text{(verzögerte } \Delta y, \Delta x \text{ Terme)} + \varepsilon_t, \qquad \lambda<0$$ (Gl. 23.12)
Granger-Darstellungssatz. Das Ergebnis, dass Kointegration und eine Fehlerkorrekturdarstellung gleichwertig sind: Ist eine Menge von I(1)-Reihen kointegriert, dann existiert ein ECM, das ihre Dynamik beschreibt, und umgekehrt. Diese Gleichwertigkeit erlaubt es, das ECM zu schätzen, sobald die Kointegration festgestellt ist.
Johansen-Verfahren. Eine Maximum-Likelihood-Methode für das ganze System, die prüft, wie viele unabhängige Kointegrationsbeziehungen zwischen einer Menge von Reihen bestehen, also den Kointegrationsrang. Wo Engle und Granger eine einzelne Beziehung in zwei Schritten behandeln, behandelt Johansen mehrere gemeinsam. Der Rang ist die Anzahl der langfristigen Gleichgewichtsbeziehungen, die das System zusammenhalten.

Ist der Fehlerkorrekturkoeffizient $\lambda = -0.2$, dann wird pro Periode ein Fünftel jeder Abweichung von der langfristigen Beziehung rückgängig gemacht. Ist $\lambda = 0$, gibt es keinen Zug und die Reihen sind überhaupt nicht kointegriert, die „Beziehung“ ist also die scheinbare aus §23.3. Der Granger-Darstellungssatz macht das Vorgehen erst seriös. Er garantiert, dass immer dann, wenn echte Kointegration vorliegt, ein Fehlerkorrekturmodell die richtige Art ist, die Dynamik aufzuschreiben. Das ECM zu schätzen ist damit die implizierte Form. Das Johansen-Verfahren dehnt den Gedanken auf Systeme mit möglicherweise mehreren Gleichgewichtsbeziehungen aus, und sein Rangtest beantwortet, wie viele es sind.

Intuition

Warum das wichtig ist: Zwei Reihen können jede für sich ewig wandern und driften doch nie voneinander weg, wie zwei Hunde, die frei laufen, aber durch eine kurze Leine miteinander verbunden sind. Jeder geht, wohin er will. Die Leine hält den Abstand zwischen ihnen davon ab, größer zu werden. Dieser gemeinsame Abstand ist das langfristige Gleichgewicht, und die Geschwindigkeit, mit der die Leine sie zurückreißt, wenn sie ausscheren, ist der Fehlerkorrekturterm. Das ist das genaue Gegenteil der Falle der Scheinregression: Dort sahen zwei wandernde Reihen nur so aus, als hingen sie zusammen, hier tun sie es wirklich, und die Leine ist echt. Vorhin lautete der Rat, wandernde Reihen vor der Untersuchung zu differenzieren, aber wer ein aneinander gebundenes Paar differenziert, zerschneidet die Leine und verliert genau das, was er sehen wollte. Die Abbildung unten hat die Straffheit der Leine als Schieberegler. Bei null treiben die beiden Reihen frei auseinander, und je straffer Sie sie ziehen, desto mehr legt sich der Abstand zwischen ihnen in ein stationäres Band, während jede Reihe weiter wandert.

Abbildung 23.5. Kointegrations- und Fehlerkorrektur-Simulator. Zwei I(1)-Reihen teilen sich einen gemeinsamen stochastischen Trend. Das zweite Panel zeigt ihren Abstand $z_t = y_t - x_t$. Bei $|\lambda|=0$ wandert der Abstand selbst (keine Kointegration). Steigt $|\lambda|$, werden Abweichungen von der langfristigen Beziehung schneller zurückgeholt und der Abstand wird sichtbar stationär. Ziehen Sie die Fehlerkorrekturgeschwindigkeit und ziehen Sie die Schocks neu.

Bisher wurde durchweg der bedingte Mittelwert modelliert, also wohin die Reihe voraussichtlich geht. Die bedingte Varianz, also wie unsicher diese Erwartung ist, hat einen vorhersagbaren Rhythmus ganz eigener Art.

23.7 ARCH- und GARCH-Modellierung der Volatilität

Betrachten Sie eine lange Reihe täglicher Aktienrenditen, und ein Merkmal springt heraus, bevor irgendein Modell geschätzt ist. Die Turbulenz kommt in Klumpen. Ruhige Strecken mit kleinen Bewegungen werden von stürmischen Strecken mit großen Bewegungen unterbrochen, und die Stürme halten tage- oder wochenlang an, bevor die Ruhe zurückkehrt. Die Renditen selbst sind nahezu unvorhersagbar, wie man es von effizienten Märkten erwarten würde, ihre Größe aber nicht. Große Bewegungen treten in Clustern mit anderen großen Bewegungen auf. Das ist das stilisierte Faktum, das die Modelle der bedingten Varianz abbilden sollen.

Volatilitätscluster. Die empirische Regelmäßigkeit, dass auf große Änderungen große und auf kleine kleine folgen, unabhängig vom Vorzeichen. Perioden hoher und niedriger Volatilität bilden zeitliche Cluster: das bestimmende stilisierte Faktum von Reihen finanzieller Renditen.
Bedingte vs. unbedingte Varianz. Die bedingte Varianz ist die Varianz der nächsten Beobachtung bei gegebener Vergangenheit. Sie ändert sich über die Zeit, wenn neue Information eintrifft. Die unbedingte Varianz ist der langfristige Durchschnitt, eine einzige Zahl. ARCH- und GARCH-Modelle machen die bedingte Varianz zeitveränderlich und prognostizierbar und lassen die unbedingte Varianz konstant.
ARCH-Prozess. Ein Modell (Engle, 1982), in dem die bedingte Varianz eine Funktion vergangener quadrierter Schocks ist: Ein großer Schock hebt die Varianz der nächsten Periode, weshalb die Turbulenz anhält. Ein ARCH(q) verwendet $q$ verzögerte quadrierte Schocks.
$$\sigma_t^2 = \omega + \sum_{i=1}^{q}\alpha_i\,\varepsilon_{t-i}^2$$ (Gl. 23.13)
GARCH-Prozess. Die Verallgemeinerung (Bollerslev, 1986), in der die bedingte Varianz auch von ihren eigenen vergangenen Werten abhängt, also von vergangenen bedingten Varianzen ebenso wie von vergangenen quadrierten Schocks. Ein GARCH(1,1) bildet die Persistenz der Volatilität mit nur zwei Parametern ab und ist das Standardmodell der empirischen Finanzmarktforschung.
$$\sigma_t^2 = \omega + \sum_{i=1}^{q}\alpha_i\,\varepsilon_{t-i}^2 + \sum_{j=1}^{p}\beta_j\,\sigma_{t-j}^2$$ (Gl. 23.14)

Im GARCH(1,1) misst die Summe $\alpha_1+\beta_1$ die Persistenz und bestimmt, wie langsam ein Volatilitätsschock abklingt, und die unbedingte Varianz $\omega/(1-\alpha_1-\beta_1)$ ist nur dann endlich, wenn $\alpha_1+\beta_1<1$. Auf täglichen Aktiendaten liegt die geschätzte Persistenz regelmäßig über $0.95$, was bedeutet, dass eine Volatilitätsspitze lange braucht, um abzuklingen: der empirische Gehalt des Satzes „Turbulenz hält an“.

Intuition

Warum das wichtig ist: Ob die Rendite von morgen steigt oder fällt, können Sie nicht vorhersagen. Könnten Sie es, wäre das Geschäft längst gemacht und die Vorhersage ausgelöscht. Sie können aber vorhersagen, ob morgen ein ruhiger oder ein wilder Tag wird, denn Ruhe und Wildheit bilden Cluster. Ein stürmischer Markt heute sagt Ihnen, dass die nächsten Tage wahrscheinlich ebenfalls stürmisch werden, auch wenn er Ihnen nichts darüber sagt, in welche Richtung die Bewegungen gehen. Diese Trennung, bei der die Richtung einer Überraschung unprognostizierbar, die Größe einer Überraschung aber prognostizierbar ist, steckt hinter diesen Modellen. Eine große Bewegung heute speist einen größeren erwarteten Ausschlag morgen. Ist diese Rückkopplung stark, dauern Stürme lange, sobald sie einmal begonnen haben. Die Abbildung unten hat einen Persistenzregler. Drehen Sie ihn nach unten, und die Reihe sieht aus wie strukturloses Rauschen. Drehen Sie ihn Richtung eins, und Sie sehen zu, wie sich ruhige und turbulente Strecken zu langen Wellen ordnen, ganz ohne Gleichung.

Abbildung 23.6. GARCH(1,1)-Simulator für Volatilitätscluster. Das obere Panel zeigt eine simulierte Renditereihe, das untere die bedingte Volatilität $\sigma_t$, die sie erzeugt hat. Bei geringer Persistenz (kleines $\alpha_1+\beta_1$) sehen die Renditen aus wie homoskedastisches Rauschen. Nähert sich die Persistenz eins, bilden ruhige und turbulente Perioden scharfe Cluster und die Volatilitätslinie zeigt lange Wellen. Der Schieberegler sichert $\alpha_1+\beta_1<1$. Ziehen Sie die Persistenzregler und ziehen Sie die Schocks neu.

Das einfache GARCH(1,1) hat eine Familie von Verfeinerungen hervorgebracht, deren jede eine Beschränkung behebt. EGARCH modelliert den Logarithmus der Varianz und lässt damit den „Leverage-Effekt“ zu, bei dem schlechte Nachrichten die Volatilität stärker heben als gleich große gute, und es braucht keine Nichtnegativitätsbedingungen. GARCH-in-Mean (GARCH-M) lässt die bedingte Varianz direkt in die Renditegleichung eingehen und formalisiert damit den Gedanken, dass Anleger bei hohem Risiko höhere erwartete Renditen verlangen. Diese Erweiterungen sind das Arbeitsvokabular der empirischen Finanzmarktforschung. Dahinter liegt eine Front aus bayesianischer Schätzung hochdimensionaler VARs und maschinellen Lernverfahren für Volatilität und Vorhersage.

Zusammenfassung

  1. Stationarität ist die Voraussetzung. Die übliche Zeitreiheninferenz unterstellt Kovarianzstationarität, also einen stabilen Mittelwert, eine stabile Varianz und eine stabile Autokovarianzstruktur. Der Lag-Operator verdichtet die Algebra, und die Wold-Zerlegung garantiert, dass jede stationäre Reihe ein gefilterter Strom von Schocks weißen Rauschens ist.
  2. ARMA-Modelle und ihre Fingerabdrücke. AR-Prozesse schreiben ihre eigene Vergangenheit fort, MA-Prozesse tragen das Echo vergangener Schocks. ACF und PACF diagnostizieren sie (die PACF bricht bei AR ab, die ACF bei MA), und das ist die Identifikationsregel von Box und Jenkins.
  3. Einheitswurzeln setzen das Instrumentarium außer Kraft. Ein Random Walk hat eine Einheitswurzel, Schocks sind deshalb dauerhaft, die Reihe wandert, und zwei unabhängig wandernde Reihen regressieren zum Schein aufeinander (hohes $R^2$, signifikantes $t$). Der ADF-Test spürt eine Einheitswurzel auf, und die Differenzenbildung ist der Ausweg.
  4. Das VAR modelliert das System. Jede Variable auf den Lags aller Variablen ergibt eine gemeinsame Prognose und einen Test auf Granger-Kausalität, doch Granger-Kausalität ist nur prognostisch, und die Schocks der reduzierten Form sind korreliert.
  5. Die strukturelle Identifikation ist der Angelpunkt. Ökonomische Schocks aus den Residuen der reduzierten Form zurückzugewinnen verlangt eine Restriktion, die die Daten nicht liefern können. Cholesky-Anordnung, Langfristrestriktionen und Vorzeichenrestriktionen erlegen jeweils eine auf, und die Wahl verändert die Impulsantwort. Das ist die Spannung um Sims' atheoretisches VAR.
  6. Die Kointegration holt langfristige Beziehungen zurück. I(1)-Reihen können eine stationäre Linearkombination teilen, ein langfristiges Gleichgewicht, das die Differenzenbildung zerstören würde. Das Fehlerkorrekturmodell misst die Geschwindigkeit der Rückkehr, und der Granger-Darstellungssatz macht Kointegration und ECM gleichwertig.
  7. ARCH und GARCH modellieren zeitveränderliches Risiko. Volatilität bildet Cluster, deshalb ist die Größe der Renditen prognostizierbar, auch wenn ihre Richtung es nicht ist. GARCH lässt die bedingte Varianz von vergangenen quadrierten Schocks und von vergangener Varianz abhängen. Eine Persistenz nahe eins bedeutet, dass Turbulenz anhält.

Wichtige Gleichungen

BezeichnungGleichungBeschreibung
Gl. 23.1$E[y_t]=\mu,\ \mathrm{Var}(y_t)=\sigma^2,\ \mathrm{Cov}(y_t,y_{t-k})=\gamma_k$Bedingungen der Kovarianzstationarität
Gl. 23.2$y_t = \sum_{j=0}^{\infty}\psi_j\varepsilon_{t-j} + \kappa_t$Wold-Zerlegung
Gl. 23.3$y_t = c + \phi_1 y_{t-1} + \cdots + \phi_p y_{t-p} + \varepsilon_t$AR(p)
Gl. 23.4$y_t = \mu + \varepsilon_t + \theta_1\varepsilon_{t-1} + \cdots + \theta_q\varepsilon_{t-q}$MA(q)
Gl. 23.5$\phi(L)y_t = \theta(L)\varepsilon_t$ARMA in Lag-Operator-Form
Gl. 23.6$y_t = y_{t-1} + \varepsilon_t$ (mit Drift $+\,\delta$)Random Walk
Gl. 23.7$\Delta y_t = \alpha + \gamma y_{t-1} + \sum\delta_i\Delta y_{t-i} + \varepsilon_t;\ H_0:\gamma=0$Erweiterte Dickey-Fuller-Regression
Gl. 23.8$\mathbf{y}_t = \mathbf{c} + A_1\mathbf{y}_{t-1} + \cdots + A_p\mathbf{y}_{t-p} + \mathbf{u}_t$VAR(p) in reduzierter Form
Gl. 23.9$\mathbf{u}_t = B\boldsymbol{\varepsilon}_t,\ E[\boldsymbol{\varepsilon}_t\boldsymbol{\varepsilon}_t']=I$SVAR-Abbildung (Residuen auf strukturelle Schocks)
Gl. 23.10$\Theta_h = A^h B$Strukturelle Impulsantwort im Horizont $h$
Gl. 23.11$z_t = y_t - \beta x_t \sim I(0)$, wenn $y_t, x_t \sim I(1)$Kointegrationsbeziehung
Gl. 23.12$\Delta y_t = \lambda(y_{t-1}-\beta x_{t-1}) + \cdots + \varepsilon_t$Fehlerkorrekturmodell
Gl. 23.13$\sigma_t^2 = \omega + \sum_{i=1}^{q}\alpha_i\varepsilon_{t-i}^2$Bedingte Varianz im ARCH(q)
Gl. 23.14$\sigma_t^2 = \omega + \sum\alpha_i\varepsilon_{t-i}^2 + \sum\beta_j\sigma_{t-j}^2$Bedingte Varianz im GARCH(p,q)

Übung

  1. Ihnen wird eine ACF gezeigt, die geometrisch abklingt, und eine PACF, die bei Lag 1 eine Spitze hat und danach im Wesentlichen null ist. Bestimmen Sie den Prozess (AR, MA oder ARMA) und seine Ordnung und nennen Sie die Regel von Box und Jenkins, die Sie verwendet haben.
  2. Klassifizieren Sie jede Reihe anhand ihrer Diagnoseplots als AR oder MA: (a) Die ACF bricht nach Lag 2 ab, die PACF klingt ab; (b) die PACF bricht nach Lag 1 ab, die ACF klingt ab. Geben Sie jeweils die Ordnung an.
  3. Eine Reihe muss zweimal differenziert werden, bevor ihre ACF und ein ADF-Test Stationarität anzeigen. Ihre erste Differenz zeigt weiterhin eine ACF nahe einer Einheitswurzel. Geben Sie ihre Integrationsordnung an und begründen Sie die Einordnung.
  4. Eine Tabelle zur Granger-Kausalität berichtet, dass der Leitzins Granger-kausal für die Inflation ist ($p<0.01$), die Inflation aber nicht Granger-kausal für den Zins ($p=0.4$). Deuten Sie das Ergebnis und sagen Sie genau, was es über Kausalität belegt und was nicht.

Anwendung

  1. Schreiben Sie die ADF-Regression auf, mit der geprüft wird, ob das logarithmierte reale BIP eine Einheitswurzel hat. Geben Sie die Nullhypothese über den Koeffizienten des verzögerten Niveaus an, erklären Sie, warum die gewöhnlichen kritischen $t$-Werte unangemessen sind, und erklären Sie das Risiko der Scheinregression, das dazu Anlass gibt, den Test durchzuführen, bevor man das BIP auf eine andere trendbehaftete Reihe regressiert.
  2. Spezifizieren Sie ein VAR(1) aus zwei Variablen in Produktionswachstum und Inflation. Schreiben Sie die beiden Gleichungen ausdrücklich auf, geben Sie an, welche Regressoren jede enthält, und erklären Sie in Worten, was eine Impulsantwort der Inflation auf einen Schock im Produktionswachstum zeigen würde.
  3. Erklären Sie am geldpolitischen VAR aus Abbildung 23.4, warum die Cholesky-Anordnung eine Rolle spielt. Geben Sie die zeitgleiche Restriktion an, die jede Anordnung auferlegt, und erklären Sie, warum ein Umordnen der Variablen die Impulsantworten verändert, obwohl die geschätzte reduzierte Form unverändert bleibt.
  4. Ein an tägliche Renditen angepasstes GARCH(1,1) liefert $\alpha_1 = 0.08$, $\beta_1 = 0.90$. Berechnen Sie die Persistenz, geben Sie an, ob die unbedingte Varianz endlich ist, und deuten Sie, was der Persistenzwert darüber besagt, wie lange eine Volatilitätsspitze anhält.

Herausforderung

  1. Leiten Sie mit dem Lag-Operator die MA($\infty$)-Darstellung eines stationären AR(1) her, $y_t = \phi y_{t-1} + \varepsilon_t$ mit $|\phi|<1$. Zeigen Sie, dass die Moving-Average-Gewichte $\psi_j = \phi^j$ sind, und erklären Sie, was die Bedingung $|\phi|<1$ über die Gewichte und die Varianz sicherstellt.
  2. Skizzieren Sie die Argumentation von Granger und Newbold dafür, warum OLS über zwei unabhängige Random Walks eine scheinbare Signifikanz erzeugt. Erklären Sie, warum die Residuen der Regression unter der Nullhypothese ohne Beziehung nichtstationär sind und warum das die gewohnte $t$-Verteilung für den Steigungskoeffizienten ungültig macht.
  3. Formulieren Sie die Gleichwertigkeit des Granger-Darstellungssatzes (Kointegration genau dann, wenn eine ECM-Darstellung existiert) in Worten. Schreiben Sie dann für einen Kointegrationsvektor, der die langfristige Beziehung $y = 1.5\,x$ impliziert, das Fehlerkorrekturmodell für $\Delta y_t$ auf, bezeichnen Sie den Fehlerkorrekturterm und deuten Sie Vorzeichen und Größe seines Koeffizienten.
  4. Stellen Sie Sims' Skepsis gegenüber der strukturellen Identifikation einem vertretbaren Schema struktureller Identifikation gegenüber (etwa einer Cholesky-Anordnung, die mit der Annahme begründet wird, dass die Geldpolitik erst mit Verzögerung auf die Inflation reagiert). Geben Sie an, was jeder Zugang einbringt und was er aufgibt, und erklären Sie, warum die Frage, welcher der richtige ist, offen bleibt.

Quellen

Genannte Literatur: Wold (1938); Box & Jenkins (1970); Granger & Newbold (1974); Dickey & Fuller (1979); Engle (1982); Sims (1980); Bollerslev (1986); Engle & Granger (1987); Blanchard & Quah (1989); Johansen (1991); Hamilton (1994); Stock & Watson (2001); Enders (2014).