Я пытаюсь произвести реконструкцию предков, используя пакет ape and phytools в Rstudio. Моя проблема в том, что в моем филогенетическом дереве ярлыки на подсказках / названия видов переполнены и неразборчивы. В настоящее время в моем дереве есть набор данных из 262 видов.
Здесь можно найти пример файла нексуса с данными.
Дерево реконструкции предков, которое я создал, находится здесь: https://i.imgur.com/WFoEu7S.png.
Каждый вид имеет состояние символа 0 или 1 и имеет метки узлов и наконечников, адресованные каждому состоянию. В конце концов, я хотел бы раскрасить ветки соответствующим состоянием символа (которое у меня было красным или черным).
В идеале я хочу создать не ультраметрическое дерево, подобное предыдущему вопросу о переполнении стека по этой ссылке здесь.
Я безуспешно пытался реализовать код R из этой ссылки для своего собственного дерева.
Ниже мой код на R. Я все еще изучаю R, не знаком с некоторыми методами построения графиков и подозреваю, что это может быть проблемой здесь:
tree = read.nexus("test_nexus")
dichot_tree = multi2di(tree)
dichot_tree$edge.length<-runif(n=nrow(dichot_tree$edge),min=0,max=1)
dichot_tree$edge.length[dichot_tree$edge.length<1]<-1
domest = read.nexus.data("test_nexus")
aceDISCRETE<-ace(as.numeric(domest), dichot_tree, type="discrete")
plot(dichot_tree, cex=0.5, label.offset=1, no.margin=TRUE)
tiplabels(pch=22, bg=as.numeric(domest),cex=1, adj=1)
nodelabels(pie=aceDISCRETE$lik.anc, piecol=c("black", "red"), cex=0.25)