Remove batch effects

In the previous exercise, you demonstrated that the effect of batch had a larger impact on the variation compared to the effect of treatment. Fortunately, the olfactory stem cell study was perfectly balanced, i.e. each treatment was included in all 4 of the batches. Thus you can remove the variation introduced by the batch processing to increase the signal-to-noise ratio.

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 olfactory stem cell data has been loaded in your workspace.

Use removeBatchEffect to remove the effect of the 4 batches from the data.
Use plotMDS to plot the principal components. Label the samples by the treatment they received.
Re-visualize the principal components, labeling the samples by their batch.

Esercizio pratico interattivo

Prova a risolvere questo esercizio completando il codice di esempio.

# Load package
library(limma)

# Remove the batch effect
exprs(eset) <- ___(eset, batch = ___)

# Plot principal components labeled by treatment
___(eset, labels = ___, gene.selection = ___)

# Plot principal components labeled by batch
___(eset, labels = ___, gene.selection = ___)

Modifica ed esegui il codice