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 run LIMMA without the error term. How can I include the error term in LIMMA? I greatly appreciate your suggestions. Thanks. Best regards, [[alternative HTML version deleted]]