arbre phylogénétique

Code R de l'exemple ggplot2

library(ggplot2)
library(ape)
library(ggtree)

data(bird.orders)
## The root splits the tree in two: ratites and fowl, then all other orders.
fowl <- c("Struthioniformes", "Tinamiformes", "Craciformes", "Galliformes", "Anseriformes")
tree <- groupOTU(bird.orders, list(fowl = fowl))

p <- ggtree(tree, aes(colour = group), linewidth = 0.4) +
  geom_tiplab(size = 2.3, offset = 0.3, family = "Quicksand") +
  annotate("text", x = 1, y = 0, label = "Ratites and fowl", colour = "#9b2367",
           hjust = 0, size = 3) +
  annotate("text", x = 1, y = 24, label = "The other 18 orders", colour = "#2e7891",
           hjust = 0, size = 3) +
  scale_colour_manual(values = c(`0` = "#2e7891", fowl = "#9b2367")) +
  scale_x_continuous(expand = expansion(mult = c(0.02, 0.32))) +
  labs(title = "The root splits five bird orders from eighteen",
       subtitle = "Orders of birds, branch lengths from how much their DNA differs",
       caption = "Data: Sibley & Ahlquist 1990, via the ape R package") +
  theme_lex_viz(grid = "", arrows = "") +
  theme(legend.position = "none", axis.text = element_blank(), axis.title = element_blank(),
        axis.line = element_blank(), axis.ticks = element_blank())