Il est parfois intéressant de calculer l'écart-relatif moyen qui existe entre des points et la droite de régression linéaire.
RÉALISER LE PREMIER GRAPHIQUE (Etalonnage)
# Absorbance à 615 nm (bleu dextran)BD=c(0,0,0,0.100,0.200,0.50,0.5,0.04,0.03,0.02,0.05,0.06,0.04,0.03,0.01,0,0,0,0,0.10) # Absorbance à 505 nm (vitamine B12)B12=c(0,0,0,0,0,0,0,0,0,0,0,0.1,0.1,0.2,0.3,0.2,0.2,0.1,0,0.1)2. Tracer la 1ere courbe
# Volumes cumulés des fractions : nombre_de_fraction = length(BD)fractions=rep(1:nombre_de_fraction)fractions=fractions*2.5# YLIM doit être changé.plot(fractions, BD, pch=16, axes=F, ylim=c(0,0.8), xlab="", ylab="", type="o",col="blue", main="Graphique à 2 axes et données lissées")axis(2, ylim=c(0,0.8),col="blue")mtext("Absorbance à XXX nm (Détection du Bleu Dextran)",side=2,line=2.5)box()3. Tracer la 2ème courbe
# Tracer les pointspar(new=T)par(mar=c(4,4,3,4))plot(fractions, B12, pch=15, xlab="", ylab="", ylim=c(0,0.8), axes=F, type="o", col="red")4. Légender
# Légende de l'axe Y de droitemtext("Absorbance à XXX nm (Détection de la vitamine B12)",side=4,col="red",line=2.5)axis(4, ylim=c(0,0.5), col="red",col.axis="red")# AXE Xaxis(1,pretty(range(fractions),10))mtext("Volume d'élution (mL)",side=1,col="black",line=2.5)# Legendlegend(x="topleft",legend=c("Bleu Dextran (A à XXXnm)","Vitamine B12 (A à XXXnm)"),text.col=c("blue","red"),pch=c(16,15),col=c("blue","red"))RÉALISER LE DEUXIÈME GRAPHIQUE (Caractérisation d'un mélange glucides/protides)
# Absorbance à 260 nmA260 = c(0,0,0,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,.1,0,0,0,0,0,.1,.1,.1,.1,.1,0)# Absorbance à 280 nmA280 = c(0,0,0,.596,.8,.9,.273,.209,.144,.08,.275,.338,.209,.145,.015,0,0,0,0,.017,.047,.060,.531,.209,0)protides = 1.55*A280-0.76*A260 # Convertit les absorbances à 260 nm et 280 nm en concentrations protéiques en mg/mL#Teneur en glucides (mesurée en fonction de la taille du précipité)glucides=c(0,0,0,0,0,0,0,0,0,0,0,1,1,2,3,2,2,1,0,0,0,0,2. Tracer les barres
# Volumes cumulés des fractionsnombre_de_fraction = length(protides)fractions=rep(1:nombre_de_fraction)volume=fractions*2.5# 1- Ouvrir une nouvelle fenêtre graphique#plot.new()# 2- Programmer des marges larges pour l'ajout ultérieur des titres des axespar(mar=c(4,4,3,4))# 3- Tracer des barres - NB : il est possible de tracer la courbe en 1erbarplot(glucides,col="orange",angle=0,axes=F,ylab="",xlab="",main="Teneur en glucides et protides\ndans des fractions obtenues par filtration sur gel",space=0, ylim=c(0,5))axis(2,col="red",col.axis="red")3. Tracer la courbe superposée aux barres
# 4- Superposer la courbepar(new=TRUE,mar=c(4,4,3,4))plot(c(0,volume),c(0,protides),col="black",type="o",lwd=2,pch=16,axes=F,ylab="",xlab="", xlim=c(0,max(volume)+2.5))axis(4,col.axis="black",col="black")axis(1,col.axis="black",col="black")box()4. Légender
# Titre de l'axe de gauchemtext("Précipités de glucose (UA)",side=2,line=2,cex=1.2,col="red")# Titre de l'axe de droitemtext("Teneur en protides (mg/mL)",side=4,col="black",line=2,cex=1.2)# AXE Xmtext("Volume cumulé des fractions (mL)",side=1,col="black",line=2.5)Mise à jour de la page nov 2018