Question: Setting up my design matrix
2
gravatar for marion.ryan
2.6 years ago by
marion.ryan20
marion.ryan20 wrote:

I think I am having an issue with my contrast matrix, in terms of comparing concentrate versus grass fed animals (RMA normalised = eset). I get an output of DE genes, but none of them are differently expressed (Adj P value) and the gene list is different from that of a previous and I believe correct analysis, by an experienced bioinformatition who has now left. Any insight in to this would be greatly appreciated, code below. 1st 7 CEL files =Concentrate 2nd 7 CEL files =Grass

Matrix design

###################################################
design <- model.matrix(~ 0+factor(c(1,1,1,1,1,1,1,2,2,2,2,2,2,2)))
colnames(design) <- c("Concentrate", "Grass")
fit <- lmFit(eset, design)**

###################################################
###Contrast Matrix-new 
###################################################
contrast.matrix <- makeContrasts(Concentrate - Grass, levels=design)
fit2 <- contrasts.fit(fit, contrast.matrix)
fit2 <- eBayes(fit2) 

###################################################
###Output –differentially expressed genes 
###################################################
topTable(fit2,coef=NULL,number=500,genelist=fit2$genes,adjust.method="BH",sort.by="B",resort.by=NULL,p.value=1,lfc=0)
tt=topTable(fit2, number = 150)
write.table(tt, "ConcentratevGrass.txt", sep = "\t")
R • 1.4k views
ADD COMMENTlink modified 2.5 years ago by Biostar ♦♦ 20 • written 2.6 years ago by marion.ryan20
1

It all looks fine to me, aside from the '**' after lmFit. Was there any blocking used in setting up the experiment?

ADD REPLYlink written 2.5 years ago by russhh4.1k
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 2.3.0
Traffic: 1340 users visited in the last hour