Nothing
y <- c(1, 3, 2, 5, 6, 8, 10, 1)
w <- seq_along(y)/length(y)
test_that("smooth_wSG works", {
halfwin = 2
d = 2
s1 <- smooth_wSG(y, halfwin, d, w)
s2 <- smooth_SG(y, halfwin, d)
MAE = 1.5
expect_true(mean(abs(y - s1)) < MAE)
expect_true(mean(abs(y - s2)) < MAE)
})
test_that("movmean works", {
s <- smooth_wSG(y, 1, 1)
m <- movmean(y, 1)
m2 <- movmean2(y, 1, 1)
expect_true(all.equal(s[2:7], m[2:7]))
expect_true(all.equal(m, m2))
})
# test_that("movmean_R works", {
# arr = rnorm(5 * 5 * 10) %>% array(dim = c(5, 5, 10))
# r = movmean_R(arr, 1) %>% array_3dTo2d()
# # only suit for 2d array
# mat = array_3dTo2d(arr)
# r2 = movmean2_mat(mat, 1, 1)
# expect_true(max(rowMeans2(abs(r - r2))) < 1e-5)
# }
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.