setwd('D:/R Stuff') library(coda) library(ape) library(phytools) library(BAMMtools) Bats<-read.nexus(file.choose()) SpAvg <- read.csv("FinalData.csv", stringsAsFactors = FALSE) ### Tree Check is.ultrametric(Bats) is.binary.tree(Bats) min(Bats$edge.length) ### Data Prep ### setwd('D:/Bamm') #Centroid count<-1 while(count