sim.Stasis: Simulate Stasis time-series for trait evolution

In paleoTS: Analyze Paleontological Time-Series

Simulate Stasis time-series for trait evolution

Description

Simulate Stasis time-series for trait evolution

Usage

sim.Stasis(
  ns = 20,
  theta = 0,
  omega = 0,
  vp = 1,
  nn = rep(20, ns),
  tt = 0:(ns - 1)
)

Arguments

ns	number of populations in the sequence
theta	mean of populations
omega	variance among populations
vp	phenotypic variance within populations
nn	vector of population sample sizes
tt	vector of population times (ages)

Value

a paleoTS object

See Also

sim.GRW, sim.OU, as.paleoTS

Examples

x <- sim.Stasis(omega = 0.5, vp = 0.1)
w.sta <- fitSimple(x, model = "Stasis")
w.ss  <- fitSimple(x, model = "StrictStasis")
compareModels(w.sta, w.ss)

paleoTS documentation built on Aug. 9, 2022, 1:06 a.m.