multiscale.rips: Multiscale Rips Persistent Homology

Description Usage Arguments Value Author(s) Examples

Description

A multiscale approach to compute Rips persistent homology.

Usage

1
multiscale.rips.ipca( X, epsGMRA,  maxD, single = FALSE )

Arguments

X

Data set. Each row is a data point

epsGMRA

Accuracy of GMRA approximation

maxD

Highest dimensionality of simplicies to insert into the filtration

single

Compute single scale TRUE or FALSE

Value

Persistent diagram, matrix with first column birth, second column death times and a third colum at which scale the death happened.

Author(s)

Samuel Gerber

Examples

1
2
3
4
5
6
7
8
9
##Noisy circle example
phi <- runif(100)*2*pi
arc <- cbind(cos(phi), sin(phi)) * (1+rnorm(length(phi)) * 0.1)

dgm <- multiscale.rips.ipca(X = arc, epsGMRA = 0, maxD = 2);

plot(x = dgm[,1], y = dgm[,2], xlim = c(0, max(dgm[,1:2]) ), ylim = c(0, max(dgm[,1:2])) ,
    asp=1, pch=19)
abline(0, 1)

mph documentation built on May 29, 2017, 12:03 p.m.

Related to multiscale.rips in mph...