# NOT RUN {
# load muscle fiber data
data(muscle_fiber_dat)
# input variables
Xlim <- c(-50,350)
Ylim <- c(-50,250)
lim <- cbind(Xlim, Ylim)
by <- 6
spseq <- seq(2,40,length.out = 9)
# compute persistence terrace
muscle_fiber_pt=computept(muscle_fiber_dat,sp=spseq,lim=lim,by=by)
# }
# NOT RUN {
# compute persistence terrace with parallel option
spseq <- seq(2,40,length.out = 30)
two_circle_density_pt <- computept(muscle_fiber_dat,sp=spseq,lim=lim,by=by,par=TRUE)
# }
# NOT RUN {
# draw terrace area plot
terracearea(muscle_fiber_pt,dimension=1,maxheight=20)
# draw persistence terrace
plotpt(muscle_fiber_pt,cmax=12,dimension=1)
# }
Run the code above in your browser using DataLab