## ----setup, include = FALSE----------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = ">"
)
library("TDAstats")

## ----install-package, eval = FALSE---------------------------------------
#  # install development version of TDAstats - advanced users
#  devtools::install_github("rrrlw/TDAstats")
#  
#  # install TDAstats from CRAN
#  install.packages("TDAstats")
#  
#  # load TDAstats for use
#  library("TDAstats")

## ----load-data-----------------------------------------------------------
# load dataset
data(circle2d)

# look at the dimensions and class of circle2d
class(circle2d)
nrow(circle2d)
ncol(circle2d)

# take a peek at first 6 rows
head(circle2d)

## ----plot-circle2d, fig.width = 4, fig.height = 4.5----------------------
# scatterplot of circle2d
plot(circle2d, xlab = "x", ylab = "y", main = "Point cloud in circle2d dataset")

## ----calc-hom------------------------------------------------------------
# calculate persistent homology
circle.phom <- calculate_homology(circle2d)

# print first 6 features (ordered by dimension and birth)
head(circle.phom)

# print last 6 features (ordered by dimension and birth)
tail(circle.phom)

## ----plot-barcode, fig.height = 4.5, fig.width = 6-----------------------
# plot topological barcode
plot_barcode(circle.phom)

## ----plot-persist, fig.height = 4.5, fig.width = 6-----------------------
# plot persistence diagram
plot_persist(circle.phom)