Hello,

I have patch-seq data (n = 16 samples). I am interested in using ZINB-WaVE along with DESeq2 to extract DE genes. I am comparing the results if I run DESeq2 only and ZINB-WaVE plus DESeq2. I don't see any difference if Observation weights = FALSE. The number of DE significant genes (p-adj < 0.01 and L2FC > 1.5) are 123 genes, for both the cases.

However, if I do set Observation weights = TRUE, the number of DE genes post-ZINB-WaVE plus DESeq2 goes up to 1515 genes. Can someone please explain how DESeq2 is using the observation weights for analysis?

Thanks in advance.