# HG changeset patch # User yhoogstrate # Date 1400767803 14400 # Node ID 00b9ceac43322728bdd5a6b775cc7b2f3b5b7c8e # Parent 185effd0163378f3352bf6889b0fb4a42dd3c299 Uploaded diff -r 185effd01633 -r 00b9ceac4332 edgeR_Differential_Gene_Expression.xml --- a/edgeR_Differential_Gene_Expression.xml Wed May 21 08:50:36 2014 -0400 +++ b/edgeR_Differential_Gene_Expression.xml Thu May 22 10:10:03 2014 -0400 @@ -277,10 +277,17 @@ ## change e into epsilon somehow dev.off() } + + if(output_heatmap_plot != "/dev/null") { + pdf(output_heatmap_plot) + etable2 <- topTags(lrt, n=100)\$table + order <- rownames(etable2) + heatmap(t(cpm(dge,normalized.lib.sizes=TRUE),sep="\t",row.names=TRUE,col.names=NA),log=TRUE)) + dev.off() + } } ##output_hierarchical_clustering_plot = args[13] - ##output_heatmap_plot = args[14] if(output_RData_obj != "/dev/null") { save.image(output_RData_obj)