Aller au contenu
Télécharger l'application

Accueil › ECG › Maths approfondies

ECG2 · Maths approfondies · semaine 33

Travaux pratiques, estimation ponctuelle et par intervalle de confiance

Méthode de Monte-Carlo, principe et garanties d'approximation. Comparaison de différents estimateurs ponctuels d'un paramètre. Comparaison des intervalles de confiance obtenus par différentes méthodes, durée indicative de 6 heures.

Sommaire

Le cours

Résumé en cinq points

  1. Monte-Carlo estime une espérance par une moyenne empirique de simulations, et une probabilité par une proportion, en s’appuyant sur la loi faible des grands nombres.
  2. La garantie d’approximation vient du théorème limite central : la demi-largeur est uσ/nu\sigma/\sqrt{n}, au plus u/(2n)u/(2\sqrt{n}) pour une probabilité, donc une précision en 1/n1/\sqrt{n}.
  3. On compare deux estimateurs en simulant NN jeux de données de taille nn, puis en comparant moyennes, variances et histogrammes des valeurs obtenues.
  4. Un intervalle de confiance se juge sur deux critères, sa demi-largeur moyenne et son niveau de confiance effectif mesuré par simulation.
  5. Pour une loi de Bernoulli, l’intervalle de Bienaymé-Tchebychev est très prudent et large ; l’intervalle asymptotique est plus court et atteint son niveau annoncé dès quelques dizaines d’observations, ce que la simulation permet de constater sans le démontrer.
Pourquoi ce chapitre

Le cours a établi des résultats asymptotiques :

  • loi faible des grands nombres ;
  • théorème limite central ;
  • convergence d’estimateurs ;
  • intervalles de confiance.

Ces énoncés portent sur des limites, ils ne disent rien sur ce qui se passe pour une taille d’échantillon donnée.

L’ordinateur répond à cette question. On peut rejouer dix mille fois la même expérience statistique :

  • observer la loi d’un estimateur ;
  • mesurer le taux réel de recouvrement d’un intervalle de confiance ;
  • comparer deux méthodes concurrentes sur des critères chiffrés.

Le thème comporte trois volets :

  • la méthode de Monte-Carlo, qui utilise la simulation pour approcher une quantité inconnue ;
  • la comparaison d’estimateurs ponctuels ;
  • la comparaison d’intervalles de confiance obtenus par différentes méthodes.

Fin de l'aperçu gratuit

Encore 11 parties de cours et 9 corrigés.

Essayer 7 jours gratuitement Sans engagement · cours, exercices et fiches complets

Méthodes et exercices

Exercices classiques corrigés

Les exercices se suivent : garde leur ordre. Difficulté et durée indicatives.

Exercice 1. Approcher π ★★

Classiques de khôlle et concours · 25 min · Niveau khôlle · Estimation par simulation

UU et VV sont indépendantes uniformes sur [0,1][0, 1].

Montrer que P(U2+V2≤1)=π4,P(U^2 + V^2 \leq 1) = \frac{\pi}{4}, écrire un programme estimant π\pi avec sa garantie à 95 pour cent, puis estimer le nombre de simulations pour trois décimales.

Voir l'indice

Technique indicatrice et garantie de Monte-Carlo.

Voir le corrigé offert
  1. Le couple (U,V)(U, V) est uniforme sur le carré unité, et l'événement est le quart de disque, d'aire π/4\pi/4.
  2. Le programme s'écrit U, V = rd.random(n), rd.random(n) puis print(4*np.mean(U**2 + V**2 <= 1)).
  3. La variable simulée 4⋅1A4\cdot\mathbb{1}_A a pour variance 16p(1−p)16p(1 - p) avec p=π/4p = \pi/4, d'où un écart type d'environ 1,641{,}64.
  4. La demi-largeur vaut environ 3,22/n3{,}22/\sqrt n. Pour qu'elle soit au plus 0,0010{,}001, il faut n≥32202n \geq 3220^2, soit environ dix millions.

piège : Oublier le facteur 44 dans la variance, ce qui sous-estime nn d'un facteur seize.

Exercice 2. Fonction de répartition d'un produit ★★

Classiques de khôlle et concours · 30 min · Niveau khôlle · Estimation par simulation

UU et VV sont indépendantes uniformes sur [0,1][0, 1].

Montrer que P(UV≤t)=t−tln⁡tP(UV \leq t) = t - t \ln t pour t∈]0,1]t \in ]0, 1], puis estimer cette valeur en t=12t = \frac{1}{2} par simulation.

Voir le début du corrigé
  1. UU et VV sont presque sûrement strictement positives. …

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 3. Tchebychev contre théorème limite central ★★

Classiques de khôlle et concours · 20 min · Niveau khôlle · Dimensionnement des simulations

Combien de simulations faut-il pour estimer une probabilité pp à 0,010{,}01 près avec un niveau d'au moins 95 pour cent, selon Bienaymé-Tchebychev puis selon l'approche asymptotique ?

