This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to extract common genes that are DE in each comparisons?

I have made a Contrasts matrix as 

design <- model.matrix(~0+f)
colnames(design) <- levels(f)
design
fit <- lmFit(eset,design)
names(fit)
cont.matrix <- makeContrasts(C9T9="Control9-Treated9",C24T24="Control24-Treated24",........................................,levels=design)
cont.matrix
fit2  <- contrasts.fit(fit, cont.matrix)
fit2  <- eBayes(fit2)
colnames(fit2)
topTable(fit2,coef=1)
topTable(fit2,coef="C9T9")

I have total of 8 comparisons and successfully extracted up and down-regulated genes from each comparisons using topTable, beside also identified the total number of up and down-regulated genes from each comparisons using

summary(decideTests(fit2))

but I am facing problem to identify/extract common genes that are DE in each of comparisons. Therefore, It is my kind request to all of you, please provides any suggestions for extracting of common genes that are DE in each of the comparisons.

Thanking you.

Kind regards,

limma bioiconductor r microarray data analysis
• 2,822 views • link
updated by GenoMax GenoMax • written by rkpathakbt • rkpathakbt

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • 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 …
  • Is it possible to apply limma on metabric data?
    written by DareDevil DareDevil
    I have downloaded the microarray expression data afrom metabric and followed below mentioned steps. Is it the right way to perform the DEG analysis? library(limma) …
  • LIMMA: ~ continuous predictor by group. How to get within-group results?
    written by sindrelee • sindrelee
    Cross-posted: https://support.bioconductor.org/p/133284/ Please see the example below. set.seed(1) pheno=cbind.data.frame(IDs=as.factor(c(1:20)), gr=as.factor(c(rep("A",10),rep("B",10))), oc=rnorm(20) ) dat=rbind.data.frame(varA=rnorm(20), varB=rnorm(20), varC=rnorm(20), varD=rnorm(20) ) colnames(dat)=IDs rownames(pheno)=IDs form=as.formula(paste0("~oc*gr")) design=model.matrix(form,data=pheno) colnames(design)=make.names(colnames(design)) colnames(design) [1] "X.Intercept." …
  • 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 …
  • Limma package, how to correct by age and sex?
    written by ellen2270 • ellen2270
    Hello everyone, I am analyzing microarray data with limma package and I have a couple of doubts as it is the first time I use …
  • Batch correction and DE analysis of miRNA Seq data
    written by georgians • georgians
    Hello everyone, I want to do the batch correction and DE analysis. I have miRNA seq read count for 448 samples. Total no of batches …
  • 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) …
  • Setting up my design matrix
    written by marion.ryan • marion.ryan
    I think I am having an issue with my contrast matrix, in terms of comparing concentrate versus grass fed animals (RMA normalised = eset). I …
  • Heatmap for Differential expressed genes
    written by Vasu Vasu
    Hi, I want to generate heatmap for differential expressed genes after using limma. But I got the following error. Group <- c("GrpA","GrpA","GrpA","GrpB","GrpB","GrpB") design <- model.matrix(~factor(Group)) …
  • limma 2 factor contrast matrix design
    written by Ming Tommy Tang Ming Tommy Tang
    Hi biostars, I have a question to ask for the limma two factor design. my experiment design is like this: ``` FileName Dox KnockDown L …