Generates species trait data under models that deviate from Brownian Motion. Supports OU (stabilising selection), regime shifts (clade-specific optima), and non-linear inter-trait correlations.
Usage
simulate_non_bm(
tree,
n_traits = 4,
scenario = c("OU", "regime_shift", "nonlinear"),
alpha = 2,
theta = 0,
sigma = 1,
shift_magnitude = 2,
seed = NULL
)Arguments
- tree
Object of class
"phylo".- n_traits
Integer. Number of traits to simulate (default 4).
- scenario
Character. One of
"OU","regime_shift","nonlinear".- alpha
Numeric. OU pull strength (only for
"OU").- theta
Numeric. OU optimum (only for
"OU").- sigma
Numeric. BM diffusion rate.
- shift_magnitude
Numeric. Regime shift size in SD units (only for
"regime_shift").- seed
Integer seed or
NULL.