Contrasts matrix for group-means

Since you're using the group-means parametrization, there is now an extra step in the limma pipeline. The coefficients now represent the mean expression level for each of the two groups of samples, so you need to specify a custom contrast to test for differences between the progressive and stable leukemias.

Diese Übung ist Teil des Kurses

Differential Expression Analysis with limma in R

Kurs anzeigen

Anleitung zur Übung

The ExpressionSet object eset with the leukemia data and the design matrix you just created (design) have been loaded in your workspace.

Use makeContrasts to define the contrast status, which is the difference in the mean expression level between the progressive and stable cancers.
Pass the design matrix design to the argument levels. This allows the column names of design (which correspond to the coefficients) to be used in the definitions of the contrasts.

Interaktive Übung

Versuche dich an dieser Übung, indem du diesen Beispielcode vervollständigst.

# Load package
library(limma)

# Create a contrasts matrix
cm <- ___(status = ___ - ___,
                    levels = ___)

# View the contrasts matrix
cm

Code bearbeiten und ausführen