MALDIquant (version 1.13)

estimateNoise-methods: Estimates the noise of a MassSpectrum object.

Description

This method estimates the noise of mass spectrometry data (represented by a MassSpectrum object).

Usage

## S3 method for class 'MassSpectrum':
estimateNoise(object,
  method=c("MAD", "SuperSmoother"),
  ...)

Arguments

object

MassSpectrum object

method

used noise estimation method, one of "MAD" or "SuperSmoother".

...

arguments to be passed to method.

Value

• Returns a two column matrix (first column: mass, second column: intensity) of the estimated noise.

Details

[object Object],[object Object]

See Also

MassSpectrum, detectPeaks,MassSpectrum-method, mad, supsmu

Website: http://strimmerlab.org/software/maldiquant/

Examples

## 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))