File: plot_me.R

package info (click to toggle)
trinityrnaseq 2.11.0%2Bdfsg-6
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 417,528 kB
  • sloc: perl: 48,420; cpp: 17,749; java: 12,695; python: 3,124; sh: 1,030; ansic: 983; makefile: 688; xml: 62
file content (9 lines) | stat: -rw-r--r-- 423 bytes parent folder | download | duplicates (5) 
1
2
3
4
5
6
7
8
9
 
data_all = read.table("kmer_histo.all.txt")
data_norm = read.table("kmer_histo.NormMaxKCov50.txt")

log_data_all = cbind(data_all[,1], log(data_all[,2]+1))
log_data_norm = cbind(data_norm[,1], log(data_norm[,2]+1))

plot(log_data_norm, col='green', xlim=c(0,200), xlab="kmer occurrence count", ylab="log(number of unique kmers)", main="Kmer composition and \nin silico read normalization")
points(log_data_all, col='red')