limma paired design
1
0
Entering edit mode
Michael • 0
@d4da273d
Last seen 8 months ago
Canada

enter image description here

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)
DifferentialExpression limma • 639 views
ADD COMMENT
0
Entering edit mode
ATpoint ★ 4.1k
@atpoint-13662
Last seen 1 hour ago
Germany

This design is not full rank:

a) condition and time is the same, drop one of it

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.

ADD COMMENT
0
Entering edit mode

Apologies, I just realized I didn't upload the data correctly and it would be moreso in this format: enter image description here

ADD REPLY
0
Entering edit mode

As ATpoint already noted, you have nested patient in batch, so the model should be ~Patient + Time*Condition.

ADD REPLY

Login before adding your answer.

Traffic: 476 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6