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

Gamme étalon avec régression linéaire avec le logiciel R project

# 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