Learn R Programming

phytools (version 0.2-40)

sim.ratebystate: Conduct simulation of state dependent rate variation

Description

This function attempts to simulate two characters under a model in which the rate of evolution for the second (y) depends on the states for the first x. See ratebystate for more details.

Usage

sim.ratebystate(tree, sig2x=1, sig2y=1, beta=c(0,1), ...)

Arguments

tree
phylogenetic tree.
sig2x
variance of the Brownian process of evolution for x.
sig2y
variance of the Brownian process of evolution for y when x-min(x)==1 (for logarithm=FALSE) or x==0 (for logarithm=TRUE).
beta
intercept and slope of the relationship between the value of x and the Brownian rate in y.
...
optional arguments which include method ("by.node" or "by.branch" indicating whether to assume the rate varies as a function of the node state or the mean branch state); plot, a logical value indicating

Value

  • This function returns a matrix.

See Also

fastBM, ratebystate