生物信息学：对比Limma如何与3组（病例/对照，年龄和性别）

Question

我是学生，我怀疑如何使用 Limma 来对比评估我的三个变量的差异基因表达：疾病（病例/对照），性别（男/女）和年龄。

我的数据是这样的:(三栏中有性别，年龄和疾病）

  

head（pData（gse））

数据：

我知道在变异疾病上建立我的实验：

    design= model.matrix(~ pData(gse)[,"disease"])
    colnames(design)= c("control", "case")

    fit= lmFit(gse, fesign)
    fit1= eBayes(fit)

    topTable(fit1, coef=2, adjust="BH", number= as.numeric(dim(gse)[1]))

但是，如何用我的3个变量做对比Limma？谢谢你的帮助！

age= pData(gse)[,"age"]
gender= pData(gse)[,"gender"]
disease= pData(gse)[,"disease"]

age.group= ifelse(age>40,"over40", "under40")

design = model.matrix(~0 + disease + gender + gender:age)

contr.matrix <- makeContrasts(  **???**