data(anoletree)
anc<-findMRCA(anoletree,c("cristatellus","cooki",
"gundlachi"))
plotTree(anoletree,type="fan",fsize=0.7,lwd=1)
nodelabels(node=anc,frame="circle",pch=21,cex=1.5,
bg="blue")
legend("topleft","most recent common ancestor\nof Puerto Rican TG anoles",
pch=21,pt.cex=1.5,pt.bg="blue",cex=0.7,bty="n")
par(mar=c(5.1,4.1,4.1,2.1)) ## reset margin to default
Run the code above in your browser using DataLab