gaus.kernel: Geodesic Gaussian Kernel (GGK)

In kernelTDA: Statistical Learning with Kernel for Persistence Diagrams

Description

Computes the Geodesic Gaussian Kernel (GGK) between persistence diagrams.

Usage

gaus.kernel(d1, d2 = NULL, h, dimension)

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).

Value

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

Author(s)

Tullia Padellini

References

\insertRef

padellini2017supervisedkernelTDA

Examples

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)
gaus.kernel(d1 = diag1, d2 = diag2, h = 1, dimension = 1)

kernelTDA documentation built on April 19, 2020, 3:56 p.m.