[Rd] cutree (PR#1067)

maechler@stat.math.ethz.ch maechler@stat.math.ethz.ch
Fri, 24 Aug 2001 17:36:23 +0200 (MET DST)


>>>>> "Anja" == Anja von Heydebreck <heydebre@molgen.mpg.de> writes:

    Anja> Full_Name: Anja von Heydebreck Version: 1.3.0 OS: Alpha Unix
    Anja> Submission from: (NULL) (141.14.19.61)

Thank you

    Anja> I repeatedly obtained meaningless results from the function
    Anja> 'cutree' in the 'mva' package, when the argument 'h' was greater
    Anja> or equal to the maximum height occuring:

Here is an executable version of Anja's example:

 library('mva')
 y <-
 rbind(c( 0,  1, -1, 1),
       c( 0, -1,  1, 0),
       c( 0, -1,  1, 1),
       c(-1, -2,  1, 0))
 d <- dist(y, method="manhattan")
 hc <- hclust(d, method = "complete")
 hc$height
 ## [1] 1 3 7
 ( cluster <- cutree(hc, h=max(hc$height)))

  ## Anja got
  ## [1] 1108612352          1 1074517880          1

  ## Martin: I get
  ## [1] 0 0 0 0

even though I don't get random results, it's still wrong.
I'll look at a fix.

Martin Maechler <maechler@stat.math.ethz.ch>	http://stat.ethz.ch/~maechler/
Seminar fuer Statistik, ETH-Zentrum  LEO D10	Leonhardstr. 27
ETH (Federal Inst. Technology)	8092 Zurich	SWITZERLAND
phone: x-41-1-632-3408		fax: ...-1228			<><

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._