pd: Calculate Faith's Phylogenetic Diversity

Description Usage Arguments Value Warning Note Author(s) References See Also Examples


Calculate the sum of the total phylogenetic branch length for one or multiple samples.


pd(samp, tree, include.root=TRUE)



Community data matrix


A phylo tree object


Should the root node be included in all PD calculations (default = TRUE)


Returns a dataframe of the PD and species richness (SR) values for all samples


If the root is to be included in all calculations (include.root=TRUE), the PD of all samples will include the branch length connecting taxa in those samples and the root node of the supplied tree. The root of the supplied tree may not be spanned by any taxa in the sample. If you want the root of your tree to correspond to the most recent ancestor of the taxa actually present in your sample, you should prune the tree before running pd:

prunedTree <- prune.sample(sample,tree)


The data sets need not be species-community data sets but may be any sample data set with an associated phylogeny. PD is not statistically independent of species richness, it positively correlates with species richness across samples. The function ses.pd compares observed PD to the values expected under various randomizations and allows a way to standardize for unequal richness across samples.

If the root is to be included in all calculations of PD (include.root=TRUE), the tree must be rooted. Single-species samples will be assigned a PD value equal to the distance from the root to the present.

If the root is not included in all calculations by default (include.root=FALSE), the tree need not rooted, but in the case of single-species samples the PD will be equal to NA and a warning will be issued.


Matthew Helmus mrhelmus@gmail.com, Jonathan Davies davies@nceas.ucsb.edu, Steven Kembel steve.kembel@gmail.com


Faith D.P. (1992) Conservation evaluation and phylogenetic diversity. Biological Conservation, 61, 1-10.

See Also

psr, ses.pd



picante documentation built on July 1, 2020, 10:57 p.m.