ich sitze seit länger Zeit an einer Ausarbeitung zum Thema der Rezessions Vorhersage und muss dazu das Thema auch empirisch bearbeiten. Nun habe ich bereits verschiedene R-Kurse belegt und bin eigentlich mit den Grundlagen vertraut, allerdings komme ich in diesem Fall einfach nicht weiter, vllt liegt es auch daran dass ich mich schon zu viel mit dem Thema beschäftigt habe und nun irgendwie den Faden verloren habe. Deshalb dachte ich einmal ich versuche es auf diese weise, evtl. kann mir ja jemand weiter helfen.
Ich möchte mit Hilfe eines Probit-Modells über die Zinsstruktur(langfristige Zinsen - kurzfristige Zinsen) Rezessionen vorhersagen. Das wurde auch schon viel in Fachzeitschriften etc. gemacht, leider finde ich nichts was mir bei meinem Problem weiter hilft. Im optimal Fall soll bei einer negativen Zinsstrukutr eine Hohe Wahrscheinlichkeit für eine Rezession herauskommen.
Ich habe also einen Datensatz mit einer binären Variable (1=Rezession, 0=keine Rezession) und den dazugehörigen Wert der oben genannten Differenz.
Mein erster Versuch war einmal, ein ganz "normal" Probit Regression durchzuführen, was auch soweit funktioniert habe. (nber steht für die Zeitreihe der binären Variable und ts_ts für den Term Spread)
Code: Alles auswählen
glm(formula=ts_nber~ts_ts, family = binomial(link="probit"))
summary(prob1)
predict(prob1, ts_nber, type="response")
Call:
glm(formula = ts_nber ~ ts_ts1, family = binomial(link = "probit"))
Deviance Residuals:
Min 1Q Median 3Q Max
-1.3024 -0.5671 -0.4450 -0.3556 2.4854
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -1.00803 0.06734 -14.969 < 2e-16 ***
ts_ts1 -0.20465 0.03879 -5.276 1.32e-07 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 496.84 on 650 degrees of freedom
Residual deviance: 468.61 on 649 degrees of freedom
AIC: 472.61
Number of Fisher Scoring iterations: 5
summary(pred1)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.04105 0.07424 0.10800 0.12750 0.16168 0.57180
Wie lässt sich das nun genau interpretieren was R macht? Mein Datensatz sind monatle Daten von 1965-01 bis 2019-03. Bedeutet dieses Ergebniss einfach nur, dass in dem Monat wo die Differenz besonders niedrig ist, eine Hohe Wahrscheinlichkeit auf eine Rezession besteht?
Wie könnte ich nun in R fortfahren und verschiedene Lags einbauen, damit ich in etwa herausbekomme bspw. für eine 10% Rezessionswahrscheinlichkeit muss die Zinsstruktur -0.1 erreichen und dies zeitlich versetzt untersuchen, also bspw. ein Quartal oder vier Monate im vorraus?
Schonmal Vielen Dank ich hoffe ich hab den Rahmen dieses Forums nicht ein wenig gesprengt bei der Menge die ich geschrieben habe.
Würde mich wirklich sehr über jede Bemühung, mir bei meinem Problem weiterzuhelfen, freuen.
Edit: bigben hat code-Tags stat i-Tags gesetzt.