Mercurial > repos > yufei-luo > differential_expression_analysis_pipeline_for_rnaseq_data

diff DiffExpAnal/DESeqTools/clusterPlot.R @ 0:63799b789162 draft
Uploaded
author: yufei-luo
date: Tue, 22 Jan 2013 10:07:03 -0500
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/DiffExpAnal/DESeqTools/clusterPlot.R	Tue Jan 22 10:07:03 2013 -0500
@@ -0,0 +1,27 @@
+# clusterPlot
+# dendrogram of sample clustering
+
+# input : counts, outputName, type of data (raw or norm)
+# output : dendrogram (jpeg)
+
+# created Sept 13th, 2012
+# modified Oct 30th, 2012
+# Marie-Agnes Dillies
+
+
+clusterPlot <- function( cds, OUT_clusterPlot, type = "raw", out = TRUE ){
+
+  if (out) png( file=OUT_clusterPlot )
+
+  if (type == "norm"){
+    cdsblind <- estimateDispersions( cds, method="blind" )
+    vsd <- getVarianceStabilizedData( cdsblind )
+  }
+  else {
+    vsd <- counts(cds)
+  }
+  hc <- hclust( dist(t(vsd)), method="ward" )
+  plot( hc, xlab = "Euclidean distance, Ward criterion", main=paste("Cluster Dendrogram, ", type, " data", sep="") )
+
+  if (out) dev.off()
+}