make_relprob_txtmatrix_sp1: Convert a observed-speciation transition matrix to an...

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

View source: R/BioGeoBEARS_basics_v1.R

Description

Convert a cladogenesis/speciation transition matrix (specifying the probability of each Left/Right descendant range pair, conditional on each ancestral state) of dimensions numstates by numstates^2 to a square transition matrix of dimensions numstates by numstates, representing the probability of a transition when only one daughter survives in the tree.

Usage

1
  make_relprob_txtmatrix_sp1(probmat, split = "\\|")

Arguments

probmat

A matrix of text, describing each of the allowed range-inheritance events. Assumes that column names are in the "A|B" format.

split

The value to split Left/Right pairs on (e.g., "A|B" –> "A", "B")

Details

This matrix could be used to quantify the probability of range-change along a branch due to unobserved speciation events; all that would be required would be an estimate of the number of unobserved speciation events on the branch, and treating this as a Poisson process. (Note: this assumes that the probability of either branch surviving is identical, which might not be the case. See the GeoSSE (Goldberg et al. (2011)) and ClaSSE ("Goldberg et al. (2012)) for the beginnings of work on this, with 2 and 3 geographic areas, respectively.

Value

newmat A new square matrix.

Note

Go BEARS!

Author(s)

Nicholas J. Matzke matzke@berkeley.edu

References

http://phylo.wikidot.com/matzke-2013-international-biogeography-society-poster https://code.google.com/p/lagrange/ http://tigger.uic.edu/~eeg/code/code.html

Goldberg_etal_2011_GeoSSE

Goldberg_Igic_2012_ClaSSE

Matzke_2012_IBS

ReeSmith2008

See Also

make_relprob_matrix_bi, make_relprob_nummatrix_sp1

Examples

1
2
3
4
5
6
7
8
testval=1
probmat = make_relprob_matrix_bi(states_list=list("_", c("A"),
c("B"), c("C"), c("A","B"), c("B","C"), c("A","C"), c("A","B","C")),
split_ABC=FALSE, splitval="", code_for_overlapping_subsets=NA, printwarn=1)
probmat

newmat = make_relprob_txtmatrix_sp1(probmat=probmat, split="\\\\|")
newmat

BioGeoBEARS documentation built on May 29, 2017, 8:36 p.m.