Emirical Bayesian shrinkage for dispertion/LFC estimation
1
1
Entering edit mode
sotogarcia.r ▴ 10
@sotogarciar-22881
Last seen 2.3 years ago

Dear all,

I am quite new in bioinformatics/data science and I need some help with this issue. While I understand the intuition behind shrinking coefficient estimates and why this is applied to RNA sequencing experiments (I've studied Lasso and regularisation before), I don't understand the Bayesian approach to shrink the dispertion and LFC parameters. A good tutorial or a brief explanation here will do.

Thank you!

rnaseq deseq2 • 285 views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 1 day ago
United States

For one, you can read the DESeq2 paper and the apeglm paper.

Secondly, there was a similar question on Biostars a while back, and I said, first, please read the papers, but also offered this concise explanation:

https://www.biostars.org/p/127756/#127941

ADD COMMENT
0
Entering edit mode

I have read the paper, it's self-explanatory! The only part I don't understand is how Bayesian shrinkage is performed. I have been following your bioinformatics repository on Github, it has been really useful. What I understand from Bayesian shrinkage is that we start from a distribution of MLE LFC estimates. We then fit a normal distribution of mean 0 to these estimates that we use as prior, right? What I don't understand is how the posterior distribution is calculated and how the normal distribution with mean 0 helps shrink the LFC estimate.

ADD REPLY
1
Entering edit mode

Those details are in the two methods papers (note we use a Cauchy now, not Normal... see paper 2).

The posterior calculation is all spelled out for Normal in the DESeq2 paper.

In apeglm we approximate the posterior using the Hessian at the MAP.

For motivation on why do hierarchical modeling, there is a nice paper from Hongkai Ji and Shirley Liu. Just google their names + hierarchical model.

ADD REPLY
0
Entering edit mode

Thanks, I will have a look. Do you recommend other resources besides your GitHub repository to get me trained up as quickly as possible?

Thanks!

Dr. Ramon Soto

ADD REPLY
0
Entering edit mode

Sorry I’m pretty busy and can’t reply more here. Ive listed plenty of references.

ADD REPLY

Login before adding your answer.

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