install.packages('ape', repos='http://cran.wustl.edu/')
library(ape)
PrimateTrees <- read.tree("PrimateTrees.txt")
ConsensusPrimate <- consensus(PrimateTrees, p=.999)
ordering <- as.numeric(ConsensusPrimate$tip.label)
names = c("Tarsius_syrichta_1",
          "Lemur_catta_2",
          "Homo_sapiens_3",
          "Pan_4",
          "Gorilla_5",
          "Pongo_6",
          "Hylobates_7",
          "Macaca_fuscata_8",
          "M_mulatta_9",
          "M_fascicularis_10",
          "M_sylvanus_11",
          "Saimiri_sciureus_12")
plot(ConsensusPrimate)
tiplabels(names[ordering], c(1:12))