DecomposeR (version 1.0.6)

check.emd: Check an EMD object

Description

Provides an ensemble of check on the quality of a decomposition presented as an emd object (see as.emd for more information)

Usage

check.emd(emd, xy = NULL, timelimit = 15)

Arguments

emd

an amd object to test

xy

the original signal that was decomposed: this parameter is simply to insure that you are indeed comparing the decomposition to the original signal, and not cheating by providing the sum of your decomposition

timelimit

a time limit for the computation of the greatest common rational divisor. A too long time may be indicative of a problem, typically depth/time values that are not rounded adequately.

Examples

set.seed(50)

h <- rnorm(n = 1000)

dt <- seq_len(length(h))

alpha <- 0.95

for(i in dt[-1]) h[i] <- alpha *  h[i-1] + h[i]

set.seed(42)

em <- extricate(h, dt, nimf = 7, repl = 1, comb = 100, sifting = 4,
                factor_noise = 20, unit_noise = "native", speak = TRUE)

if (FALSE) {
plot_emd(em, adapt.axis = TRUE)}

check.emd(em, h)