CommencezCommencez gratuitement

Gérer la tendance et l'hétéroscédasticité

Ici, nous allons rendre des données non stationnaires stationnaires en calculant le rendement ou le taux de croissance comme suit.

Souvent, les séries chronologiques sont générées selon $$X_t = (1 + p_t) X_{t-1}$$ ce qui veut dire que la valeur de la série observée au temps \(t\) est égale à la valeur observée au temps \(t-1\) multipliée par une petite variation en pourcentage \(p_t\) au temps \(t\).

Un exemple déterministe simple est de déposer de l'argent dans une banque avec un intérêt fixe \(p\). Dans ce cas, \(X_t\) est la valeur du compte à la période \(t\) avec un dépôt initial \(X_0\).

Typiquement, on désigne \(p_t\) comme le rendement ou le taux de croissance d'une série chronologique, et ce processus est souvent stable.

Pour des raisons qui dépassent la portée de ce cours, on peut montrer que le taux de croissance \(p_t\) peut être approximé par $$Y_t = \log X_t - \log X_{t-1} \approx p_t.$$

En R, on calcule souvent \(p_t\) avec diff(log(x)) et on peut le tracer en une seule ligne : plot(diff(log(x))).

Cette activité fait partie du cours

Modèles ARIMA en R

Voir le cours

Instructions de l’exercice

  • Comme précédemment, les forfaits astsa et xts sont préchargés.
  • Générez un graphique multifenêtre pour (1) tracer les données trimestrielles du PIB des É.-U. (gnp) et constater qu'elles ne sont pas stationnaires, puis (2) tracer le taux de croissance approximatif du PIB américain en utilisant diff() et log().
  • Utilisez un graphique multifenêtre pour (1) tracer les clôtures quotidiennes du DJIA (djia$Close) et constater qu'elles ne sont pas stationnaires. Les données sont un objet xts. Ensuite, (2) tracez les rendements approximatifs du DJIA en utilisant diff() et log(). Comment cela se compare-t-il au taux de croissance du PIB ?

Exercice interactif pratique

Essayez cet exercice en complétant ce code d’exemple.

# astsa and xts are preloaded 

# Plot GNP series (gnp) and its growth rate
par(mfrow = c(2,1))
plot(gnp)


# Plot DJIA closings (djia$Close) and its returns
par(mfrow = c(2,1))
plot(djia$Close)


 
Modifier et exécuter le code