Question: DeSeq RNASeq differential expression between region samples
0
gravatar for ecg1g15
18 days ago by
ecg1g150
ecg1g150 wrote:

Hi,

I have a dataset of RNASeq samples from 7 different regions. I am trying to see which genes are differentially expressed over all the regions.

When using the differential expression analysis of DESeq2 seems like it compares the differential expression between two conditions, whereas if I have 7 regions (conditions) A, B, C, D, E, F, G, I would like to know if there are significant expression changes on the individual gene transcriptions overall.

dds <- DESeq(dds) 
res <- results(dds)

log2 fold change (MLE): location A vs G 
Wald test p-value: location A vs G

I am using a vst normalised correction. And I have been told I should use the adjusted p-value for this. Is there a way where I can obtain a list of all the genes that are significant differentially expression across all regions?

Thanks

deseq2 R deg • 60 views
ADD COMMENTlink modified 18 days ago by Michael Love25k • written 18 days ago by ecg1g150
Answer: DeSeq RNASeq differential expression between region samples
1
gravatar for Michael Love
18 days ago by
Michael Love25k
United States
Michael Love25k wrote:

It sounds like you should be using an LRT of ~region vs ~1. Take a look at the vignette section on the Likelihood Ratio Test, and also see ?DESeq and ?results.

ADD COMMENTlink written 18 days ago by Michael Love25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 232 users visited in the last hour