## 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
n <- estimateNoise(s)
## draw noise on the plot
lines(n, col="red")Run the code above in your browser using DataLab