ses.mpd {picante}R Documentation

Standardized effect size of mpd

Description

Standardized effect size of mean pairwise distances in communities. When used with a phylogenetic distance matrix, equivalent to -1 times the Nearest Relative Index (NRI).

Usage

ses.mpd(samp, dis, null.model = c("taxa.labels", "sample.pool",
    "phylogeny.pool", "weighted.sample.pool"), runs = 99)

Arguments

samp Community data matrix
dis Distance matrix (generally a phylogenetic distance matrix)
null.model Null model to use
taxa.labels
Shuffle distance matrix labels (across all taxa included in distance matrix)
sample.pool
Randomize community data matrix by drawing species from pool of species occurring in at least one community (sample pool) with equal probability
phylogeny.pool
Randomize community data matrix by drawing species from pool of species occurring in the distance matrix (phylogeny pool) with equal probability
weighted.sample.pool
Randomize community data matrix by drawing species from sample pool with probability weighted by occurrence frequency
runs Number of randomizations

Value

A data frame of results for each community

ntaxa Number of taxa in community
mpd.obs Observed mpd in community
mpd.rand.mean Mean mpd in null communities
mpd.rand.sd Standard deviation of mpd in null communities
mpd.obs.rank Rank of observed mpd vs. null communities
mpd.obs.z Standardized effect size of mpd vs. null communities (= (mpd.obs - mpd.rand.mean) / mpd.rand.sd, equivalent to -NRI)
mpd.obs.p P-value (quantile) of observed mpd vs. null communities (= mpd.obs.rank / runs + 1)
runs Number of randomizations

Warning

Weighted null model currently only works with presence-absence data. Convert your data to presence-absence before using this null model (e.g. decostand(x,method="pa"))

Author(s)

Steve Kembel <skembel@berkeley.edu>

References

Webb, C.O., Ackerly, D.D., and Kembel, S.W. 2008. Phylocom: software for the analysis of phylogenetic community structure and trait evolution. Version 4.0.1. http://www.phylodiversity.net/phylocom/.

See Also

mpd

Examples

data(phylocom)
ses.mpd(phylocom$sample, cophenetic(phylocom$phylo),null.model="taxa.labels")

[Package picante version 0.2-0 Index]