spb/gsf6053_h25 Public

1% =========================================================================2% GSF6053 — Économétrie Financière (Hiver 2025)3% Université Laval — Faculté des sciences de l'administration4% Auteur  : Simon-Pierre Boucher5% Contact : contact@spboucher.ai6% =========================================================================7\documentclass[14pt]{extarticle} % Utilisation de extarticle pour supporter 14pt89% ------------------------------------------------------------------------10% Packages indispensables et recommandés11% ------------------------------------------------------------------------12\usepackage[utf8]{inputenc}     % Encodage des caractères13\usepackage[T1]{fontenc}        % Encodage de la police14\usepackage[french]{babel}      % Support de la langue française15\usepackage{amsmath, amssymb}   % Environnements mathématiques enrichis16\usepackage{amsthm}             % Environnements théorèmes17\usepackage{graphicx}           % Inclusion d'images18\usepackage{hyperref}           % Liens hypertextes19\usepackage{geometry}           % Gestion des marges20\usepackage{titlesec}           % Personnalisation des titres21\usepackage{setspace}           % Gestion de l'interligne22\usepackage{xcolor}             % Gestion des couleurs23\usepackage{booktabs}           % Pour améliorer les tableaux24\usepackage{fancyhdr}           % Personnalisation des en-têtes et pieds de page25\usepackage{cleveref}           % Références intelligentes26\usepackage{caption}            % Pour personnaliser les légendes27\usepackage{enumitem}           % Pour personnaliser les listes28\usepackage{pdflscape}          % Pour placer les tableaux en paysage29\usepackage[backend=biber,style=apa]{biblatex} % Gestion de la bibliographie30\addbibresource{references.bib} % Fichier de bibliographie3132% ------------------------------------------------------------------------33% Configuration de la mise en page34% ------------------------------------------------------------------------35\geometry{36    a4paper,37    margin=25mm38}3940% ------------------------------------------------------------------------41% Définition des environnements théorèmes, définitions, etc.42% ------------------------------------------------------------------------43\theoremstyle{definition}44\newtheorem{definition}{Définition}[section]4546\theoremstyle{plain}47\newtheorem{proposition}{Proposition}[section]48\newtheorem{theorem}{Théorème}[section]4950% ------------------------------------------------------------------------51% Personnalisation des sections52% ------------------------------------------------------------------------53\titleformat{\section}{\Large\bfseries}{\thesection}{1em}{}54\titleformat{\subsection}{\large\bfseries}{\thesubsection}{1em}{}55\titleformat{\subsubsection}{\normalsize\bfseries}{\thesubsubsection}{1em}{}5657% ------------------------------------------------------------------------58% Commande pour mettre en avant les livres en bleu59% ------------------------------------------------------------------------60\newcommand{\livre}[1]{\textcolor{blue}{#1}}6162% ------------------------------------------------------------------------63% Configuration des hyperliens64% ------------------------------------------------------------------------65\hypersetup{66    colorlinks=true,          % Les liens sont colorés67    linkcolor=blue,           % Couleur des liens internes (TOC, références, etc.)68    urlcolor=blue,            % Couleur des URL69    citecolor=blue,           % Couleur des citations70    filecolor=blue,           % Couleur des liens vers des fichiers71    pdfborder={0 0 0},        % Pas de bordure autour des liens72    breaklinks=true            % Permet les sauts de ligne dans les liens73}7475% ------------------------------------------------------------------------76% Configuration des en-têtes et pieds de page77% ------------------------------------------------------------------------78\pagestyle{fancy}79\fancyhf{}80\fancyhead[L]{Autocorrélation des erreurs}81\fancyhead[R]{Hiver 2025}82\fancyfoot[C]{\thepage}83\fancyfoot[R]{\includegraphics[height=1cm]{logo_universite_laval.png}} % Ajout du logo ERS en bas à droite8485% ------------------------------------------------------------------------86% Informations sur le document87% ------------------------------------------------------------------------88\title{\textbf{Autocorrélation des erreurs}}89\author{GSF-6053}90\date{Hiver 2025}9192% ------------------------------------------------------------------------93% Début du document94% ------------------------------------------------------------------------95\begin{document}96\maketitle9798\tableofcontents99100\onehalfspacing101\newpage102103\section{Introduction}104105\textbf{\textcolor{blue}{\cite{gujarati2010}, \cite{wooldridge2010}}}106107L’autocorrélation ou la corrélation sérielle des erreurs dans les séries chronologiques est assez fréquente. Elle découle souvent du fait qu’il y a une certaine inertie dans les données économiques et financières, c’est-à-dire que les observations passées se reflètent souvent dans les observations présentes et futures. Cela peut amener des problèmes de corrélation sérielle des erreurs si une variable autocorrélée est omise du modèle, par exemple.108109Il est important de s’assurer que lorsque vous détectez de l’autocorrélation, elle n’est pas due à la non-stationnarité de $Y$ et $X$. En effet, si ces deux quantités sont non stationnaires, les erreurs le seront possiblement aussi et seront possiblement autocorrélées.110111Même si le modèle est bien spécifié, on peut penser que les erreurs (les chocs) affectant les marchés boursiers aujourd’hui ont des chances d’influencer la magnitude et le signe des chocs affectant les marchés boursiers demain aussi. Dans la plupart des cas en finance, l’autocorrélation sera positive, mais il est possible en théorie d’avoir de l’autocorrélation négative.112113\section{Modèle de Base avec Autocorrélation}114115\textbf{\textcolor{blue}{\cite{gujarati2010}, \cite{wooldridge2010}}}116117Soit le modèle de base :118\begin{equation}119y = X\beta + u120\end{equation}121122\begin{equation}123Y_t = \beta_0 + \beta_1 X_{1t} + \beta_2 X_{2t} + \dots + \beta_{K-1} X_{K-1,t} + u_t124\end{equation}125126\begin{equation}127u_t = \rho u_{t-1} + \epsilon_t128\end{equation}129130Où $\epsilon_t$ est un bruit blanc de moyenne nulle et de variance constante dans le temps.131132\textbf{Différence avec le modèle linéaire de base :}133134La matrice de variance-covariance des erreurs sera affectée, ce qui impliquera, comme dans le cas hétéroscédastique, que l’estimateur OLS standard ne sera pas BLUE et que la variance des estimateurs OLS sera incorrecte, bien que l’estimateur soit sans biais. La démonstration de ceci est du même ressort que celle faite pour l’hétéroscédasticité. En effet, bien que l’autocorrélation implique des erreurs homoscédastiques (la variance est constante), les termes hors diagonale ne sont pas nuls comme dans le cas de base du modèle linéaire.135136\subsection{Matrice de Variance-Covariance des Erreurs}137138\[139\mathbf{E}(uu') = 140\begin{bmatrix}141\mathbb{E}(u_1^2) & \mathbb{E}(u_1 u_2) & \mathbb{E}(u_1 u_3) & \dots & \mathbb{E}(u_1 u_T) \\142\mathbb{E}(u_2 u_1) & \mathbb{E}(u_2^2) & \mathbb{E}(u_2 u_3) & \dots & \mathbb{E}(u_2 u_T) \\143\mathbb{E}(u_3 u_1) & \mathbb{E}(u_3 u_2) & \mathbb{E}(u_3^2) & \dots & \mathbb{E}(u_3 u_T) \\144\vdots & \vdots & \vdots & \ddots & \vdots \\145\mathbb{E}(u_T u_1) & \mathbb{E}(u_T u_2) & \mathbb{E}(u_T u_3) & \dots & \mathbb{E}(u_T^2)146\end{bmatrix}147\]148149Avec :150151\[152\mathbf{E}(uu') = \sigma^2 153\begin{bmatrix}1541 & \rho & \rho^2 & \dots & \rho^{T-1} \\155\rho & 1 & \rho & \dots & \rho^{T-2} \\156\rho^2 & \rho & 1 & \dots & \rho^{T-3} \\157\vdots & \vdots & \vdots & \ddots & \vdots \\158\rho^{T-1} & \rho^{T-2} & \rho^{T-3} & \dots & 1159\end{bmatrix}160= \sigma^2 \Omega161\]162163Où :164165\[166\Omega = \begin{bmatrix}1671 & \rho & \rho^2 & \dots & \rho^{T-1} \\168\rho & 1 & \rho & \dots & \rho^{T-2} \\169\rho^2 & \rho & 1 & \dots & \rho^{T-3} \\170\vdots & \vdots & \vdots & \ddots & \vdots \\171\rho^{T-1} & \rho^{T-2} & \rho^{T-3} & \dots & 1172\end{bmatrix}173\]174175\subsection{Preuve de la Forme de $\Omega$}176177\textbf{Hypothèses de départ :}178179\begin{enumerate}180    \item $\epsilon_t$ est un bruit blanc, c'est-à-dire :181    \[182    \mathbb{E}(\epsilon_t) = 0, \quad \mathbb{V}(\epsilon_t) = \sigma^2, \quad \forall t \quad \text{et} \quad \mathbb{E}(\epsilon_t \epsilon_s) = 0 \quad \forall t \neq s183    \]184    185    \item $|\rho| < 1$ (absence de racine unitaire).186    187    Ceci implique que $u_t$ est un processus stationnaire, c'est-à-dire :188    \begin{enumerate}189        \item $\mathbb{E}(u_t) = 0$190        \item $\mathbb{V}(u_t) = \frac{\sigma^2}{1 - \rho^2}$, \quad $\forall t$191        \item $\mathbb{E}(u_t u_{t-j}) = \frac{\rho^j \sigma^2}{1 - \rho^2} \quad \forall j \geq 1$192    \end{enumerate}193    194    \item $\epsilon_t$ n’est pas corrélé avec les $u_s$ dont l’indice $s$ est antérieur à $t$, c'est-à-dire :195    \[196    \epsilon_t \perp u_{t-s}, \quad \forall s \geq 1197    \]198\end{enumerate}199200\textbf{Développement :}201202\begin{align}203u_t &= \rho u_{t-1} + \epsilon_t \nonumber \\204&= \rho (\rho u_{t-2} + \epsilon_{t-1}) + \epsilon_t \nonumber \\205&= \rho^2 u_{t-2} + \rho \epsilon_{t-1} + \epsilon_t \nonumber \\206&= \dots \nonumber \\207u_t &= \sum_{s=0}^{\infty} \rho^s \epsilon_{t-s} \label{eq:autoregressive}208\end{align}209210\textbf{Calcul de la Variance de $u_t$ :}211212\begin{align}213\mathbb{V}(u_t) &= \mathbb{V}\left( \sum_{s=0}^{\infty} \rho^s \epsilon_{t-s} \right) \nonumber \\214&= \sum_{s=0}^{\infty} (\rho^s)^2 \mathbb{V}(\epsilon_{t-s}) \nonumber \\215&= \sigma^2 \sum_{s=0}^{\infty} \rho^{2s} \nonumber \\216&= \sigma^2 \left( \frac{1}{1 - \rho^2} \right) \quad \text{(série géométrique convergente pour } |\rho| < 1) \label{eq:variance_ut}217\end{align}218219Ainsi :220221\[222\mathbb{V}(u_t) = \frac{\sigma^2}{1 - \rho^2}223\]224225\textbf{Calcul des Covariances :}226227\begin{align}228\mathbb{E}(u_t u_{t-j}) &= \mathbb{E}\left( \sum_{s=0}^{\infty} \rho^s \epsilon_{t-s} \sum_{k=0}^{\infty} \rho^k \epsilon_{t-j-k} \right) \nonumber \\229&= \sum_{s=0}^{\infty} \sum_{k=0}^{\infty} \rho^{s+k} \mathbb{E}(\epsilon_{t-s} \epsilon_{t-j-k}) \nonumber \\230&= \sum_{s=0}^{\infty} \rho^{s+j+s} \mathbb{E}(\epsilon_{t-s} \epsilon_{t-j-s}) \nonumber \\231&= \sum_{s=0}^{\infty} \rho^{2s+j} \mathbb{E}(\epsilon_{t-s} \epsilon_{t-j-s}) \nonumber \\232&= \rho^j \sigma^2 \sum_{s=0}^{\infty} \rho^{2s} \nonumber \\233&= \rho^j \sigma^2 \left( \frac{1}{1 - \rho^2} \right) \nonumber \\234&= \frac{\rho^j \sigma^2}{1 - \rho^2} \label{eq:covariance_ut_uj}235\end{align}236237Ainsi, pour toute autocorrélation d'ordre $j$, on a :238239\[240\mathbb{E}(u_t u_{t-j}) = \frac{\rho^j \sigma^2}{1 - \rho^2}241\]242243\section{La Transformation de Prais-Winsten}244245\textbf{\textcolor{blue}{\cite{gujarati2010}, \cite{wooldridge2010}}}246247Comme on l’a fait dans le cas de l’hétéroscédasticité, on peut utiliser un modèle transformé (donc un GLS ou un FGLS) pour utiliser les OLS en présence d’autocorrélation des erreurs. La transformation de Prais-Winsten consiste à multiplier les éléments de la régression par la matrice $P^{-1}$.248249\[250P^{-1} \begin{bmatrix}251Y_1 \\252Y_2 \\253Y_3 \\254\vdots \\255Y_{T-1} \\256Y_T257\end{bmatrix}258=259\begin{bmatrix}260\sqrt{1 - \rho^2} Y_1 \\261-\rho Y_1 + Y_2 \\262-\rho Y_2 + Y_3 \\263\vdots \\264-\rho Y_{T-1} + Y_T265\end{bmatrix}266\]267268La transformation sera la même pour chaque régresseur $X$. Le modèle transformé s’écrit :269270\begin{align}271\sqrt{1 - \rho^2} Y_1 &= \sqrt{1 - \rho^2} \beta_0 + \beta_1 \sqrt{1 - \rho^2} X_{11} + \beta_2 \sqrt{1 - \rho^2} X_{21} + \dots + \sqrt{1 - \rho^2} u_1 \label{eq:prais_winsten1} \\272Y_t - \rho Y_{t-1} &= (1 - \rho) \beta_0 + \beta_1 (X_{1t} - \rho X_{1,t-1}) + \beta_2 (X_{2t} - \rho X_{2,t-1}) + \dots \nonumber \\273&\quad + \beta_{K-1} (X_{K-1,t} - \rho X_{K-1,t-1}) + u_t - \rho u_{t-1} \label{eq:prais_winsten2}274\end{align}275276Il faut porter attention à la constante dans ce modèle qui n’est plus la constante d’origine. Si $\rho$ est inconnu, ce modèle n’est plus linéaire, mais on peut toujours l’estimer par OLS.277278\section{La Transformation de Cochran-Orcutt}279280\textbf{\textcolor{blue}{\cite{gujarati2010}, \cite{wooldridge2010}}}281282Une autre façon de transformer le modèle est de laisser tomber la première observation et de transformer les $Y$ comme suit :283284\[285Y_t - \rho Y_{t-1}, \quad t = 2, \dots, T286\]287288Les $X$ sont transformés de la même façon.289290Le modèle transformé de Cochran-Orcutt s’écrit de la même façon que Prais-Winsten, mais sans la première observation :291292\begin{align}293Y_t - \rho Y_{t-1} &= (1 - \rho) \beta_0 + \beta_1 (X_{1t} - \rho X_{1,t-1}) + \beta_2 (X_{2t} - \rho X_{2,t-1}) + \dots \nonumber \\294&\quad + \beta_{K-1} (X_{K-1,t} - \rho X_{K-1,t-1}) + u_t - \rho u_{t-1} \label{eq:cochran_orcutt}295\end{align}296297\section{Le Problème de $\rho$}298299$\rho$ est inconnu. Il faut donc trouver un estimateur. Plusieurs sont suggérés dans la littérature économétrique.300301\subsection{Estimateur Usuel (\textcolor{blue}{\cite{gujarati2010}})}302303\[304\hat{\rho} = \frac{\sum_{t=2}^{T} \hat{u}_t \hat{u}_{t-1}}{\sum_{t=2}^{T} \hat{u}_{t-1}^2}305\]306307Où $\hat{u}_t$ est le résidu OLS de la régression de $Y$ sur $X$. On peut montrer que $\text{PLIM} \ \hat{\rho} = \rho$.308309\subsection{Estimateur de Hildreth-Lu}310311Cet estimateur ne passe pas par les OLS. En balayant l’intervalle $(-1, 1)$, on peut choisir la valeur de $\rho$ qui minimise la somme des carrés des erreurs du modèle transformé. Les étapes sont :312313\begin{enumerate}314    \item Choisir $\rho^{(1)} = -0.99999$.315    \item Obtenir le modèle transformé correspondant.316    \item Appliquer les OLS et sauvegarder $\hat{u}'\hat{u}_{\rho^{(1)}}$.317    \item Choisir une nouvelle valeur : $\rho^{(2)} = \rho^{(1)} + \text{step}$.318    \item Obtenir le modèle transformé correspondant.319    \item Appliquer les OLS et sauvegarder $\hat{u}'\hat{u}_{\rho^{(2)}}$.320    \item Continuer jusqu’à ce que $(-1, 1)$ soit couvert.321    \item Choisir la valeur $\rho^{(s)}$ telle que322    \[323    \hat{u}'\hat{u}_{\rho^{(s)}} = \min \left[ \hat{u}'\hat{u}_{\rho^{(i)}} \right]_i324    \]325\end{enumerate}326327L’estimateur FGLS obtenu à partir d’un estimateur de $\rho$ sera convergent. Donc, lorsque la matrice d’information $P$ dépend de paramètres inconnus qu’il faut estimer, on perd la propriété BLUE. L’estimateur qui est un FGLS au lieu de GLS reste convergent si l’estimateur des paramètres inconnus sur lequel il est fondé est convergent.328329\section{Diagnostic de l’Autocorrélation}330331\textbf{\textcolor{blue}{\cite{gujarati2010}}}332333On peut suspecter l’autocorrélation en faisant un graphique entre les résidus aujourd’hui et hier. Ce n’est pas un test statistique. Il existe des tests d’autocorrélation des erreurs.334335\subsection{Test de Durbin-Watson}336337C’est un test à bornes. Vous avez deux points critiques à regarder dans la table. À l’origine, ce test était à bornes, car on ne connaissait pas les vrais points critiques seulement des bornes sous l’hypothèse de normalité. Maintenant, on a trouvé les vrais points critiques.338339On teste :340341\[342H_0 : \rho = 0 \quad \text{contre} \quad H_A : \rho \neq 0343\]344345La statistique de test est :346347\begin{equation}348d = \frac{\sum_{t=2}^{T} ( \hat{u}_t - \hat{u}_{t-1} )^2}{\sum_{t=1}^{T} \hat{u}_t^2} \approx 2(1 - \hat{\rho})349\end{equation}350351\textbf{Hypothèses du test de Durbin-Watson :}352\begin{enumerate}353    \item Le modèle de régression inclut une constante pour pouvoir calculer RSS.354    \item Les variables explicatives sont non stochastiques.355    \item L’autocorrélation est d’ordre 1. On ne peut pas utiliser Durbin-Watson avec des ordres supérieurs.356    \item On ne peut pas utiliser ce test dans le cas de modèles autorégressifs.357    \item Le test n’accommode pas les observations manquantes.358\end{enumerate}359360Pour faire le test de Durbin-Watson, il faut effectuer une régression OLS et obtenir les résidus. Ensuite, on calcule la statistique de Durbin-Watson avec ces résidus. On compare la statistique avec les points critiques de la distribution pour les statistiques de Durbin-Watson pour une taille d’échantillon donnée et un $K$ donné. On suit la règle de décision suivante :361362\[363\begin{cases}364d < d_L & \rightarrow \text{rejeter } H_0 \\365d_L < d < d_U & \rightarrow \text{test non concluant} \\366d_U < d < 4 - d_U & \rightarrow \text{ne pas rejeter } H_0 \text{ et } H_0^* \\3674 - d_U < d < 4 - d_L & \rightarrow \text{test non concluant} \\368d > 4 - d_L & \rightarrow \text{rejeter } H_0^*369\end{cases}370\]371372La statistique de Durbin-Watson est une « institution » en économétrie parce qu’il s’agit d’un des premiers tests à bornes dérivés, mais les hypothèses sous-jacentes sont assez contraignantes.373374\subsection{Le Test de Breusch-Godfrey}375376Ce test est plus général et plus populaire que Durbin-Watson. Il permet de tester pour des formes d’autocorrélation avec plusieurs retards ($q$). L’intérêt du test est qu’il se généralise facilement pour $q$ lags. Pour $q$ retards, la régression artificielle est de $\hat{u}_t$ sur une constante, $X_{1t}, \dots, X_{kt}$ et $\hat{u}_{t-1}, \dots, \hat{u}_{t-q}$.377378\[379\hat{u}_t = \alpha_1 + \alpha_2 X_t + \dots + \rho_1 \hat{u}_{t-1} + \dots + \rho_q \hat{u}_{t-q} + \epsilon_t380\]381382\textbf{Hypothèse nulle :}383\[384H_0 : \rho_1 = \rho_2 = \dots = \rho_q = 0 \quad \text{contre l’alternative qu’au moins une des autocorrélations est non nulle.}385\]386387La statistique de test est :388389\[390BG = (T - q) R^2 \sim \chi^2(q)391\]392393\subsection{Estimateur Robuste de Newey-West}394395\textbf{\textcolor{blue}{\cite{gujarati2010}, \cite{greene2018}}}396397Il existe une forme de correcteur robuste pour les problèmes d’autocorrélation des erreurs qui est une alternative au FGLS lorsque le nombre d’observations est assez élevé. En effet, cette correction est asymptotique. C’est l’estimateur de la variance de Newey-West. L’intuition est la même que le correcteur de White, mais appliquée aux problèmes d’autocorrélation des erreurs. En fait, le charme de l’estimateur de Newey-West est qu’il corrige pour l’autocorrélation ET l’hétéroscédasticité des erreurs (c’est ce qu’on appelle l’estimateur « HAC 1 »). C’est un correcteur très populaire en finance.398399Plus généralement, la classe d’estimateurs corrige les déviations dans la matrice de variance-covariance des erreurs par rapport au cas simple.400401Nous n’allons pas faire la preuve ou la dérivation, car elle est assez avancée. Pour ceux qui veulent aller plus loin, elle est présentée dans la section 20.5.2 de \textcolor{blue}{\cite{greene2018}}. Le correcteur de Newey-West est déjà programmé dans les logiciels d’analyse de données comme Stata.402403\section{Conclusion}404405L'autocorrélation des erreurs est une problématique courante dans les modèles de régression linéaire appliqués aux séries chronologiques. Les tests diagnostiques tels que Durbin-Watson et Breusch-Godfrey permettent de détecter sa présence et de choisir les méthodes appropriées pour la corriger, assurant ainsi la fiabilité des estimations des paramètres du modèle. Des méthodes de correction comme Prais-Winsten, Cochran-Orcutt et l'estimateur de Newey-West offrent des solutions robustes pour traiter l'autocorrélation et améliorer les performances des modèles économétriques.406407% ------------------------------------------------------------------------408% Bibliographie409% ------------------------------------------------------------------------410\printbibliography411412\end{document}413