data(anoletree)
anole_tree<-as.phylo(anoletree)
data(anole.data)
anole_data<-cbind(phyl.resid(anole_tree,
x=as.matrix(anole.data[,"SVL",drop=FALSE]),
Y=as.matrix(anole.data[,c(6,4,2)]))$resid,
exp(anole.data[,"SVL",drop=FALSE]))
plotTree.lollipop(anole_tree,anole_data,
ylab=c("relative (TL)","relative (FLL)",
"relative (HL)","SVL"))
par(mar=c(5.1,4.1,4.1,2.1))
Run the code above in your browser using DataLab