Calculates measures of community phylogenetic structure (quantile regression between co-occurrence and phylogenetic distance) to patterns expected under various null models
comm.phylo.qr(samp, phylo, metric = c("cij", "checkerboard", "jaccard", "doij"),
null.model = c("sample.taxa.labels", "pool.taxa.labels",
"frequency", "richness", "independentswap","trialswap"),
quant = 0.75, runs = 999, show.plot = FALSE, ...)Community data matrix
Phylogenetic tree
Metric of co-occurrence to use (see species.dist)
Null model to use (see Details section for description)
Quantile of slope to be fit (using rq)
Number of runs (randomizations)
Option to display a plot of co-occurrence versus phylogenetic distance with quantile regression slope fit
Additional arguments to randomizeMatrix
A list with elements:
Observed co-occurrence/phylogenetic distance quantile regression intercept
Observed co-occurrence/phylogenetic distance quantile regression slope
P-value of observed quantile regression slope significance versus null model (calculated based on comparison with randomizations)
Rank of observed quantile regression slope vs. random
Number of runs (randomizations)
A vector of quantile regression slopes calculated for each randomization
This function fits a quantile regression of co-occurrence versus phylogenetic distances separating species, and compares observed patterns to the patterns expected under some null model. The quantile regressions are fit using the rq function from the quantreg package.
Currently implemented null models (arguments to null.model):
Shuffle phylogeny tip labels (only within set of taxa present in community data)
Shuffle phylogeny tip labels (across all taxa included in phylogenetic tree)
Randomize community data matrix abundances within species (maintains species occurence frequency)
Randomize community data matrix abundances within samples (maintains sample species richness)
Randomize community data matrix maintaining species occurrence frequency and site richnessing using independent swap
Randomize community data matrix maintaining species occurrence frequency and site richnessing using trial swap
Cavender-Bares J., D.A. Ackerly, D. Baum and F.A. Bazzaz. 2004. Phylogenetic overdispersion in Floridian oak communities, American Naturalist, 163(6):823-843. Slingsby, J. A. and G. A. Verboom. 2006. Phylogenetic relatedness limits coexistence at fine spatial scales: evidence from the schoenoid sedges (Cyperaceae: Schoeneae) of the Cape Floristic Region, South Africa. The American Naturalist 168:14-27.
# NOT RUN {
data(phylocom)
comm.phylo.qr(phylocom$sample, phylocom$phylo, metric="cij",
null.model="sample.taxa.labels", runs=99)
# }
Run the code above in your browser using DataLab