/
ltt.Rd
64 lines (59 loc) · 3.27 KB
/
ltt.Rd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
\name{ltt}
\alias{ltt}
\alias{gtt}
\alias{mccr}
\title{Creates lineage-through-time plot (including extinct lineages)}
\usage{
ltt(tree, plot=TRUE, drop.extinct=FALSE, log.lineages=TRUE, gamma=TRUE, ...)
gtt(tree, n=100, ...)
mccr(obj, rho=1, nsim=100, ...)
}
\arguments{
\item{tree}{is a phylogenetic tree in \code{"phylo"} format, or an object of class \code{"multiPhylo"} containing a list of phylogenetic trees.}
\item{plot}{a logical value indicating whether or not to create LTT plot.}
\item{drop.extinct}{logical value indicating whether or not to drop extinct tips from the tree.}
\item{log.lineages}{logical value indicating whether LTT plot should be on log-linear (default) or linear-linear scale.}
\item{gamma}{logical value indicating whether or not to compute eqn{gamma} from Pybus & Harvey (2000; \emph{Proc. Roy. Soc. B}).}
\item{n}{for \code{gtt} the number of time intervals to use to track \eqn{\gamma} through time.}
\item{obj}{for \code{mccr} an object of class \code{"ltt"}.}
\item{rho}{for \code{mccr} sampling fraction.}
\item{nsim}{for \code{mccr} number of simulations to use for the MCCR test.}
\item{...}{other arguments to be passed to plotting methods. See \code{\link{plot.default}}.}
}
\description{
The function \code{ltt} computes LTT plot with extant and extinct lineages, and optionally conducts \eqn{\gamma}-test of Pybus & Harvey (2000). The object returned by \code{ltt} can be plotted or re-plotted using \code{\link{plot}}.
The function \code{gtt} computes the value of Pybus & Harvey's \eqn{\gamma} statistic through time by slice the tree at various points - by default in even intervals from the time above the root at which \emph{N} = 3 to the present day.
The function \code{mccr} performs the MCCR test of Pybus & Harvey (2000) which takes into account incomplete taxon sampling in computing a P-value of the \eqn{\gamma} statistic.
}
\details{
Although it is calculated here, it's unclear how to interpret the \eqn{\gamma}-statistic if not all the tips in the tree are contemporaneous.
}
\value{
\code{ltt} returns an object of class \code{"ltt"} which includes the following components:
\item{times}{a vector of branching times.}
\item{ltt}{a vector of linages.}
\item{gamma}{optionally, a value for the gamma-statistic.}
\item{p}{two-tailed P-value for the gamma-test.}
If \code{tree} is an object of class \code{"multiPhylo"}, then an object of class \code{"multiLtt"} is returned consisting of a list of object of class \code{"ltt"}.
\code{gtt} returns an object of class \code{"gtt"}.
\code{mccr} returns of object of class \code{"mccr"}.
}
\references{
Pybus, O. G., and P. H. Harvey (2000) Testing macro-evolutionary models using incomplete molecular phylogenies. \emph{Proc. R. Soc. Lond. B}, \bold{267}, 2267-2272.
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). \emph{Methods Ecol. Evol.}, \bold{3}, 217-223.
}
\author{Liam Revell \email{liam.revell@umb.edu}}
\seealso{
\code{\link{gammatest}}, \code{\link{ltt95}}
}
\examples{
trees<-pbtree(n=100,scale=100,nsim=10)
obj<-ltt(trees,plot=FALSE)
plot(obj,log="y",log.lineages=FALSE,main="lineage through time plots")
tree<-pbtree(b=1,d=0.25,t=4)
obj<-ltt(tree,gamma=FALSE)
obj
}
\keyword{phylogenetics}
\keyword{plotting}
\keyword{diversification analysis}