I have been using limma to perform linear regression to find genes differentially expressed with age, also include sex in the model.
Here is my model.
design <- model.matrix(~0 + age + sex)
My questions are:
1. I just want to make sure whether the model is make sense?
2. Does the logFC got from topTable mean the fold change of gene expression for each year of age?