This is a test version of Biostars. For the public version, visit https://www.biostars.org.
why my volcano plot looks strange?

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
• 4,601 views • link
updated by cpad0112 cpad0112 • written by hellocita • hellocita

is it log10p or -log10p on y-axis?

• 0 views • link
written by cpad0112 cpad0112

1 answer

Reverse and log scale the y axis for starters.

• 0 views • link
written by swbarnes2 swbarnes2

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

• 0 views • link
written by hellocita • hellocita

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?

• 0 views • link
written by ATpoint ATpoint

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

• 0 views • link
written by hellocita • hellocita

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.

• 0 views • link
written by ATpoint ATpoint

Log in to answer this question.

More posts like this