Le cours
Résumé en cinq points
- 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.
- La garantie d’approximation vient du théorème limite central : la demi-largeur est , au plus pour une probabilité, donc une précision en .
- On compare deux estimateurs en simulant jeux de données de taille , puis en comparant moyennes, variances et histogrammes des valeurs obtenues.
- Un intervalle de confiance se juge sur deux critères, sa demi-largeur moyenne et son niveau de confiance effectif mesuré par simulation.
- 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.
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
et sont indépendantes uniformes sur .
Montrer que écrire un programme estimant 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
- Le couple est uniforme sur le carré unité, et l'événement est le quart de disque, d'aire .
- Le programme s'écrit
U, V = rd.random(n), rd.random(n)puisprint(4*np.mean(U**2 + V**2 <= 1)). - La variable simulée a pour variance avec , d'où un écart type d'environ .
- La demi-largeur vaut environ . Pour qu'elle soit au plus , il faut , soit environ dix millions.
piège : Oublier le facteur dans la variance, ce qui sous-estime 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
et sont indépendantes uniformes sur .
Montrer que pour , puis estimer cette valeur en par simulation.
Voir le début du corrigé
- et 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é à 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é
- Bienaymé-Tchebychev donne . …
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
sont indépendantes de loi uniforme sur .
On pose et .
- Déterminer la loi de .
- vérifier que est sans biais.
- calculer et .
- puis proposer une simulation de contrôle.
Voir le début du corrigé
- Pour , par indépendance, d'où la densité sur .
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 pour .
Calculer son risque quadratique, trouver optimal, comparer à et à .
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 est un intervalle de confiance de niveau .
Écrire une fonction Python de paramètres , et qui renvoie les demi-largeurs et les niveaux effectifs des deux intervalles du cours à 95 pour cent.
Voir le début du corrigé
- 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 , de 5 à 200, pour puis .
Décrire et expliquer les courbes.
Voir le début du corrigé
- On appelle la fonction de l'exercice 6 pour chaque 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
sont indépendantes de loi avec connu.
Montrer que pour tout .
Simuler ensuite l'intervalle en remplaçant par l'écart type empirique, pour , et .
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 . Écrire 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é
- 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 avec uniforme sur .
Tracer les moyennes empiriques en fonction de et expliquer.
Voir le début du corrigé
- La fonction de répartition de la loi de Cauchy est , 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.