Dear All,
I am using the LIMMA package to create 2 contrasts for my data and then
calculating the vennCounts of the decideTests from the contrast.fit to be able
to create venn Diagrams.
The code works fine but the summary(results) shows zeros for all i.e. no gene
were up regulated or downregulated. This is not true for my data since toptable
output shows Log fold change greater than > 2.
I am certain it is a small glitch somewhere at my end that i get zero counts
for my summary(decideTests). Please find below my code and i would really
appreciate any help here at all.
Thanks,
Ekta
## R Script ###
> numGenes <- nrow(eset)
> library(limma)
> samples <- c("Un","Un","DMSO10","DMSO10","DMSO5","DMSO5");
> fl <- as.factor(samples)
> design <- model.matrix(~ 0+ fl)
> colnames(design) <- levels(fl)
> fit <- lmFit(eset, design)
> cont.matrix <- makeContrasts(DMSO10-Un, DMSO5-Un, levels=design)
> fit2 <- contrasts.fit(fit, cont.matrix)
> fit2 <- eBayes(fit2)
> tTUni<- topTable(fitUni, adjust="fdr", sort.by="B", number=numGenes)
> results <- decideTests(fit2)
> vennDiagram(results,include=c("up","down"),counts.col=c("red","green"))
## ### I get zero genes for upregulation and downregulation #####
Senior Research Associate
Bioinformatics Department
Jubilant Biosys Pvt Ltd,
#96, Industrial Suburb, 2nd Stage
Yeshwantpur, Bangalore 560 022
Ph No : +91-80-66628346
The information contained in this electronic message and in any attachments to
this message is confidential, legally privileged and intended only for use by
the person or entity to which this electronic message is addressed. If you are
not the intended recipient, and have received this message in error, please
notify the sender and system manager by return email and delete the message and
its attachments and also you are hereby notified that any distribution,
copying, review, retransmission, dissemination or other use of this electronic
transmission or the information contained in it is strictly prohibited. Please
note that any views or opinions presented in this email are solely those of the
author and may not represent those of the Company or bind the Company. Any
commitments made over e-mail are not financially binding on the company unless
accompanied or followed by a valid purchase order. This message has been
scanned for viruses and dangerous content by Mail Scanner, a!
nd is believed to be clean. The Company accepts no liability for any damage
caused by any virus transmitted by this email.
www.jubl.com
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
