Design matrix for group-means model
In the previous chapter, you tested the leukemia data for differential expression using the traditional treatment-contrasts parametrization. As a first step to learning the more flexible group-means parametrization, you will re-test the leukemia data to confirm you obtain the same results.
Questo esercizio fa parte del corso
Differential Expression Analysis with limma in R
Istruzioni dell'esercizio
The ExpressionSet object eset with the leukemia data has been loaded in your workspace.
- Use
model.matrixto create a design matrix with no intercept. Recall that the variable of interest for this study (progressive vs. stable cancers) is in the columnDiseaseof the phenotype data frame.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# Create design matrix with no intercept
design <- ___(~___ + ___, data = ___(eset))
# Count the number of samples modeled by each coefficient
colSums(design)