pflamelet (version 0.1.1)

flamelet.plot: Plot Persistence Flamelet

Description

Plot the Persistence Flamelet in its original 3-dimensional form, or in a 2-dimensional projection.

Usage

flamelet.plot(
  flamelet,
  scale.param,
  tseq,
  flat = FALSE,
  scale.name = "Bandwidth",
  band = NULL
)

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 has been evaluated.

tseq

a vector of length k containing of values at which the Flamelet function is evaluated for a fixed scale level.

flat

a logical denoting whether the plot should be a 2-d projection of the Flamelet (TRUE) or a 3-d object (FALSE).

scale.name

name of the scale parameter.

band

scalar representing the confidence band for Persistence Flamelet. Only available when flat = FALSE.

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))
flamelet.plot(foo.flamelet, scale.param = seq(0.01, 1, length.out = 40),
             tseq = seq(0, .75, length.out = 500) )


foo.band = flamelet.band(X = xx, B = 10, alpha = 0.05,
                   tseq = seq(0, .75, length.out = 500), diag.fun = kde,
                   h.grid = seq(0.01, 1, length.out = 40), lim = lim, by = by)

flamelet.plot(foo.flamelet, scale.param = seq(0.01, 1, length.out = 40),
             tseq = seq(0, .75, length.out = 500), band = foo.band)
# }