MALDIquant (version 1.6)

movingAverage: Moving Average

Description

This function runs a simple 2-side moving average.

Usage

movingAverage(y, halfWindowSize=2)

Arguments

y

intensity values

halfWindowSize

half window size. The resulting window reaches from mass[currentIndex-halfWindowSize] to mass[currentIndex+halfWindowSize] (window size is 2*halfWindowSize+1).

Value

• a smoothed intensity vector

Details

The movingAverage function is defined as: movingAverage <- function(y, halfWindowSize = 2) { windowSize <- halfWindowSize * 2 + 1 y <- filter(x=y, filter=rep(1, windowSize)/windowSize, sides=2) return(y) }

See Also

filter, MassSpectrum, transformIntensity,AbstractMassObject-method

Examples

## load library
library("MALDIquant")

## 5 point moving average
movingAverage(rep(1:5, times=2), halfWindowSize=2)
# NA NA  3  3  3  3  3  3 NA NA

## 9 point moving average
movingAverage(rep(1:5, times=2), halfWindowSize=4)
# NA NA NA NA 2.777778 3.222222 NA NA NA NA

## load example data
data("fiedler2009subset", package="MALDIquant")

## smooth spectra
s <- transformIntensity(fiedler2009subset, movingAverage, halfWindowSize=2)