fix counting matrix header name

likelet · Mar 28, 2018 · 98d0a60 · 98d0a60
1 parent 4290288
commit 98d0a60
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/bin/get_kallisto_matrix.R b/bin/get_kallisto_matrix.R
@@ -44,10 +44,12 @@ map.file2=unique(map.file2)
 row.names(map.file2)=map.file2[,1]
 typelist=map.file2[levels(as.factor(genelist)),2]
 trans.kallisto.count.merge <- cbind(trans.kallisto.coun.out[,1],typelist,trans.kallisto.coun.out[,-1])
+colnames(trans.kallisto.count.merge)[1:2]=c("ID","type")
 fwrite(trans.kallisto.count.merge ,file="kallisto.count.txt",sep="\t",quote=F,row.names=F)
 
 #write tpm  matrix of  protein coding and lncRNA
-trans.kallisto.TPM.merge <- cbind(trans.kallisto.coun.out[,1],typelist,trans.kallisto.TPM.out[,-1])
+trans.kallisto.TPM.merge <- cbind(trans.kallisto.TPM.out[,1],typelist,trans.kallisto.TPM.out[,-1])
+colnames(trans.kallisto.TPM.merge)[1:2]=c("ID","type")
 #annoted gene type
 fwrite(trans.kallisto.TPM.merge ,file="kallisto.tpm.txt",sep="\t",quote=F,row.names=F)