ggdensitree + revts

[thijsjanzen]

I would like to make a plot as hinted here: #87
e.g. a densitree with a timescale with 0 at the tips, and -crown_age at the root.

However, the following code does not work as expected, and only reverts one of the two tree plots.

trees <- list()
trees[[1]] <- rtree(n = 10)
trees[[2]] <- rtree(n = 10)
class(trees) <- "multiPhylo"
p <- ggdensitree(trees) + theme_tree2()
p2 <- revts(p)
multiplot(p, p2, ncol=2)

[brj1]

I have built a workaround for this in 5a70208. ggdensitree will set the maximum tip date to 0 when align.tips = TRUE (the default). So,

trees <- list()
trees[[1]] <- rtree(n = 10)
trees[[2]] <- rtree(n = 10)
class(trees) <- "multiPhylo"
p <- ggdensitree(trees) + theme_tree2()
p

will produce the plot you want, with no revts required.

[thijsjanzen]

Works like a charm! Also, works for all densitrees