IniziaInizia gratis

Specify a linear model to compare 2 groups

To identify differentially expressed genes for the leukemia experiment, you need to translate the following linear model to R:

where \(X_{1}\) is equal to 1 for progressive cancers and 0 for stable cancers (note: R automatically chooses the base condition by alphabetical order).

Questo esercizio fa parte del corso

Differential Expression Analysis with limma in R

Visualizza il corso

Istruzioni dell'esercizio

The ExpressionSet object eset with the leukemia data has been loaded in your workspace.

  • Use model.matrix to construct a design matrix with an intercept coefficient and a coefficient that indicates the disease status.

Esercizio pratico interattivo

Prova a risolvere questo esercizio completando il codice di esempio.

# Create design matrix for leukemia study
design <- ___(~___, data = ___(eset))

# Count the number of samples modeled by each coefficient
colSums(design)
Modifica ed esegui il codice