context("estimateNoise-methods")
s <- createMassSpectrum(mass=1:20, intensity=rep(10:1, 2))
test_that("estimateNoise throws errors", {
expect_error(estimateNoise(s, method="foobar"),
".*arg.* should be one of .*MAD.*, .*SuperSmoother.*")
})
test_that("estimateNoise shows warnings", {
expect_warning(estimateNoise(
createMassSpectrum(mass=double(), intensity=double())),
"empty")
expect_identical(suppressWarnings(estimateNoise(
createMassSpectrum(mass=double(), intensity=double()))),
0L)
})
test_that(".estimateNoiseMad", {
n <- rep(stats::mad(intensity(s)), length(s))
m <- cbind(mass=mass(s), intensity=n)
expect_identical(MALDIquant:::.estimateNoiseMad(1:20, intensity(s)), n)
expect_identical(MALDIquant:::.estimateNoise(mass(s), intensity(s)), n)
expect_identical(MALDIquant:::.estimateNoise(mass(s), intensity(s),
method="MAD"), n)
expect_identical(estimateNoise(s), m)
expect_identical(estimateNoise(s, method="MAD"), m)
})
test_that(".estimateNoiseSuperSmoother", {
n <- stats::supsmu(x=1:20, y=intensity(s))$y
m <- cbind(mass=mass(s), intensity=n)
expect_identical(MALDIquant:::.estimateNoiseSuperSmoother(mass(s),
intensity(s)), n)
expect_identical(MALDIquant:::.estimateNoise(mass(s), intensity(s),
method="SuperSmoother"), n)
expect_identical(estimateNoise(s, method="SuperSmoother"), m)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.