Question: Interpreting Fold Change Calculation in Limma to Non-Statisticians: Sanity Check
0
gravatar for JMallory
4 weeks ago by
JMallory0
JMallory0 wrote:

I have seen this question posed elsewhere, but wanted to double check for some of my colleagues who want to make sure. It is the following: when defining contrasts in Limma, how is the syntax to be interpreted? I've included some code. If clarification is required, I can provide that as needed.

#### edata is expression matrix of aptamer probes; design is the design matrix; Used a repeated measures design so estimated the patient-wise correlation in expression over visits where StudyID defines this pairing ; pheno_table details sample metadata 

fit<-lmFit(edata,design,
               block=pheno_table$StudyID,
               correlation=corfit$consensus)

cm <- makeContrasts(
    `Diabetpost.PUFA.Met-Diabetpre.PUFA.Met` = Diabetpost.PUFA.Met-Diabetpre.PUFA.Met,
    `Diabetpost.MUFA.Met-Diabetpre.MUFA.Met` = Diabetpost.MUFA.Met-Diabetpre.MUFA.Met,
        `Diabetpost.PUFA.NoMet-Diabetpre.PUFA.NoMet` = Diabetpost.PUFA.NoMet-Diabetpre.PUFA.NoMet,
    `Diabetpost.MUFA.NoMet-Diabetpre.MUFA.NoMet` = Diabetpost.MUFA.NoMet-Diabetpre.MUFA.NoMet,
     levels=design)

fit2 <- contrasts.fit(fit, cm)
fit2 <- eBayes(fit2,robust=TRUE)
T<-topTable(fit2,number=100000)


In plain English, are the above comparisons made "between post with respect to pre" or "between pre with respect to post"? That is, would a positive logFC indicate post > pre or pre > post? Again, I believe I know the answer, but I would like confirmation for my colleagues. Thanks.

proteomics limma • 125 views
ADD COMMENTlink modified 25 days ago by LKal0 • written 4 weeks ago by JMallory0
Answer: Interpreting Fold Change Calculation in Limma to Non-Statisticians: Sanity Check
1
gravatar for Aaron Lun
4 weeks ago by
Aaron Lun24k
Cambridge, United Kingdom
Aaron Lun24k wrote:

It's easy to understand in terms of how you called makeContrasts. Take your first comparison:

    `Diabetpost.PUFA.Met-Diabetpre.PUFA.Met` = Diabetpost.PUFA.Met-Diabetpre.PUFA.Met,

The log-fold change will be computed as Diabetpost.PUFA.Met-Diabetpre.PUFA.Met, i.e., the difference between the two coefficients. So if the log-fold change is positive, you've basically got:

Diabetpost.PUFA.Met - Diabetpre.PUFA.Met > 0 # i.e.,
Diabetpost.PUFA.Met > Diabetpre.PUFA.Met

So a positive log-fold change would represent an increase in "Post" relative to "Pre". And so forth.

ADD COMMENTlink modified 4 weeks ago • written 4 weeks ago by Aaron Lun24k

Thanks Aaron. This is exactly what I was looking for.

ADD REPLYlink written 4 weeks ago by JMallory0
Answer: Interpreting Fold Change Calculation in Limma to Non-Statisticians: Sanity Check
0
gravatar for LKal
25 days ago by
LKal0
LKal0 wrote:

I have also wondered this question. How I came to my solution was to check a single gene that had a log-fold change and then go look at the CPM values for the same gene. It is very easy to tell by this method and there is basically no chance to misinterpret the answer.

ADD COMMENTlink written 25 days ago by LKal0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 276 users visited in the last hour