silhouette: The Persistence Silhouette Function

Description

This function computes the silhouette function corresponding to a given persistence diagram.

Usage

silhouette(Diag, p = 1, dimension = 1, 
           tseq = seq(min(Diag[, 2:3]), max(Diag[, 2:3]), length = 500))

Arguments

Diag

an object of class diagram or a $P$ by $3$ matrix, storing a persistence diagram with columns: "dimension", "Birth", "Death".

number: the power of the weights of the silhouette function. See the definition of silhouette function, Section 5 in the reference.

dimension

the dimension of the topological features under consideration. Default is 1 (loops).

tseq

a vector of values at which the silhouette function is evaluated.

Value

Returns a numeric vector of the same length of tseq, with the values of the silhouette function evaluated at the points of tseq.

References

Frederic Chazal, Brittany T. Fasy, Fabrizio Lecci, Alessandro Rinaldo, and Larry Wasserman, (2014), "Stochastic Convergence of Persistence Landscapes and Silhouettes", Proceedings of the 30th Symposium of Computational Geometry (SoCG). (arXiv:1312.0308)

Examples

Run this code

Diag=matrix(c(0,0,10,1,0,3,1,3,8), ncol=3, byrow=TRUE)
DiagLim=10
colnames(Diag)=c("dim", "Birth", "Death")

#persistence silhouette
tseq=seq(0,DiagLim, length=1000)
Sil=silhouette(Diag, p=1,  dimension=1, tseq)

par(mfrow=c(1,2))
plot.diagram(Diag)
plot(tseq, Sil, type="l", xlab="t", ylab="silhouette", asp=1)