Fwd: Re: Fwd: Re: [BioC] Heatmap function
Marcus
marcusb at biotech.kth.se
Thu Oct 30 09:31:21 MET 2003
>Hello again... I still have some questions. If you have the time Sean (or
>someone else) I would be very grateful.
>
>Yes, the M-values are log, but I have checked and I do not have any
>infinite or NaN values in the matrix.
>
>I have som questions regarding hclust. It says that the function
>calculates the distances row-wise in the helpfile. But what I do not
>understand is that I have chips with M-values that are seperated
>column-wise. Am I supposed to have them by row instead?
>
>Another question.
>The function looks like this hclust(as.dist(myobject), method ="complete")
>The as.dist have euclidean distance by default and the hclust have
>"complete" as default. What does that mean, that It calculates two??
>distances or that on distancemeasure is "shut off"?
>
>And the final question.
>If I use hclust to create an object to use in the heatmap function, is
>that the object that I give to the heatmapfunction both as colv and rowv?
>heatmap(x=mymatrixwithcorrelations, Rowv =
>distanceofmymatrixwithcorrelations,
>Colv=distanceofmymatrixwithcorrelations, distfun=dist, hclustfun = hclust
>... etc) ?
>
>Or is there a way to fist cluster by row in hclust and then by column?
>
>Thanks for the help
>
>/ A very confused person
>
>
>
>
>>>I assume that the M values are a log, so you might check to be sure that you
>>>do not have infinite or NaN results from taking logs of zero and the like.
>>>
>>>As for dist, use ?dist to get the help file for dist. There are other
>>>distances that you can use (by changing the arguments to dist). However,
>>>one trick that I have found useful is to do your clustering as you want it
>>>done before calling heatmap and then call heatmap with Rowv=rowclust and
>>>Colv=colclust where rowclust and colclust are "cluster" objects that you can
>>>obtain from calling hclust. In this way, you can make your clusters "look
>>>good" before spending the time having heatmap do it for you each time.
>>>
>>>Hope that helps. Oh, by the way, I am absolutely not an expert, so check
>>>with other folks, also.
>>>
>>>Sean
>>*******************************************************************************************
>>
*******************************************************************************************
Marcus Gry Björklund
Royal Institute of Technology
AlbaNova University Center
Stockholm Center for Physics, Astronomy and Biotechnology
Department of Molecular Biotechnology
106 91 Stockholm, Sweden
Phone (office): +46 8 553 783 39
Fax: + 46 8 553 784 81
Visiting adress: Roslagstullsbacken 21, Floor 3
Delivery adress: Roslagsvägen 30B
More information about the Bioconductor
mailing list