Guide de Distribution Binomiale en R : Un Guide pour Débutants

Salut à toi, futur·e programmeur·euse en R ! Aujourd'hui, nous allons plonger dans le monde fascinant des distributions binomiales en R. Ne t'inquiète pas si tu n'as jamais codé auparavant - je serai ton guide amical tout au long de ce voyage. À la fin de ce tutoriel, tu manipuleras les distributions binomiales comme un·e pro !

R - Binomial Distribution

Qu'est-ce qu'une Distribution Binomiale ?

Avant de nous lancer dans le code, comprenons ce qu'est une distribution binomiale. Imagine que tu lances une pièce de monnaie 10 fois. La distribution binomiale t'aiderait à calculer la probabilité d'obtenir un certain nombre de faces. Il s'agit de calculer la probabilité de succès dans un nombre fixe de trials indépendantes.

Maintenant, explorons les quatre principales fonctions que R nous offre pour travailler avec les distributions binomiales.

Les Quatre Fantastiques de la Distribution Binomiale en R

R nous offre quatre fonctions puissantes pour travailler avec les distributions binomiales. Rencontrons-les :

Fonction But
dbinom() Calcule la densité de probabilité
pbinom() Calcule la probabilité cumulative
qbinom() Calcule le quantile
rbinom() Génére des nombres aléatoires

Explorons chacun de ces détails.

dbinom() : Le Maître de la Probabilité

Qu'est-ce que dbinom() fait ?

La fonction dbinom() calcule la probabilité d'obtenir exactement k succès dans n trials. C'est comme demander : " Quelle est la chance d'obtenir exactement 3 faces lorsque je lance une pièce de monnaie 5 fois ?"

Syntaxe et Paramètres

dbinom(x, size, prob)
  • x : le nombre de succès que nous intéressons
  • size : le nombre de trials
  • prob : la probabilité de succès à chaque trial

Exemple : Lancers de Pièce

Calculons la probabilité d'obtenir exactement 3 faces en 5 lancers de pièce :

probability <- dbinom(3, size = 5, prob = 0.5)
print(probability)

Lorsque vous exécutez ce code, vous verrez :

[1] 0.3125

Cela signifie qu'il y a une chance de 31.25% d'obtenir exactement 3 faces en 5 lancers de pièce.

pbinom() : Le Calculateur de Probabilité Cumulative

Qu'est-ce que pbinom() fait ?

La fonction pbinom() calcule la probabilité cumulative - la probabilité d'obtenir jusqu'à k succès dans n trials. C'est comme demander : " Quelle est la chance d'obtenir 3 faces ou moins lorsque je lance une pièce de monnaie 5 fois ?"

Syntaxe et Paramètres

pbinom(q, size, prob, lower.tail = TRUE)
  • q : le nombre de succès que nous intéressons
  • size : le nombre de trials
  • prob : la probabilité de succès à chaque trial
  • lower.tail : si TRUE (par défaut), calcule P(X ≤ x) ; si FALSE, calcule P(X > x)

Exemple : Scores d'Examen

Imaginons un examen à choix multiples avec 10 questions, chacune avec 4 options. Quelle est la probabilité d'obtenir 6 questions ou moins juste en devinant ?

probability <- pbinom(6, size = 10, prob = 0.25)
print(probability)

Lorsque vous exécutez ce code, vous obtenez :

[1] 0.9803073

Cela signifie qu'il y a une chance de 98.03% d'obtenir 6 questions ou moins juste en devinant !

qbinom() : La Quête du Quantile

Qu'est-ce que qbinom() fait ?

La fonction qbinom() est comme l'inverse de pbinom(). Elle trouve le nombre minimum de succès pour une probabilité cumulative donnée.

Syntaxe et Paramètres

qbinom(p, size, prob, lower.tail = TRUE)
  • p : la probabilité cumulative
  • size : le nombre de trials
  • prob : la probabilité de succès à chaque trial
  • lower.tail : si TRUE (par défaut), utilise P(X ≤ x) ; si FALSE, utilise P(X > x)

Exemple : Contrôle de Qualité

Une usine produit des ampoules. Ils veulent savoir le nombre maximum d'ampoules défectueuses qu'ils peuvent avoir dans un lot de 100 pour maintenir un standard de qualité à 95%.

max_defects <- qbinom(0.05, size = 100, prob = 0.03, lower.tail = FALSE)
print(max_defects)

Ce code affichera :

[1] 6

Cela signifie qu'ils peuvent avoir au maximum 6 ampoules défectueuses pour maintenir leur standard de qualité à 95%.

rbinom() : Générateur de Nombres Aléatoires

Qu'est-ce que rbinom() fait ?

La fonction rbinom() génère des nombres aléatoires à partir d'une distribution binomiale. C'est comme simuler des trials réelles !

Syntaxe et Paramètres

rbinom(n, size, prob)
  • n : le nombre de valeurs aléatoires à générer
  • size : le nombre de trials
  • prob : la probabilité de succès à chaque trial

Exemple : Simulation de Lancers de Pièce

Simulons le lancer d'une pièce 10 fois, et nous ferons cette expérience 5 fois :

simulations <- rbinom(5, size = 10, prob = 0.5)
print(simulations)

Vous pourriez obtenir un résultat comme :

[1] 4 6 5 3 7

Chaque nombre représente le nombre de faces dans un ensemble de 10 lancers de pièce. Comme vous pouvez le voir, c'est aléatoire à chaque fois !

Mettre Tout Ensemble

Maintenant que nous avons exploré chaque fonction, utilisons-les dans une situation pratique. Imaginons que vous êtes un·e prévisionniste météo essayant de prédire les jours de pluie.

# Probabilité de exactement 3 jours de pluie dans une semaine de 7 jours
exactly_three <- dbinom(3, size = 7, prob = 0.3)

# Probabilité de 3 jours de pluie ou moins
three_or_fewer <- pbinom(3, size = 7, prob = 0.3)

# Nombre de jours de pluie avec une probabilité de 80%
days_with_80_percent <- qbinom(0.8, size = 7, prob = 0.3)

# Simuler 10 semaines de pluie
rain_simulations <- rbinom(10, size = 7, prob = 0.3)

print(paste("Probabilité de exactement 3 jours de pluie :", exactly_three))
print(paste("Probabilité de 3 jours de pluie ou moins :", three_or_fewer))
print(paste("Jours de pluie avec une probabilité de 80% :", days_with_80_percent))
print("Simulations de jours de pluie pour 10 semaines:")
print(rain_simulations)

Cet exemple complet montre comment ces fonctions travaillent ensemble pour analyser et prédire les modèles de jours de pluie.

Conclusion

Félicitations ! Vous avez fait vos premiers pas dans le monde des distributions binomiales en R. Souvenez-vous, la pratique rend parfait. Essayez de modifier les chiffres dans ces exemples et voyez ce qui se passe. Bientôt, vous utiliserez ces fonctions comme un·e data scientist chevronné·e !

Bonne programmation, et puissent les probabilités toujours être en votre faveur !

Credits: Image by storyset