DEseq2 Analysis over multiple tissues and conditions only gives 16 coefficients
1
0
Entering edit mode
@3cbfcf40
Last seen 5 months ago
Germany

Hi, I am trying to perform a RNAseq analysis using DEseq2. My Experimental setup includes 8 different tissues with 2 genotypes each. I need to compare conditions within tissues and also tissues within genotypes. I tried using a ~genotype + tissue + genotype:tissue design but resultsNames(dds) only shows 16 coefficients and not all possible combinations, thus I can't select the comparisons I need for the results. In various posts, I read the suggestion to instead create a surrogate variable gentoype_tissue for each sample and use a simplified design formula ~genotype_tissue. I can now select every comparison I need using the contrasts. However, resultsNames(dds) still only gives me 16 coefficients and not all possible combinations. I would like to perform l2fc shrinkage for my results using the apeglm method, for which I need to use the coefficients which are not included for every possible combination.

Thanks a lot for your help, best Niko

DESeq2 RNASeq • 1.1k views
ADD COMMENT
0
Entering edit mode

Check contrasts in the vignette.

ADD REPLY
0
Entering edit mode

Thanks for your reply. Of course, I read the vignette, multiple times actually. Maybe I just don't understand how GLM work correctly... If I use a single combined group variable for the design ~genotype_tissue I don't understand why the coefficients don't yield all the different combinations. I must be missing something here... Thanks so much for your help, Niko

ADD REPLY
0
Entering edit mode
@mikelove
Last seen 1 minute ago
United States

only shows 16 coefficients and not all possible combinations

I'd recommend to collaborate with a local statistician or someone familiar with linear model coefficients in R. Due to time constrains I can only respond to software related questions on the support site.

ADD COMMENT

Login before adding your answer.

Traffic: 691 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