Displaying 1 result from an estimated 1 matches for "dmso5".
Did you mean:
dmso
2012 May 04
0
LIMMA decideTests result zero from contrast matrix
...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)
>...