Mercurial > repos > elixir-it > vinyl_score_heatmap

comparison tool_heat_score.R @ 0:a7aa56155a91 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author elixir-it
date Thu, 17 Sep 2020 09:40:37 +0000
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:a7aa56155a91
1 #!/usr/bin/env Rscript
2 args <- commandArgs(trailingOnly = TRUE)
3 file1=args[1]
4 cutoff=as.numeric(args[2])
5 ofile=args[3]
6 data=read.table(file1,header=T,sep="\t")
7 rownames(data)=paste(data[,1], data[,2], data[,3], data[,4] ,sep="_")
8 png(ofile,width=1600,height=1800)
9 heatmap(as.matrix(data[data$ScoreT>=cutoff,9:21]),scale="row",col=colorRampPalette(c("white","purple","magenta"))(100),mar=c(12,16),cexRow=1.2,cexCol=2.5)
10 dev.off()
11
12