I'm working with data from a crossover disease study to run limma on. The data would look something like this:
If I wanted to create a design matrix where I cared about how gene expression changes with condition interacting with time, would it look something like:
model.matrix(~Person + Condition + Time + Condition:Time) or model.matrix(~Condition*Time)?
