Présentation
On souhaite simuler le lancer d’un dé équilibré.
Obtenir un échantillon
Créons la liste des faces :
faces<-1:6 faces
[1] 1 2 3 4 5 6
Choisissons une taille d’échantillon :
tailleEchantillon<-10000 taillEchantillon
[1] 10000
Construisons cet échantillon avec la fonction sample :
echantillon<-sample(faces,tailleEchantillon,replace=TRUE)
Voici ces 100 premières valeurs :
head(echantillon,100)
[1] 3 2 2 5 3 1 2 3 3 1 1 1 1 3 5 1 1 5 6 2 5 1 1 1 5 3 4 5 1 6 1 6 1 2 3 2 5 4 1 3 5 1 2 6 3
[46] 1 6 6 5 5 1 2 1 3 5 4 2 4 1 5 1 3 1 5 6 4 4 6 6 6 1 1 5 6 5 4 1 1 2 2 4 5 5 1 1 6 2 1 6 5
[91] 1 2 3 5 4 3 1 5 4 3
Analyse de l’échantillon
Pour obtenir le tableau des effectifs :
tabEffectifs<-table(echantillon) tabEffectifs
echantillon
1 2 3 4 5 6
1687 1684 1651 1693 1672 1613
Pour obtenir la tableau des fréquences :
tabFrequences<-prop.table(tabEffectifs)
echantillon
1 2 3 4 5 6
0.1687 0.1684 0.1651 0.1693 0.1672 0.1613
On peut alors les représenter par un diagramme en bâtons.
Le diagramme basique est le suivant :
barplot(tabEffectifs)
Utilisation d’un script
Pour plus de facilité, il suffit de sauvegarder le code suivant dans un fichier script et de le lancer. On pourra si besoin changer la valeur de la taille de l’échantillon.
###################################################################################
###################################################################################
# Simulation du lancer d'un dé équilibré
###################################################################################
###################################################################################
# Liste des faces
faces<-1:6
print("Liste des faces :")
print(faces)
# Taille de l'échantillon
tailleEchantillon<-1000 #-----------------------------------A modidier
print("Taille de l'échantillon")
print(tailleEchantillon)
# Création de l'échantillon
echantillon<-sample(faces,tailleEchantillon,replace=TRUE)
# Tableau des fréquences
tabEffectifs<-table(echantillon)
print("Tableau des effectifs")
print(tabEffectifs)
# diagramme en bâtons des effectifs
barplot(tabEffectifs)
tabFrequences<-prop.table(tabEffectifs)
print("Tableau des fréquences")
print(tabFrequences)
titreEffectifs<-paste("Lancer d'un dé équilibré","\n",
"Taille de l'échantillon : ",bquote(.(tailleEchantillon)))
barplot(tabEffectifs,ylim=c(0,tailleEchantillon/5),las=1,col="blue"
, xlab="Faces",ylab="Effectifs",
,main=titreEffectifs)
abline(a=tailleEchantillon/6,b=0,lwd=1.5,col="grey",lty=2)
# diagramme en bâtons des fréquences
titreFrequences<-paste("Lancer d'un dé équilibré","\n",
"Taille de l'échantillon : ",bquote(.(tailleEchantillon)))
barplot(tabFrequences,ylim=c(0,0.2),las=1,col="brown"
, xlab="Faces",ylab="Fréquences",
,main=titreFrequences)
abline(a=1/6,b=0,lwd=1.5,col="grey",lty=2)
On obtient dans la console :
[1] "Liste des faces :"
[1] 1 2 3 4 5 6
[1] "Taille de l'échantillon"
[1] 1000
[1] "Tableau des effectifs"
echantillon
1 2 3 4 5 6
187 153 170 162 144 184
[1] "Tableau des fréquences"
echantillon
1 2 3 4 5 6
0.187 0.153 0.170 0.162 0.144 0.184
On obtient les deux graphiques suivants.
