Re: [Rd] error in plot.dendrogram (PR#7300)

From: <>
Date: Fri 22 Oct 2004 - 00:52:41 EST


First. If you should not do it like you write here. You will get an error loading the rda file. It is a binary format.

Error in load(dFile) : input has been corrupted, with LF replaced by CR

Hence, you should specify mode="wb" for downloading binary formats.

Try this code.


dFile <- paste(getwd(),"/hres.rda",sep="") if(!file.exists(dFile))

    download.file(url ="", dest= dFile,mode="wb") load(dFile)


In one think you are right, the X axis. But there is an ERROR.

On my Machine (see previous mail) not the complete dendrogram is drawn. It just draws the first 100 leavs out of 380 when I run Second. Even if I increase the size of the window maximally the dendrogram is not plotted.
And I checked it just again on a linux box. And it does not work there either.

How the dendrogram should look like you can figure out looking at the following.

hclustObj <- paste(getwd(),"/hress.rda",sep="") if(!file.exists(dFile))

    download.file(url ="", dest= hclustObj ,mode="wb") load(dFile)



Third. I have expected this comment about if it is meaningfull or not.  First for what I need it is it meaningfull. Second it is a valid dendrogram generated by as.dendrogram from a vaild hclust object. I do not need a plot routine which teaches me what she thinks is meaningfull or not.


Martin Maechler wrote:

>>>>>>"Eryk" == Eryk Wolski <>
>>>>>> on Thu, 21 Oct 2004 13:41:29 +0200 (CEST) writes:
> Eryk> Hi,
> Eryk> hres <- hclust(smatr,method="single")
> Eryk> hresd<-as.dendrogram(hres)
> Eryk> as.dendrogram(hres)
> Eryk> `dendrogram' with 2 branches and 380 members total, at height 2514.513
> Eryk> plot(hresd,leaflab="none") #<-error here.
>definitely no error here.. maybe your graphic window is too
>small or otherwise unable to show all the leaf labels?
> Eryk> #the plotted dendrogram is incomplete. The x axis is not drawn.
>ha! and why should this be a bug????
>Have you RTFHP and looked at its example??
>There's never an x-axis in such a plot!
>[You really don't want an x-axis overlayed over all the labels]
> Eryk> #The interested reader can download the
> Eryk> save(hresd,file="hres.rda")
> Eryk> #from the following loacation
> Eryk>
>If you send a bug report (and please rather don't..),
>it should be reproducible, i.e., I've just wasted my time for
>dFile <- "/u/maechler/R/MM/Pkg-ex/stats/wolski-hres.rda"
> download.file(url ="", dest= dFile)
>If you look at this plot I hope you rather see that "single" has
>been an extremly unuseful clustering method for this data / dissimilarities,
>and you'd rather tried other methods than to which for an
>If you really want one (just to see that it doesn't make sense),
>you can always add
> axis(1, fg = "red")

Dipl. bio-chem. Witold Eryk Wolski
MPI-Moleculare Genetic
Ihnestrasse 63-73 14195 Berlin
tel: 0049-30-83875219                 __("<    _      \__/    'v'    ||    /   \
mail:    ^^     m m

______________________________________________ mailing list
Received on Fri Oct 22 05:34:34 2004

This archive was generated by hypermail 2.1.8 : Wed 03 Nov 2004 - 22:45:23 EST