# NOT RUN {
set.seed(42)
n <- 600
t <- seq_len(n)
p1 <- 30
p2 <- 240
xy <- (1 + 0.6 * sin(t*2*pi/p2)) * sin(t*2*pi/p1) + 2 * sin(t*2*pi/p2) +
rnorm(n, sd = 0.5)
xy <- xy - mean(xy)
inter_dt <- round(runif(length(xy), min = 0.5, max = 1.5),1)
dt <- cumsum(inter_dt)
dec <- extricate(xy, dt, nimf = 7, repl = 1, comb = 40, factor_noise = 10,
speak = TRUE)
integrity(xy, dec)
parsimony(dec)
n.extrema(dec$m, dec$mode)
plot_emd(dec, select = c(6,8,9), pdf = FALSE, adapt.axis = TRUE)
# }
# NOT RUN {
plot_emd(dec, li = list(v = 0), adapt.axis = TRUE, dir = tempdir())
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab