Hi,
I am analyzing Microarray data using R Bioconductor. I am using the following script. But the problem all the log2fc are in -ve (mostly). I am not sure what's going wrong. Can someone help me to modify it.
library(limma)
targets <- readTargets("affy_targets.txt")
library(affy)
data <- ReadAffy(filenames=targets$FileName)
eset <- rma(data)
pData(eset)
write.exprs(eset, file="(affy_salt_all.txt")
design <- model.matrix(~ -1+factor(c(1,1,2,2,3,3)))
colnames(design) <- c("group1", "group2", "group3")
fit <- lmFit(eset, design)
contrast.matrix <- makeContrasts(group2-group1, group3-group2, group3-group1, levels=design)
fit2 <- contrasts.fit(fit, contrast.matrix)
fit2 <- eBayes(fit2)
topTable(fit2, coef=1, adjust="fdr", sort.by="B", number=10)
write.table(topTable(fit2, coef=1, adjust="fdr", sort.by="B", number=50000), file="limma_salt_complete.xls", row.names=F, sep="\t")
Output
"ID" "logFC" "AveExpr" "t" "P.Value" "adj.P.Val" "B"
"249897_at" 0.602237312171529 5.64509582185098 6.73094814125857 0.00162896558995517 0.876788085929335 -3.5957825791825
"249188_at" 0.44578148174508 7.25973883239111 5.68295138543025 0.00326382351423818 0.876788085929335 -3.65128570595124
.... .... ......... ...... ........... ..........
1 answer
The log2FC are in the column marked logFC and are positive for the two example probesets you have shown.
From your code, I am not convinced that anything is wrong. It is certainly possible to have more down-regulated genes than up-regulated.
thanks alot for your reply. I agree that there could be more chances of occurring down-regulated genes in data as compared to up- regulated.
Now the problem is I am using cel files from one of the experiment published by some group and they have demonstrated the results using Affymetrix Microarray Suite 5.0. The analysis is based on two time period i.e 10 and 23 dpi and has 12 .cel file (1,2,3,7,8,9) for 10 dpi and (4,5,6,10,11,12) for 23 dpi. So in their analysis they have shown Upregulation >20 fold at 23 vs. 10 days postinfection and Downregulation <0.05 fold at 23 vs. 10 days postinfection. I am not sure how can I implement my code because I am using this approach for all .cel files (1,2,3,7,8,9) for 10 dpi in one run and (4,5,6,10,11,12) for 23 dpi in another run. Then substraction of 23 values vs 10 values (logfc) for each probeset. I wondering if I am doing the things in right way.
But after analysis I am getting the values like (for 1 probeset only) "247718_at" -2.60385457687538 7.36250717411517 -1.33369969784976 0.246361557614825 0.84824112853189 -4.53046067905237 (23 dpi)
"247718_at" -0.0743259153248443 5.00612574664665 -0.667284645810624 0.53718062033035 0.886783970000752 -4.75126236806355 (10 dpi)
and in article they mentioned like "247718_at" AT5G59310 663.576923 (Upregulation >20 fold at 23 vs. 10 dpi)
I am not sure I can I proceed for such output.
You may want to contact the authors if there are discrepancies between your analysis and theirs. I don't think we can help you much here.
Another upvote for Biostars!
Log in to answer this question.
-
DE analysis in paired samples with limma
written by dv_f •
-
topTable doesn't have LogFc and B column
written by xxxxxxxx •
-
Why different DEGs between my result and the author's result?
written by fernardo
-
Setting up my design matrix
written by marion.ryan •
-
Why Could I Only Get Two Different Expressed Genes By Limma?
written by chenyongzi77 •
-
What Are Usual Logfc Values In Differential Ge Analysis?
written by munch •
-
Affymetrix Data Using Mas5 Giving Large Fc Values
written by Bade •
-
Microarray Relative Expression Disease V Controls
written by alicia.berkenwood •
-
Creating Heatmaps Of Diffrentially Expressed Genes In Agilent Two-Color Array Analysis
written by Ly •
-
Creating Heatmaps Of Diffrentially Expressed Genes In Angilent Single-Color Array Analysis
written by Tonig
