Calculates mean pairwise distance separating taxa in each pair of communities in a given community matrix.
bmpd(comm, pd, abundance.weighted = TRUE, na.zero = TRUE,
time.output = FALSE, unit.sum = NULL)
matrix or data.frame, community data matrix, rownames are sample names, colnames are OTU ids.
matrix, pairwise phylogenetic distance matrix.
logic, whether weighted by species abundance, default is TRUE, means weighted.
logic. when the sum of a row (a sample) is zero in community data matrix, the relative abundance will be NAN. Sometimes, to avoid some problem in following calculation, this kind of NAN value need be set as zero. Defalt is TRUE.
logic, whether to count calculation time, default is FALSE.
When a beta diversity index is calculated for a bin, the taxa abundances will be divided by unit.sum to calculate the relative abundances. usually, unit.sum can be set as the sequencing depth in each sample. Default setting is NULL, means not to do this special transformation.
Output is a distance object of pairwise betaMPD between samples.
beta mean pairwise distance.
Webb CO, Ackerly DD, and Kembel SW. 2008. Phylocom: software for the analysis of phylogenetic community structure and trait evolution. Bioinformatics 18:2098-2100
# NOT RUN {
data("example.data")
comm=example.data$comm
pd=example.data$pd
bmpd.wt=bmpd(comm, pd, abundance.weighted = TRUE)
# }
Run the code above in your browser using DataLab