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 points
par(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 droite
mtext("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 X
axis(1,pretty(range(fractions),10))
mtext("Volume d'élution (mL)",side=1,col="black",line=2.5)
# Legend
legend(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 nm
A260 = 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 nm
A280 = 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 fractions
nombre_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 axes
par(mar=c(4,4,3,4))
# 3- Tracer des barres - NB : il est possible de tracer la courbe en 1er
barplot(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 courbe
par(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 gauche
mtext("Précipités de glucose (UA)",side=2,line=2,cex=1.2,col="red")
# Titre de l'axe de droite
mtext("Teneur en protides (mg/mL)",side=4,col="black",line=2,cex=1.2)
# AXE X
mtext("Volume cumulé des fractions (mL)",side=1,col="black",line=2.5)
Mise à jour de la page nov 2018