[BioC] clustering in R

priya [guest] guest at bioconductor.org
Tue Oct 23 09:34:14 CEST 2012


I have a RMA normalized genes expression datset with 22810 rows and 9 columns( types of promoters) and a subset of the data is as follows:

    ID_REF GSM362180    GSM362181  GSM362188    GSM362189  GSM362192
    244901 5.094871713 4.626623079 4.554272515 4.748604391 4.759221647
    244902 5.194528083 4.985930299 4.817426064 5.151654407 4.838741605
    244903 5.412329253 5.352970877 5.06250609  5.305709079 8.365082403
    244904 5.529220594 5.28134657  5.467445095 5.62968933  5.458388909
    244905 5.024052699 4.714631878 4.792865831 4.843975286 4.657188246
    244906 5.786557533 5.242403911 5.060605782 5.458148567 5.890061836





 -- output of sessionInfo(): 

I want to do a clustering of the above and tried the hierarchical clustering:

    d <- dist(as.matrix(deg), method = "euclidean")
where deg is the a matrix of the differentially expressed genes ( 4300 in number ).And I get the following warning:

      Warning message:
     In dist(as.matrix(deg), method = "euclidean") : NAs introduced by coercion

 Is it allright to proceed with the clustering inspite of the warning ?


    hc <- hclust(d)
    plot(hc, hang = -0.01, cex = 0.7)

I get a dendrogram which is very dense and the labels are not clear: Also I do not know which of the 9 promoters are classified in the tree for the several genes: How would it be possible to label the tree with the promoters and also how to visualize the genes into a clearer dendrogram? There are around 4300 genes and would like to get a better dendrogram so that I could visualize it better.




--
Sent via the guest posting facility at bioconductor.org.



More information about the Bioconductor mailing list