Mercurial > repos > xuebing > sharplab_interval_analysis
annotate alignvis.r @ 24:8dd2a3f51c42 draft default tip
Deleted selected files
author | xuebing |
---|---|
date | Wed, 09 May 2012 10:53:44 -0400 |
parents | 16ba480adf96 |
children |
rev | line source |
---|---|
20 | 1 args <- commandArgs(TRUE) |
2 x <- read.table(args[1]) | |
3 pdf(args[2]) | |
4 #visualize the profile with heatmap | |
5 srt <- sort(x[,2],index.return=T) # sort by total number of reads | |
6 image(-t(log(as.matrix(x[srt$ix[1:nrow(x)],3:ncol(x)],nc=ncol(x)-2))),col=gray.colors(100)) | |
7 dev.off() | |
8 |