HTseq header
0
0
Entering edit mode
3.3 years ago
Nai ▴ 50

I have htseq count file of all samples (150) together with out header. How can I add in header in that count.mat.tsv file

count Htseq matrix • 1.2k views
ADD COMMENT
0
Entering edit mode

Please post headers you want to add and first few lines of the htseq counts

ADD REPLY
0
Entering edit mode

The Header of the column as file name which same used for count matrix. Example of

0259.bam
126.bam 
180  .bam
45.bam
136.bam

Current count matrix file is combined for 150 samples: top 5 are: (count matrix file..)

ENSG00000000003.15  0   0   0   0
ENSG00000000005.6   1   0   0   0
ENSG00000000419.14  8   14  14  14
ENSG00000000457.14  75  51  48  63
ENSG00000000460.17  231 140 100 172
ENSG00000000938.13  0   0   0   0
ENSG00000000971.16  0   0   2   0
ENSG00000001036.14  3   7   2   0
ENSG00000001084.13  2   1   3   0
ENSG00000001167.15  2   2   2   8
ADD REPLY
0
Entering edit mode

how do you link each column with sample number? do the columns correspond to samples, in same order as those in bam list? For the list you posted:

$ cat <(datamash transpose < bam_list.txt) htseq.txt  | sed -e '1s/^/gene\t/;s/\.bam//g' 

gene    0259    126     180     45  136
ENSG00000000003.15  0   0   0   0
ENSG00000000005.6   1   0   0   0
ENSG00000000419.14  8   14  14  14
ENSG00000000457.14  75  51  48  63
ENSG00000000460.17  231 140 100 172
ENSG00000000938.13  0   0   0   0
ENSG00000000971.16  0   0   2   0
ENSG00000001036.14  3   7   2   0
ENSG00000001084.13  2   1   3   0
ENSG00000001167.15  2   2   2   8
ADD REPLY

Login before adding your answer.

Traffic: 1964 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6