الرئيسية

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

مسابقة تخصص · 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, 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 — Loi uniforme, exponentielle et loi gamma : sommes partielles

#probability#uniform-distribution#exponential-distribution#gamma-distribution#beta-distribution

Soit X1,,XnX_1, \ldots, X_n une suite de variables aléatoires indépendantes et identiquement distribuées (iid) de loi de probabilité continue FF.

  1. Montrer que F(X1),,F(Xn)F(X_1), \ldots, F(X_n) sont iid et uniformément distribuées sur [0,1][0, 1].
  2. Montrer que E1=logF(X1),,En=logF(Xn)E_1 = -\log F(X_1), \ldots, E_n = -\log F(X_n) sont iid exponentielles de paramètre 1.
  3. Quelle est, pour chaque ii, la loi de probabilité de la variable aléatoire Si=E1++EiS_i = E_1 + \ldots + E_i ?
  4. Soit Z1,Z2Z_1, Z_2 deux variables aléatoires indépendantes de lois gamma de paramètres (α,θ)(\alpha, \theta) et (β,θ)(\beta, \theta) respectivement. Indiquer (sans le démontrer) quelle est la loi de la variable aléatoire Y=Z1/(Z1+Z2)Y = Z_1/(Z_1 + Z_2) ?
  5. Déduire, de la question 4, la loi de probabilité de la variable aléatoire Vi=Si/SnV_i = S_i/S_n, i=1,,ni = 1, \ldots, n. Indication : écrire SnS_n sous forme de la somme de deux termes.
الحل

1.

Si XX a pour fonction de répartition FF continue, alors U=F(X)Unif(0,1)U = F(X) \sim \text{Unif}(0,1) (transformation intégrale de probabilité). Par iid, les F(Xi)F(X_i) sont iid Unif(0,1)\text{Unif}(0,1).

2.

Si UUnif(0,1)U \sim \text{Unif}(0,1), alors logUExp(1)-\log U \sim \text{Exp}(1). Donc les EiE_i sont iid Exp(1)\text{Exp}(1).

3.

Si=E1++EiS_i = E_1 + \cdots + E_i est somme de ii v.a. iid Exp(1)\text{Exp}(1).

SiGamma(i,1)\boxed{S_i \sim \text{Gamma}(i, 1)}

4.

Si Z1Gamma(α,θ)Z_1 \sim \text{Gamma}(\alpha, \theta) et Z2Gamma(β,θ)Z_2 \sim \text{Gamma}(\beta, \theta) indépendantes, alors

Y=Z1Z1+Z2Beta(α,β)\boxed{Y = \frac{Z_1}{Z_1+Z_2} \sim \text{Beta}(\alpha, \beta)}

5.

Sn=Si+(SnSi)S_n = S_i + (S_n - S_i)SiGamma(i,1)S_i \sim \text{Gamma}(i,1) et SnSiGamma(ni,1)S_n - S_i \sim \text{Gamma}(n-i, 1) sont indépendantes. Par la question 4 :

Vi=Si/SnBeta(i,ni)\boxed{V_i = S_i/S_n \sim \text{Beta}(i, n-i)}

التمرين 2

Exercice 2 — Loi de Rayleigh : EMV et propriétés

#statistics#rayleigh-distribution#maximum-likelihood#bias#convergence

Soit XX une variable aléatoire de densité de probabilité ff définie sur R\mathbb{R} par

