sim.char {geiger} | R Documentation |
Simulates evolution of discrete or continuous characters on a phylogenetic tree
sim.char(phy, model.matrix, nsims = 1, model = "brownian", root.state = 1)
phy |
Phylogenetic tree |
model.matrix |
Matrix describing model: either vcv matrix or q matrix |
nsims |
Number of simulations to run |
model |
"brownian", "speciational", or "discrete" |
root.state |
Starting state at root |
Array of simulated data
Luke J, Harmon
data(geospiza) attach(geospiza) # Multivariage continuous drop.tip(geospiza.tree, "olivacea")->g.tree s<-ic.sigma(g.tree, geospiza.data) sims<-sim.char(g.tree, s, 100) # Discrete character evolution q<-list(rbind(c(-.01, .01), c(.01, -.01))) sims<-sim.char(geospiza.tree, q, model="discrete", n=100)