Is it possible to bind two simmaps together? When I attempt to run the
following code, the call to bind.tree doesn't throw any errors but I get an
error when I attempt to plot the resulting tree:
library(ape)
library(phytools)
t1 <- pbtree(n=100, scale=5)
t2 <- pbtree(n=100, scale=5)
Q1 <- matrix(c(-0.1,0.1,0.1,-0.1), 2,2)
Q2 <- matrix(c(-1,1,1,-1), 2,2)
s1 <- sim.history(t1, Q1)
s2 <- sim.history(t2, Q2)
out <- bind.tree(s1, s2)
plot(out)
> Error in tree$mapped.edge[ii, ] : subscript out of bounds
Thanks,
- Christie Ziegler
[[alternative HTML version deleted]]
_______________________________________________
R-sig-phylo mailing list - [email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/[email protected]/
