Median normalisation log2 transformed intensities from proteomics data
2
1
Entering edit mode
@sandramurphy-22841
Last seen 3.6 years ago

Hello,

I have a question about some proteomics data that I am analysing in R. I have log2 transformed the intensity data (from mass spec) and would like to perform median normalisation on this dataset. I came across the function "normalizeMedianValues(x)" in Limma. In my case my matrix x contains the log transformed intensities where each column is a different sample (in this case patients) and each row is a different protein.

Is this function okay to use in this situation? It seems like it was originally designed for normalising between different arrays.

Any advice is much appreciated!

Thanks, Sandra

median normalisation intensity data Limma • 2.4k views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 1 hour ago
United States

You could, but you need to 'unlog' first, as normalizeBetweenArrays (which is the function you should be using, rather than a helper function like normalizeMedianValues). Assuming your data are log2 transformed, you could do

x <- 2^x
normx <- normalizeBetweenArrays(x, "scale")

And go on from there. Note that normx will now contain log2 transformed data.

ADD COMMENT
0
Entering edit mode

Great, thanks so much for your detailed response! That's super helpful.

ADD REPLY
0
Entering edit mode
edammer • 0
@edammer-19863
Last seen 5 months ago
Atlanta, GA, United States

You could try median polish of ratio, tunable by which samples you choose for central tendency. The TAMPOR package has been published by Johnson ECB, et al, Nat Med, (2020) and is available on GitHub, designed for proteomic data with global internal standard (GIS) samples, but it works even for biological replicates of a particular treatment or disease (or control) group present in each batch.

ADD COMMENT
1
Entering edit mode

Okay great, sounds interesting. I'll definitely check it out. Thank you!

ADD REPLY

Login before adding your answer.

Traffic: 842 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