I have integrated multiple sc datasets using scvi. Now I’m working on annotating the data. The first step of manual annotation is to conduct differential expression analysis between cluster. It seems that SCVI.differential_expression by default conducts differential expression analysis on the original expression matrix (X), though the clusters are assigned according to the latent space embedding. I wonder if this makes sense since the batch effect that lies in X hasn’t been removed. Should I instead force SCVI.differential_expression to conduct differential expression analysis on the corrected expression matrix (Expected
frequency fw(zn, sn) in the paper)?