Discussion :

[pzorba75]

Bonjour,
j'écris des programmes Python, niveau lycée, en m'appuyant sur des sujets proposés dans les ouvrages scolaires.
Dans u programme de probabilité, niveau 1ère, on demande de simuler le nombre de filles dans des famille de deux enfants, sachant qu'en moyenne, il naît 105 garçons pour 100 filles (pas de jumeaux dans ces familles).
Le code est pré-mâché sur 4 ou 5 lignes, sans placer plusieurs instructions par ligne.
J'ai écrit la fonction ci-dessous :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
from random import randint
def famille():
    N=0
    if randint(1,205)<=100:
        N=N+1
    if randint(1,205)<=100:
        N=N+1
    return(N)

Python peut-il permettre une écriture plus concise?
J'ai cherché quelques exemples, sans succès.
Merci d'avance pour les conseils "Pythonic".

[wiztricks]

Salut,

Python peut-il permettre une écriture plus concise?

Est-ce que vous êtes satisfait par la traduction faite de ces probabilités?

- W

[pzorba75]

En faisant un grand nombre de simulations (99999) de la fonction famille() et en, calculant le nombre moyen de filles, la moyenne obtenue du nombre de filles est 0.9779997799978 pour des familles ayant 2 enfants, ce qui me convient tout à fait pour l'utilisation de randint(1,205) et simuler 100 naissances de filles sur 205 naissances.

[marco056]

Tu as 2 fois la même instruction, c'est normal ?

[pzorba75]

À mon sens oui, un tirage aléatoire pour la première naissance et un second tirage pour la deuxième et dernière naissance dans cette famille.

[Sve@r]

Bonjour

À mon sens oui, un tirage aléatoire pour la première naissance et un second tirage pour la deuxième et dernière naissance dans cette famille.

La question de marco06 était faite pour t'amener à penser "boucle".

Code c :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
def famille():
    N=0
    for i in range(2):
        if randint(1,205) <= 100:
            N=N+1
    return(N)

Bien entendu ce n'est pas pour remplacer la solution de wiztricks qui utilises les possibilités de random.choice, c'est pour te montrer la façon de faire "n" choses identiques dans un cadre plus général. Ecrire ça sous forme de boucle, c'est

• factoriser les instructions (le jour où les probas passent à 106 filles pour 102 garçons tu ne modifies qu'une seule ligne)
• offrir la possibilité facile de pouvoir passer à 3, 4, 5 naissances sans changer quoi que ce soit au traitement (juste le nombre d'itérations)

Et puisque tu aimes les truc "pythonic", tu peux utilises les listes en comprehension disponibles dans Python pour adapter la solution de tyrtamos => return tuple(int(randint(1,205) > 100) for i in range(2))

[wiztricks]

Salut,

Dans tous les cas, lorsque vous partez à écrire une fonction essayez de voir si c'est un problème classique.
La plupart du temps vous allez avoir de quoi faciliter votre travail dans les bibliothèques (ici, c'est random et les fonctions choices, sample,...)
Mais lorsque ce n'est pas trop compliqué à écrire, coder son truc va bien aussi.


- W

[tyrtamos]

Bonjour,


Petit problème amusant => petite proposition:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
from random import randint
 
def famille():
 
    enfant1 = int(randint(1,205)>100)
 
    enfant2 = int(randint(1,205)>100)
 
    return enfant1, enfant2

Explication:
Pour chaque famille de deux enfants, on tire 2 fois au hasard un garçon ou une fille. Pour ça, on tire un nombre entier avec randint entre 1 et 205 (bornes comprises). Si c'est >100, c'est un garçon (plage de 101 à 205), sinon, c'est une fille. (plage de 1 à 100) Avec int(), on convertit le booléen en entier (int(True) => 1 et int(False) => 0) . Donc, à chaque tirage on connaitra le sexe des 2 enfants sous forme d'un tuple de type (1, 0). La fonction sum permet d'avoir simplement le nombre de garçon de chaque famille.

On pourrait d'ailleurs faire le calcul en une seule ligne, mais il vaut mieux faciliter la lecture pour faciliter la compréhension:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

return int(randint(1,205)>100),  int(randint(1,205)>100)

Test: on fait 1 million de tirage (1 millions de familles de 2 enfants):

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
n = 1000000
ng = 0
for i in range(0, n):
    ng += sum(famille())
 
print(ng, ng/(n*2))   
1024524 0.512262

Sur 1 million de familles, on devrait avoir 2 millions d'enfants, dont 2000000/205*105=1024390 garçons et 2000000/205*100=975610 filles soit un rapport de garçons de 1024390/2000000=0.512195 très proche du résultat obtenu ici.