[BioC] split-plot factorial design with limma?
Valeria Edefonti
valeria.edefonti at unimi.it
Thu Jan 11 15:49:44 CET 2007
Hi all,
I'm analyzing an experiment using 10 Affymetrix chips.
I have 5 subjects, who are either healthy or having two different
illnesses, TTD or CS.
Each subject is considered twice, at different levels of growth, P or
D, of some of his cells.
Pairs are A and G, B and H and so on:
>targets
GROWTHLEVEL ILLNESS
A P H
B P H
C P H
D P TTD
F P CS
G D H
H D H
I D H
L D TTD
N D CS
I adopted a factorial design with a blocking variable, blocco, for
identifying the five subjects.
The code follows:
SS<-paste(targets$GROWTHLEVEL, targets$ILLNESS, sep=".")
SS<-factor(SS, levels=c("P.H","P.TTD","P.CS","D.H","D.TTD","D.CS"))
design.ch<-model.matrix(~0+SS)
colnames(design.ch)<-levels(SS)
rownames(design.ch)<-rownames(targets)
blocco<-c(1,2,3,4,5,1,2,3,4,5)
corfit<-duplicateCorrelation(MA, design=design.ch,block=blocco)
fit<-lmFit(MA, design.ch, block=blocco, cor=corfit$consensus)
contrast.matrix<-makeContrasts(H.PvsD=P.H-D.H, TTD.PvsD=P.TTD-D.TTD,
CS.PvsD=P.CS-D.CS, P.TTDvsH=P.TTD-P.H, P.CSvsH=P.CS-P.H,
D.TTDvsH=D.TTD-D.H, D.CSvsH=D.CS-D.H,levels=design.ch)
fit2<-contrasts.fit(fit, contrast.matrix)
fit3<-eBayes(fit2)
I guess so far so good, right?
However, a two treatment split-plot factorial design may be slightly
better for my case, though I'm definitely not an expert in this field.
I see three groups of blocks:
Healthy: 3 blocks
TTD: TTD block
CS: CS block,
each one corresponding to a different level of factor "illness".
Is it a good idea?
Is there a way to perform such an analysis with limma?
If not, can I still stick to my current solution?
I know I can use lme() or aov(), but I should modify them to moderate
the random effects across genes.
Thank you very much for your help.
Best wishes,
Valeria
More information about the Bioconductor
mailing list