For paired data analysis using limma with multiple factor, if I was interested in the interaction effect between Condition and Time would the structure look like this?
mm <- model.matrix(~Patient + Condition*Time + Batch)
b) if you're doing a paired analysis then you're anyway comparing within patient so batch is not needed as the two samples per patient are always in the same batch.
Do ~Patient+Condition to test condition, accounting for patient differences.
Apologies, I just realized I didn't upload the data correctly and it would be moreso in this format:![enter image description here](/media/images/6808e793-38fd-4c01-888d-445e255c)
As ATpoint already noted, you have nested patient in batch, so the model should be
~Patient + Time*Condition
.