data(pwed_pd)
data(spcolmap)
pwed_pd <- matrix2list.2(pwed_pd)
cladeI <- spcolmap$species[spcolmap$host %in% "M.buchanani"]
#We just want to look at distance between LM1 and LM3 in for dorsal anchor
boxplotSort(lapply(pwed_pd, function(k) k[,which(colnames(k)=="D1_3")]),
italic=TRUE, col=c("dodgerblue","violetred"), clade=cladeI,
ylab=expression(paste("Length ", "(", italic(mu),"m", ")")))
#Separation of two lineage seems possible at 15 micrometers
abline(h=15)
Run the code above in your browser using DataLab