[R] identify.hclust() doesn't cut tree at the vertical position of the mouse pointer
Uwe Ligges
ligges at statistik.tu-dortmund.de
Wed Jul 11 12:24:05 CEST 2012
On 10.07.2012 10:48, WATSON Mick wrote:
> Dear All
>
> According to the identify.hclust documentation the function "cuts the tree at the vertical position of the pointer and highlights the cluster containing the horizontal position of the pointer".
>
> When I carry out this, the tree isn't cut where I click - in fact, there seems to be a limit below which I cannot go.
>
> Consider the following code:
>
> mat <- matrix(rnorm(5000), ncol=5)
> hc <- hclust(dist(mat))
> plot(hc)
> identify(hc)
>
> No matter where I click on the tree, I cannot cut below around about 5. I can cut above that value, but not below.
So time to read the help page:
?identify.hclust
tells you about an argument MAXCLUSTER
and hence
identify(hc, MAXCLUSTER=500)
will probably due the trick.
Best,
Uwe Ligges
> Any help is much appreciated.
>
>> sessionInfo()
> R version 2.15.1 (2012-06-22)
> Platform: x86_64-pc-mingw32/x64 (64-bit)
>
> locale:
> [1] LC_COLLATE=English_United Kingdom.1252 LC_CTYPE=English_United Kingdom.1252 LC_MONETARY=English_United Kingdom.1252 LC_NUMERIC=C
> [5] LC_TIME=English_United Kingdom.1252
>
> attached base packages:
> [1] grid stats graphics grDevices utils datasets methods base
>
> other attached packages:
> [1] geneplotter_1.34.0 lattice_0.20-6 annotate_1.34.1 AnnotationDbi_1.18.1 Biobase_2.16.0 BiocGenerics_0.2.0 BiocInstaller_1.4.7 gplots_2.11.0 MASS_7.3-18
> [10] KernSmooth_2.23-7 caTools_1.13 bitops_1.0-4.1 gdata_2.11.0 gtools_2.7.0
>
> loaded via a namespace (and not attached):
> [1] DBI_0.2-5 IRanges_1.14.4 RColorBrewer_1.0-5 RSQLite_0.11.1 stats4_2.15.1 tools_2.15.1 XML_3.9-4.1 xtable_1.7-0
>
> Thanks
> Mick
>
More information about the R-help
mailing list