File: plot_comparison.Rscript

package info (click to toggle)
trinityrnaseq 2.15.2%2Bdfsg-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 468,004 kB
  • sloc: perl: 49,905; cpp: 17,993; java: 12,489; python: 3,282; sh: 1,989; ansic: 985; makefile: 717; xml: 62
file content (17 lines) | stat: -rwxr-xr-x 360 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
#!/usr/bin/env Rscript

charlotte = read.table("tximport.counts.dat")
mine = read.table("salmon.gene.counts.matrix")

mine = mine[rownames(charlotte),] # ensure in same order

samples = colnames(mine)

pdf("compare.pdf")

for (sample in samples) {

	c_sample = paste0('c_', sample)
    plot(log2(mine[,sample]+1), log2(charlotte[,c_sample]+1), main=sample)
}