library(paran) library(MASS) library(ape) library(splits) #Single threshold GMYC model tr <- read.nexus("cox1.nexus") result <- gmyc(tr) #Summary results summary(result) plot(result) #Node support GMYC_support <- gmyc.support(result) is.na(GMYC_support[GMYC_support ==0]) <- TRUE plot(result, cex=.4, no.margin = TRUE) nodelabels(round(GMYC_support, 2), cex= .6) #Species-sample list spec.list(result)