## Load the Beck & Lee 2014 data
data(BeckLee_mat50) ; data(BeckLee_tree)
## Run a simple disparity through time analysis (with three time bins)
result <- dispRity.through.time(BeckLee_mat50, BeckLee_tree, 3)
summary(result) ; plot(result)
## This is equivalent to run the following decomposed code
dispRity(boot.matrix(chrono.subsets(BeckLee_mat50, BeckLee_tree, time = 3, method = "discrete"),
bootstraps = 100),
metric = c(median, centroids))
Run the code above in your browser using DataLab