###Code for pic example in R

mytree="((((Homo:0.21,Pongo:0.21):0.28,Macaca:0.49):0.13,Ateles:0.62):0.38,Galago:1.00);"
tree.primates=read.tree(text=mytree)
body=c(4.09434, 3.61092, 2.37024, 2.02815, -1.46968)
longevity=c(4.74493,3.3322,3.3673,2.89037,2.30259)
par(mfrow=c(1,2))
plot(body,longevity,pch=15)

names(body)=names(longevity)=c("Homo","Pongo","Macaca","Ateles","Galago")
pic.body=pic(body,tree.primates)
pic.longevity=pic(longevity,tree.primates)
pic.body
pic.longevity
plot(pic.body,pic.longevity,pch=15)

plot(tree.primates)
nodelabels(round(pic.body,3),adj=c(0,-0.5),frame="n")
nodelabels(round(pic.longevity,3),adj=c(0,1),frame="n")

