pflamelet (version 0.1.1)

flamelet.clean: Clean Persistence Flamelet

Description

Remove all the values of the Persistence Flamelet which are not significantly different than 0

Usage

flamelet.clean(flamelet, band)

Arguments

flamelet

a kxm matrix corresponding to the Persistence Flamelet.

band

a scalar representing the confidence band.

Value

the value of the bandwidth corresponding to the most persistent feature of the Flamelet.

References

T. Padellini and P. Brutti (2017) Persistence Flamelets: multiscale Persistent Homology for kernel density exploration https://arxiv.org/abs/1709.07097

Examples

# NOT RUN {
library(TDA)
xx = rbind(circleUnif(50, 1), circleUnif(50, 1.5) + 3)
Xlim = c(-1, 5);  Ylim = c(-1, 5);  by = 0.05
lim = cbind(Xlim, Ylim)
foo.flamelet = build.flamelet(X = xx, h.grid = seq(0.01, 1, length.out = 40),
base.type = "landscape", dimension = 1,base.param = 1, lim = lim, by = by,
                            tseq = seq(0, .75, length.out = 500))
foo.band = flamelet.band(X = xx, B = 10, alpha = 0.95,
                   tseq = seq(0, .75, length.out = 500), diag.fun = kde,
                   h.grid = seq(0.01, 1, length.out = 40), lim = lim, by = by)
new.flamelet = flamelet.clean(foo.flamelet, foo.band)

# }