Plot a Traitgram for Continuous Traits
plotTraitgram plots a traitgram showing the evolution of a continuous trait.
If node values are not given (i.e. the data is empirical data collected from tips,
rather than simulated data), maximum-likelihood ancestral trait estimation is used
to calculate node values. (Ackerly, 2009) given a tree and a set of continuous trait
A vector of continuous trait values. If the length of
Main title of traitgram plot.
The line width used for branches in the figure.
By default, this function will use
ace from the library
reconstruct ancestral traits and confidence intervals using the PIC method, if internal
node values (i.e. ancestral node values) are not given.
As with many functions in the paleotree library, absolute time is always decreasing, i.e. the present day is zero.
Return no value, just plot the traitgram.
One should probably never do ancestral trait estimation without looking at the confidence intervals, as these reconstructed estimates tend to be very uncertain.
David W. Bapst
Ackerly, D. 2009 Conservatism and diversification of plant functional traits: Evolutionary rates versus phylogenetic signal. Proceedings of the National Academy of Sciences 106(Supplement 2):19699–19706.
Also see the functions
traitgram in the library picante and
phenogram in the library phytools.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
set.seed(444) tree <- rtree(10) trait <- rTraitCont(tree) #first, traitgram without conf intervals plotTraitgram(trait,tree,conf.int=FALSE) #now, with plotTraitgram(trait,tree) #not much confidence, eh? # plotting simulated data # with values for ancestral nodes as input trait <- rTraitCont(tree, ancestor=TRUE) plotTraitgram(tree=tree,trait=trait)