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
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 utilisantdiff()etlog(). - 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 objetxts. Ensuite, (2) tracez les rendements approximatifs du DJIA en utilisantdiff()etlog(). 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)