Mercurial > repos > xuebing > sharplab_seq_motif
diff mytools/alignvis.r @ 0:39217fa39ff2
Uploaded
author | xuebing |
---|---|
date | Tue, 13 Mar 2012 23:34:52 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/mytools/alignvis.r Tue Mar 13 23:34:52 2012 -0400 @@ -0,0 +1,8 @@ +args <- commandArgs(TRUE) +x <- read.table(args[1]) +pdf(args[2]) +#visualize the profile with heatmap +srt <- sort(x[,2],index.return=T) # sort by total number of reads +image(-t(log(as.matrix(x[srt$ix[1:nrow(x)],3:ncol(x)],nc=ncol(x)-2))),col=gray.colors(100)) +dev.off() +