Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
MEME multiple testing (Read 2366 times)
matias
YaBB Newbies
*
Offline


Curious HyPhy user

Posts: 6
MEME multiple testing
May 15th, 2013 at 10:19am
 
Hello,

I am running MEME on 2 thousand genes, so I will be searching for selection in that number of genes. My question is, which is the best way to correct for multiple testing on MEME results.

Thanks in advance for your response
Back to top
 
 
IP Logged
 
konrad
Junior Member
**
Offline


I love YaBB 1G - SP1!

Posts: 53
Re: MEME multiple testing
Reply #1 - May 15th, 2013 at 4:16pm
 
Dear matias,

It depends on what you want to report. If you want a list for which you are confident (at a certain significance level) that there are no false positives in the list, you need to do a family-wise correction, e.g. Holm-Bonferroni. We (and others in similar contexts) have argued that this is typically not what one wants when screening multiple genes - instead, we suggest reporting the uncorrected list as is along with the q-values (see Multimedia File Viewing and Clickable Links are available for Registered Members only!!  You need to Login Login), as we did in the MEME paper (tables S4-S19). The interpretation is then that the list is expected to contain false positives, and each individual q-value can be interpreted (roughly; see Multimedia File Viewing and Clickable Links are available for Registered Members only!!  You need to Login Login and Multimedia File Viewing and Clickable Links are available for Registered Members only!!  You need to Login Login for a detailed description) as the probability that that particular result is a false positive. The sum of the q-values gives the expected number of false positives in the list.

So we advocate a philosophy under which we are willing to report lists that almost certainly contain false positive results, along with numbers from which we can estimate how many false positives there are likely to be, but not which ones they are. Of course, it is then important to make sure this is clear when reporting the list.

Hope this helps,
Konrad
Back to top
 
WWW WWW  
IP Logged