Etudier des rendements

En langage R

L'étude des rendements permet d'estimer les pertes à chaque étape d'un processus ou d'une filière.

Graphique d'étude des rendements d'un processus - logiciel R
  1. Voici un vecteur matiere contenant l'évolution de la masse sur 3 étapes partant d'une masse initiale de 100%
# Valeurs dans une liste matiere
matiere = c(100,97.5,94,88)
  1. On peut ainsi calculer les rendements intermédiaires
rendement1 = matiere[2]/matiere[1]*100 ; cat("Rendement intermédiaire de l'étape 1 : ",rendement1,"%.\n")
rendement2 = matiere[3]/matiere[2]*100 ; cat("Rendement intermédiaire de l'étape 2 : ",rendement2,"%.\n")
rendement3 = matiere[4]/matiere[3]*100 ; cat("Rendement intermédiaire de l'étape 3 : ",rendement3,"%.\n")
  1. Tracer le graphique non encore annoté
plot(matiere, 
    tck=-0.02,  # extériorise les encoches des graduations des axes
    type = "s",         # Type de courbe en marches
    ylim= c(85,102),
    xlim =c(1,4.8),
    yaxp = c(88,100,6), # Définit les graduations à afficher et le nombre d'intervalles
    yaxs="i", # Fait partir la première valeur de l'axe y directement du bas
    xaxs = "i",
    las = 1, # Ecriture des valeurs de graduations de l'axe des x et y horizontalement
    font.lab=2, # écriture des titres des axes (1= normal, 2=gras, 3 = italique...)
    xaxt = "n", # Supprime les graduations de l'axe x
    lwd = 2,
    xlab = "Etapes",
    ylab = "Quantité de matière (%)"
)
  1. Afficher l'axe des abscisses
axis(
  side=1, # X axis
  at=c(2,3,4),
  labels=c("Etape 1","Etape 2","Etape 3")
)
  1. Afficher l'axe des ordonnées avec une cassure entre 0 et 88
# Afficher la valeur 0
axis(side=2, # Y axis
  las=0, # text horizontal
  at=c(85), # position du 0
  labels=c("0") # Le 0 qui s'affiche là où il aurait eu un 85 (configuration manuelle)
)
# Afficher l'axe Y avec une coupure entre la première valeur et la valeur 0
install.packages("plotrix") # si le package plotrix n'est pas installé
library(plotrix)
axis.break(axis=2,breakpos=86.5,style="slash") # break the bottom X axis
  1. Afficher les rendements en configurant manuellement la position des traits et légendes
# rendement intermédiaire 1
arrows(2.1,97.7,2.1,99.8,
    angle=90,
    lwd=3,lty=1,
    length=0.1,
    code = 3 # flèche dans les 2 directions
)
text(2.5, 98.75, "Ri = 97.5%")
# rendement intermédiaire 2
arrows(3.1,94.2,3.1,97.3,
    angle=90,
    lwd=3,lty=1,
    length=0.1,
    code = 3 # flèche dans les 2 directions
)
text(3.5, 95.5, "Ri = 96.4%")
# rendement intermédiaire 3
arrows(4.1,93.8,4.1,88.2,
    angle=90,
    lwd=3,lty=1,
    length=0.1,
    code = 3 # flèche dans les 2 directions
)
text(4.5, 91, "Ri = 93.6%")