mtag icon indicating copy to clipboard operation
mtag copied to clipboard

Change genome-wide P threshold to control FDR

Open seanlzheng opened this issue 1 year ago • 1 comments

Thank you for developing and maintaining this excellent tool for the genomics community!

My understanding is that the computed maxFDR refers to the FDR at P 5E-8 threshold. If I am correct, is it possible to adjust the genome-wide P threshold to account for this e.g. calculate how much to reduce the P threshold by to control the FDR to a given %?

Thanks in advance

seanlzheng avatar Dec 13 '23 11:12 seanlzheng

You are right. It's the maximum potential FDR given the Omega and Sigma matrices for a threshold of 5e-8.

That's a really interesting idea you have though, but not one that we incorporated into the MTAG software. To do this, you'd probably need to dig into the hardcoded thresholds of the MTAG software and make those changes.

On Wed, Dec 13, 2023 at 6:10 AM seanlzheng @.***> wrote:

Thank you for developing and maintaining this excellent tool for the genomics community!

My understanding is that the computed maxFDR refers to the FDR at P 5E-8 threshold. If I am correct, is it possible to adjust the genome-wide P threshold to account for this e.g. calculate how much to reduce the P threshold by to control the FDR to a given %?

Thanks in advance

— Reply to this email directly, view it on GitHub https://github.com/JonJala/mtag/issues/198, or unsubscribe https://github.com/notifications/unsubscribe-auth/AFBUB5OUJNRQQBKT2HY44CLYJGEK5AVCNFSM6AAAAABAS7CM46VHI2DSMVQWIX3LMV43ASLTON2WKOZSGAZTSNBYGA4TCMQ . You are receiving this because you are subscribed to this thread.Message ID: @.***>

paturley avatar Dec 13 '23 15:12 paturley