الرئيسية

مسابقة دكتوراه 2016Université Mohamed Khider - Biskra — الموضوع 13

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

Concours d'accès à la formation de 3ème Cycle LMD — Mathématiques Appliquées (Option Probabilités et Statistique) — Épreuve de Statistique (variante 2), Université Mohamed Khider - Biskra, Faculté des Sciences Exactes et des Sciences de la Nature et de la Vie, Département de Mathématiques — 17/10/2015 — Durée 01h30.

التمرين 1

Exercice 1 — Perturbation d'une densité : mélange gaussien

#statistics#mixture-distribution#gaussian#expectation#variance

Soit f(x)f(x) une fonction de densité définie sur R\mathbb{R}. On définit une perturbation de f(x)f(x) par

fϵ,b(x)=(1ϵ)f(x)+ϵbf(xb),f_{\epsilon,b}(x) = (1 - \epsilon)f(x) + \frac{\epsilon}{b} f\left(\frac{x}{b}\right),

0<ϵ0,050 \lt \epsilon \leq 0{,}05 et b>1b \gt 1.

  1. Montrer que fϵ,b(x)f_{\epsilon,b}(x) est une fonction de densité.
  2. Si f(x)f(x) est la densité d'une distribution normale d'espérance 0 et de variance 1, calculer E(X)\mathbb{E}(X) et Var(X)\text{Var}(X), où Xfϵ,b(x)X \sim f_{\epsilon,b}(x).
  3. Si b=3b = 3 et ϵ=0,05\epsilon = 0{,}05, interpréter fϵ,b(x)f_{\epsilon,b}(x) de la question 2 et la comparer avec f(x)f(x).
الحل

1.

fϵ,b(x)0f_{\epsilon,b}(x) \geq 0 car combinaison convexe de densités positives. fϵ,bdx=(1ϵ)fdx+ϵbf(x/b)dx=(1ϵ)+ϵ=1\int f_{\epsilon,b} dx = (1-\epsilon)\int f dx + \frac{\epsilon}{b}\int f(x/b)dx = (1-\epsilon) + \epsilon = 1.

2.

E(X)=(1ϵ)0+ϵ0=0\mathbb{E}(X) = (1-\epsilon) \cdot 0 + \epsilon \cdot 0 = 0 (par symétrie des deux composantes).

E(X2)=(1ϵ)1+ϵb2=1ϵ+ϵb2\mathbb{E}(X^2) = (1-\epsilon) \cdot 1 + \epsilon \cdot b^2 = 1 - \epsilon + \epsilon b^2. Donc Var(X)=1+ϵ(b21)\text{Var}(X) = 1 + \epsilon(b^2 - 1).

E(X)=0,Var(X)=1+ϵ(b21)\boxed{\mathbb{E}(X) = 0, \quad \text{Var}(X) = 1 + \epsilon(b^2-1)}

3.

Avec b=3,ϵ=0,05b=3, \epsilon=0{,}05 : Var(X)=1+0,05×8=1,4\text{Var}(X) = 1 + 0{,}05 \times 8 = 1{,}4. C'est un mélange de 95% d'une N(0,1)\mathcal{N}(0,1) et 5% d'une N(0,9)\mathcal{N}(0,9). La contamination par la composante à grande variance crée des queues plus épaisses que la gaussienne standard.

التمرين 2

Exercice 2 — Moments de la moyenne empirique et coefficients d'asymétrie

#statistics#moments#skewness#kurtosis#sample-mean

Soit Xˉ\bar{X} la moyenne empirique d'un échantillon (X1,X2)(X_1, X_2) d'une variable aléatoire XX centrée de kk-ème moment μk\mu_k (k1k \geq 1). On désigne par μk(Xˉ)\mu_k(\bar{X}) les moments d'ordre k1k \geq 1 de Xˉ\bar{X}.

  1. Exprimer μk(Xˉ)\mu_k(\bar{X}) pour k=1,,4k = 1, \ldots, 4 en termes des μk\mu_k.
  2. Exprimer les coefficients de dissymétrie γ1(Xˉ)\gamma_1(\bar{X}) et d'aplatissement γ2(Xˉ)\gamma_2(\bar{X}) de Xˉ\bar{X} en termes de ceux de XX.
  3. On suppose maintenant que XX est symétrique et on désigne par MkM_k le moment empirique d'ordre kk (k1k \geq 1). Calculer Cov(M1,M2)\text{Cov}(M_1, M_2).
