Diagramme de Pareto avec R
Il est possible de réaliser un diagramme de Pareto avec la librairie qcc destinée à la qualité (Méthode 2). On peut aussi faire appel au code donné dans la méthode 1 plus adaptée pour la mise en forme.
1- Réaliser un diagramme de Pareto par la méthode 1
Télécharger le script complet ici pour réaliser un tel diagramme de Pareto sinon suivre le code ci-dessous.
Suggestion d'amélioration : permettre d'afficher les données non-converties en fréquence sur l'axe des ordonnées.
Diagramme de Pareto réalisé avec la fonction R ci-dessous.
- Voici une commande permettant de réaliser un diagramme de Pareto : cette commande peut être modifiée ou tapée directement dans R
# Fonction pareto
pareto = function(x, bar.col="cyan", line.col="red", pch=16, h=80, h.lty=3,main="",xlab="Défauts",ylab="Fréquence (%)", names.arg=c(), ylab2="Cumul",mar=c(5,4,3,4)) {
if (length(names.arg)>0) {names.arg=names.arg[order(x, decreasing = TRUE)]}
x = sort(x,decreasing=T); x = x*100/sum(x);
cumul = (cumsum(x)/sum(x))*100
simulation = barplot(x,col=bar.col) ; simulation
plot.new()
par(mar=mar)
barplot(x,col=bar.col,axes=F,ylim=c(0,100),main=main,xlab=xlab,ylab="",names.arg=names.arg)
#par(new=TRUE)
points(simulation,cumul,pch=pch,col=line.col,xlab="",ylab="",type="o")
abline(h=h,lty=h.lty) ; box()
axis(2) ; axis(4,c(0,20,40,60,80,100),col.axis=line.col,col=line.col)
mtext(ylab,side=2,line=2,cex=1.2) ; mtext(ylab2,side=4,col="red",line=2,cex=1.2)
result = c(x , cumul) ; result = matrix(result,nc=length(x), byrow=T)
if (length(names.arg)>0) {colnames(result) = names.arg }
rownames(result) = c("frequency","cumul")
return(result)}
2. Préparer les données à traiter. Voici un vecteur contenant des valeurs qui estiment le nombre de problèmes (ou le pourcentage) apparus liés à une cause.
# Valeurs dans une liste valeurs
valeurs = c(20,10,12,5,2,80)
description = c("Para1","Para2","Para3","Para4","Para5","Para6")
3. Tracer le diagramme de Pareto et récupérer les résultats d'analyse
pareto(valeurs,main="Diagramme de Pareto", names.arg=description, bar.col="blue")
2- Réaliser un diagramme de Pareto par la méthode 2 (librairie qcc)
Diagramme de Pareto réalisé avec la librairie qcc de R.
- Installer la librairie qcc et l'exécuter
install.packages("qcc")
library(qcc)
2. Voici un vecteur x contenant des valeurs qui estiment le nombre de problèmes (ou le pourcentage) apparu lié à une cause
# Valeurs dans une liste x
x = c(20,10,12,5,2,80)
3. Voici une commande permettant de réaliser un diagramme de Pareto : cette commande vient de la librairie qcc
pareto.chart(x, cumperc = seq(0, 100, by = 5), ylab = "Fréquence", ylab2 = "Cumul")
?pareto.chart # accéder à l'aide