📚 الرئيسية

مسابقة دكتوراه 2022Centre Universitaire Abdelhafid Boussouf - Mila — الموضوع 01

مسابقة تخصص · Probabilités & Statistiques · المدة: 2سا

Concours d'accès à la formation de doctorat 3ème cycle, Spécialité : Probabilités et statistiques, Matière : Séries chronologiques (Sujet 1), Centre Universitaire Abdelhafid Boussouf - Mila, Institut des Sciences et Technologie, Département de Mathématiques et Informatique, Laboratoire de Mathématiques et leurs Interactions, 05 mars 2022, durée 02 heures.

التمرين 1

Exercice 1 — Série chronologique : tendance linéaire et exponentielle

#time-series#trend-estimation#least-squares#additive-model

Considérons la série chronologique présentée dans le tableau suivant :

t12345678910xt2.005.736.0210.3610.0613.0315.0116.3918.2220.46\begin{array}{|c|c|c|c|c|c|c|c|c|c|c|} \hline t & 1 & 2 & 3 & 4 & 5 & 6 & 7 & 8 & 9 & 10 \\\\ \hline x_t & 2.00 & 5.73 & 6.02 & 10.36 & 10.06 & 13.03 & 15.01 & 16.39 & 18.22 & 20.46 \\\\ \hline \end{array}

  1. Représenter graphiquement cette série puis proposer un modèle parmi la famille des modèles additifs Yt=f(t)+St+XtY_t = f(t) + S_t + X_t ? Justifier votre réponse.
  2. En supposant que la série n'a pas de composante saisonnière et que la tendance est linéaire (f(t)=at+bf(t) = at + b), estimer cette tendance puis calculer les valeurs de la variable centrée.
  3. En supposant que la série n'a pas de composante saisonnière et que la tendance est exponentielle (f(t)=ect+df(t) = e^{ct + d}), estimer cette tendance puis calculer les valeurs de la variable centrée.
  4. Quel est le meilleur modèle pour cette série ?
الحل

1.

Le nuage de points croît de façon quasi linéaire sans oscillation périodique : on propose un modèle additif à tendance seule (pas de saisonnalité), Yt=f(t)+XtY_t = f(t) + X_t.

2.

Moindres carrés : tˉ=5.5\bar t = 5.5, xˉ=11.728\bar x = 11.728, tx=806.21\sum t x = 806.21, t2=385\sum t^2 = 385.

a=txntˉxˉt2ntˉ2=806.21645.04385302.51.953,b=xˉatˉ0.988.a = \frac{\sum tx - n\bar t\bar x}{\sum t^2 - n\bar t^2} = \frac{806.21 - 645.04}{385 - 302.5} \approx 1.953, \quad b = \bar x - a\bar t \approx 0.988.

f(t)1.953t+0.988\boxed{f(t) \approx 1.953\,t + 0.988}

La variable centrée est Xt=xtf(t)X_t = x_t - f(t) (résidus).

3.

On linéarise : lnxt=ct+d\ln x_t = ct + d. Une régression linéaire de lnxt\ln x_t sur tt donne c^\hat c et d^\hat d, puis f(t)=ec^t+d^f(t) = e^{\hat c t + \hat d} et Xt=xtf(t)X_t = x_t - f(t).

4.

