Hi all,
Does anyone know of a good implementation of an LTT plot that can draw a
Confidence Interval or HPD interval from a set of trees? I've seen things like
`ltt` in phytools that can draw one line for each tree in the sample. However,
this can look a bit messy, and I'd ideally love to just
Hi Simon,
one solution is to extract a vector of branching times and sort this vector
for each tree (using branching.times in ape).
From that you can create a matrix with the branching events (ordered in
time) as columns and the corresponding time in each tree as rows.
Using the quantile function
The function multiDiv in my package 'paleotree' plots the 95% quantile
for a diversity curve / LTT as a gray shaded region, among other
options.
-Dave Bapst
On Tue, Apr 16, 2013 at 12:33 AM, Simon Greenhill si...@simon.net.nz wrote:
Hi all,
Does anyone know of a good implementation of an LTT