data(SpikePos)
plot(SpikePos$annotation[,c("rt", "mz")],
xlab = "Time (s)", ylab = "m/z",
main = "Positive ionization mode")
points(pos.markers[!is.na(pos.markers$feature.nr), c("rt", "mz")],
pch = 19, col = 2)
data(SpikeNeg)
plot(SpikeNeg$annotation[,c("rt", "mz")],
xlab = "Time (s)", ylab = "m/z",
main = "Negative ionization mode")
points(neg.markers[!is.na(neg.markers$feature.nr), c("rt", "mz")],
pch = 19, col = 2)
Run the code above in your browser using DataLab