## load library
library("MALDIquant");
## load example data
data("sA1", package="MALDIquant");
## show some information
sA1
## running typical workflow
## transform intensities
t <- lapply(sA1, function(x)transformIntensity(x, sqrt));
## smoothing function
movingAvg <- function(y) {
ma <- rep(1, 5)/5;
return(filter(y, ma, sides=2));
}
## smooth spectra
s <- lapply(t, function(x)transformIntensity(x, movingAvg));
## baseline correction
b <- lapply(s, function(x)removeBaseline(x));
## detect peaks
p <- lapply(b, function(x)detectPeaks(x));
## plot results
par(mfcol=c(3, length(sA1)));
for (i in seq(along=sA1)) {
plot(sA1[[i]], main="raw spectrum");
plot(s[[i]], main="transformed (sqrt) and smoothed (ma5)");
plot(b[[i]], main="baseline removed and peaks detected");
points(p[[i]], pch=4, col=2);
}
par(mfrow=c(1, 1));Run the code above in your browser using DataLab