windex.sim.test: Tests for presence of convergent evolution (over and above...

In windex: Analysing Convergent Evolution using the Wheatsheaf Index

Tests for presence of convergent evolution (over and above that expected from Brownian motion)

Description

The P-value returned is for the null hypothesis that the calculated Wheatsheaf index is no higher than expected for traits evolving under Brownian motion (parameterised with rates of evolution and trait covariances estimated from the original traits).

Usage

windex.sim.test(dat, tree, traits, focal = dat[, 2], SE = TRUE, Nsims,
plot = TRUE, fossil = FALSE, ...)

Arguments

dat	A dataframe containing a column of 0s and 1s to denote non-focal and focal taxa respectively, and columns of trait data which you which to plot. The first column must be named 'species' and contain species names that correspond to those in the phylogenetic tree.
tree	Phylogenetic tree of class 'phylo' containing branch lengths. The tree should also be ultrametric.
traits	Column numbers (or names) for the traits for which you want to calculate a Wheatsheaf index.
focal	Column in the dataframe containing the focal designations.
SE	A logical specifying whether to standardise the traits by their standard error across species, default is SE=TRUE.
Nsims	Number of simulations on which to base the P-value.
plot	A logical indicating whether to plot the simulated distribution. If TRUE, a histogram is plotted with the calculated Wheatsheaf index and its 95% confidence interval overlayed on the histogram as a solid and dashed lines (respectively).
fossil	A logical specifying whether tree is ultrametric (for instance to to presence of fossil tips), in which case a different phylogenetic distance penalty based on shared branch length rather than shared time before divergence will be used instead.
...	Additional arguments passed to hist() to customise the histogram (when plot=TRUE).

Value

w	Calculated Wheatsheaf index
low95	Lower bound of 95% confidence interval for the Wheatsheaf index obtained by jackkniving
up95	Upper bound of 95% confidence interval for the Wheatsheaf index obtained by jackkniving
P	P-value obtained from comparing observed Wheatsheaf index to simulations under Brownian motion on the phylogenetic tree
sim.dist	Wheatsheaf indices of simulated datasets used to calculate P-value

Author(s)

Kevin Arbuckle

References

Arbuckle, K., Bennett, C.M. and Speed, M.P. 2014. A simple measure of the strength of convergent evolution. Methods in Ecology and Evolution 5:685 - 693.

See Also

windex test.windex

Examples

## Not run: 
data(sample.data)
data(sample.tree)
windex.sim.test(sample.data,sample.tree,traits=c("bm1","bm2"),focal=sample.data[,2],Nsims=1000,
plot=TRUE,col="light grey")

## End(Not run)

windex documentation built on Oct. 11, 2023, 5:16 p.m.