why my volcano plot looks strange?
1
0
Entering edit mode
29 days ago
hellocita ▴ 30

Hi all, I am dealing with an LC-MS Metabolomics data. I perform log transform, quantile normalization and use linear model to regress out batch. And use wilcox test to figure out which metabolite is significantly changed. However, the pvalue and log2 fold change of metabolites are weird. There are lots of dots near the middle line. Please help and give me some advice! the volcano plot volcano plot after the y-axis is reverse and loged

wilcox.test RNAseq Metabolomics LCMS • 447 views
ADD COMMENT
0
Entering edit mode

is it log10p or -log10p on y-axis?

ADD REPLY
0
Entering edit mode
29 days ago
swbarnes2 9.9k

Reverse and log scale the y axis for starters.

ADD COMMENT
0
Entering edit mode

Thanks for your help! Now it looks better, and there are still many dots near the middle line, was this acceptable?

ADD REPLY
0
Entering edit mode

First thing to do would be to not use custom statistics but a proper framework such as limma. There are quite some peptides (or whatever you measure) with tiny effect sizes but still small pvalues. Is this corrected for multiple testing?

ADD REPLY
0
Entering edit mode

The smallest p value is 0.02, so I did not correct it

ADD REPLY
2
Entering edit mode

Then your analysis is basically meaningless, FDR is not optional, it always needs to be applied to avoid a large amount of false-positives. I strongly recommend to either work with a statistician or use dedicated statistical software such as limma which does all that under the hood.

ADD REPLY

Login before adding your answer.

Traffic: 2475 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6