Avec la fonction bernoulli(p) du cours, écrire une fonction estimation(N) qui simule N fois trois essais indépendants de paramètre 0,7 et renvoie la fréquence des séries comportant exactement deux succès. Quelle valeur attend-on pour N grand ?
Voir la correction
chaque série se simule par trois appels indépendants à bernoulli(0.7), dont la somme s est le nombre de succès ; on compte les séries où s vaut 2 :
def estimation(N):
favorables = 0
for repetition in range(N):
s = bernoulli(0.7) + bernoulli(0.7) + bernoulli(0.7)
if s == 2:
favorables = favorables + 1
return favorables / NValeur théorique : les chemins avec exactement deux succès sont SSE, SES et ESS, chacun de probabilité 0,72 × 0,3 = 0,147 ; la probabilité vaut 3 × 0,147 = 0,441.
Pour N grand, la fréquence renvoyée fluctue autour de 0,441 : par exemple 0,4378 ou 0,4452 pour N = 10 000. Deux exécutions donnent des résultats différents, et la fluctuation diminue quand N augmente. La simulation confirme le calcul mais ne le remplace pas : elle donne une valeur approchée, l'arbre donne la valeur exacte.