Mémoire de Séries Temporelles
Forfatter:
LaurentLesage
Sidst opdateret:
9 år siden
Licens:
Creative Commons CC BY 4.0
Resumé:
Séries temporelles
\begin
Opdag hvorfor 18 millioner mennesker verden rundt stoler på Overleaf med deres arbejde.
Séries temporelles
\begin
Opdag hvorfor 18 millioner mennesker verden rundt stoler på Overleaf med deres arbejde.
\documentclass[a4paper,12pt]{article}
%%%%%%%%%%%%%%%%%%%%%
% pour écrire avec les accents directement dans le texte
\usepackage[frenchb,english]{babel}
\usepackage[T1]{fontenc}
\usepackage[colorlinks=false, urlcolor=black, breaklinks, pagebackref, citebordercolor={0 0 0}, filebordercolor={0 0 0}, linkbordercolor={0 0 0}, pagebordercolor={0 0 0}, runbordercolor={0 0 0}, urlbordercolor={0 0 0}, pdfborder={0 0 0}]{hyperref}
%%%%%%%%%%%%%%%%%%%%%
\usepackage{float}
\usepackage{amsmath}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{geometry}
\usepackage{graphicx}
\usepackage{hyperref}
\usepackage{icomma}
\usepackage{latexsym}
\usepackage[numbers]{natbib}
\usepackage{textcomp}
\usepackage[all]{xy}
\usepackage{cancel}
\usepackage{color}
\usepackage{dsfont}
\usepackage{subfig}
\usepackage{caption}
\usepackage{dsfont}
\usepackage{extarrows}
\usepackage{shorttoc}
\usepackage{listings}
\usepackage[utf8]{inputenc}
\usepackage{lmodern}
\usepackage{xcolor}
\usepackage{moreverb}
\lstset{language=scilab}
\lstset{
basicstyle=\ttfamily,
emphstyle=\color{blue},
columns=flexible,
breaklines=true,
keywordstyle=\color{blue},
commentstyle=\color{red},
}
\newcommand{\bigO}[1]{\ensuremath{\mathop{}\mathopen{}O\mathopen{}\left(#1\right)}}
\newcommand{\smallO}[1]{\ensuremath{\mathop{}\mathopen{}o\mathopen{}\left(#1\right)}}
\captionsetup{font=footnotesize}
\setcounter{tocdepth}{3}
\makeatletter
\def\clap#1{\hbox to 0pt{\hss #1\hss}}%
\def\ligne#1{%
\hbox to \hsize{%
\vbox{\centering #1}}}%
\def\haut#1#2#3{%
\hbox to \hsize{%
\rlap{\vtop{\raggedright #1}}%
\hss
\clap{\vtop{\centering #2}}%
\hss
\llap{\vtop{\raggedleft #3}}}}%
\def\bas#1#2#3{%
\hbox to \hsize{%
\rlap{\vbox{\raggedright #1}}%
\hss
\clap{\vbox{\centering #2}}%
\hss
\llap{\vbox{\raggedleft #3}}}}%
\def\maketitle{%
\thispagestyle{empty}\vbox to \vsize{%
\haut{}{\@blurb}{}
\vfill
\vspace{0.5cm}
\hrule height 2pt
\vspace{0.2cm}
\par
\begin{center}
\fontseries{bx}
\huge \@title
\end{center}
\vspace{0.2cm}
\par
\hrule height 2pt
\par
\vspace{2cm}
\begin{center}
\Large
\textit{Auteurs:}
\end{center}
\begin{center}
\Large
\@author
\par
\end{center}
\vspace{1cm}
\begin{center}
\Large
\textit{Responsable:}
\end{center}
\begin{center}
\Large
Christian Robert
\par
\end{center}
\vfill
\vfill
\bas{}{\begin{flushright}
\large
\@location
\end{flushright}}{}
}%
\cleardoublepage
}
\def\date#1{\def\@date{#1}}
\def\author#1{\def\@author{#1}}
\def\title#1{\def\@title{#1}}
\def\location#1{\def\@location{#1}}
\def\blurb#1{\def\@blurb{#1}}
\date{\today}
\author{}
\title{}
\location{Lyon}\blurb{}
\makeatother
\title{Mémoire de Séries Temporelles}
\author{Lesage Laurent - Maillart Arthur - Montagne Pierre}
\location{1 avril 2016}
\blurb{%
\begin{center}
\fbox{\includegraphics[width=6em]{Images/isfa.jpg}}
\end{center}
%Pour\\
%\textbf{type de stage}\\[1em]
%Maître de stage : moralles\\
}%
\begin{document}
\maketitle
\shorttoc{Plan}{1}
\newpage
\newpage
\part*{Introduction}\label{sec:name}
\addcontentsline{toc}{part}{Introduction}
\textcolor{blue}{En bleu : Laurent}
\vspace{0.5cm}
\hspace{1.5cm} L'objectif de ce mémoire est de mettre en oeuvre les techniques d'estimation de tendance, de saisonnalité et de composantes aléatoires, pour ensuite réaliser des prédictions sur deux jeux de données. Les données concernent un historique de la production hydraulique d'électricité et de consommation d'énergie, de 1981 jusque 2015. Les deux prédictions porent sur l'année 2016. \\\\
\textcolor{red}{Proposition modification plan.}
I-Série 1\\
a) méthode 1\\
b) méthode 2\\
II-Série 2\\
a) méthode 1\\
b) méthode 2\\
Permet de bien dérouler les méthodes.
\section{Modélisation de la série de la consommations réelle}
Nous commençons par quelques remarques sur le chronogramme de la série de la consommation réelle. Sur la représentation graphique ci-dessous de la série nous pouvons observer un trend qui a priori est quadratique ainsi qu'une saisonnalité d'ordre 12. Notons que la variance semble augmenter après 2006. \textbf{Mettre un graphique propre}
\begin{center}
\includegraphics[scale=0.5]{Images/ConsommationReelle.jpeg}
\end{center}
Dans un premier temps, nous avons utilisé la méthode de Buys-Ballot \textcolor{red}{(Attention aux spec du modèle)} sur le logarithme de la série (pour atténuer les effets de la variance après 2006). Nous avons effectué un décomposition avec un trend quadratique et une saisonnalité d'ordre 12. Cependant, nous n'avons pas réussi à modéliser la composante aléatoire de manière satisfaisante. En observant bien les ACF et PACF obtenus, nous avons préféré éliminer uniquement le trend et modéliser avec un SARIMA la saisonnalité et le bruit (a reformuler). Nous détaillons maintenant notre approche.
\subsection{Buys-Ballot}
\begin{itemize}
\item \textbf{Objectif de la méthode :}
\\
La méthode de Buys-Ballot généralisée consiste à proposer une décomposition du type tendance + saisonnalité + composante aléatoire.
\begin{center}
$X_t = m_t + S_t + Y_t$ avec t=1,..., T.
\end{center}
\item \textbf{Conditions :}
\\
Les conditions à respecter sont les suivantes :
\begin{itemize}
\item La tendance $m_t$ se décompose de la sorte :\\
\begin{center}
$m_{t} = \sum_{i=1}^{k} b_{i} m_{t}^{(i)}$
\end{center}
- les $b_{i}$ des constantes à déterminer \\
- k l'ordre de la tendance.
\item La tendance $S_t$ se décompose de la sorte :\\
\begin{center}
$S_t = \sum_{j=1}^{l} c_{j} S_{t}^{(j)}$
\end{center}
- avec les $c_{i}$ des constantes à déterminer \\
- k la saisonnalité à modéliser.
\item La composante aléatoire $Y_{t}$ est caractérisée par :\\
- Son Espérance : $\mathbb{E}$[$Y_{t}$] = 0 \\
- Sa Variance : V[$Y_{t}$] = constante
\begin{center}
\end{center}
\end{itemize}
\item \textbf{Intérêt principal :}
\\
Pouvoir faire des prévisions en utilisant la tendance et la saisonnalité (la composante aléatoire est d'espérance nulle).
\item \textbf{Résumé des éléments :}
\\
Décrire le tableau avec les éléments
\end{itemize}
\subsection{Modéliser un SARIMA sur la série entière}
\textcolor{blue}{Grand schéma synthétique avec les différentes étapes successives : méthodologie de Box et Jenkins}
\textcolor{blue}{Ajout de sous titres correspondant aux différentes étapes}
On commence par représenter l'autocorrélogramme de la série brute. Cet autocorrélogramme montre une décroissance lente ainsi qu'une saisonnalité d'ordre 12. Pour stationnariser la série nous allons donc différencier la série une fois. Ceci a pour effet de faire disparaître le trend mais pas la saisonnalité. Maintenant, nous faisons une différence avec un lag de 12 suivie d'une différenciation classique. Cette fois ci le trend et la saisonnalité on été éliminés. L'autocorrélogramme que nous obtenons semble stationnaire.
\begin{center}
\includegraphics[scale=0.5]{Images/ACF_serie_conso_ini.jpg}
\end{center}
Nous effectuons donc un test pour confirmer notre hypothèse.
\textcolor{blue}{Test bruit blanc}
\textcolor{blue}{Pas Dickey Fuller qui donne de la merde}
Les résultats du test confirment notre intuition. Nous pouvons donc continuer l'analyse.
\textcolor{blue}{Hétéroscédasticité}
Notre idée ici est d'ajuster un modèle SARIMA sur le logarithme de la série. Nous venons de déterminer en stationnarisant la série que le paramètre $d$ vaut 1, le paramètre de saisonnalité vaut 12, et que $D=1$. Reste maintenant à déterminer $p,q,P,D,Q$. Pour obtenir $p_{max}$, $q_{max}$,$P_{max}$, $Q_{max}$ nous allons observer plus en détail l'autocorrélogramme et l'autocorrélogramme partiel. \textcolor{red}{Les graphiques sont dégueux mais c'est pour illustrer}
Posons $s=12$, alors on observe un pic significatif en s puis, il n'y en a plus. Par ailleurs, il y a un pic significatif qui suit les pics en $0 \times s $ et $ 1 \times s$, comme nous pouvons le voir ci-dessous.
\begin{center}
\includegraphics[scale=0.4]{Images/ACF_serie_conso_statio.jpg}
\end{center}
Ainsi nous posons $p_{max}=1$ et $P_max=1$. Sur l'autocorrélogramme partiel, on observe 2 pics saisonniers ainsi qu'un pic significatif en dehors de ces pics saisonniers.
\begin{center}
\includegraphics[scale=0.4]{Images/PACF_serie_conso_statio.jpg}
\end{center}
Ceci nous incite à poser $Q_{max}=2$ et $q_{max}=2$
Rappeler le modèle théorique.
Désormais, nous avons tous les paramètres pour ajuster notre modèle. Après quelques essais infructueux, nous observons un modèle qui se détache.
\vspace{0.5cm}
\hspace{1cm} A développer ...
\textcolor{blue}{Comparaison des modèles acceptables avec plusieurs critères : AIC, BIC}
\newpage
\part*{Conclusion}\label{sec:name}
\addcontentsline{toc}{part}{Conclusion}
\newpage
\tableofcontents
\end{document}