Differences

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

--- mkatari-bioinformatics-august-2013-clustering [2014/12/15 12:05] – mkatari
+++ mkatari-bioinformatics-august-2013-clustering [2014/12/15 12:45] – mkatari
@@ Line 125: / Line 125: @@
+</code>
+The code below plots k-means clustering results. You simply have to provide the k-means output and the labels.
+<code>
+plotClusterCenters<-function(kmeansres,
+                             myxlab="Treatment",
+                             myylab="Expression",
+                             mymain="K-means Clusters") {
+  mycolors=c("blue","red","green","orange","pink","black")
+  centersdim = dim(kmeansres$centers)
+  plot(kmeansres$centers[1,],
+       type="b",
+       col=mycolors[1],
+       xlab=myxlab,
+       ylab=myylab,
+       main=mymain,
+       ylim=c(round(min(kmeansres$centers)),
+                               round(max(kmeansres$centers))),
+                               xaxt="n")
+  axis(1, at=c(1:centersdim[2]), labels=names(kmeansres$centers[1,]))
+  for (i in 2:centersdim[1]) {
+    lines(kmeansres$centers[i,], type="b", col=mycolors[i])
+  }
+}
+plotClusterCenters(SigGenes.kmeans.2)
 </code>