context("alignSpectra")
test_that("alignSpectra throws errors", {
expect_error(alignSpectra(list()),
".*spectra.* is no list of MALDIquant::MassSpectrum objects")
})
test_that("alignSpectra", {
data(fiedler2009subset)
spectra <- fiedler2009subset[1:8]
peaks <- detectPeaks(spectra, halfWindowSize=20, method="MAD", SNR=2)
wf <- determineWarpingFunctions(peaks, tolerance=0.002, method="lowess")
expect_equal(alignSpectra(spectra), warpMassSpectra(spectra, wf))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.