Skip to content
Snippets Groups Projects
Commit 35a91ea1 authored by Gervaise Henry's avatar Gervaise Henry :cowboy:
Browse files

Generalize count column earlier

parent b36bb08e
Branches
Tags
2 merge requests!37v0.0.1,!22Resolve "process_count"
Pipeline #6239 passed with stages
in 53 minutes and 49 seconds
......@@ -17,9 +17,10 @@ if (!("count" %in% names(opt))){
repRID <- basename(gsub(".featureCounts","",opt$count))
count <- read.delim(opt$count, comment.char="#")
count <- read.delim(opt$count, comment.char="#") # if featureCounts file changes structure, be sure to update count and Length columns below
colnames(count)[7] <- "count"
rpk <- count$Q.Y5JA.sorted.deduped.bam/count$Length/1000
rpk <- count$count/count$Length/1000
scale <- sum(rpk)/1000000
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment