Mercurial > repos > xuebing > sharplab_interval_analysis

view alignvis.r @ 24:8dd2a3f51c42 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Deleted selected files
author xuebing
date Wed, 09 May 2012 10:53:44 -0400
parents 16ba480adf96
children
line wrap: on
line source

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()