Entering edit mode
Guest User
★
13k
@guest-user-4897
Last seen 10.5 years ago
Hi!
I would like to use Limma to compare gene expression between two
treatment groups (PO vs C). In this analysis I need to adjust for
differences in a quantitative covariate (age) between samples. Is the
following setup appropriate for this analysis? Do I accurately adjust
for age in the final analysis?
Ingrid
-- output of sessionInfo():
eset<-readExpressionSet("eset.txt","target.txt",header=TRUE)
GROUP <- factor(target$GROUP, levels=c("C","PO"))
AGE <- factor(target$AGE)
design <- model.matrix(~0+GROUP+AGE)
colnames(design) <- c("C","PO")
fit <- lmFit(eset,design)
cont.matrix <- makeContrasts(CvsPO=C-PO,levels=design)
fit2 <- contrasts.fit(fit, cont.matrix)
fit2 <- eBayes(fit2)
topTable(fit2, n=100, coef="CvsPO", adjust="BH")
--
Sent via the guest posting facility at bioconductor.org.