compute.brtime

0th

Percentile

Compute and Set Branching Times

This function computes the branch lengths of a tree giving its branching times (aka node ages or heights).

Keywords
manip
Usage
compute.brtime(phy, method = "coalescent", force.positive = NULL)
Arguments
phy

an object of class "phylo".

method

either "coalescent" (the default), or a numeric vector giving the branching times.

force.positive

a logical value (see details).

Details

By default, a set of random branching times is generated from a simple coalescent, and the option force.positive is set to TRUE so that no branch length is negative.

If a numeric vector is passed to method, it is taken as the branching times of the nodes with respect to their numbers (i.e., the first element of method is the branching time of the node numbered \(n + 1\) [= the root], the second element of the node numbered \(n + 2\), and so on), so force.positive is set to FALSE. This may result in negative branch lengths. To avoid this, one should use force.positive = TRUE in which case the branching times are eventually reordered.

Value

An object of class "phylo" with branch lengths and ultrametric.

See Also

compute.brlen, branching.times

Aliases
  • compute.brtime
Examples
# NOT RUN {
tr <- rtree(10)
layout(matrix(1:4, 2))
plot(compute.brtime(tr)); axisPhylo()
plot(compute.brtime(tr, force.positive = FALSE)); axisPhylo()
plot(compute.brtime(tr, 1:9)); axisPhylo() # a bit nonsense
plot(compute.brtime(tr, 1:9, TRUE)); axisPhylo()
layout(1)
# }
Documentation reproduced from package ape, version 5.0, License: GPL (>= 2)

Community examples

Looks like there are no examples yet.