inst/tinytest/test_rescale.R

In alkahest: Pre-Processing XY Data from Experimental Methods

data("XRD")

# Normalize by AUC =============================================================
XRD_area <- rescale_area(XRD)
expect_equal(integrate_rectangle(XRD_area), 1)

# Normalize by total intensity =================================================
XRD_total <- rescale_total(XRD, total = 1)
expect_equal(sum(XRD_total$y), 1)

# Normalize max to 1 ===========================================================
expect_error(rescale_max(XRD, max = 1), "must be lower than")

XRD_max <- rescale_max(XRD, max = 35000)
expect_equal(max(XRD_max$y), 35000)

# Normalize min to 0 ===========================================================
XRD_min <- rescale_min(XRD, min = 0)
expect_equal(min(XRD_min$y), 0)

# Normalize 0-1 ================================================================
XRD_range <- rescale_range(XRD, min = 0, max = 1)
expect_equal(min(XRD_range$y), 0)
expect_equal(max(XRD_range$y), 1)

# Transform ====================================================================
XRD_trans <- rescale_transform(XRD, f = sqrt)
expect_equal(max(XRD_trans$y), sqrt(max(XRD$count)))

# SNV ==========================================================================
data("Raman")
Raman_snv <- rescale_snv(Raman)
expect_equal(mean(Raman_snv$y), 0)
expect_equal(sd(Raman_snv$y), 1)