library(gplots)
library(circular)
data(va_mean)
data(estimated_ancestral_va)
data(spcolmap)
cladeII <- spcolmap$species[spcolmap$host %in% "L.subviridis"]
shapeEvo(va_mean, estimated_ancestral_va, col.lab="dodgerblue",
clade=cladeII, exfac=2, tit="Ventral anchors")
#Some journals want the title to be left-adjusted, so set tit="" and then:
#title("a)", adj=0)
Run the code above in your browser using DataLab