similar to: Can hclust() generate numerical results?

Hello, I would like to cut a hclust tree into several groups at a specific similarity. I assume this can be achieved by specifying the "h" argument with the specified similarity, e.g.: clust<-hclust(dist,"average") cut<-cutree(clust,h=0.65) Now, I would like to draw rectangles around the branches of the dendrogram highlighting the corresponding clusters, as is done by

I've been trying to get plot.hclust to work with a hclust object I created and have not had much success. It seems that there is some "hidden" characteristic of a hclust object that I can't see. This is most easily seen in the following example, where plot.hclust works on one object, but when this object is "dumped" and then re-read, plot.hclust no longer works. Is

This worked for me with your example: source("dumpdata.R") storage.mode(x.hc$merge) <- "integer" plot(x.hc) (R-1.4.1 compiled from source on WinNT4.) Andy > -----Original Message----- > From: Hugh Chipman [mailto:hachipma at icarus.math.uwaterloo.ca] > Sent: Wednesday, February 20, 2002 5:32 PM > To: andy_liaw at merck.com > Cc: r-help at stat.math.ethz.ch

Having some problems with using hclust.identify on my data, I revert to trying it out on the example in the help manual. My result is still the same (as with my own data): data(iris) hci <- hclust(dist(iris[,1:4])) plot(hci) testx<-identify.hclust(hci) Error in rect.hclust(x, k = k, x = x$x, cluster = cluster[, k - 1], border = "red") : k must be between 2 and 0

Hello all, I'm now working with hclust objects and was hoping to perform some basic editing on them like: - Joining = the merging of two hclust objects (so they will share one root) - Splicing = So to cut/extract a branch out of an hclust object - that by itself will be an hclust object. I noticed I could extract one element of an hclust object by turning it into a dendrogram,

Hi When I use plot(hclust(dist..)...)...) etc to create a dendrogram of a hierarchial cluster analysis, I end up with a vertical tree. What do I need to do to get a horizontal tree? Also, my users are used to seeing trees who's leaves all "end" at the same place (eg. Like in minitab). Is this possible in R? Thanks Mick Michael Watson Head of Informatics Institute for Animal

I have tried to use rect.hclust() to draw a rectangle around a set of leaves, but am running into trouble. The rect.hclust() is drawing two rects instead of one, and of the wrong size: -------------------- scoreClusterObj <- hclust(scoreDistanceObj, method=clustMethod) order <- scoreClusterObj$order orderedLabels <- rep(0, length(order)) for (orderIndex in 1:length(order)) { # this

There seems to be a bug in hclust (package mva) for clustering methods 'median' and 'centroid'. I have written a clustering program in C and discovered that the results for 'median' differ from those of hclust in R. I used a third program, written by someone else in Pascal, and that program agrees with the output of my program. I found yet another clustering program that

Good afternoon, I ran heatmap and hclust on the same matrix x (strictly, I ran heatmap(x), and hclust(dist(t(x))), and realized that the two dendrograms were slightly different, in that the left-right arrangement of one pair of subclusters (columns) was reversed in the two functions (but all individual columns were grouped correctly). Looking through the code for heatmap as a most definite

Thanks to Andy and Thomas, Reading help(hclust) more carefully would have done it but sometimes you do not see the wood for the trees... So hc$merge does exactly what I want. I have never been aware of the command str to get the structure of an R-object. It seems pretty useful to me. Thanks, Arne > -----Original Message----- > From: Liaw, Andy [mailto:andy_liaw at merck.com] >

I am making heatmaps for a dataset (~ 300*600 matrix) with the following R script (I am not familiar with R and this is the first time I am using it). library("gplots") library("Cairo") mydata <- read.csv(file="data.csv", header=TRUE, sep=",") rownames(mydata)=mydata$Name mydata <- mydata[,2:297] mydatamatrix <- data.matrix(mydata) mydatascale

Hello, I'm having trouble figuring out how to see resulting groups (clusters) from my hclust() output. I have a very large matrix of 4371 plots and 29 species, so simply looking at the graph is impossible. There must be a way to 'print' the results to a table that shows which plots were in what group, correct? I've attached the matrix I'm working with (the whole thing

Hi Sorry if this is in the help :-S I've looked at example(dendrogram) and though it gives some indication of what I want, it doesn't do all. OK, so here is what I want to do: draw a tree, and then have an action, on user-click, to either draw a sub tree or a plot of the data. I also want users to be able to search for a particular label and have it highlighted on the tree, say in

Hi All Not sure if this a bioconductor question or general R mailing list so apologies if this has gone to the wrong one................. When plotting dendrograms created by hclust you can "identify" clusters by clicking on the graphics and returning a list of what is contained in each cluster. However I'd like to be able to "zoom in" on specific clusters and plot

Dear R users, I have noted a difference in the merge distances given by hclust using centroid method. For the following data: x<-c(1009.9,1012.5,1011.1,1011.8,1009.3,1010.6) and using Euclidean distance, hclust using centroid method gives the following results: > x.dist<-dist(x) > x.aah<-hclust(x.dist,method="centroid") > x.aah$merge [,1] [,2] [1,] -3 -6

Hello, The heatmap function conveniently has a "reorder.dendrogram" function so that clusters follow a certain logic. It seems that the hclust function doesn't have such feature. I can use the "reorder" function on the dendrogram obtained from hclust, but this does not modify the hclust object itself. I understand that the answer should be within the "heatmap"

Dear group, I am currently struggling with the following problem for a while: I want to create a list whose elements consists of lists themselves. More concise: The list elements are HCLUST objects. However, when I try to append the HCLUST objects to my list via: cluster_list <- append(cluster_list, HCLUSTobject) the HCLUST object is appended - but not as an object but as its components. So

Hi all, it appears that plot.hclust disregards the argument 'cex'. Up until 2.15.1 (and probably all 2.15.x), specifying cex in plot.hclust changed the size of the plotted labels on each leaf. In 3.0.0 the character size remains the same irrespective of the cex setting. Since the help text for plot.hclust does not mention the cex argument, it is possible that an "undocumentde

I have been attempting to do some work using hclust, and have run into a (possibly subtle) problem. The background is that I constructed a dissimilarity matrix ``d1'' (it involved something called the ``Jaccard similarity coefficient''; I won't go into the details unless requested). I then did d2 <- as.dist(d1) try <- hclust(d2,method=ward)

I am having difficulty getting the dynamic tree cut package to work. Given the data table "myddtable" LengthPlaceColorAge5HRed224ABlue205WGreen243GRed222GBlue236WGreen255ARed194H Blue23 I created a similarity matrix using DAISY and Gower metric and specified Place and Color columns as characters (since they are categorical variables) > dd.daisy<-daisy(myddtable, metric =