R/plot.barcode.R

In BenjaminDelory/topotree: Topological Analysis of Tree Architectures using Persistent Homology

plot.barcode<-function(x, xlab=NULL, ylab=NULL, ...){
  if (is.null(xlab)==TRUE){xlab<-"Geodesic distance"}
  if (is.null(ylab)==TRUE){ylab<-expression(H[0])}
  xcoord<-c(min(x[,2:3]), max(x[,2:3]))
  ycoord<-c(0, nrow(x))
  plot(xcoord, ycoord, type="n", xlab=xlab, ylab=ylab, ...)
  segments(x0=x[,"birth"], y0=c(1:nrow(x)), x1=x[,"death"], y1=c(1:nrow(x)), ...)}