Learn R Programming

xcms (version 1.48.0)

specNoise: Calculate noise for a sparse continuum mass spectrum

Description

Given a sparse continuum mass spectrum, determine regions where no signal is present, substituting half of the minimum intensity for those regions. Calculate the noise level as the weighted mean of the regions with signal and the regions without signal. If there is only one raw peak, return zero.

Usage

specNoise(spec, gap = quantile(diff(spec[, "mz"]), 0.9))

Arguments

spec
matrix with named columns mz and intensity
gap
threshold above which to data points are considerd to be separated by a blank region and not bridged by an interpolating line

Value

A numeric noise level

Details

The default gap value is determined from the 90th percentile of the pair-wise differences between adjacent mass values.

See Also

getSpec, specPeaks