`wsmooth` weighted whittaker smoothing.
raw_GCxGC or preproc_GCxGC object with name and mod_time slots.
an integer of the order of the penalty. Only penalty of first (penalty = 1) and second (penalty = 2) order are allowed. By default, the smooth function is performed with first penalty order.
smoothing parameter: larger values lead to more smoothing.
minimum intensity value. The smoothing routine usually creates low intensity artifacts which can obscure other compounds signals. The min intensity value replace signals bellow the given value with 0. For quadrupole mass detector this artifacts may range from 0-100, while for TOF mass analyzers it can be 0-1e3.
This function takes a raw two-dimensional chromatogram and performs the weighted wittaker smoothing. It smooths the signal with linear or quadratic penalty, depending on the provided penalty, along side the first dimension, based on Whittaker smoother (Eilers 2003) .
chrom_name <- system.file("extdata", "08GB.cdf", package = "RGCxGC") chrom_2D <- read_chrom(chrom_name, 5L) #> Warning: The last 51 signals will be omitted #> Warning: data length  is not a sub-multiple or multiple of the number of rows  #> Retention time ranges: #> 1D (min): 7.98 18.16 #> 2D (sec): 0 5 #> Acquisition rate: 100 chrom_smooth <- wsmooth(chrom_2D, penalty = 1, lambda = 1e1) plot(chrom_smooth, nlevels = 150, color.palette = colorRamps::matlab.like, main = expression(paste(lambda, "= 10, penalty = 1")) ) # Remove intensities bellow 1.75e5 (too high) chrom_smooth2 <- wsmooth(chrom_2D, penalty = 1, lambda = 1e1, min_int = 1.75e5) plot(chrom_smooth2, nlevels = 150, color.palette = colorRamps::matlab.like, main = expression(paste(lambda, "= 10, penalty = 1, min_int = 1.75e5")) )