Basic question about Limma
2
0
Entering edit mode
Guido Leoni ▴ 200
@guido-leoni-3328
Last seen 9.4 years ago
European Union
Dear List I've a very basic question about the limma output (and I'm newbie of MA analysis). I'm looking some analysis output .Usually in the M column (if I well understand ) I should expect to find the log2-fold change of my conditions.Instead I find the difference between the normalized expressions... Some one could suggest to me how to transform these values in log2 fold change? Thank you very much and sorry for my silly question Guido Leoni [[alternative HTML version deleted]]
limma limma • 1.3k views
ADD COMMENT
0
Entering edit mode
Neel Aluru ▴ 460
@neel-aluru-3760
Last seen 7.3 years ago
United States
Hi Guido, What do you mean by "difference between the normalized expressions...." Neel On Apr 27, 2010, at 1:05 PM, Guido Leoni wrote: > Dear List I've a very basic question about the limma output (and I'm newbie > of MA analysis). > I'm looking some analysis output .Usually in the M column (if I well > understand ) I should expect to find the log2-fold change of my > conditions.Instead I find the difference between the normalized > expressions... > Some one could suggest to me how to transform these values in log2 fold > change? > Thank you very much and sorry for my silly question > Guido Leoni > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor > Neel Aluru Postdoctoral Scholar Biology Department Woods Hole Oceanographic Institution Woods Hole, MA 02543 USA 508-289-3607
ADD COMMENT
0
Entering edit mode
Heidi Dvinge ★ 2.0k
@heidi-dvinge-2195
Last seen 9.6 years ago
Hello Guido, the M-values are indeed the log2 fold change, so log2(R/G) for the red and green channel. However, after normalisation the R and G intensities are usually reported on a log2 scale (from 0-16), and since log2(R/G)=log2(R)-log2(G) that's probably what you see. HTH \Heidi > Dear List I've a very basic question about the limma output (and I'm > newbie > of MA analysis). > I'm looking some analysis output .Usually in the M column (if I well > understand ) I should expect to find the log2-fold change of my > conditions.Instead I find the difference between the normalized > expressions... > Some one could suggest to me how to transform these values in log2 fold > change? > Thank you very much and sorry for my silly question > Guido Leoni > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor >
ADD COMMENT
0
Entering edit mode
Thank you Heidi and Neel I'll give you more details about my data: The experiments are performed with an affymetrix custom 3' array here is an example for one probeset : Normalized expression with gcRMA control control control case case case 219070_s_at 679,15 734,98 777,32 667,32 507,76 645,37 This is the corrisponden row from limma output for the same probeset AFFy ID EG Symbol M A t P.value B 219070_s_at 64598 MOSPD3 -124.008563236060 666.759701175327 -2.27688122341071 0.0489596557364987 -4.19921794672031 2010/4/27 Heidi Dvinge <heidi@ebi.ac.uk> > Hello Guido, > > the M-values are indeed the log2 fold change, so log2(R/G) for the red and > green channel. However, after normalisation the R and G intensities are > usually reported on a log2 scale (from 0-16), and since > log2(R/G)=log2(R)-log2(G) that's probably what you see. > > HTH > \Heidi > > > Dear List I've a very basic question about the limma output (and I'm > > newbie > > of MA analysis). > > I'm looking some analysis output .Usually in the M column (if I well > > understand ) I should expect to find the log2-fold change of my > > conditions.Instead I find the difference between the normalized > > expressions... > > Some one could suggest to me how to transform these values in log2 fold > > change? > > Thank you very much and sorry for my silly question > > Guido Leoni > > > > [[alternative HTML version deleted]] > > > > _______________________________________________ > > Bioconductor mailing list > > Bioconductor@stat.math.ethz.ch > > https://stat.ethz.ch/mailman/listinfo/bioconductor > > Search the archives: > > http://news.gmane.org/gmane.science.biology.informatics.conductor > > > > > -- Guido Leoni National Research Institute on Food and Nutrition (I.N.R.A.N.) via Ardeatina 546 00178 Rome Italy tel + 39 06 51 49 41 (operator) + 39 06 51 49 4519 (direct) [[alternative HTML version deleted]]
ADD REPLY
0
Entering edit mode
On 28 Apr 2010, at 11:11, Guido Leoni wrote: > Thank you Heidi and Neel > I'll give you more details about my data: > The experiments are performed with an affymetrix custom 3' array > here is an example for one probeset : > Normalized expression with gcRMA > control > control control case case case > > 219070_s_at 679,15 734,98 777,32 667,32 507,76 645,37 > So, you data is definitely not on a log2 scale, since it'd then be from 0-16 (rather than 2^0-2^16). limma expects the MAList (or the other inputs to lmFit) to contain log2 values, c.f. the documentation for ?lmFit lmFit (object,design=NULL,ndups=1,spacing=1,block=NULL,correlation,weights=N UL L,method="ls",...) Arguments: object: object of class ?numeric?, ?matrix?, ?MAList?, ?EList?, ?marrayNorm?, ?ExpressionSet? or ?PLMset? containing log-ratios or log-values of expression for a series of microarrays Without knowing exactly what you did, i.e. seeing your actual code, I'm afraid this list can't give you any more specific help. But this is the source of the unexpected behaviour you see for your log2 fold changes. HTH \Heidi > This is the corrisponden row from limma output for the same probeset > AFFy ID EG Symbol M > A t > P.value B > 219070_s_at 64598 MOSPD3 -124.008563236060 666.759701175327 > -2.27688122341071 0.0489596557364987 -4.19921794672031 > > > > 2010/4/27 Heidi Dvinge <heidi at="" ebi.ac.uk=""> > Hello Guido, > > the M-values are indeed the log2 fold change, so log2(R/G) for the > red and > green channel. However, after normalisation the R and G intensities > are > usually reported on a log2 scale (from 0-16), and since > log2(R/G)=log2(R)-log2(G) that's probably what you see. > > HTH > \Heidi > > > Dear List I've a very basic question about the limma output (and > I'm > > newbie > > of MA analysis). > > I'm looking some analysis output .Usually in the M column (if I well > > understand ) I should expect to find the log2-fold change of my > > conditions.Instead I find the difference between the normalized > > expressions... > > Some one could suggest to me how to transform these values in > log2 fold > > change? > > Thank you very much and sorry for my silly question > > Guido Leoni > > > > [[alternative HTML version deleted]] > > > > _______________________________________________ > > Bioconductor mailing list > > Bioconductor at stat.math.ethz.ch > > https://stat.ethz.ch/mailman/listinfo/bioconductor > > Search the archives: > > http://news.gmane.org/gmane.science.biology.informatics.conductor > > > > > > > > -- > Guido Leoni > National Research Institute on Food and Nutrition > (I.N.R.A.N.) > via Ardeatina 546 > 00178 Rome > Italy > > tel + 39 06 51 49 41 (operator) > + 39 06 51 49 4519 (direct)
ADD REPLY

Login before adding your answer.

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