This is a test version of Biostars. For the public version, visit https://www.biostars.org.
different results with normal fit and contrast matrix fit

I have a data set with loge base expression values of 90 normal and 10 tumor sample.

mydesign<- cbind(Tumor=c(rep(1,90),rep(0,10)),Control=c(rep(0,90),rep(1,10)))

Now, I proceeded with two way:

fit <- lmFit(data, design=mydesign)
fit <- eBayes(fit)
topTable(fit)

numGenes <- nrow(data)
case1.list <- topTable(fit,number=numGenes)

or

cont.matrix <- makeContrasts(TumorvsControl=Tumor-Control, levels=mydesign)
fit2 <- lmFit(data, design=mydesign)
fit2 <- contrasts.fit(fit2, cont.matrix)
fit2 <- eBayes(fit2)

case2.list <- topTable(fit2,number=numGenes)

However, when I searched for a particular gene, the results are quite contrasting. e.g.

case1.list["ENSG1",]
                  Tumor  Control AveExpr        F      P.Value    adj.P.Val
ENSG1 1.53268 5.134342 1.87246 1802.497 1.115791e-83 4.074213e-83

case2.list ["ENSG1",]
                    logFC AveExpr         t      P.Value    adj.P.Val        B
ENSG1 -3.601662 1.87246 -29.42565 2.183333e-53 3.353327e-49 110.5354

Could anyone please advice that with given tumor and control estimation, how log FC were obtained?

r limma
• 4,570 views • link
updated by Ram Ram • written by Ice P • Ice P

"contrast" in limma means something similar to "What conditions I want to compare".

• 0 views • link
updated by Ram Ram • written by Antonio R. Franco Antonio R. Franco

2 answers

This ends up not being a limma question, but more of a question of what a model matrix is and how it's used.

Case 1 is asking, "Is the signal in the 'Tumour' group greater than 0?" That is unlikely to be a biologically interesting question to ask.

Case 2 is asking, "Is there a difference between 'Tumour' and 'Control'?" That is a much more biologically interesting question to ask.

I would strongly discourage you from manually making model matrices until you're more familiar with them (and even then it's rarely more convenient to do so).

• 0 views • link
written by Devon Ryan Devon Ryan

This should be the accepted answer.

• 0 views • link
written by russhh russhh

Please read the documentation for limma. Start with Chapter 9 on page 40.

• 0 views • link
written by Kamil Kamil

Log in to answer this question.

More posts like this
  • Limma for DEG analysis
    written by rahma.khoualdia • rahma.khoualdia
    Hello everyone ! this is my first time using the limma package for running DEG analysis and I'm a bit lost cause when I compare …
  • Single sample vs. a control group in limma
    written by druggable • druggable
    Hi All, I was reading a paper that performed differential expression analysis for each sample vs. samples from a control group. Here was the code …
  • topTable doesn't have LogFc and B column
    written by xxxxxxxx • xxxxxxxx
    Following this [workflow][1]. topTable contains these data- topTable(fit2, coef=1) logFC AveExpr t P.Value adj.P.Val B XYZ_11111 -0.5288943 3.348961 -31.490438 1.904000e-20 7.676928e-17 33.060675 XYZ_54843 -0.6151814 3.428361 …
  • LIMMA time-dependent analysis
    written by anoopta88 • anoopta88
    Hi, I want to analyse time-series data. My design is as pasted below, sample group H24_Untreated_R1 H24_Untreated H24_Untreated_R2 H24_Untreated H24_Untreated_R3 H24_Untreated H24_Treated_R1 H24_Treated H24_Treated_R2 H24_Treated …
  • limma paired t test
    written by xaimi0813 • xaimi0813
    Hello guys, One question I have met recently is that when i handle the GEO data(GSE100186) with limma package, the result by t~test I could …
  • gene symbol in r
    written by miravet65 • miravet65
    I am having a problem in getting gene symbol after analysis GEO data. I can get logFC, P.value but in gene.symbol only show `<NA>` gr<-as.factor(gr) …
  • Tutorial: Importance of Array Quality Control - arrayQualityMetrics (PART I)
    written by majuang66 majuang66
    Today's topic is the importance of array quality control. Many articles and websites related to microarray analysis have demonstrated that array quality control is privotal …
  • Error in R commands
    written by 786 786
    I'm working with **GSE12643** and want to extract differentially expressed gene.Can any one help me in removing the error... study4 = "GSE12643" gse12643dat = getGEO('GSE12643') …
  • Microarray Data Analysis Using Bioconductor
    written by H@rry • H@rry
    <p>Hi,</p> <p>I am analyzing Microarray data using R Bioconductor. I am using the following script. But the problem all the log2fc are in -ve (mostly). …
  • What Are Usual Logfc Values In Differential Ge Analysis?
    written by munch • munch
    <p>Hello, I want compare GE from experiment A (three replicates) with experiment B (three replicates)</p> <p>Here is my workflow: </p> <pre><code>library(limma) celFiles &lt;- list.celfiles(path_to_cel_files,full.names=TRUE) affyExpressionFS …