I am very new to limma and trying to build a logistic regression like model using limma where my goal is to understand associations between the cohort status and gene expression while accounting for variables like age and sex. My gene expression is quantified using RSEM with TPM counts. I have normalized my data using TMM normalization. The following is my code snippet for limma.
design <- model.matrix(~0+Cohort+sex+age, data=PhenoType) #Cohort:1/0 v <- voomWithQualityWeights(myNormalized_data, design=design, normalization="none", plot=TRUE) fit1 <- lmFit(v,design) fit2 <- eBayes(fit1)
Is this design correct? Again, I am not interested in group comparison but predicting cohort status using expression values.