
print(santatraits)
heft transport first.appearance
Coca.cola.Santa fat reindeer 1900
Department.Store.Santa fat reindeer 1900
Salvation.Army.Santa fat reindeer 1890
Kriss.Kringle fat reindeer 1800
Santa.Claus fat reindeer 1700
Pere.Noel thin donkey 1400
Father.Christmas fat foot 1400
Sinterklaas thin horse 400
Saint.Nicholas thin foot 300
Ded.Moroz thin foot 1937
We are fortunate to also have sequence data (from DNA on milk glasses and lost beard hairs), so we can add additional information about relatedness amongst these species.
plot(xmas.tree, type = "c", FALSE, edge.color="darkgreen", edge.lty=1, edge.width=18, label.offset = 1, direction="downward", font=3, tip.color="darkred")

The phylogeny shows that there seems to be an early divergence between European and North American santas. Indeed, there is a group of North American santas (Mall Santa, Coca-cola Santa, Salvation Army Santa) which are closely related (and also appear to share very similar traits, based on the table above). (Note that branch lengths in this phylogeny show nucleotide substitutions, and it is not time-calibrated, due to the absence of santa fossils).
One approach is to identify a few traits and plot them on the phylogeny to compare how traits vary among santas. Let's start with anatomical characteristics:
#Plot traits (fatness) against Santa co1 = c("blue", "purple") tiplabels(pch = 12, col = co1[as.factor(heft)], cex = 3.5, adj=c(0.5, 0), lwd=2) #Let's see the transportation mode trait too: co2 <- c("yellow", "gold", "darkorange", "red") tiplabels(pch = 8, col = co2[as.factor(transport)], cex = 2, adj=c(0.5, 0), lwd=2) #Legends legend("topleft", legend=c("reindeer", "donkey", "foot", "horse"), fill=rev(c("yellow", "gold", "darkorange", "red"))) legend("topright", legend=c("fat", "thin"), fill=c("blue", "purple"))

For a future study, we could ask whether the apparent correlation between fatness and reindeer usage is significant, once the underlying phylogenetic relationships were controlled for.
We can also reconstruct santa traits (here, we look at the form of transportation) to explore what form of transportation ancestral santas likely used:
#reconstruct ancestral state cc = ace(transport, xmas.tree, type="discrete") co2 = c("yellow", "gold", "darkorange", "red") nodelabels(pie = cc$lik.anc, piecol = co2, cex = c(1.5, rep(1, 8)))
The markers at each node show the probability that this ancestral taxa used each of the four possible types of transportation. It seems that the North American santas and their ancestors have long relied on a trusty reindeer mutualism.
Finally, we can look at the geography of all the various santas:
phylo.to.map(xmas.tree, locales)

To run this caRd yourself, follow the link to the R code: https://gist.github.com/cmtucker/8e5677bdd5c409d70738









