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 taux de croissance comme suit.
Souvent, les séries temporelles sont générées selon $$X_t = (1 + p_t) X_{t-1}$$ ce qui signifie que la valeur observée à l’instant \(t\) est égale à la valeur observée à l’instant \(t-1\) plus une petite variation en pourcentage \(p_t\) à l’instant \(t\).
Un exemple déterministe simple consiste à placer de l’argent à la 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\).
En général, \(p_t\) est appelé le rendement ou le taux de croissance d’une série temporelle, et ce processus est souvent stable.
Pour des raisons qui dépassent le cadre de ce cours, on peut montrer que le taux de croissance \(p_t\) peut être approché par $$Y_t = \log X_t - \log X_{t-1} \approx p_t.$$
En R, on calcule souvent \(p_t\) via diff(log(x)) et on peut le tracer en une seule ligne : plot(diff(log(x))).
Cet exercice fait partie du cours
Modèles ARIMA en R
Instructions
- Comme précédemment, les packages astsa et xts sont préchargés.
- Générez un graphique multifigure pour (1) tracer les données trimestrielles du PIB des États-Unis (
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 multifigure pour (1) tracer les cours de clôture quotidiens du DJIA (
djia$Close) et constater qu’ils ne sont pas stationnaires. Les données sont un objetxts. Puis (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 cet exemple de code.
# 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)