الحل

1.

Xˉ=(X1+X2)/2\bar{X} = (X_1+X_2)/2. μ1(Xˉ)=0\mu_1(\bar{X}) = 0. μ2(Xˉ)=μ2/2\mu_2(\bar{X}) = \mu_2/2. μ3(Xˉ)=μ3/4\mu_3(\bar{X}) = \mu_3/4. μ4(Xˉ)=(μ4+3μ22)/8+μ223/4\mu_4(\bar{X}) = (\mu_4 + 3\mu_2^2)/8 + \mu_2^2 \cdot 3/4... Par calcul direct des moments centrés.

2.

γ1(Xˉ)=μ3(Xˉ)/μ2(Xˉ)3/2=γ1(X)/2\gamma_1(\bar{X}) = \mu_3(\bar{X})/\mu_2(\bar{X})^{3/2} = \gamma_1(X)/\sqrt{2}.

γ2(Xˉ)=μ4(Xˉ)/μ2(Xˉ)23\gamma_2(\bar{X}) = \mu_4(\bar{X})/\mu_2(\bar{X})^2 - 3. Par calcul, γ2(Xˉ)=γ2(X)/2\gamma_2(\bar{X}) = \gamma_2(X)/2.

3.

Si XX est symétrique, μ3=0\mu_3 = 0. M1=XˉM_1 = \bar{X}, M2=12Xi2M_2 = \frac{1}{2}\sum X_i^2.

Cov(M1,M2)=E[M1M2]E[M1]E[M2]\text{Cov}(M_1, M_2) = E[M_1 M_2] - E[M_1]E[M_2]. Par symétrie et calcul : Cov(M1,M2)=μ3/n=0\text{Cov}(M_1, M_2) = \mu_3/n = 0.

Cov(M1,M2)=0\boxed{\text{Cov}(M_1, M_2) = 0}

التمرين 3

Exercice 3 — Loi Gamma(2,λ) : EMV et intervalle de confiance

#statistics#gamma-distribution#maximum-likelihood#confidence-interval

Soient X1,,XnX_1, \ldots, X_n une suite de variables aléatoires indépendantes et identiquement distribuées de densité

fλ(x)={λ2xexp(λx)si x00si x<0,  λ>0.f_\lambda(x) = \begin{cases} \lambda^2 x \exp(-\lambda x) & \text{si } x \geq 0 \\\\ 0 & \text{si } x \lt 0 \end{cases} \quad , \; \lambda \gt 0.
  1. Calculer μ=E(X1)\mu = \mathbb{E}(X_1) puis exprimer Var(X1)\text{Var}(X_1) en fonction de μ\mu.
  2. Construire un intervalle de confiance approché au niveau 95% pour μ\mu basé sur i=1nXi\sum_{i=1}^n X_i.
الحل

1.

XGamma(2,λ)X \sim \text{Gamma}(2, \lambda). μ=E(X)=2/λ\mu = \mathbb{E}(X) = 2/\lambda. Var(X)=2/λ2=μ2/2\text{Var}(X) = 2/\lambda^2 = \mu^2/2.

μ=2/λ,Var(X)=μ2/2\boxed{\mu = 2/\lambda, \quad \text{Var}(X) = \mu^2/2}

2.

Par le TCL, XˉN(μ,μ2/(2n))\bar{X} \approx \mathcal{N}(\mu, \mu^2/(2n)). On remplace μ\mu par Xˉ\bar{X} dans la variance : IC95%=Xˉ±1,96Xˉ2n\text{IC}_{95\%} = \bar{X} \pm 1{,}96 \cdot \frac{\bar{X}}{\sqrt{2n}}.

IC=[Xˉ(11,962n),  Xˉ(1+1,962n)]\boxed{\text{IC} = \left[\bar{X}\left(1 - \frac{1{,}96}{\sqrt{2n}}\right), \; \bar{X}\left(1 + \frac{1{,}96}{\sqrt{2n}}\right)\right]}