Entering edit mode
Sharon Anbu
▴
480
@sharon-anbu-1524
Last seen 10.2 years ago
Hi,
I have 2-factors (treatment and task) with repeated measures design.
Each
subject has to preform 2 different tasks before and after treatment.
We have
5 subjects, so it is 20 arrays in total. For each probeset, I have
tried
with
aov.res <- aov( expr[1,] ~ (treat * task) + Error( subject / (treat *
task)
), rma_norm )
Now, I would like to translate the above model into LIMMA. I was bale
to do
run LIMMA without the error term. How can I include the error term in
LIMMA?
I greatly appreciate your suggestions.
Thanks.
Best regards,
Sharon
[[alternative HTML version deleted]]