Barplot et courbe pourcentage cumulé

**did10** · 19/08/2009, 11h26

Bonjour,

J'ai tracé un barplot et je souhaiterais y ajouter une courbe de fréquence cumulée.

Les effectifs sont classé dans une variable x par ordre décroissant, j'ai donc calculé d'abord les fréquences de chacun:

Code :

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

1
2
3
h=list()
for(i in 1:length(x)){
h[i]=x[i]*100/sum(x)}

Puis les fréquences cumulées:

Code :

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

1
2
3
4
h2=list()
h2[1]=h[[1]]
for(i in 2:length(h2)){
h2[i]=h2[[i]]+h2[[i-1]]}

Ensuite j'ai mis les valeurs à l'échelle (celle de gauche):

Code :

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

1
2
for(i in 1:length(h2)){
h2[i]=h2[[i]]*max(x)/100}

Et j'ai tracé mon barplot de cette façon:

Code :

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

1
2
3
4
5
6
par(mfrow=c(1,1),....)
for(i in 1:length(h2)){
bar=barplot(x,...)
points(h2[[i]],type='b',lty=1,pch=10,col=2)
axis(1,bar,labels=l,tick=T)
axis(side=4,at=c(0,max(x)*0.25,max(x)*0.5,max(x)*0.75,max(x)),labels=c(0.0,25.0,50.0,75.0,100.0))

Problème: ma courbe de fréquence ne s'affiche pas ...

Merci de votre aide
Did

**pitipoisson** · 19/08/2009, 12h41

Bonjour,

Je pense qu'il y a une erreur dans ton calcul de h2... et surtout que tu te complique beaucoup la vie. Tu devrais essayer autant que possible de vectoriser les calculs et éviter les boucles.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
x <- sort(rnorm(20, 2.5, 1), decreasing=TRUE)
x
 
h <- x * 100 / sum(x)
h
 
h2 <- cumsum(h)
h2
 
x2 <- barplot(x)
par(new=TRUE)
plot(x2, h2, xaxt="n", yaxt="n", bty="n")
axis(side=4)

**did10** · 19/08/2009, 14h13

Merci, c'est vari que c'est beaucoup plus simple comme ça.
Derniers points:

- comment faire pour l'axe de droite part de 0 ?
- comment virer les noms h2 et bar sur le graphe ? (cf. figure)

Merci

**pitipoisson** · 19/08/2009, 14h21

Code :

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

1
2
3
4
x2 <- barplot(x)
par(new=TRUE)
plot(x2, h2, xaxt="n", yaxt="n", bty="n", ylab="", xlab="")
axis(side=4, at=c(seq(0, 1, length=5) * max(h2)), labels=paste(seq(0, 1, length=5) * max(h2), "%"), xpd=NA)

**did10** · 19/08/2009, 14h56

Ok, il s'avère maintenant que pour un autre pareto j'ai été obligé de rajouter un autre axis() pour l'axe des abscisses. Auparavant les labels étaient bien centrés au milieu de chaque bande horizontale, mais depuis que j'ai ajouté axis(side=4), ça n'est plus centré, voici mon code:

Code :

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

1
2
3
4
5
6
7
8
h21=x*100/sum(x)
h22=cumsum(h21)
 
bar2=barplot(x,las=2,ylim=c(0,max(x)+2),col=4,border="darkblue",main="Out Of Control EWMA_MEAN",ylab="Count",cex.names=0.6)
par(new=TRUE,bty='l',cex.axis=0.6,cex.main=1.5,cex.lab=1.5,las=2,adj=0,col.main="darkblue",yaxs="r")
plot(bar2,h22,type='b',lty=1,pch='19',col=2,xaxt="n",yaxt="n",bty="n",ylab="",xlab="")
axis(1,bar2,labels=l,tick=T)
axis(side=4,at=c(seq(0,1,length=5)*max(h22)),labels=paste(seq(0,1,length=5)*max(h22),"%"),xpd=NA)