الرئيسية

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

مسابقة تخصص · 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 octobre 2015, durée 01h30.

التمرين 1

Exercice 1 — Variables aléatoires iid et lois dérivées

#probability#uniform-distribution#exponential-distribution#gamma-distribution#order-statistics

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.

Puisque FF est continue et strictement croissante, U=F(X)U[0,1]U = F(X) \sim \mathcal{U}[0,1]. Les F(Xi)F(X_i) sont iid par indépendance des XiX_i.

2.

Ei=logUiE_i = -\log U_i avec UiU[0,1]U_i \sim \mathcal{U}[0,1]. On a P(Ei>t)=P(Ui<et)=etP(E_i \gt t) = P(U_i \lt e^{-t}) = e^{-t} pour t0t \geq 0, donc EiE(1)E_i \sim \mathcal{E}(1).

3.

SiS_i est somme de ii variables iid E(1)\mathcal{E}(1), donc SiΓ(i,1)S_i \sim \Gamma(i, 1).

4.

Si Z1Γ(α,θ)Z_1 \sim \Gamma(\alpha,\theta) et Z2Γ(β,θ)Z_2 \sim \Gamma(\beta,\theta) indépendantes, alors

Y=Z1Z1+Z2Beta(α,β).Y = \frac{Z_1}{Z_1+Z_2} \sim \mathcal{B}eta(\alpha, \beta).

5.

Sn=Si+(SnSi)S_n = S_i + (S_n - S_i) avec SiΓ(i,1)S_i \sim \Gamma(i,1) et SnSiΓ(ni,1)S_n - S_i \sim \Gamma(n-i,1) indépendants. Par la question 4 :

Vi=Si/SnBeta(i,ni).\boxed{V_i = S_i/S_n \sim \mathcal{B}eta(i,\, n-i).}

التمرين 2

Exercice 2 — Densité de Rayleigh, estimateur du maximum de vraisemblance

#statistics#mle#rayleigh-distribution#unbiasedness#consistency

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} \dfrac{x}{\theta} \exp\!\left\{-\dfrac{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, l'estimateur θ^\hat{\theta} du maximum de vraisemblance 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.

Pour y>0y \gt 0 : P(Yy)=P(X2y)=1ey/θP(Y \leq y) = P(X \leq \sqrt{2y}) = 1 - e^{-y/\theta}. Donc

fY(y)=1θey/θ,y>0.f_Y(y) = \frac{1}{\theta}e^{-y/\theta},\quad y \gt 0.

C'est la loi exponentielle E(1/θ)\mathcal{E}(1/\theta).

2.

La log-vraisemblance (θ)=nlogθ12θxi2\ell(\theta) = -n\log\theta - \frac{1}{2\theta}\sum x_i^2. En annulant (θ)\ell'(\theta) :

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

3.

Yi=Xi2/2E(1/θ)Y_i = X_i^2/2 \sim \mathcal{E}(1/\theta), donc E(Yi)=θE(Y_i)=\theta et Var(Yi)=θ2\mathrm{Var}(Y_i)=\theta^2. Ainsi E(θ^)=θE(\hat{\theta})=\theta (sans biais) et Var(θ^)=θ2/n0\mathrm{Var}(\hat{\theta}) = \theta^2/n \to 0, donc θ^\hat{\theta} est convergent en MSE.

التمرين 3

Exercice 3 — Test d'hypothèse unilatéral, fonction puissance

#statistics#hypothesis-testing#power-function#normal-distribution#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 nulle 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.840.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)=Φ ⁣(n(2μ)1)\beta(\mu) = P_\mu(\bar{X} \lt 2 - 1/\sqrt{n}) = \Phi\!\left(\sqrt{n}(2-\mu) - 1\right).

β\beta est décroissante en μ\mu, tend vers 1 quand μ\mu \to -\infty et vers 0 quand μ+\mu \to +\infty.

2.

Risque de 1ère espèce =supμ2β(μ)=β(2)=Φ(1)=1Φ(1)=10.84=0.16= \sup_{\mu \geq 2}\beta(\mu) = \beta(2) = \Phi(-1) = 1 - \Phi(1) = 1 - 0.84 = 0.16.

3.

Le test est sans biais si β(μ)α\beta(\mu) \leq \alpha pour μH0\mu \in H_0 et β(μ)α\beta(\mu) \geq \alpha pour μH1\mu \in H_1. Puisque β\beta est décroissante, supH0β=β(2)=0.16\sup_{H_0}\beta = \beta(2) = 0.16 et infH1β=β(2)=0.16\inf_{H_1}\beta = \beta(2) = 0.16, donc

β(μ)0.16 pour tout μ<2: le test est sans biais.\boxed{\beta(\mu) \geq 0.16 \text{ pour tout } \mu \lt 2 : \text{ le test est sans biais.}}