Only a certain percentage of most intense MS2 peaks are taken into
account to use the most significant signal. Default value is 10% (see
percentage
argument). The difference between peaks is then
computed for all individual spectra and their distribution is plotted
as a histogram where single bars represent 1 m/z differences. Delta
m/z between 40 and 200 are plotted by default, to encompass the
residue masses of all amino acids and several common contaminants,
although this can be changes with the xlim
argument.
In addition to the processing described above, isobaric reporter tag
peaks (see the reporters
argument) and the precursor peak (see
the precMz
argument) can also be removed from the MS2 spectrum,
to avoid interence with the fragment peaks.
Note that figures in Foster et al 2011 have been produced and
optimised for centroided data. Application of the plot as is for data
in profile mode has not been tested thoroughly, although the example
below suggest that it might work.
The methods make use the ggplot2
system. An object of class
ggplot
is returned invisibly.
Most of the code for plotMzDelta has kindly been contributed by Guangchuang Yu.
MSnExp
or mzRramp
(from the mzR
package) containing MS2 spectra. "ReporterIons"
that defines which reporter ion
peaks to set to 0. The default value NULL
leaves the spectra
as they are. object
's MS2 spectra. numeric
of length one or NULL
default. In the latter (and preferred) case, the precursor m/z
values are extracted from the individual MS2 spectra using the
precursorMz
method. numeric
of length 1 that specifies the
width around the precursor m/z where peaks are set to 0. Default is
2. numeric
specifying the bandwith to be used to bin
the delta m/z value to plot the histogram. Default if 1. See
geom_histogram
for more details. numeric
of length 2 specifying the range of
delta m/z to plot on the histogram. Default is c(40,200)
. logical
defining if amino acid residue
labels are plotted on the figure. Default is TRUE
. numeric
of length 1 specifying the font size of
amino acids lables. Default is 2.5.logical
of length 1 that defines whether the
figure should be plotted on the active device. Default is
TRUE
. Note that the ggplot
object is always returned
invisibly. logical
of length 1 specifying whether
textual output and a progress bar illustration the progress of data
processing should be printed. Default is TRUE
signature(object = "MSnExp", ...)
plotDensity
and plot2d
methods for
other QC plots.
mzdplot <- plotMzDelta(itraqdata,
subset = 0.5,
reporters = iTRAQ4,
verbose = FALSE, plot = FALSE)
## let's retrieve peptide sequence information
## and get a table of amino acids
peps <- as.character(fData(itraqdata)$PeptideSequence)
aas <- unlist(strsplit(peps,""))
## table of aas
table(aas)
## mzDelta plot
print(mzdplot)
Run the code above in your browser using DataLab