📚 الرئيسية

مسابقة دكتوراه 2019Université Badji Mokhtar - Annaba — الموضوع 02

مسابقة تخصص · Probabilités & Statistiques

Concours d'entrée en formation doctorale 2018–2019 — Épreuve de statistiques, Université Badji-Mokhtar - Annaba, Département de Mathématiques — Durée : 02h.

التمرين 1

Exercice 1 — Loi de Pareto : fiabilité, estimation et efficacité

#pareto-distribution#maximum-likelihood#gamma-distribution#bias#reliability

On rappelle que la fonction gamma est définie par : Γ(x)=0tx1etdt\Gamma(x) = \int_0^{\infty} t^{x-1} e^{-t} dt et que Γ(n)=(n1)!\Gamma(n) = (n-1)!.

La densité d'une loi Gamma γ(θ,α)\gamma(\theta, \alpha) est f(x)=θαΓ(α)xα1exp(θx)f(x) = \frac{\theta^\alpha}{\Gamma(\alpha)} x^{\alpha-1} \exp(-\theta x) pour x0x \geq 0.

Sa fonction caractéristique est 1(1itθ)α\frac{1}{\left(1 - \frac{it}{\theta}\right)^\alpha}.


On considère un élément dont la durée de vie est une variable aléatoire TT de loi de Pareto (notée par la suite TP(a)T \to P(a)) et définie par sa densité :

f(t)=ata+1si t1,a>0f(t) = \frac{a}{t^{a+1}} \quad \text{si } t \geq 1, \quad a > 0

aa est un réel strictement positif inconnu.

  1. Déterminer la fonction de fiabilité et le taux de panne d'un tel élément.
  2. On souhaite estimer aa, pour cela on observe nn éléments identiques de durées de vie respectives T1,T2,,TnT_1, T_2, \ldots, T_n. Montrer que l'estimateur par maximum de vraisemblance de aa s'écrit :

a^=nk=1nlntk\hat{a} = \frac{n}{\sum_{k=1}^{n} \ln t_k}

  1. Montrer que la variable Z=lnTZ = \ln T suit une loi gamma dont on déterminera les paramètres. En déduire que la loi de la variable aléatoire U=k=1nZkU = \sum_{k=1}^{n} Z_k est une loi Γ(a,n)\Gamma(a, n).
  2. Calculer l'espérance et la variance de la variable 1a^\frac{1}{\hat{a}}.
  3. Déterminer le biais de l'estimateur a^\hat{a}. En déduire un estimateur non biaisé de aa noté a~\tilde{a} et déterminer sa variance.
  4. L'estimateur a~\tilde{a} est-il l'estimateur efficace de aa ?
  5. On monte un système SS dont les éléments identiques T1,T2,,TnT_1, T_2, \ldots, T_n sont en série. Calculer la fiabilité RS(t)R_S(t) du système SS.
الحل

1.

La fonction de répartition pour t1t \geq 1 :

F(t)=11taF(t) = 1 - \frac{1}{t^a}

La fiabilité :

R(t)=1F(t)=1taR(t) = 1 - F(t) = \frac{1}{t^a}

Le taux de panne :

h(t)=f(t)R(t)=ath(t) = \frac{f(t)}{R(t)} = \frac{a}{t}

2.

La vraisemblance : L(a)=k=1natka+1=antk(a+1)L(a) = \prod_{k=1}^n \frac{a}{t_k^{a+1}} = a^n \prod t_k^{-(a+1)}. Log-vraisemblance : (a)=nlna(a+1)lntk\ell(a) = n \ln a - (a+1)\sum \ln t_k. En annulant la dérivée :

a^=nk=1nlntk\boxed{\hat{a} = \frac{n}{\sum_{k=1}^n \ln t_k}}

3.

Posons Z=lnTZ = \ln T. Pour z0z \geq 0, P(Zz)=P(Tez)=1eazP(Z \leq z) = P(T \leq e^z) = 1 - e^{-az}. Donc ZExp(a)=γ(a,1)Z \sim \text{Exp}(a) = \gamma(a,1). Par propriété de la loi gamma, U=ZkΓ(a,n)U = \sum Z_k \sim \Gamma(a, n).

4.

On a 1a^=Un\frac{1}{\hat{a}} = \frac{U}{n} avec UΓ(a,n)U \sim \Gamma(a,n). Donc :

E(1a^)=1a,V(1a^)=1na2E\left(\frac{1}{\hat{a}}\right) = \frac{1}{a}, \quad V\left(\frac{1}{\hat{a}}\right) = \frac{1}{na^2}

5.

E(a^)=nan1E(\hat{a}) = \frac{na}{n-1}, biais =an1= \frac{a}{n-1}. Estimateur non biaisé :

a~=n1na^=n1lntk\boxed{\tilde{a} = \frac{n-1}{n}\hat{a} = \frac{n-1}{\sum \ln t_k}}