Comme les données sont presque alignées (croissance additive constante d'environ 22 par pas), le modèle linéaire ajuste mieux (résidus plus faibles) que l'exponentiel.

Modeˋle lineˊaire f(t)=at+b\boxed{\text{Modèle linéaire } f(t) = at + b}

التمرين 2

Exercice 2 — Processus AR(1) : causalité, stationnarité et estimation

#ar-process#stationarity#maximum-likelihood#time-series

Soit (Xt)(X_t) le processus donné par

Xt=ϕXt1+εt,X_t = \phi X_{t-1} + \varepsilon_t,

ϕ<1|\phi| \lt 1 et les εt\varepsilon_t sont des variables aléatoires i.i.d. de loi normale N(0,σ2)N(0, \sigma^2).

  1. Montrer que le processus XtX_t est causal.
  2. Le processus XtX_t est-il stationnaire au second ordre ? Justifier votre réponse.
  3. Écrire la fonction de vraisemblance de (X1,X2,,XN)(X_1, X_2, \ldots, X_N).
  4. En déduire les estimateurs du maximum de vraisemblance des paramètres inconnus ϕ\phi et σ2\sigma^2.
الحل

1.

Comme ϕ<1|\phi| \lt 1, la série Xt=j0ϕjεtjX_t = \sum_{j \geq 0} \phi^j \varepsilon_{t-j} converge (dans L2L^2) et n'utilise que le passé du bruit : le processus est causal.

2.

E(Xt)=0\mathbb{E}(X_t) = 0, Var(Xt)=σ21ϕ2\mathrm{Var}(X_t) = \dfrac{\sigma^2}{1 - \phi^2}, et γ(h)=σ2ϕh1ϕ2\gamma(h) = \dfrac{\sigma^2 \phi^{|h|}}{1 - \phi^2} ne dépend que du retard hh. Le processus est donc stationnaire au second ordre.

3.

En conditionnant, XtXt1N(ϕXt1,σ2)X_t \mid X_{t-1} \sim N(\phi X_{t-1}, \sigma^2) et X1N ⁣(0,σ21ϕ2)X_1 \sim N\!\left(0, \frac{\sigma^2}{1-\phi^2}\right) :

L=1ϕ22πσ2e(1ϕ2)x122σ2t=2N12πσ2exp ⁣((xtϕxt1)22σ2).L = \sqrt{\frac{1-\phi^2}{2\pi\sigma^2}}\,e^{-\frac{(1-\phi^2)x_1^2}{2\sigma^2}} \prod_{t=2}^{N}\frac{1}{\sqrt{2\pi\sigma^2}}\exp\!\left(-\frac{(x_t - \phi x_{t-1})^2}{2\sigma^2}\right).

4.

En maximisant la vraisemblance conditionnelle :

ϕ^=t=2Nxtxt1t=2Nxt12,σ^2=1N1t=2N(xtϕ^xt1)2\boxed{\hat\phi = \frac{\sum_{t=2}^N x_t x_{t-1}}{\sum_{t=2}^N x_{t-1}^2}, \qquad \hat\sigma^2 = \frac{1}{N-1}\sum_{t=2}^N (x_t - \hat\phi x_{t-1})^2}

التمرين 3

Exercice 3 — Espérance conditionnelle et martingale des sommes centrées

#conditional-expectation#martingales#independence#filtration

Soit (Ω,F,P)(\Omega, \mathcal{F}, P) un espace de probabilité, G\mathcal{G} sous-tribu de F\mathcal{F} et XX une variable aléatoire.

  1. Montrer que l'application φ:XE(XG)\varphi : X \longmapsto \mathbb{E}(X \mid \mathcal{G}) est linéaire.
  2. Si XX est une variable aléatoire de L1(Ω,F,P)\mathbb{L}^1(\Omega, \mathcal{F}, P) indépendante de G\mathcal{G}, montrer que E(XG)=E(X)\mathbb{E}(X \mid \mathcal{G}) = \mathbb{E}(X).
  3. Montrer que si XX est G\mathcal{G}-mesurable, alors E(XG)=X\mathbb{E}(X \mid \mathcal{G}) = X p.s.
  4. Soient X0,X1,,XnX_0, X_1, \ldots, X_n une suite de variables aléatoires réelles indépendantes et centrées et notons Fn=σ(X0,X1,,Xn)\mathcal{F}_n = \sigma(X_0, X_1, \ldots, X_n). Montrer que Sn=k=0nXkS_n = \sum_{k=0}^n X_k est une martingale par rapport à la filtration Fn\mathcal{F}_n.
الحل

1.

Pour a,bRa, b \in \mathbb{R}, E(aX+bYG)=aE(XG)+bE(YG)\mathbb{E}(aX + bY \mid \mathcal{G}) = a\mathbb{E}(X\mid\mathcal{G}) + b\mathbb{E}(Y\mid\mathcal{G}) p.s. (linéarité de l'espérance conditionnelle).

2.

Pour tout AGA \in \mathcal{G}, par indépendance E[1AX]=P(A)E[X]=E[1AE(X)]\mathbb{E}[\mathbf{1}_A X] = P(A)\mathbb{E}[X] = \mathbb{E}[\mathbf{1}_A \mathbb{E}(X)]. Comme E(X)\mathbb{E}(X) est G\mathcal{G}-mesurable (constante), E(XG)=E(X)\mathbb{E}(X\mid\mathcal{G}) = \mathbb{E}(X) p.s.

3.

Si XX est G\mathcal{G}-mesurable, elle vérifie la propriété caractéristique E[1AX]=E[1AX]\mathbb{E}[\mathbf{1}_A X] = \mathbb{E}[\mathbf{1}_A X] pour tout AGA \in \mathcal{G}, donc E(XG)=X\mathbb{E}(X\mid\mathcal{G}) = X p.s.

4.

SnS_n est Fn\mathcal{F}_n-adaptée et intégrable. De plus

E(Sn+1Fn)=Sn+E(Xn+1Fn)=Sn+E(Xn+1)=Sn,\mathbb{E}(S_{n+1}\mid\mathcal{F}_n) = S_n + \mathbb{E}(X_{n+1}\mid\mathcal{F}_n) = S_n + \mathbb{E}(X_{n+1}) = S_n,

car Xn+1X_{n+1} est indépendante de Fn\mathcal{F}_n et centrée. Donc (Sn)(S_n) est une martingale.

(Sn) est une Fn-martingale\boxed{(S_n) \text{ est une } \mathcal{F}_n\text{-martingale}}