phytools: bind.tip – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

bind.tip

Attaches a new tip to a tree

Description

Functions adds a new tip to the tree. If the tree is ultrametric and no branch length is specified, then edge.length is scaled so that the tree remains ultrametric after the new tip is added.

Usage

bind.tip(tree, tip.label, edge.length=NULL, where=NULL, position=0,
   interactive=FALSE, ...)

Arguments

`tree`	receptor tree.
`tip.label`	a string containing the species name for the new tip.
`edge.length`	edge length for the new tip (a scalar).
`where`	node number to attach new tip. If `position>0` then then tip will be attached rootward of the specified node. Node numbers can also be tips, in which case the new tip will be added along the terminal edge. To find out the tip number for given species with name "species" type: `which(tree$tip.label=="species")`.
`position`	distance rootward of the node to add the new tip.
`interactive`	logical value indicating whether or not the species should be added interactively. (Defaults to `FALSE`.)
`...`	arguments to be passed to `plotTree` (for `interactive=TRUE`.)

Details

Wrapper function for 'ape' bind.tree. Note that interactive=TRUE works only for right-facing phylograms.

Value

An object of class "phylo".

Author(s)

Liam Revell liam.revell@umb.edu

References

Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.

Examples

set.seed(123)
## generate tree
tree<-pbtree(b=0.1, n=10)

## plot original tree
plot(tree)
axisPhylo()

## add an extant tip ("t_extant") sister to taxon 't5' 
## with divergence time of 4.5 Ma
node <- which(tree$tip.label=="t5")
tree <- bind.tip(tree, tip.label="t_extant", 
    where=node, position=4.5)
# plot to see the result
plot(tree)
axisPhylo()

## add an extinct tip ("t_extinct") sister to 't2' with 
## divergence time of 7.8 Ma and duration (edge length) of 
## 3.3 Ma
node <- which(tree$tip.label=="t2") 
tree <- bind.tip(tree, tip.label="t_extinct", where=node, 
    position=7.8, edge.length=3.3)

## plot to see the result
plot(tree)
axisPhylo()

phytools

Phylogenetic Tools for Comparative Biology (and Other Things)

v0.7-70

GPL (>= 2)

Authors

Liam J. Revell

Initial release

2020-9-19