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)
Run the code above in your browser using DataLab