Plots of a Divisive Hierarchical Clustering
Creates plots for visualizing a diana object.
## S3 method for class 'diana'
plot(x, ask = FALSE, which.plots = NULL, main = NULL,
sub = paste("Divisive Coefficient = ", round(x$dc, digits = 2)),
adj = 0, nmax.lab = 35, max.strlen = 5, xax.pretty = TRUE, ...)x |
an object of class |
ask |
logical; if true and |
which.plots |
integer vector or NULL (default), the latter
producing both plots. Otherwise, |
main, sub |
main and sub title for the plot, each with a convenient
default. See documentation for these arguments in
|
adj |
for label adjustment in |
nmax.lab |
integer indicating the number of labels which is considered too large for single-name labelling the banner plot. |
max.strlen |
positive integer giving the length to which strings are truncated in banner plot labeling. |
xax.pretty |
logical or integer indicating if
|
... |
graphical parameters (see |
When ask = TRUE, rather than producing each plot sequentially,
plot.diana displays a menu listing all the plots that can be produced.
If the menu is not desired but a pause between plots is still wanted
one must set par(ask= TRUE) before invoking the plot command.
The banner displays the hierarchy of clusters, and is equivalent to a tree.
See Rousseeuw (1986) or chapter 6 of Kaufman and Rousseeuw (1990).
The banner plots the diameter of each cluster being splitted.
The observations are listed in the order found by the diana
algorithm, and the numbers in the height vector are represented
as bars between the observations.
The leaves of the clustering tree are the original observations. A branch splits up at the diameter of the cluster being splitted.
An appropriate plot is produced on the current graphics device. This can
be one or both of the following choices:
Banner
Clustering tree
In the banner plot,
observation labels are only printed when the number of observations is
limited less than nmax.lab (35, by default), for readability.
Moreover, observation labels are truncated to maximally
max.strlen (5) characters.
see those in plot.agnes.
example(diana)# -> dv <- diana(....)
plot(dv, which = 1, nmax.lab = 100)
## wider labels :
op <- par(mar = par("mar") + c(0, 2, 0,0))
plot(dv, which = 1, nmax.lab = 100, max.strlen = 12)
par(op)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.