This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Limma voom with missing values from mass spectrometry

Hi everyone. Firstly thank you in advance for any help you can give, I am new to working with mass spec data and biostars has been immensely helpful.

I data from three technical repeats and three biological repeats (9 Controls, 9 Treated - 18 Total) Each time the mass spec was performed some genes have an N/A value. I have performed a differential protein analysis using limma which worked really well for each technical repeat however when combining all the biological repeats together the number of total of genes reduces dramatically as only genes that have no N/A values across all samples can be used.

Is there any way to work with what is there to expand the gene list out (i.e. still calculate p.value for those genes with only a few missing values)

Thanks again

differential protein missing t-test. limma
• 2,329 views • link
updated by Gordon Smyth Gordon Smyth • written by harrydolan.dc • harrydolan.dc

Why did you delete the question, harrydolan.dc?

• 0 views • link
written by Ram Ram

1 answer

You can use limma for mass spec data, but not voom. voom is designed specifically for count data and is neither needed for appropriate for mass spec data. If you use non-voom limma, then there is no need to remove NAs, so the problem you mention will go away. I would input log-expression values to limma then set trend=TRUE when you run eBayes.

• 0 views • link
written by Gordon Smyth Gordon Smyth

Log in to answer this question.

More posts like this