📚 الرئيسية

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

مسابقة تخصص · 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 (2ème version), Université Mohamed Khider - Biskra, Faculté des Sciences Exactes et des Sciences de la Nature et de la Vie, Département de Mathématiques, 17 octobre 2015, durée 01h30.

التمرين 1

Exercice 1 — Perturbation d'une densité : vérification, moments et interprétation

#density-function#normal-distribution#moments#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_{\varepsilon,b}(x) = (1-\varepsilon)f(x) + \frac{\varepsilon}{b}f\!\left(\frac{x}{b}\right),

0<ε0.050 \lt \varepsilon \leq 0.05 et b>1b \gt 1.

  1. Montrer que fε,b(x)f_{\varepsilon,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)\mathbf{E}(X) et Var(X)\mathbf{Var}(X), où Xfε,b(x)X \sim f_{\varepsilon,b}(x).
  3. Si b=3b = 3 et ε=0.05\varepsilon = 0.05, interpréter fε,b(x)f_{\varepsilon,b}(x) de la question 2 et la comparer avec f(x)f(x).
الحل

1.

fε,b(x)0f_{\varepsilon,b}(x) \geq 0 car f0f \geq 0, ε(0,1)\varepsilon \in (0,1). Pour l'intégrale : fε,b=(1ε)f+εbf(x/b)dx=(1ε)+ε=1\int f_{\varepsilon,b} = (1-\varepsilon)\int f + \frac{\varepsilon}{b}\int f(x/b)dx = (1-\varepsilon) + \varepsilon = 1 (changement de variable u=x/bu = x/b). ✓

2.

Si Xfε,bX \sim f_{\varepsilon,b}, alors E[X]=(1ε)0+εb0=0E[X] = (1-\varepsilon)\cdot 0 + \varepsilon \cdot b \cdot 0 = 0.

E[X2]=(1ε)1+εb21=1ε+εb2E[X^2] = (1-\varepsilon)\cdot 1 + \varepsilon \cdot b^2 \cdot 1 = 1 - \varepsilon + \varepsilon b^2.

Var(X)=1ε+εb2=1+ε(b21)\boxed{\text{Var}(X) = 1 - \varepsilon + \varepsilon b^2 = 1 + \varepsilon(b^2-1)}

3.

Pour b=3b = 3, ε=0.05\varepsilon = 0.05 : Var(X)=1+0.05×8=1.4\text{Var}(X) = 1 + 0.05 \times 8 = 1.4. La distribution est un mélange de N(0,1)\mathcal{N}(0,1) (avec poids 0.95) et de N(0,9)\mathcal{N}(0,9) (avec poids 0.05) — une distribution plus lourde de queue que la normale standard, modélisant des observations aberrantes rares.

التمرين 2

Exercice 2 — Moments empiriques de la moyenne d'un échantillon de taille 2

#empirical-moments#sample-mean#skewness#kurtosis

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

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

1.

X=(X1+X2)/2\overline{X} = (X_1+X_2)/2 avec X1,X2X_1, X_2 iid de loi XX. μ1(X)=0\mu_1(\overline{X}) = 0. μ2(X)=μ2/2\mu_2(\overline{X}) = \mu_2/2. μ3(X)=μ3/4\mu_3(\overline{X}) = \mu_3/4. μ4(X)=(μ4+3μ22)/43(μ2/2)2+3(μ2/2)2\mu_4(\overline{X}) = (\mu_4 + 3\mu_2^2)/4 - 3(\mu_2/2)^2 + 3(\mu_2/2)^2... En développant : μ4(X)=μ4/8+3μ22/8\mu_4(\overline{X}) = \mu_4/8 + 3\mu_2^2/8.

2.

γ1(X)=μ3(X)/[μ2(X)]3/2=(μ3/4)/(μ2/2)3/2=γ1/2\gamma_1(\overline{X}) = \mu_3(\overline{X})/[\mu_2(\overline{X})]^{3/2} = (\mu_3/4)/(\mu_2/2)^{3/2} = \gamma_1/\sqrt{2}. γ2(X)=γ2/2\gamma_2(\overline{X}) = \gamma_2/2.

3.

Mk=(X1k+X2k)/2M_k = (X_1^k + X_2^k)/2. Comme XX est symétrique, μ2j+1=0\mu_{2j+1} = 0. Cov(M1,M2)=14[E(X1X12)E(X1)E(X12)]=14(μ30)=0\text{Cov}(M_1, M_2) = \frac{1}{4}[E(X_1 \cdot X_1^2) - E(X_1)E(X_1^2)] = \frac{1}{4}(\mu_3 - 0) = \boxed{0} (car μ3=0\mu_3 = 0).

التمرين 3

Exercice 3 — Densité gamma, espérance, variance et intervalle de confiance

#gamma-distribution#confidence-interval#method-of-moments#central-limit-theorem

Soient X1,...,XnX_1, ..., 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 = \mathbf{E}(X_1) puis exprimer Var(X1)\mathbf{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.

On reconnaît X1Γ(2,λ)X_1 \sim \Gamma(2, \lambda), donc μ=E[X1]=2/λ\mu = E[X_1] = 2/\lambda et Var(X1)=2/λ2=μ2/2\text{Var}(X_1) = 2/\lambda^2 = \mu^2/2.

2.

Par le TCL, XˉnN(μ,σ2/n)\bar{X}_n \approx \mathcal{N}(\mu, \sigma^2/n) avec σ2=μ2/2\sigma^2 = \mu^2/2. Un IC approché à 95% est :

[Xˉn1.96Xˉn2n,  Xˉn+1.96Xˉn2n]\boxed{\left[\bar{X}_n - \frac{1.96\,\bar{X}_n}{\sqrt{2n}},\; \bar{X}_n + \frac{1.96\,\bar{X}_n}{\sqrt{2n}}\right]}

(en remplaçant σ=μ/2\sigma = \mu/\sqrt{2} par Xˉn/2\bar{X}_n/\sqrt{2}).