The goal of sgolay is to offer efficient and vectorized Savitzky-Golay filters.
You can install the CRAN version with
# install.packages("sgolay")Or you can install the development version of sgolay from GitHub with:
# install.packages("remotes")
remotes::install_github("zeehio/sgolay")sgolay is faster than signal, especially on
either larger filter lengths or when applied on matrices, since it uses
the Fast Fourier Transform and avoids several memory copies and extra
allocations.
library(sgolay)
x <- matrix(runif(1000*1000), nrow = 1000, ncol = 1000)
filt <- signal::sgolay(p = 2, n = 51)
timing <- bench::mark(
signal = apply(x, 2L, function(s) signal::sgolayfilt(s, filt)),
sgolay = sgolay::sgolayfilt(x, filt),
min_iterations = 50
)
#> Warning: Some expressions had a GC in every iteration; so filtering is disabled.
plot(timing, type = 'ridge')
#> Loading required namespace: tidyr
#> Picking joint bandwidth of 0.0135