Biostar Beta. Not for public use.
Mann-Whitney test -- is it applicable to compare log2 fold change in a subset of genes to all genes?
0
Entering edit mode
20 months ago
biostart • 290
Germany

Subj. Mann-Whitney test -- is it applicable to compare log2 fold change in a subset of genes to all genes?

(See the history of this question here)

RNA-Seq • 1.2k views
ADD COMMENTlink
0
Entering edit mode

How big are the subsets and full sets? You likely want GSEA instead...

ADD REPLYlink
0
Entering edit mode

The full set is 7,000 genes, and the subsets are like 1,000, 2,000 genes

ADD REPLYlink
1
Entering edit mode

N.B., you'll need abs(log2 fold change) for Mann-Whitney to make sense.

With a subset of that size the question becomes how you want to think about the comparison. Given your background question, I suspect that GSEA is answering the question you want to ask (i.e., "given two groups, is the subset 'perturbed' in one group versus the other in respect to the other genes").

ADD REPLYlink
0
Entering edit mode

Hmm, if I take abs(log2 fold change), then how do I distinguish between up- and downregulated? This will be a problem... But what if I just use the two sample t-test, would that be problematic?

ADD REPLYlink
0
Entering edit mode

You're comparing a subset versus the whole. The direction of change of individual genes is completely irrelevant to that question.

ADD REPLYlink
0
Entering edit mode

I actually don't understand the argument about "abs". Why negative values are a problem for Mann-Whitney?

ADD REPLYlink
0
Entering edit mode

I just need one quick calculation, would not want install GSEA for that. Do you know, which statistical test would be used in GSEA in this case?

ADD REPLYlink
1
Entering edit mode

It's largely a more complicated version of ks.test() in R.

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1