given_a_starting_state_simulate_split

An integer
  index value, between 0 and <code>(numstates-1)</code>, which
  specifies what state is the starting point for the
  branch.

index_Qmat_0based_of_starting_state

A speciation/cladogenesis
  transition matrix, in COO-like form, as produced by
  <code><a rd-options="cladoRcpp" href="/link/rcpp_calc_anclikes_sp_COOweights_faster?package=BioGeoBEARS&version=0.2.1&to=cladoRcpp" data-mini-rdoc="cladoRcpp::rcpp_calc_anclikes_sp_COOweights_faster">rcpp_calc_anclikes_sp_COOweights_faster</a></code>.

COO_probs_columnar

numstates


  This function simulates a biogeographical history during
  a speciation/cladogenesis range inheritance event, given
  a cladogenesis probability transition matrix and a
  starting state.


BioGeoBEARS allows probabilistic inference of both historical biogeography (ancestral geographic ranges on a phylogeny) as well as comparison of different models of range evolution.  It reproduces the model available in LAGRANGE (Ree and Smith 2008), as well as making available numerous additional models. For example, LAGRANGE as typically run has two free parameters, d (dispersal rate, i.e. the rate of range addition along a phylogenetic branch) and e (extinction rate, really the rate of local range loss along a phylogenetic branch). LAGRANGE also has a fixed cladogenic model which gives equal probability to a number of allowed range inheritance events, e.g.: (1) vicariance, (2) a new species starts in a subset of the ancestral range, (3) the ancestral range is copied to both species; in all cases, at least one species must have a starting range of size 1.  LAGRANGE assigns equal probability to each of these events, and zero probability to other events.  BioGeoBEARS adds an additional cladogenic event: founder-event speciation (the new species jumps to a range outside of the ancestral range), and also allows the relative weighting of the different sorts of events to be made into free parameters, allowing optimization and standard model choice procedures to pick the best model. The relative probability of different descendent range sizes is also parameterized and thus can also be specified or estimated.  The flexibility available in BioGeoBEARS also enables the natural incorporation of (1) imperfect detection of geographic ranges in the tips, and (2) inclusion of fossil geographic range data, when the fossils are tips on the phylogeny. Bayesian analysis has been implemented through use of the "LaplacesDemon" package, however this package is now maintained off of CRAN, so its usage is not formally included in BioGeoBEARS at the current time.  CITATION INFO: This package is the result of my Ph.D. research, please cite the package if you use it!  Type: citation(package="BioGeoBEARS") to get the citation information.

given_a_starting_state_simulate_split: Given the state just below a node, simulate the states after speciation

Description

Usage

Arguments

Value

References

See Also

Examples