Hi all,
I'm trying to use geiger's rjmcmc.bm function for some continuous trait
evolution analyses, and running into an error. Here's code to simulate a
tree with trait data, then run the rjmcmc.bm analysis:
phy<-sim.bdtree(b=1,d=0, n=100, t=20, extinct=TRUE)
phy$tip.label=paste("sp",1:100,sep="
Hello,
sorry, I forgot to refer where the problem was. On the beginning, I didn't
check all over thousand trees manually. I found trees on which it was failing
by simple for loop (could be simpler and more straightforward, I know):
Y <- trees.rooted.nexus
class(Y) <- "multiPhylo"
X <- list()
tree