Nothing
if (at_home()) {
noise <- c(-0.002, 0, -0.014, -0.001, -0.008, -0.001, -0.021, 0.009, 0.016,
-0.006, -0.013, -0.01, -0.017, 0, -0.012, 0.01, 0.005, 0.003,
0.012, 0.002, -0.022, -0.004, 0.007, 0.002, -0.006, -0.014, 0.012,
0.015, -0.01, 0.001, -0.007, 0.01, 0.002, -0.018, -0.009, 0,
-0.009, -0.008, -0.018, -0.007, -0.013, -0.003, -0.001, -0.002,
0.002, -0.02, 0.012, -0.008, 0.018, -0.001, 0.009, -0.012, -0.003,
0.012, -0.001, -0.016, 0.013, -0.015, -0.01, 0.014, 0.013, -0.003,
0.007, 0.005, 0.008, 0.011, -0.007, 0.004, -0.005, 0.006, 0.006,
0.003, 0.003, -0.024, 0.014, 0.001, 0.01, -0.01, 0.009, -0.003,
-0.004, 0.004, -0.006, 0.008, -0.009, 0.006, -0.017, 0.001, -0.004,
-0.01, -0.001, 0, 0.014, 0.007, 0.009, 0.012, 0.013, -0.007,
0.017, -0.007)
x <- seq(-4, 4, length = 100)
y <- dnorm(x) + noise
z <- data.frame(x, y)
# Rectangular smoothing ======================================================
unweighted <- smooth_rectangular(z, m = 3)
expect_equal_to_reference(unweighted, file = "_snaps/smooth_rectangular.rds")
# plot(z, type = "l", xlab = "", ylab = "")
# lines(unweighted, type = "l", col = "red")
# Triangular smoothing =======================================================
weighted <- smooth_triangular(z, m = 3)
expect_equal_to_reference(weighted, file = "_snaps/smooth_triangular.rds")
# plot(z, type = "l", xlab = "", ylab = "")
# lines(weighted, type = "l", col = "red")
# Loess smoothing ============================================================
loess <- smooth_loess(z, span = 0.2)
expect_equal_to_reference(loess, file = "_snaps/smooth_loess.rds")
# plot(z, type = "l", xlab = "", ylab = "")
# lines(loess, type = "l", col = "red")
# Savitzky–Golay filter ======================================================
savitzky <- smooth_savitzky(z, m = 21, p = 2)
expect_equal_to_reference(savitzky, file = "_snaps/smooth_savitzky.rds")
# plot(z, type = "l", xlab = "", ylab = "")
# lines(savitzky, type = "l", col = "red")
# Whittaker smoothing ========================================================
whittaker <- smooth_whittaker(z, lambda = 1000, d = 3, sparse = FALSE)
expect_equal_to_reference(whittaker, file = "_snaps/smooth_whittaker.rds")
# plot(z, type = "l", xlab = "", ylab = "")
# lines(whittaker, type = "l", col = "red")
if (requireNamespace("Matrix", quietly = TRUE)) {
whittaker_sparse <- smooth_whittaker(z, lambda = 1000, d = 3, sparse = TRUE)
expect_equal(whittaker, whittaker_sparse)
# plot(z, type = "l", xlab = "", ylab = "")
# lines(whittaker_sparse, type = "l", col = "red")
}
}
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.