Mercurial > repos > elixir-it > vinyl_score_heatmap

annotate 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
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
a7aa56155a91 Uploaded
elixir-it
parents:
diff changeset
1 #!/usr/bin/env Rscript
a7aa56155a91 Uploaded
elixir-it
parents:
diff changeset
2 args <- commandArgs(trailingOnly = TRUE)
a7aa56155a91 Uploaded
elixir-it
parents:
diff changeset
3 file1=args[1]
a7aa56155a91 Uploaded
elixir-it
parents:
diff changeset
4 cutoff=as.numeric(args[2])
a7aa56155a91 Uploaded
elixir-it
parents:
diff changeset
5 ofile=args[3]
a7aa56155a91 Uploaded
elixir-it
parents:
diff changeset
6 data=read.table(file1,header=T,sep="\t")
a7aa56155a91 Uploaded
elixir-it
parents:
diff changeset
7 rownames(data)=paste(data[,1], data[,2], data[,3], data[,4] ,sep="_")
a7aa56155a91 Uploaded
elixir-it
parents:
diff changeset
8 png(ofile,width=1600,height=1800)
a7aa56155a91 Uploaded
elixir-it
parents:
diff changeset
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)
a7aa56155a91 Uploaded
elixir-it
parents:
diff changeset
10 dev.off()
a7aa56155a91 Uploaded
elixir-it
parents:
diff changeset
11
a7aa56155a91 Uploaded
elixir-it
parents:
diff changeset
12