data(AA)
AA
AA.lm<-lm(AA$Monoisotopic ~ AA$Average)
plot(AA$Monoisotopic, AA$Average);
abline(AA.lm, col='grey')
text(AA$Monoisotopic, AA$Average, AA$letter1, pos=3)
plot(AA$Average-AA$Monoisotopic)
axis(3,1:20,AA$letter1);
abline(v=1:20,col='grey')
Run the code above in your browser using DataLab