Question: How to get fitted per sample means and dispersion from DESeq2?
0
gravatar for kjkjindal
3 months ago by
kjkjindal0
kjkjindal0 wrote:

Hi,

DESeq uses the following model for differential expression analysis:

K_ij ~ NB(mu_ij, alpha_i)

mu_ij = s_j q_ij

log2(q_ij) = x_j. beta_i

Is there a way to access the specific 'mu_ij' and 'alpha_i' values?

I have been looking around and found assays(dds)[['mu']] giving something like mu_ij but am not sure if those are indeed the values that I am looking for.

Thanks!

sequencing rna-seq next-gen R • 132 views
ADD COMMENTlink modified 3 months ago by Kevin Blighe26k • written 3 months ago by kjkjindal0
2
gravatar for Kevin Blighe
3 months ago by
Kevin Blighe26k
USA / Europe / Brazil
Kevin Blighe26k wrote:

Take a look at the Access to all calculated values of the tutorial Analyzing RNA-seq data with DESeq2.

Dispersion is calculated on a gene-wise basis and the final modelled dispersion can be obtained with mcols(dds)$dispersion or dispersions(dds). The estimated gene dispersions that go into the Bayesian model as priors can be obtained with mcols(dds)$dispGeneEst

For the mean, it states:

Screen_Shot_2018_05_16_at_16_00_51

ADD COMMENTlink written 3 months ago by Kevin Blighe26k

Hi, that worked great. Thanks a lot Kevin!

ADD REPLYlink written 3 months ago by kjkjindal0

You're welcome - no problem.

ADD REPLYlink written 3 months ago by Kevin Blighe26k
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 2.3.0
Traffic: 1498 users visited in the last hour