Réaliser une gamme étalon avec R
Une gamme étalon est fondamentale pour réaliser de nombreux dosages en biochimie ou en chimie.
Il est possible de tracer une telle courbe ou droite avec le logiciel R project en réalisant une régression linéaire passant par l'origine ou au contraire coupant l'axe des ordonnées.
Exemple : Gamme étalon permettant un dosage X à 505 nm
Charger des données (une liste pour les coordonnées en x (abscisses) et une liste pour pour les coordonnées en y (ordonnées))
# Exemple - chargement de données
absorbance = c(0.256,0.347,0.522,0.657,0.804)
masse = c(3,6,9,12,15)
2. Rédiger les titres des axes (cette étape permet de les mettre en forme s'ils présentent du texte à mettre en indice ou en exposant) (pour en savoir plus sur la mise en forme d'expressions)
# Pour rédiger des expression mathématiques avec du texte en exposant ou en indice (fonction expression)
x_lab = expression(Concentration~~en~~mol.L^-1) # ~~ pour les espaces ici mol.L-1 voit -1 mis en exposant
y_lab = expression(A[505]) # ici la longueur d'onde de 505 nm est mise en indice pour avoir A505
3. Tracer le graphique
plot(masse,absorbance,xlim=c(0,15),ylim=c(0,0.9),pch=4,cex=1.5,font.lab=4,ylab=y_lab,xlab=x_lab)
# pch : type de points (ici des croix avec pch=4)
# cex=1.5 (taille de la police)
# font.lab : indique si les titres des axes sont normaux (1) , en gras, italique ou les 2 (4).
# Cette commande qui suit permet de quadriller le graphique
panel.first = grid()
4. Afficher une régression linéaire passant par l'origine
# Calculer la régression
regression = lm(absorbance~masse +0)$coefficients
# Tracer la régression
abline(0,(regression),col="orange",lwd=2)
# Afficher l'équation de la droite
text(10,0.2,paste("Equation de la droite\ny =",round(regression[1],3),"x"),col="orange")
# round ici arrondi les valeurs, le 3 indique le nombre de décimales, à paramétrer
5. Afficher une régression linéaire ne passant pas par l'origine
# Calculer la régression
regression = lm(absorbance~masse)$coefficients
# Tracer la régression
abline(regression,col="red",lwd=2)
# Afficher l'équation de la droite
text(10,0.1,paste("Equation de la droite\ny =",round(regression[2],3),"x+",round(regression[1],3)),col="red")
# round ici arrondi les valeurs, le 3 indique le nombre de décimales, à paramétrer
Mise à jour de la page avril 2024