[Rd] cluster - clusplot.default (PR#1249)
Mon, 7 Jan 2002 19:42:44 +0100 (MET)
The following code in clusplot.default (package cluster) is in error:
x1 <- cmdscale(x, k = 2, eig = TRUE)
var.dec <- sum(x1$eig)/sum(diag(x1$x))
if (var.dec < 0)
var.dec <- 0
if (var.dec > 1)
var.dec <- 1
x1 <- x1$points
x1 has components with names "points" and "eig", not "x", so
sum(diag(x1$x)) returns 0, the division gives Inf which is later replaced
So in the plot it is reported (always) that "These two components explain
of the variability".
Besides, is it reasonable that sum(NULL) returns 0 without at least a
Another small point about the cluster package: it loads automatically
mva, but that is not mentioned in the Depends field in the description file.
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: firstname.lastname@example.org