Can I Use ComBat-Seq to adjust multigroups RNAseq?

zhangyuqing / ComBat-seq

Batch effect adjustment based on negative binomial regression for RNA sequencing count data

154 stars 39 forks source link

Can I Use ComBat-Seq to adjust multigroups RNAseq? #28

Open yirenheihei opened 2 years ago

yirenheihei commented 2 years ago

hi,I am confused to use the combat_seq to deal with my RNAseq data .I have two batchs and I have thress groups such as groupA vs groupB,groupA vs groupC.How I use the combat_seq? below is my code . group=c(rep(0,1,2),6) batch=c(rep(0,1),9) adjusted_counts=combat_seq(counts,batch=batch,group=group) is this right? thank you very much.