if(require("TDAstats"))
{
# create a persistence diagram from a 2D Gaussian
df = data.frame(x = rnorm(n = 20,mean = 0,sd = 1),y = rnorm(n = 20,mean = 0,sd = 1))
# compute persistence diagram with calculate_homology from package TDAstats
phom_TDAstats = TDAstats::calculate_homology(mat = df,dim = 0,threshold = 1)
# convert to data frame
phom_TDAstats_df = diagram_to_df(d = phom_TDAstats)
}
Run the code above in your browser using DataLab