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

Reall remove unnecessary cols from tpm table

parent 9924e777
Branches
Tags
2 merge requests!41Develop,!40Remove.tpm.col
Pipeline #7944 passed with stages
in 1 hour, 33 minutes, and 51 seconds
......@@ -10,9 +10,6 @@ rm(option_list)
countTable <- read.csv(paste0(opt$repRID,".countData.countTable.csv"), stringsAsFactors=FALSE)
geneID <- read.delim("geneID.tsv", header=FALSE, stringsAsFactors=FALSE)
Entrez <- read.delim("Entrez.tsv", header=FALSE, stringsAsFactors=FALSE)
convert <- merge(x=Entrez,y=geneID,by.x="V1",by.y="V1",all.x=TRUE)
convert <- data.frame(geneID=countTable$Geneid)
convert <- merge(x=convert,y=geneID[,1:2],by.x="geneID",by.y="V2",all.x=TRUE)
......@@ -22,6 +19,6 @@ convert <- convert[,-1]
colnames(convert) <- c("GeneID","EntrezID")
convert <- unique(convert)
output <- merge(x=convert,y=countTable,by.x="GeneID",by.y="Geneid",all.x=TRUE)
output <- merge(x=convert,y=countTable[,c("Geneid","count","tpm")],by.x="GeneID",by.y="Geneid",all.x=TRUE)
write.table(output,file=paste0(opt$repRID,".tpmTable.csv"),sep=",",row.names=FALSE,quote=FALSE)
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