if (FALSE) {
if (requireNamespace("adephylo", quietly = TRUE) & requireNamespace("ape", quietly = TRUE)) {
data(carni70)
carni70.phy <- newick2phylog(carni70$tre)
plot(carni70.phy)
size <- scalewt(log(carni70$tab))[,1]
names(size) <- row.names(carni70$tab)
symbols.phylog(carni70.phy,size)
tre <- ape::read.tree(text = carni70$tre)
adephylo::orthogram(size, tre = tre)
yrange <- scalewt(carni70$tab[,2])
names(yrange) <- row.names(carni70$tab)
symbols.phylog(carni70.phy,yrange)
adephylo::orthogram(as.vector(yrange), tre = tre)
if(adegraphicsLoaded()) {
g1 <- s.label(cbind.data.frame(size, yrange), plabel.cex = 0)
g2 <- addhist(g1)
} else {
s.hist(cbind.data.frame(size, yrange), clabel = 0)
}
}}
Run the code above in your browser using DataLab