Question: microarray data analysis

0

au.rinki.bio •

**10**wrote:I am working on microarray data analysis. In Bonferroni correction α is devided by total sample number of microarray experiment or gene number.

ADD COMMENT
• link
•
modified 23 months ago
by
WouterDeCoster ♦

**42k**• written 23 months ago by au.rinki.bio •**10**
I'd strongly urge you to use Benjamini-Hochberg, if it is available to you, rather than the unforgiving Bonferroni approach. Strictly, in Bonferroni, alpha should be divided by the number of hypotheses tested, which is typically (the number of contrasts X the number of features), so in your terms, if you're only testing a single contrast it would be the number of genes

4.8kthank you so much for giving suggestion. actually i am new in this field i have one more question i am working on microarray data, i have a GSE14007 for the analysis. i have normalyzed this data and tried to find out differentially expressed genes. GSE14007 has 12 sample for normal condition and 12 sample for treated condition. after the preprocessing of this data their arround 22000 genes. to find out DE gene i have apply t-statistic and filtered the gene which p-vale<0.05. from litrature review i have found p-adjusted value and fdr based selection are more appropriate. now first i want to go with benferoni correction,alpha should be divided by the number of hypotheses tested but i am confuse here what is tested hypothesis? in my case is total sample number 24? or total gene number 22000? i hope you understand much better now. kindly give me your suggestion.

10I understood you to begin with, which is why I stated that you had to divide alpha by the "number of hypotheses tested". You only have one contrast (treated vs normal), and because this defines a hypothesis for each of your genes, you should divide by the number of genes.

But please be aware that the Bonferroni correction is too stringent, and is rarely used in microarray analysis. You should use the BH correction that I mentioned.

4.8kThank you so much for insightful reply.

10What is your question?

46kI interpreted this:

42ki am working on microarray data, i have a GSE14007 for the analysis. i have normalyzed this data and tried to find out differentially expressed genes. GSE14007 has 12 sample for normal condition and 12 sample for treated condition. after the preprocessing of this data their arround 22000 genes. to find out DE gene i have apply t-statistic and filtered the gene which p-vale<0.05. from litrature review i have found p-adjusted value and fdr based selection are more appropriate. now first i want to go with benferoni correction,alpha should be divided by the number of hypotheses tested but i am confuse here what is tested hypothesis? in my case is total sample number 24? or total gene number 22000? i hope you understand much better now. kindly give me your suggestion.

10