Discussion :

[enicnath]

Voici le programme suivant :

Code :

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

1
2
3
4
A=rbind(c(1,0,0,0,0,1),c(0,0,0,2,1,2))
B=rbind(c(1,NA,NA,NA,NA,0.9),c(NA,NA,NA,1.8,0.8,1.7))
barplot(A,col=c('blue','red'))			#barres verticales bleues et rouges
barplot(B,col=c('green','cyan'),add=TRUE)	#barres verticales vertes et cyan

Lorsqu'on lance le 1er barplot, les barres bleues et rouges s'affichent bien. Elles ont les hauteur qui correspondent à leur valeur. 1 et 1 en bleu aux position 1 et 6, 2-1-2 en rouge aux positions 3, 4, 5.

Lorsqu'on lance le 2ème barplot, qui est en superposition du 1er, les barres vertes s'affichent, en position 1 et 6, et avec une hauteur de 1 et 0.9.
Jusque là, rien à dire.

En revanche, une seule barre cyan s'affiche, en position 6, et sur une hauteur correcte de 1.7.

On a donc deux barres verticales de couleur cyan qui n'apparaissent pas.
Elles devraient être en position 4 et 5, et avoir les hauteurs de 1.8 et 0.8.

Maintenant, si on remplace les NA dans le vecteur B par des 0, tout s'affiche correctement.

Est-ce que quelqu'un saurait me dire ce qui se passe ?

[mgdondon]

Bonjour,

Pour comprendre ce qui se passe, commencez par réaliser le deuxième barplot sans le superposer au premier.

Code :

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

1
2
3
4
5
> B
     [,1] [,2] [,3] [,4] [,5] [,6]
[1,]    1   NA   NA   NA   NA  0.9
[2,]   NA   NA   NA  1.8  0.8  1.7
> barplot(B,col=c('green','cyan'))

R ne dessine pas les barres correspondant à B[2,4] et B[2,5] parce qu'il ne sait pas à quel niveau les faire commencer.

Cordialement,

[enicnath]

Ok, mais alors pourquoi arrive-t-il à dessiner les barres vertes, qui correspondent à B[1,] ?

[mgdondon]

Je pense que c'est parce que c'est la première donc on est sûr qu'elle commence à 0.

En inversant l'ordre des vecteurs :

Code :

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

1
2
3
4
5
6
> C=rbind(c(NA,NA,NA,1.8,0.8,1.7),c(1,NA,NA,NA,NA,0.9))
> C
     [,1] [,2] [,3] [,4] [,5] [,6]
[1,]   NA   NA   NA  1.8  0.8  1.7
[2,]    1   NA   NA   NA   NA  0.9
> barplot(C,col=c('green','cyan'))

Cordialement,

[faubry]

En effet, barplot commence par exploiter la première ligne. Quand la fonction trouve des NA, elle interprète cela comme étant des positions pour lesquelles on n'a rien à dire, c'est-à-dire des positions à ne pas prendre en compte dans le graphique. Quand il y a zéro, cela signifie que la mesure est nulle. Ce comportement est logique car NA signifie formellement qu'il existe une valeur mais qu'on ne sait pas laquelle. Donc assimiler NA à zéro est formellement faux et, comme le dit l'adage, dans le doute abstiens-toi. C'est ce que fait barplot.

Preuve

Code :

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

1
2
B=rbind(c(1,NA,NA,NA,NA,0.9),c(NA,NA,NA,1.8,0.8,1.7), cbind( 2, 2, 2, 2, 2, 2))
barplot(B,col=c('green','cyan'))

En position 1, la troisième ligne ne sera pas tracée.

[enicnath]

Ok, J'éviterais donc les NA dorénavant lorsque j'aurais des doutes

Merci à vous.