fpath <- system.file("extdata", "ecospat.testTree.tre", package="ecospat")
library(ape)
tree <-ape::read.tree(fpath)
data <- ecospat.testData[9:52]
pd <- ecospat.calculate.pd(tree, data, method = "spanning", type = "species", root = FALSE,
average = FALSE, verbose = TRUE )
plot(pd)
Run the code above in your browser using DataLab