Voir le début du corrigé
  1. Bienaymé-Tchebychev donne P(∣Xˉn−p∣≥ε)≤1/(4nε2)P(|\bar X_n - p| \geq \varepsilon) \leq 1/(4n\varepsilon^2). …

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 4. Deux estimateurs pour la loi uniforme ★★

Classiques de khôlle et concours · 45 min · Niveau écrit Ecricome · Comparaison d'estimateurs

X1,…,XnX_1, \dots, X_n sont indépendantes de loi uniforme sur [0,θ][0, \theta].

On pose Tn=2X‾nT_n = 2\overline{X}_n et Mn=max⁡XiM_n = \max X_i.

  • Déterminer la loi de MnM_n.
  • vérifier que Un=(n+1)MnnU_n = \frac{(n + 1) M_n}{n} est sans biais.
  • calculer V(Tn)V(T_n) et V(Un)V(U_n).
  • puis proposer une simulation de contrôle.
Voir le début du corrigé
  1. Pour t∈[0,θ]t \in [0, \theta], P(Mn≤t)=(t/θ)nP(M_n \leq t) = (t/\theta)^n par indépendance, d'où la densité ntn−1/θnn t^{n-1}/\theta^n sur [0,θ][0, \theta].

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 5. Un estimateur biaisé meilleur ★★

Classiques de khôlle et concours · 30 min · Niveau écrit HEC ou ESSEC · Comparaison d'estimateurs

Avec les notations de l'exercice 4, on considère cMncM_n pour c>0c > 0.

Calculer son risque quadratique, trouver cc optimal, comparer à UnU_n et à MnM_n.

Corrigé détaillé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 6. Deux intervalles pour une loi de Bernoulli ★★

Classiques de khôlle et concours · 30 min · Niveau écrit EDHEC · Intervalles de confiance

Démontrer que X‾n±12nα\overline{X}_n \pm \frac{1}{2\sqrt{n\alpha}} est un intervalle de confiance de niveau 1−α1 - \alpha.

Écrire une fonction Python de paramètres nn, pp et NN qui renvoie les demi-largeurs et les niveaux effectifs des deux intervalles du cours à 95 pour cent.

Voir le début du corrigé
  1. On a E(X̄ₙ) = p et V(X̄ₙ) = p(1 − p)/n ≤ 1/(4n). …

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 7. À partir de quel rang ? ★★

Classiques de khôlle et concours · 25 min · Niveau écrit EDHEC · Intervalles de confiance

Tracer le niveau effectif de l'intervalle asymptotique en fonction de nn, de 5 à 200, pour p=0,5p = 0{,}5 puis p=0,05p = 0{,}05.

Décrire et expliquer les courbes.

Voir le début du corrigé
  1. On appelle la fonction de l'exercice 6 pour chaque nn et on trace la liste des niveaux avec plt.plot.

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 8. Loi normale, σ connu puis estimé ★★

Classiques de khôlle et concours · 20 min · Niveau écrit Ecricome · Intervalles de confiance

X1,…,XnX_1, \dots, X_n sont indépendantes de loi N(m,σ2)\mathcal{N}(m, \sigma^2) avec σ\sigma connu.

Montrer que P(∣X‾n−m∣≤1,96 σn)=2Φ(1,96)−1≈0,95\begin{aligned} &P\left(|\overline{X}_n - m| \leq 1{,}96\,\frac{\sigma}{\sqrt{n}}\right) \\ &= 2\Phi(1{,}96) - 1 \\ &\approx 0{,}95 \end{aligned} pour tout nn.

Simuler ensuite l'intervalle en remplaçant σ\sigma par l'écart type empirique, pour n=5n = 5, 1010 et 5050.

Corrigé détaillé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 9. Une espérance avec garantie ★★

Classiques de khôlle et concours · 30 min · Niveau écrit HEC ou ESSEC · Estimation par simulation

Soit I=∫0+∞e−x1+x dxI = \int_0^{+\infty} \frac{e^{-x}}{1 + x}\,dx. Écrire II comme une espérance, proposer un estimateur de Monte-Carlo, montrer sa convergence et donner une garantie sans calculer de variance.

Voir le début du corrigé
  1. Si X suit la loi exponentielle de paramètre 1, le théorème de transfert donne I = E(Y) avec Y = 1/(1 + X). …

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 10. Quand Monte-Carlo s'effondre ★★

Classiques de khôlle et concours · 20 min · Type oral HEC · Limites du Monte-Carlo

Simuler une variable de Cauchy par tan⁡(π(U−12))\tan\left(\pi\left(U - \frac{1}{2}\right)\right) avec UU uniforme sur ]0,1[]0, 1[.

Tracer les moyennes empiriques en fonction de nn et expliquer.

Voir le début du corrigé
  1. La fonction de répartition de la loi de Cauchy est F(x)=1/2+arctan⁡(x)/πF(x) = 1/2 + \arctan(x)/\pi, et sa réciproque donne la formule.

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Fiche de cours rédigée à partir du programme officiel de la voie EC, mathématiques approfondies de seconde année, partie « Travaux pratiques de mathématiques avec Python », thème « Estimation ponctuelle et par intervalle de confiance ». Durée indicative de six heures.