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.
- 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)
- 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")
- 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 (%)"
)
- Afficher l'axe des abscisses
axis(
side=1, # X axis
at=c(2,3,4),
labels=c("Etape 1","Etape 2","Etape 3")
)
- 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
- 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%")