Question: Mann-Whitney test -- is it applicable to compare log2 fold change in a subset of genes to all genes?
0
gravatar for biostart
3.9 years ago by
biostart350
Germany
biostart350 wrote:

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.5k views
ADD COMMENTlink written 3.9 years ago by biostart350

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

ADD REPLYlink written 3.9 years ago by Devon Ryan96k

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

ADD REPLYlink written 3.9 years ago by biostart350
1

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 written 3.9 years ago by Devon Ryan96k

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 written 3.9 years ago by biostart350

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

ADD REPLYlink written 3.9 years ago by Devon Ryan96k

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

ADD REPLYlink written 3.9 years ago by biostart350

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 written 3.9 years ago by biostart350
1

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

ADD REPLYlink written 3.9 years ago by Devon Ryan96k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1687 users visited in the last hour