pflamelet (version 0.1.1)

mpbandwidth: Maximum Persistence Bandwidth

Description

Select the bandwidth that maximises the persistence of the scale-family, highlighting the most prominent topological feature of the KDE.

Usage

mpbandwidth(flamelet, scale.param)

Arguments

flamelet

a kxm matrix corresponding to the Persistence Flamelet.

scale.param

a vector of lenght m corresponding to the values of the scale parameter at which the Flamelet X has been evaluated.

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))
mpband(flamelet = foo.flamelet, scale.param = seq(0.01, 1, length.out = 40))
# }