Differences

This shows you the differences between two versions of the page.

--- mkatari-bioinformatics-august-2013-clustering [2013/10/09 13:30] – mkatari
+++ mkatari-bioinformatics-august-2013-clustering [2013/10/11 15:01] – mkatari
@@ Line 2: / Line 2: @@
-Clustering rna-seq data, continuation from [[mkatari-bioinformatics-august-2013-deseq|DESeq]]
+====== Clustering rna-seq data ======
+continuation from [[mkatari-bioinformatics-august-2013-deseq|DESeq]]
 Get the significant genes
@@ Line 35: / Line 36: @@
 <code>
 sigGenes.hclust.k2<-cutree(sigGenes.normalized.hclust, k=2)
+</code>
+Now to get all the genes that are in cluster 2 simply type.
+<code>
+hclust.k2.cluster2=names(which(sigGenes.hclust.k2==2))
+</code>
+Now we can create a new matrix/data frame with just these genes. This new matrix can be used to plot a heatmap to make it easier to see a expression profile of the cluster (see below).
+<code>
+hclust.k2.cluster2.normalized = sigGenes.normalized[hclust.k2.cluster2,]
 </code>
@@ Line 53: / Line 66: @@
 </code>
-Heatmap
+====== Heatmap ======
 <code>
@@ Line 71: / Line 85: @@
 </code>
-Create heatmap. We can save it to a pdf file
+Create heatmap. We can save it to a pdf file. Note that sigGenes.normalized is just a matrix. Here we can provide any matrix of values, for example hclust.k2.cluster2.normalized which is the expression values of genes in cluster 2 (see above)
 <code>
@@ Line 91: / Line 105: @@
 dev.off()
 </code>