plot.phylo.extra | R Documentation |
These are extra functions to plot and annotate phylogenies, mostly calling basic graphical functions in ape.
plotBreakLongEdges(phy, n = 1, ...)
drawSupportOnEdges(value, ...)
phy |
an object of class |
n |
the numner of long branches to be broken. |
value |
the values to be printed on the internal branches of the tree. |
... |
further arguments to be passed to |
drawSupportOnEdges
assumes the tree is unrooted, so the vector
value
should have as many values than the number of internal
branches (= number of nodes - 1). If there is one additional value, it
is assumed that it relates to the root node and is dropped (see examples).
NULL
Emmanuel Paradis
plot.phylo
, edgelabels
,
boot.phylo
, plotTreeTime
tr <- rtree(10)
tr$edge.length[c(1, 18)] <- 100
op <- par(mfcol = 1:2)
plot(tr); axisPhylo()
plotBreakLongEdges(tr, 2); axisPhylo()
## from ?boot.phylo:
f <- function(x) nj(dist.dna(x))
data(woodmouse)
tw <- f(woodmouse) # NJ tree with K80 distance
set.seed(1)
## bootstrap with 100 replications:
(bp <- boot.phylo(tw, woodmouse, f, quiet = TRUE))
## the first value relates to the root node and is always 100
## it is ignored below:
plot(tw, "u")
drawSupportOnEdges(bp)
## more readable but the tree is really unrooted:
plot(tw)
drawSupportOnEdges(bp)
par(op)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.