library(ape)
tr1 <- read.tree(text ="(((a:2,(b:1,c:1):1):1,d:3):1,e:4);")
com <- matrix(c(1,0,1,1,0,0,
1,0,0,1,1,0,
1,1,1,1,1,1,
1,0,1,1,0,1,
0,0,0,1,1,0), 6, 5,
dimnames=list(paste0("g",1:6), tr1$tip.label))
collapse_range(com, tr1, n=1)
Run the code above in your browser using DataLab