Rupture structurelle de la crise : III
Vous pouvez maintenant tout rassembler pour effectuer le test de Chow.
Les données de 2005 à 2010 ont été séparées en deux DataFrames, before et after, en utilisant le 30 juin 2008 comme point de rupture structurelle (identifié dans le premier exercice de cette série). Les colonnes des deux DataFrames sont mort_del et returns pour, respectivement, les données de défaut hypothécaire et les rendements.
Vous exécuterez deux régressions OLS sur before et after, en régressant la colonne returns sur la colonne mort_del dans chaque DataFrame, puis vous en tirerez la somme des carrés des résidus.
Ensuite, vous calculerez la statistique du test de Chow comme dans la vidéo, en utilisant ssr_total (fourni par le deuxième exercice) et les résidus obtenus. La valeur critique de F à 99 % de confiance est d'environ 5,85. Quelle valeur obtenez-vous pour votre statistique de test?
Cette activité fait partie du cours
Gestion quantitative des risques en Python
Instructions de l’exercice
- Ajoutez un terme d'ordonnée à l'origine (interception) OLS à
mort_delpourbeforeetafter. - Ajustez une régression OLS de la colonne
returnssur la colonnemort_del, pourbeforeetafter. - Placez la somme des carrés des résidus dans
ssr_beforeetssr_after, pourbeforeetafter, respectivement. - Créez et affichez la statistique du test de Chow.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# Add intercept constants to each sub-period 'before' and 'after'
before_with_intercept = sm.____(before['mort_del'])
after_with_intercept = sm.____(____['mort_del'])
# Fit OLS regressions to each sub-period
r_b = sm.____(____['returns'], before_with_intercept).____
r_a = sm.____(after['returns'], after_with_intercept).____
# Get sum-of-squared residuals for both regressions
ssr_before = r_b.____
ssr_after = ____.ssr
# Compute and display the Chow test statistic
numerator = ((ssr_total - (ssr_before + ____)) / 2)
denominator = ((____ + ssr_after) / (24 - 4))
print("Chow test statistic: ", numerator / ____)