V(a~)=a2n2V(\tilde{a}) = \frac{a^2}{n-2} pour n>2n > 2.

6.

L'information de Fisher est I(a)=na2I(a) = \frac{n}{a^2}. La borne de Cramér-Rao est a2n\frac{a^2}{n}. Comme V(a~)=a2n2>a2nV(\tilde{a}) = \frac{a^2}{n-2} > \frac{a^2}{n}, a~\tilde{a} n'est pas efficace.

7.

Système série :

RS(t)=(1ta)n=1tna\boxed{R_S(t) = \left(\frac{1}{t^a}\right)^n = \frac{1}{t^{na}}}

التمرين 2

Exercice 2 — Test d'ajustement du khi-deux pour la loi binomiale

#binomial-distribution#chi-square-test#goodness-of-fit#genetics

Au cours d'essais génétiques sur la vigne, on a vu apparaître un caractère nouveau décelable uniquement à la floraison. Empiriquement, on a établi que, dans la souche ainsi créée, la fréquence d'apparition du caractère était de 0,150{,}15. Dans le but de tester les conditions de multiplication des individus présentant ce caractère, on répartit des échantillons de 5 individus dans des terrains et des terroirs différents.

  1. On considère la variable aléatoire XX donnant le nombre d'apparition du caractère étudié dans un lot de 5 individus. Quelle est la distribution de probabilité de XX ?
  2. On étudie 1000 échantillons de 5 individus. À la floraison, on détermine dans chaque échantillon le nombre d'individus présentant le caractère, on obtient les résultats suivants :
nombre d'individus ayant le caractère012345
nombre d'échantillons3994011504631

Étudier la conformité de ces résultats à la loi théorique déterminée dans la question 1).

الحل

1.

XX compte le nombre de succès (apparition du caractère) parmi 5 individus avec p=0,15p = 0{,}15, donc :

XB(5,  0,15)\boxed{X \sim \mathcal{B}(5,\; 0{,}15)}

2.

Probabilités théoriques : P(X=k)=(5k)(0,15)k(0,85)5kP(X=k) = \binom{5}{k}(0{,}15)^k(0{,}85)^{5-k}.

Effectifs théoriques (sur 1000) : 443,7443{,}7 ; 391,5391{,}5 ; 138,2138{,}2 ; 24,424{,}4 ; 2,22{,}2. On regroupe les classes k3k \geq 3.

On calcule χobs2=(OiEi)2Ei\chi^2_{\text{obs}} = \sum \frac{(O_i - E_i)^2}{E_i} et on compare à χ0,05,22=5,99\chi^2_{0{,}05,\, 2} = 5{,}99. L'ajustement est acceptable si χobs2<5,99\chi^2_{\text{obs}} < 5{,}99.

التمرين 3

Exercice 3 — Comparaison de trois files d'attente

#queueing-theory#poisson-process#exponential-distribution#performance-analysis

Dans une entreprise de service, où les usagers arrivent suivant un processus de Poisson de taux λ>0\lambda > 0, la direction a le choix entre 3 modèles suivants :

  1. Installation d'un serveur dont la durée suit une loi exponentielle de paramètre 2μ2\mu.
  2. Installation de 2 serveurs séparés et identiques, chacun admettant une durée de service suivant une loi exponentielle de paramètre μ\mu. Tout usager venant d'arriver, choisit avec la probabilité 0,50{,}5 l'un des serveurs et il n'y a pas d'échange d'usagers potentiels chez les serveurs.
  3. Installation de deux serveurs identiques non séparés. Dès que l'un des serveurs est libre, s'il y a des usagers en attente, le premier arrivé se fait servir chez le serveur libéré.

Supposons que λ2μ<1\frac{\lambda}{2\mu} < 1. Calculer le temps moyen passé dans le système (pour les trois systèmes en question). Lequel des 3 systèmes est le plus performant ?

الحل

Modèle 1 : M/M/1 avec taux 2μ2\mu

ρ1=λ2μ\rho_1 = \frac{\lambda}{2\mu}. Temps moyen :

W1=12μλ\boxed{W_1 = \frac{1}{2\mu - \lambda}}

Modèle 2 : Deux files M/M/1 indépendantes

Chaque serveur reçoit λ/2\lambda/2, taux μ\mu. Temps moyen :

W2=1μλ/2=22μλ\boxed{W_2 = \frac{1}{\mu - \lambda/2} = \frac{2}{2\mu - \lambda}}

Modèle 3 : File M/M/2

File M/M/2 avec taux λ\lambda, service μ\mu par serveur. Le temps moyen est inférieur à W1W_1.

Comparaison

W3<W1<W2W_3 < W_1 < W_2. Le modèle 3 (serveurs non séparés) est le plus performant.