Description Usage Arguments Value See Also Examples
Find proteins that possibly have multiple charges. This is achieved by detecing proteins whose mz values are nearly exact multiples of others and hence potentially represent the same protein.
1 2 |
x |
An object of class |
... |
Not used. |
event |
A character string denoting the name of the
event to register with the (embedded) event history object of the input
after processing the input data. Default: |
mz.precision |
A numeric value, used to construct the threshold
when comparing the multiples of mz values.
The default value is 0.003 because SELDI data is often assumed to have
+/- 0.3% mass drift, i.e., a peak at mass |
ncharge |
A numeric integer vector denoting the multiple charges of interest. All of its elements must be larger than 1. |
An object of class msSet
with charge estimate attached
as element "peak.charge"
.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 | if (!exists("qcset")) data("qcset", package="msProcess")
## extract several spectra from the build-in
## dataset
z <- qcset[, 1:8]
## denoising
z <- msDenoise(z, FUN="wavelet", n.level=10, thresh.scale=2)
## local noise estimation
z <- msNoise(z, FUN="mean")
## baseline subtraction
z <- msDetrend(z, FUN="monotone", attach=TRUE)
## intensity normalization
z <- msNormalize(z)
## peak detection
z <- msPeak(z, FUN="simple", use.mean=FALSE, snr=2)
## peak alignment
z <- msAlign(z, FUN="cluster", snr.thresh=10,
mz.precision=0.004)
## charge detection
z <- msCharge(z, ncharge=2:5, mz.precision=0.003)
## extract the peak.charge
z[["peak.charge"]]
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.