This is a test version of Biostars. For the public version, visit https://www.biostars.org.
50000 significant p values

I have ADNI microarray dataset (https://ida.loni.usc.edu/pages/access/geneticData.jsp#207). The paper says it is normalized data. I have used R -limma package to do the expression analysis among healthy controls and patients. It has around 50000 probe set and sample size of 700 (equal amount of control and patients). After the analysis, all 50000 probe ids showed high significance level (<0.01). Is it possible?

r microarray pvalue

What p-value correction did you apply (if any)?

Is this normal?

No.

But you'll need to add more information to this question in order to get accurate troubleshooting, such as which technology, which experimental setup, which packages, which code,...

I have ADNI microarray dataset (https://ida.loni.usc.edu/pages/access/geneticData.jsp#207). The paper says it is normalized data. I have used R -limma package to do the expression analysis among healthy controls and patients. It has around 50000 probe set and sample size of 700 (equal amount of control and patients). After the analysis, all 50000 probe ids showed high significance level (<0.01). Is it possible

Is it possible

No.

all 50000 probe ids showed high significance level (<0.01)

Are all transcripts up- or downregulated or both directions?

The paper

Which paper? Be specific.

And please use ADD REPLY to answer to earlier comments, as such this thread remains logically structured and easy to follow.

1 answer

Judging from your comment you did not apply any p-value correction. As you are doing multiple tests, you have to change your acceptance level from 0.01 to 0.01/50000 in the simplest case (Bonferroni). That then is the true alpha-value you have to compare with.

Bonferonni is not appropriate (too stringent) for gene expression analysis. In addition I would expect Limma to perform multiple testing correction.

I agree, False Discovery Rate (FDR) in limma is appropriate.

Log in to answer this question.