zhangyuqing / ComBat-seq

Batch effect adjustment based on negative binomial regression for RNA sequencing count data
154 stars 39 forks source link

Combat-seq, Pseudoaligners and DESeq2 #36

Open pintoa1-mskcc opened 1 year ago

pintoa1-mskcc commented 1 year ago

Hi

I have been working with combat-seq for a little while, where I quantify with kallisto, import with tximport, load to deseq2 to get "raw" counts, then run combat-seq on that data.

I noticed that when i use combat_seq, I need to load my data into deseq2 as a matrix, which removes all the special normalization deseq2 will do with tximported data.

Would this be the proper way to analyze pseudo counts? I feel I am losing a lot of the transcript/size factor normalization that is expected with pseudoaligners when I feed the data into combat-seq this way