kernelTDA (version 1.0.0)

wass.kernel: L_infty q-Wasserstein Kernel (WK)

Description

Computes the L_infty q-Wasserstein Kernel (WK) between persistence diagrams.

Usage

wass.kernel(d1, d2 = NULL, h, dimension, q)

Arguments

d1

A persistence diagram (matrix with 3 col where the first one is the dimension, the second is the birth-time and the third is the death-time) or a list of diagrams

d2

A persistence diagram (matrix with 3 col where the first one is the dimension, the second is the birth-time and the third is the death-time).

h

bandwidth of the kernel.

dimension

The dimension of the topological feature (0 for connected components, 1 for cycles etc).

q

order of the q-Wasserstein distance.

Value

If d1 is a list of Persistence Diagrams, this function returns a matrix whose (i,j) entry is the WK computed in (d1[[i]], d2[[j]]), otherwise it returns the value for the WK computed in (d1, d2).

Examples

# NOT RUN {
diag1 <- matrix(c(1,1,1,0,2,3,2,2.5,4), ncol = 3, byrow = FALSE)
diag2 <- matrix(c(1,1,0,1,1,2), ncol = 3, byrow = FALSE)
wass.kernel(d1 = diag1, d2 = diag2, h = 1, dimension = 1, q = 2)
# }