f(x)={xθexp{x22θ}si x>00si x0;  θ>0.f(x) = \begin{cases} \frac{x}{\theta} \exp\left\{-\frac{x^2}{2\theta}\right\} & \text{si } x \gt 0 \\\\ 0 & \text{si } x \leq 0 \end{cases} \quad ; \; \theta \gt 0.
  1. Déterminer la fonction de répartition de Y=X2/2Y = X^2/2, puis déduire sa densité de probabilité. À quelle loi de probabilité usuelle correspond-elle ?
  2. Déterminer, à partir d'un échantillon (X1,,Xn)(X_1, \ldots, X_n), de taille n1n \geq 1, de XX, l'estimateur du maximum de vraisemblance θ^\hat{\theta} du paramètre θ\theta.
  3. Montrer que θ^\hat{\theta} est sans biais et convergent en moyenne quadratique (pour θ\theta). On rappelle que l'espérance et la variance d'une variable aléatoire exponentielle, de paramètre α>0\alpha \gt 0, sont respectivement égales à 1/α1/\alpha et 1/α21/\alpha^2.
الحل

1.

Y=X2/2Y = X^2/2. FY(y)=P(X2/2y)=P(X2y)=FX(2y)F_Y(y) = P(X^2/2 \leq y) = P(X \leq \sqrt{2y}) = F_X(\sqrt{2y}). fY(y)=1θey/θf_Y(y) = \frac{1}{\theta} e^{-y/\theta} pour y>0y \gt 0.

YExp(1/θ)\boxed{Y \sim \text{Exp}(1/\theta)}

2.

L(θ)=xiθexi2/(2θ)L(\theta) = \prod \frac{x_i}{\theta} e^{-x_i^2/(2\theta)}. (θ)=lnxinlnθ12θxi2\ell(\theta) = \sum \ln x_i - n\ln\theta - \frac{1}{2\theta}\sum x_i^2.

(θ)=n/θ+12θ2xi2=0\ell'(\theta) = -n/\theta + \frac{1}{2\theta^2}\sum x_i^2 = 0 donne

θ^=12ni=1nXi2=Yˉ\boxed{\hat{\theta} = \frac{1}{2n}\sum_{i=1}^n X_i^2 = \bar{Y}}

3.

θ^=1nYi\hat{\theta} = \frac{1}{n}\sum Y_i avec Yi=Xi2/2Exp(1/θ)Y_i = X_i^2/2 \sim \text{Exp}(1/\theta). E(θ^)=θ\mathbb{E}(\hat{\theta}) = \theta (sans biais). Var(θ^)=θ2/n0\text{Var}(\hat{\theta}) = \theta^2/n \to 0. Donc convergent en moyenne quadratique.

التمرين 3

Exercice 3 — Test de puissance et risque de première espèce

#statistics#hypothesis-testing#power-function#type-I-error

Soit (X1,,Xn)(X_1, \ldots, X_n) un échantillon, de taille n1n \geq 1, d'une population normale d'espérance μ\mu et de variance 1. Pour trancher entre l'hypothèse H0:μ2H_0 : \mu \geq 2 contre l'hypothèse alternative H1:μ<2H_1 : \mu \lt 2, on propose le test dont la région de rejet est W={(x1,,xn)Rn:xˉ<21/n}W = \{(x_1, \ldots, x_n) \in \mathbb{R}^n : \bar{x} \lt 2 - 1/\sqrt{n}\}.

  1. Déterminer la fonction puissance de ce test. Étudier ses variations puis tracer (approximativement) son graphe.
  2. Sachant que la valeur de la fonction de répartition de la loi normale standard en 1 est égale à 0.84, déterminer le risque de première espèce de ce test.
  3. Vérifier que ce test est sans biais.
الحل

1.

π(μ)=Pμ(Xˉ<21/n)=P(Xˉμ1/n<21/nμ1/n)=Φ(n(21/nμ)n)\pi(\mu) = P_{\mu}(\bar{X} \lt 2 - 1/\sqrt{n}) = P\left(\frac{\bar{X}-\mu}{1/\sqrt{n}} \lt \frac{2-1/\sqrt{n}-\mu}{1/\sqrt{n}}\right) = \Phi(\sqrt{n}(2-1/\sqrt{n}-\mu) \cdot \sqrt{n}).

Simplification : π(μ)=Φ((2μ)n1)\pi(\mu) = \Phi((2-\mu)\sqrt{n} - 1).

π\pi est décroissante en μ\mu (car Φ\Phi est croissante et l'argument est décroissant).

2.

Le risque de première espèce est α=supμ2π(μ)=π(2)=Φ(1)=1Φ(1)=10.84=0.16\alpha = \sup_{\mu \geq 2} \pi(\mu) = \pi(2) = \Phi(-1) = 1 - \Phi(1) = 1 - 0.84 = 0.16.

α=0.16\boxed{\alpha = 0.16}

3.

Pour μ<2\mu \lt 2 (sous H1H_1) : (2μ)n1>1(2-\mu)\sqrt{n} - 1 \gt -1, donc π(μ)>Φ(1)=α\pi(\mu) \gt \Phi(-1) = \alpha. Le test est sans biais car la puissance sous H1H_1 est toujours >α\gt \alpha.