Voici mont problème:
Je travaille sur la mise au point d'un ELISA de VTG. Je dois tester le parallélisme entre une courbe standard type et des dilutions successives de plasma.
Alors voici mon approche.
Graph 1: J'ai fait un graphique de ma courbe standard illustrant une régression à 4 paramètres. L'axe des y correspond a mon % de binding et l'axe des x à ma concentration en VTG.
Graph 2 : J'ai fait un graphique avec mes dilutions successives de plasma.L'axe des y correspond toujours à mon % de binding, mais l'axe des x correspond à mon facteur de dilution.
Dans les deux cas mon axe des x est une échelle logarithmique.
Voici mon problème. Je dois collé les graphiques cote à cote pour en faire un seul afin de comparer l'allure des pentes. Mais dans mon graphique 1, un cm = un faccteur de 10 et dans mon graphique 2, 1cm= un facteur de 100. Je veux avoir deux échelles comparable donc de facteur identique.
J'espère que vous comprenez ce que je veux?
Voici un exemple
Code : Tout sélectionner
library(drc)
#////////////
#standard
#////////////
FACT <- c( 8000, 8000, 8000, 3200, 3200, 3200, 1280, 1280, 1280, 512, 512, 512, 204, 204, 204, 81, 81, 81, 32, 32, 32, 13, 13, 13, 5, 5, 5, 2, 2, 2)
SLOPE <-c(24,20,20,27,20,23,29,25,20,42,28,32,44,39,33,52,52,50,79,78,65,99,82,79,95,90, 89,104,102,89)
standard <- data.frame(FACT,SLOPE)
modelstandard <- multdrc(standard$SLOPE~standard$FACT)
#////////////
#Échantillon
#////////////
FACT <- c(1.0e-02, 4.0e-03, 1.6e-03, 1.6e-03, 1.0e-02, 4.0e-03, 6.4e-04, 6.4e-04, 1.0e-02, 4.0e-03, 2.5e-04, 2.5e-04, 2.5e-04, 1.0e-04, 1.0e-04, 4.2e-09, 1.7e-09, 4.0e-05, 4.0e-05, 4.2e-09, 1.7e-09, 1.6e-05, 1.6e-05, 1.6e-05, 6.8e-10, 6.5e-06, 6.5e-06, 6.5e-06, 6.8e-10, 2.7e-10, 2.6e-06, 2.6e-06, 2.6e-06, 6.8e-10, 2.7e-10, 1.0e-06, 1.0e-06, 1.0e-06, 2.6e-08, 1.0e-08, 4.1e-07, 4.19e-07, 4.1e-07, 2.6e-08, 1.0e-08, 1.6e-07, 1.6e-07, 1.6e-07, 2.6e-08, 6.7e-08, 6.7e-08, 6.7e-08)
SLOPE <- c(2.822958, 5.683305, 5.290709 , 6.861096 , 1.140400,3.103384, 5.795476, 7.534119 , 1.981679 , 4.449430,8.319312, 9.048420, 13.815666 , 6.580669 , 9.497102,101.701253, 82.912694, 14.096093, 13.759581 , 88.969901,87.231258, 25.425313, 25.481398 ,27.500467 ,101.028230,36.922789, 28.734343, 38.324921, 101.532997 ,104.112918, 44.606468, 48.252010, 54.477472 , 95.475790 ,105.402879, 51.168443, 59.469060, 67.713591 ,108.375397 ,103.832492, 74.612077, 73.322116 , 72.985605 ,99.682184 , 95.531875, 87.511684, 86.053468 ,89.923350 ,99.233502 ,94.241914, 89.306412, 90.876799)
echantillon <- data.frame(FACT,SLOPE)
modelechantillon <- multdrc(echantillon$SLOPE~echantillon$FACT)
par(mfrow=c(1,2))
plot(modelstandard, ylim=c(0,100),ylab= "Bi/B0(%)", xlab= "Vtg ng/ml",legend=FALSE)
plot(modelechantillon, conLevel= -3 ,ylim=c(0,100),ylab= "Bi/B0(%)", xlab= "Dilution factor")
par(mfrow=c(1,2))
Donc si vous regardez bien le graphique échantillon, l'échelle logarithmique est de facteur 100 contrairement de facteur 10 pour le graphique standard. Je veux avoir le même échelle afin de pouvoir comparez l'allure des deux pentes.
pour voir un exemple du type de graphique, je vous invite à voir l'article de Brion fig no 5 : Cross-reaction of whole-body...
http://www.setacjournals.org/archive/1552-8618/21/8/pdf/i1552-8618-21-8-1699.pdf
Merci beaucoup de votre aide
Bernard-Antonin Dupont Cyr