Latent factors for differential expression
1
1
Entering edit mode
@vincentcroset-14226
Last seen 23 months ago
United Kingdom

This is a copy of this post https://github.com/drisso/zinbwave/issues/58#issue-656953981 but I thought it might reach out to more people by posting it here too.

Hi everyone,

I am measuring differential expression between two conditions in 10X data, using zinbwave and DESeq2. I am using K=2 in zinbwave as there are some latent factors (e.g. library size) that I would like to infer from the data. Now, would it make sense to include W1 (as calculated with zinbwave) in the DESeq2 design model (something like design <- model.matrix(~condition + reducedDim(exp, "zinbwave")[,1])? Or is W somehow already reflected in the observational weights, in which case including it in the model would be redundant?

Similarly, as I use the LRT test in DESeq2 should I include W in the reduced formula?

Many thanks for your input! Vincent

deseq2 zinbwave 10X scRNA-seq Differential expression • 1.6k views
ADD COMMENT
2
Entering edit mode
@mikelove
Last seen 6 hours ago
United States

I'm not sure but I believe that the latent factors are not redundant with the observation weights. I'll await the zinbwave team to answer definitely.

Yes, if you included a latent factor in the full design which is controlling for unwanted variation, you should also include it in the reduced design.

ADD COMMENT
2
Entering edit mode

Mike is right. The observation weights only account for the zero inflation, so if you want to account for the latent factor in the DE model you have to include W in both the full and reduced model.

ADD REPLY
0
Entering edit mode

That makes sense. Thank you both!

ADD REPLY

Login before adding your answer.

Traffic: 763 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6