search

BUStools / bustools

Tools for working with BUS files
https://bustools.github.io/
BSD 2-Clause "Simplified" License
92 stars 23 forks source link

how to compare eq classes between multiple runs #42

Closed vivekbhr closed 4 years ago

vivekbhr commented 4 years ago

Hi authors

I wanted to compare the TCCs between two kallisto-bus runs, which shared the same set of barcodes. I want to do the gene level testing, and I want to avoid dealing with multiple tx-to-ec maps. Is there a way to subset the outputs such that I can use ONE tx-to-ec map for all batches?

Best, Vivek

vivekbhr commented 4 years ago

Update: I tried to merge the two matrices by Kallisto merge (after tweaking the barcodes to make them unique), but it doesn't seem to do that I want. The ECs that were not present in the second sample (judging from the Tx contained in the set) somehow got counts for the second sample in the merged output. So the problem remains unsolved

vivekbhr commented 4 years ago

I implemented a custom solution. closing for now.