Mercurial > repos > xuebing > sharplab_interval_analysis

annotate 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
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
20
16ba480adf96 Uploaded
xuebing
parents:
diff changeset
1 args <- commandArgs(TRUE)
16ba480adf96 Uploaded
xuebing
parents:
diff changeset
2 x <- read.table(args[1])
16ba480adf96 Uploaded
xuebing
parents:
diff changeset
3 pdf(args[2])
16ba480adf96 Uploaded
xuebing
parents:
diff changeset
4 #visualize the profile with heatmap
16ba480adf96 Uploaded
xuebing
parents:
diff changeset
5 srt <- sort(x[,2],index.return=T) # sort by total number of reads
16ba480adf96 Uploaded
xuebing
parents:
diff changeset
6 image(-t(log(as.matrix(x[srt$ix[1:nrow(x)],3:ncol(x)],nc=ncol(x)-2))),col=gray.colors(100))
16ba480adf96 Uploaded
xuebing
parents:
diff changeset
7 dev.off()
16ba480adf96 Uploaded
xuebing
parents:
diff changeset
8