[R] plot.hclust with lots of objects
Martin Maechler
maechler at stat.math.ethz.ch
Wed Jun 19 17:25:50 CEST 2002
>>>>> "Michael" == Michael Dondrup <Michael.Dondrup at Genetik.Uni-Bielefeld.DE> writes:
Michael> I found that converting my hclust object to a
Michael> dendrogram and cutting it with cut.dendrogram, does
Michael> almost what I want, showing only a branch of the
Michael> tree. But now my labels in the plot are gone,
Michael> although they appear in the dendrogram object
Michael> itself, when printed. How can I bring them back?
plot.dendrogram() is still not as good as it should be.
R-devel has an improved version -- which still isn't good
enough.
But for hierarchical clustering results (hclust being one kind;
there are agnes(), diana() and others in the "cluster" package),
you shouldn't need the detour through dendrograms at all.
cutree() is the function to cut hclust objects; returning
grouping vectors only
but the `maptree' package has a prune.clust function (andmore).
Hoping this helps, yours,
Martin Maechler <maechler at stat.math.ethz.ch> http://stat.ethz.ch/~maechler/
Seminar fuer Statistik, ETH-Zentrum LEO C16 Leonhardstr. 27
ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND
phone: x-41-1-632-3408 fax: ...-1228 <><
Michael> Laurent Gautier wrote:
>> On Sat, Jun 15, 2002 at 11:11:31PM +0200, Michael
>> Dondrup wrote: > Dear all,
>> >
>> > the default plotting method for hclust trees looks just
>> fine for few > objects like in the > example dataset. But
>> when it comes to many objects (some 1000 - I'm > trying
>> to visualize clustered microarray data) it renders a
>> tree, that > one cannot inspect, because of overlapping
>> text and lines. My question > is, is there a way or a
>> plotting parameter for plotting a tree which is > wide
>> enough to have all leaves separated and readable labels
>> even for > that many objects? > This would produce a
>> very big image, so I think scrolling is essential.
>> >
>> > I believe the answer is simple, but I'm unable to
>> figure it out.
>> >
>> > Thanks in advance > Michael
>>
>> Dear Michael,
>>
>> For the labels, you may want to try use smaller
>> characters to dsiplay them. This can be done using the
>> incantation: par(cex=0.5) # for example
>>
>> For the overlapping lines, well I think there is not much
>> to do except magnifying. Did you try to save your plot in
>> a format like postscript or pdf then display it with a
>> wiever that let you zoom in ?
>>
>> Ultimately, you can try to prune your tree and inspect
>> specific branches (I send you to R help files in this
>> case)
>>
>> Hopin' it helps,
>>
>> Laurent
>>
>> --------------------------------------------------------------
>> Laurent Gautier CBS, Building 208, DTU PhD. Student
>> DK-2800 Lyngby,Denmark tel: +45 45 25 24 85
>> http://www.cbs.dtu.dk/laurent
>> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
>> r-help 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-help-request at stat.math.ethz.ch
>> _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Michael> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
Michael> r-help mailing list -- Read
Michael> http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Michael> Send "info", "help", or "[un]subscribe" (in the
Michael> "body", not the subject !) To:
Michael> r-help-request at stat.math.ethz.ch
Michael> _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help 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-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list