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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/rotate.R
\name{rotate}
\alias{rotate}
\alias{rotate.default}
\alias{rotate.dendrogram}
\alias{rotate.hclust}
\alias{rotate.phylo}
\alias{sort.dendrogram}
\alias{sort.hclust}
\alias{sort.dendlist}
\alias{rev.hclust}
\title{Rotate a tree object}
\usage{
rotate(x, ...)
\method{rotate}{default}(x, order, ...)
\method{rotate}{dendrogram}(x, order, ...)
\method{rotate}{hclust}(x, order, ...)
\method{rotate}{phylo}(x, ..., phy)
\method{sort}{dendrogram}(x, decreasing = FALSE, type = c("labels",
"nodes"), ...)
\method{sort}{hclust}(x, decreasing = FALSE, ...)
\method{sort}{dendlist}(x, ...)
\method{rev}{hclust}(x, ...)
}
\arguments{
\item{x}{a tree object (either a \code{dendrogram} or \code{hclust})}
\item{...}{parameters passed (for example, in case of sort)}
\item{order}{Either numeric or character vector.
Is numeric: it is a numeric vector with the order of the value to be
assigned to object's label. The numbers say are just like when you use \link{order}:
which of the items on the tree-plot should be "first" (e.g: most left),
second etc. (this is relevant only to \code{rotate})
Is character: it must be a vector with the content of labels(x), in the
order we'd like to have the new tree.}
\item{phy}{a placeholder in case the user uses "phy ="}
\item{decreasing}{logical. Should the sort be increasing or decreasing? Not available for partial sorting. (relevant only to \code{sort})}
\item{type}{a character indicating how to sort. If "labels" then by lexicographic
order of the labels. If "nodes", then by using \link{ladderize} (order so that
recursively, the leftmost branch will be the smallest)}
}
\value{
A rotated tree object
}
\description{
Rotates, rev and sort the branches of a tree object (dendrogram, hclust)
based on a vector - eithor of labels order (numbers) or the labels in their
new order (character).
}
\details{
The motivation for this function came from the function
\code{\link{order.dendrogram}} NOT being very intuitive.
What \code{rotate} aims to do is give a simple tree rotation function which
is based on the order which the user would like to see the tree rotated by
(just as \code{\link{order}} works for numeric vectors).
\code{\link{rev.dendrogram}} is part of base R, and returns the tree object
after rotating it so that the order of the labels is reversed.
Here we added an S3 method for hclust objects.
The \code{sort} methods sort the labels of the tree (using \code{order})
and then attempts to rotate the tree to fit that order.
The hclust method of "\code{rotate}" works by first changing the object into
dendrogram, performing the rotation, and then changing it back to hclust.
Special care is taken in preserving some of the properties of the hclust
object.
The {ape} package has its own \code{\link[ape]{rotate}}({ape}) function
(Which is sadly not S3, so cannot be easily connected with the
current implementation). Still, there is an S3 plug that makes sure people
loading first ape and then dendextend will still be able to
use \code{rotate} without a problem.
Notice that if you will first load {ape} and only then {dendextend},
using "rotate" will fail with the error: "Error in rotate(dend, ____) :
object "phy" is not of class "phylo"" - this is because rotate in ape
is not S3 and will fail to find the rotate.dendrogram function.
In such a case simply run \code{unloadNamespace(ape)}. Or, you can run:
\code{unloadNamespace("dendextend"); attachNamespace("dendextend")}
The solution for this is that if you have {ape} installed on your machine,
It will be loaded when you load {dendextend} (but after it).
This way, \code{rotate} will work fine for both dendrogram AND phylo
objects.
}
\examples{
hc <- hclust(dist(USArrests[c(1,6,13,20, 23),]), "ave")
dend <- as.dendrogram(hc)
# For dendrogram objects:
labels_colors(dend) <- rainbow(nleaves(dend))
# let's color the labels to make the followup of the rotation easier
par(mfrow = c(1,2))
plot(dend, main = "Original tree")
plot(rotate(dend, c(2:5,1)), main =
"Rotates the left most leaf \\n into the right side of the tree")
par(mfrow = c(1,2))
plot(dend, main = "Original tree")
plot(sort(dend), main = "Sorts the labels by alphabetical order \\n
and rotates the tree to give the best fit possible")
par(mfrow = c(1,2))
plot(dend, main = "Original tree")
plot(rev(dend), main = "Reverses the order of the tree labels")
# For hclust objects:
plot(hc)
plot(rotate(hc, c(2:5,1)), main = "Rotates the left most leaf \\n
into the right side of the tree")
par(mfrow = c(1,3))
dend \%>\% plot(main = "Original tree")
dend \%>\% sort \%>\% plot(main = "labels sort")
dend \%>\% sort(type = "nodes") \%>\% plot(main = "nodes (ladderize) sort")
}
\seealso{
\code{\link{order.dendrogram}}, \code{\link{order}},
\code{\link{rev.dendrogram}}, \code{\link[ape]{rotate}} ({ape}), \link{ladderize}
}
|
