Skip to content
Snippets Groups Projects

Resolve "process_count"

Merged Gervaise Henry requested to merge 6-processCounts into develop
All threads resolved!
+ 30
0
gc()
library(optparse)
option_list=list(
make_option("--count",action="store",type='character',help="Count File")
)
opt=parse_args(OptionParser(option_list=option_list))
rm(option_list)
if (!("count" %in% names(opt))){
stop("No count file passed, exiting.")
} else if (!file.exists(opt$count)) {
user avatar user avatar user avatar +2
stop("No count file passed, exiting.")
}
repRID <- basename(gsub(".featureCounts","",opt$count))
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$count/count$Length/1000
scale <- sum(rpk)/1000000
tpm <- rpk/scale
output <- cbind(count,tpm)
colnames(output)[7] <- "count"
write.table(output,file=paste0(repRID,".countTable.csv"),sep=",",row.names=FALSE,quote=FALSE)