Entering edit mode
Jason Skelton
▴
510
@jason-skelton-135
Last seen 10.2 years ago
Hi
I'm using the hclust function from mva R v 1.8
I've plotted a dendrogram with genenames
e.g.
plot(hclustoutput, labels=genenames etc)
I'm currently using the identify.hclust command:
ClusterID <- identify(hclustoutplut, MAXCLUSTER=60)
which returns
something like:
375
376
377
378
379
380 2.7004349
381 1.8060133 1.1440009
382 0.7630247 3.1794324 2.1269562
383 0.4950094 2.4738040 1.4827588 0.7810937
384 0.9291393 2.2186614 1.1459034 1.0344517 0.6617410
385 1.4570777 2.5765989 1.7992481 1.6056936 1.5644727 1.1153774
386 0.3695838 2.8249040 1.8551363 0.5220643 0.5924702 0.8052963
1.2503008
for each cluster specified.
My question is can I specify the genenames that I used in the original
plot and return them for each of the clusters I choose rather than the
information that is returned by default.
many thanks
Jason
--
--------------------------------
Jason Skelton
Pathogen Microarrays
Wellcome Trust Sanger Institute
Hinxton
Cambridge
CB10 1SA
Tel +44(0)1223 834244 Ext 7123
Fax +44(0)1223 494919