data(geospiza)
attach(geospiza)
name.check(geospiza.data, geospiza.tree)->r
drop.tip(geospiza.tree, r[[1]])->g.tree
data<-geospiza.data[,1]
names(data)<-rownames(geospiza.data)
ape2ouch(g.tree, data)->geospiza.ouch
#Interface with ouch package
library(ouch)
brown.fit(geospiza.ouch$d, geospiza.ouch$node, geospiza.ouch$ancestor, geospiza.ouch$time)
Run the code above in your browser using DataLab