Thread: DE Seq
View Single Post
Old 09-02-2012, 03:41 AM   #10
Senior Member
Location: Stockholm, Sweden

Join Date: Feb 2008
Posts: 319

Anybody know how to do this with a BASH script??
You should be able to just loop over all files and use the 'join' command.

This of course assumes that the genes are always in the same order (which they are from HTSeq) but you may get into trouble if you have, e.g., HTSeq output using different versions of GTF files (like different ENSEMBL releases), when you will lose some genes from one of the files. That's why I think it's still better to join in R, Python etc where it is easier to check that nothing gets lost.
kopi-o is offline   Reply With Quote