Gene wise dispersion estimate Process - From the prior to the posterior
1
0
Entering edit mode
stew • 0
@d3c5f0ec
Last seen 4 days ago
Germany

In the paper of Love (2014) it says about the initial GLM:

To get a gene-wise dis- persion estimate for a gene i, we start by fitting a negative binomial GLM without an LFC prior for the design matrix X to the gene’s count data. This GLM uses a rough method-of-moments estimate of dispersion, based on the within-group variances and means.

What does it mean to fit a negative binomial GLM 'WITHOUT AN LFC PRIOR FOR THE DESIGN MATRIX X'? If we are fitting a negative binomial GLM then we would have to estimate the mean and the dispersion parameters. So where are the LFC priors taking place and what are they? And how do we get the initial fitted values \hat{\mu}_ij^0?

I would be very glad to get an explanation for this or a paper/book to read this in.

GLM prior DESeq LFC • 118 views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 2 days ago
United States

Just think, standard GLM as you would get with glm.nb.

We ended up breaking out the prior on coefficients into a separate function, so you can totally ignore the LFC prior part of the 2014 DESeq2 paper.

ADD COMMENT

Login before adding your answer.

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