Discussion :

[marjo92blabla]

Salut à tous,

J'aimerais simplement placer des barres d'erreur sur mon diagramme en barres qui affiche la longueur des sépales selon l'espèce (3 colonnes pour 3 espèces).
Avant l'ajout des barres d'erreur il n'y a pas de souci.
J'essaye ce code sans succès, j'obtiens le message d'erreur:

"Erreur : Aesthetics must be either length 1 or the same as the data (150): ymin"

Code :

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

1
2
3
4
5
 
Variance <- var(iris$Sepal.Length, na.rm = T)
Ecart <- sqrt(Variance)
 
ggplot(iris, aes(x = Species, y = iris$Sepal.Length)) + geom_col() + geom_errorbar(aes(ymax = iris$Sepal.Length+Ecart, ymin = iris$Sepal.Lenght-Ecart), data=iris)

J'ai retourné le problème dans tous les sens mais je ne comprends pas ...
J'ai pris le jeu de données "iris" sur R.

Auriez-vous une explication à me donner ... ?

[olivier.decourt]

Bonjour.
Déjà, si le jeu de données est toujours le même, il suffit de l'indiquer dans ggplot() et pas nécessaire de le rappeler ensuite.

Code :

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

1
2
3
4
5
ggplot(iris, 
       aes(x = Species, y = Sepal.Length)) + 
  geom_col() + 
  geom_errorbar(aes(ymax = Sepal.Length+Ecart, 
                    ymin = Sepal.Length-Ecart))

Reste la question de ce que tu veux obtenir.

Sepal.Length+Ecart est un vecteur de 150 valeurs (chaque valeur de Sepal.Length à laquelle on ajoute l'écart-type). Comment utiliser ces 150 valeurs comme extrémités de barres d'erreur ? En superposant 150 barres d'erreurs ?
On pourrait ajouter une fonction statistique (mean(Sepal.Length)+Ecart par exemple) mais la moyenne ne serait pas calculée séparément pour chaque barre, mais au global sur toutes les espèces.

Dans ce type de graphique, le plus simple est souvent de préparer les données en amont avec toutes les quantités déjà calculées : la hauteur des barres, les extrémités des barres d'erreur, et d'avoir autant de lignes que de bâtons dans le graphique. Avec cette organisation des données, le code s'écrit de manière beaucoup plus facile.

[marjo92blabla]

Bonjour,

Tout d'abord merci de votre réponse. J'ai pu avancer sur mon code mais il n'est pas encore abouti.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
Ecart_type <- rep(c(Sd_Setosa, Sd_Versicolor, Sd_Virginica), each=50)
Ecart_type
 
## Ajout d'une colonne avec les écarts-type
 
Tab_Complet <- mutate(iris, Ecart_type) %>%
  select(Ecart_type, Sepal.Length, Species) 
Tab_Complet
 
library(ggplot2)
ggplot(Tab_Complet, aes(x = Species, y = Sepal.Length)) +
  labs (x="Espèce", y="Longueur", title="Liaison entre espèce et longueur de sépale") +
  geom_bar(stat="identity") 
  #geom_errorbar(aes(ymin=Sepal.Length-Ecart_type,ymax=Sepal.Length+Ecart_type))

J'ai créé le vecteur "Ecart_Type" avec les 50 valeurs associées à chaque espèce.

Cependant, je me rends compte que la simple réalisation du mon diagramme en barre pose problème.
Les valeurs sur l'axe y ne correspondent pas aux données sur la longueur des sépales.

Les 3 moyennes sont les suivantes :

Code :

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

1
2
3
4
5
6
> mean(Tab_Versicolor$Sepal.Length,na.rm=T)
[1] 5.936
> mean(Tab_Setosa$Sepal.Length,na.rm=T)
[1] 5.006
> mean(Tab_Virginica$Sepal.Length,na.rm=T)
[1] 6.588

J'aurais grand besoin de votre aide ... !

[olivier.decourt]

Bonjour.
Les hauteurs des barres correspondent à la somme des longueurs de sépales. Plutôt que stat="identity" il faudrait écrire stat="summary", fun=mean pour obtenir des barres proportionnelles aux moyennes. Ou mettre en entrée, comme je le proposais dans ma réponse précédente, une table déjà agrégée (3 lignes, une par espèce), avec 4 colonnes : espèce, moyenne, borne inf et borne sup de la barre d'erreur. Dans cette configuration stat="identity" sera approprié.