estimateNoise-methods | R Documentation |
This method estimates the noise of mass spectrometry data
(represented by a MassSpectrum
object).
## S4 method for signature 'MassSpectrum'
estimateNoise(object,
method=c("MAD", "SuperSmoother"),
...)
object |
|
method |
used noise estimation method, one of
|
... |
arguments to be passed to |
"MAD"
:This function estimates the noise of mass spectrometry data
by calculating the median absolute deviation, see also
mad
.
"SuperSmoother"
:This function estimates the noise of mass spectrometry data using
Friedman's Super Smoother. Please refer
supsmu
for details and additional arguments.
Returns a two column matrix (first column: mass, second column: intensity) of the estimated noise.
Sebastian Gibb mail@sebastiangibb.de
MassSpectrum
,
detectPeaks,MassSpectrum-method
,
mad
,
supsmu
Website: https://strimmerlab.github.io/software/maldiquant/
## load package
library("MALDIquant")
## load example data
data("fiedler2009subset", package="MALDIquant")
## choose only the first mass spectrum
s <- fiedler2009subset[[1]]
## transform intensities
s <- transformIntensity(s, method="sqrt")
## remove baseline
s <- removeBaseline(s)
## plot spectrum
plot(s)
## estimate noise
nm <- estimateNoise(s, method="MAD")
nss <- estimateNoise(s, method="SuperSmoother")
## draw noise on the plot
lines(nm, col=2)
lines(nss, col=4)
## draw legend
legend(x="topright", lwd=1, legend=c("MAD", "SuperSmoother"),
col=c(2, 4))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.